The biest is the trilinear Fourier multiplier operator T(f1,f2,f3)(x)=∫ξ1<ξ2<ξ3f^1(ξ1)f^2(ξ2)f^3(ξ3)e2πix(ξ1+ξ2+ξ3) dξ1dξ2dξ3T(f_1,f_2,f_3)(x)=\int_{\xi_1<\xi_2<\xi_3}\hat f_1(\xi_1)\hat f_2(\xi_2)\hat f_3(\xi_3)e^{2\pi ix(\xi_1+\xi_2+\xi_3)}\,d\xi_1d\xi_2d\xi_3, whose symbol is the characteristic function of a frequency simplex; it emerged from the study of eigenfunction expansions of one-dimensional Schrodinger operators. The 2004 paper 'L^p estimates for the biest I. The Walsh case' proved Holder-type bounds for the dyadic (Walsh) model on an open exponent polytope DD and explicitly deferred three directions: the Fourier model, endpoint behavior on the boundary of DD, and higher-order ('triest'/'multiest') analogues; the classification of reparametrized simplex multipliers was subsequently posed as an open question by the later literature (the 2003 Christ-Kiselev counterexample paper, the multiest paper, and Kesler's papers). The Fourier case was resolved by the 2004 sequel (and extended to arbitrary product symbols by Jung in 2017); the reparametrized-multiplier question is largely delineated through degeneracy classes (degenerate tuples admit no LpL^p estimates, non-degenerate tuples a wide range, semi-degenerate tuples a sharp intermediate range); and the nn-linear simplex operators are bounded from L2×⋯×L2L^2\times\cdots\times L^2 into L2/nL^{2/n} for every nn. The surviving open core comprises: (a) endpoint and weak-type behavior of the biest exactly on the boundary of DD, in both the Walsh and Fourier models; (b) completion of the parameter delineation to all tuples and all degrees with sharp exponents; and (c) L2L^2/LpL^p estimates for the maximal Carleson-type ordering operators MnM^n with n≥3n\ge3, announced as forthcoming by the authors in 2007, never published, and identified by them as the obstruction to pointwise convergence of eigenfunction expansions for general upper-triangular systems.

Background

A multilinear Fourier multiplier operator acts on several input functions by multiplying their Fourier transforms by a fixed function of the frequency variables (the symbol) and transforming back. The emblematic case is the bilinear Hilbert transform H(f,g)(x)=p.v.∫Rf(x+y)g(x−y) dy/yH(f,g)(x)=\mathrm{p.v.}\int_{\mathbb{R}}f(x+y)g(x-y)\,dy/y, whose symbol is singular along a diagonal in frequency space; Lacey and Thiele proved that it maps Lp×LqL^p\times L^q boundedly into LrL^r whenever 1<p,q≤∞1<p,q\le\infty and 1/p+1/q=1/r1/p+1/q=1/r with 2/3<r<∞2/3<r<\infty, resolving a 1960s conjecture of A. P. Calderon. The present problem concerns the next orders of this theory.

The biest operator of Muscalu, Tao and Thiele is the trilinear Fourier multiplier

T(f1,f2,f3)(x)=∫ξ1<ξ2<ξ3f^1(ξ1)f^2(ξ2)f^3(ξ3)e2πix(ξ1+ξ2+ξ3) dξ1dξ2dξ3,T(f_1,f_2,f_3)(x)=\int_{\xi_1<\xi_2<\xi_3}\hat f_1(\xi_1)\hat f_2(\xi_2)\hat f_3(\xi_3)e^{2\pi ix(\xi_1+\xi_2+\xi_3)}\,d\xi_1d\xi_2d\xi_3,

whose symbol is the characteristic function of the frequency simplex ξ1<ξ2<ξ3\xi_1<\xi_2<\xi_3; the name (a deliberate misspelling of 'beast') signals that it was the first genuinely trilinear time-frequency operator to be estimated. It is not an isolated curiosity: as the source paper explains, multilinear singular integrals of exactly this type occur in the WKB-type eigenfunction expansions of one-dimensional Schrodinger operators (WKB, after Wentzel, Kramers and Brillouin: the high-frequency asymptotic expansion of the eigenfunctions) and, more generally, of AKNS systems (the integrable-systems framework of Ablowitz, Kaup, Newell and Segur that contains Schrodinger and Dirac operators as special cases) - iterating the expansion produces a multilinear series whose terms are multipliers supported on ordering regions of the frequencies. The bilinear member of this family is essentially the bilinear Hilbert transform, and the trilinear member is the biest. In the Schrodinger application the inputs are L2L^2 functions, so the desired bound is L2×L2×L2→L2/3L^2\times L^2\times L^2\to L^{2/3} (note 2/3<12/3<1: the target is quasi-Banach).

The source paper proved its theorem in the Walsh model. The Walsh model is a dyadic toy version of Fourier analysis in which the characters are the Walsh functions, {0,1}\{0,1\}-valued analogues of the exponentials built by repeatedly dilating and reflecting χ[0,1)\chi_{[0,1)}, so that wave packets are exactly localized on quartiles - time-frequency rectangles I×ωI\times\omega that are products of dyadic intervals with area ∣I∣ ∣ω∣=4|I|\,|\omega|=4. Exact localization isolates the combinatorial core of the problem from kernel regularity, and it is standard to prove the Walsh theorem first and transfer to the Fourier setting afterwards. The Walsh biest is decomposed into two pieces (one coupling the middle input with the first, the other with the third), each analyzed by organizing quartiles into trees and controlling them through size and energy functionals, and the LpL^p bounds are then deduced by multilinear interpolation of restricted weak-type estimates, i.e. bounds tested against indicator functions of finite-measure sets, which is how quasi-Banach targets are handled.

The shape of the result is best expressed in exponent coordinates. Writing the target exponent as p4′p_4' (so that 1/p1+1/p2+1/p3=1/p4′1/p_1+1/p_2+1/p_3=1/p_4', the same homogeneity as the pointwise product of three functions) and setting α=(1/p1,1/p2,1/p3,1−1/p4′)\alpha=(1/p_1,1/p_2,1/p_3,1-1/p_4'), the tuples for which Theorem 1.8 of the source paper asserts the bounds form an open polytope D=D′∩D′′D=D'\cap D'' inside the hyperplane α1+α2+α3+α4=1\alpha_1+\alpha_2+\alpha_3+\alpha_4=1, where D′D' and D′′D'' are the interiors of the convex hulls of the twelve explicit vertices

A1=(1,1/2,1,−3/2)A2=(1/2,1,1,−3/2)A3=(1/2,1,−3/2,1)A4=(1,1/2,−3/2,1)A5=(1,−1/2,0,1/2)A6=(1,−1/2,1/2,0)A7=(1/2,−1/2,0,1)A8=(1/2,−1/2,1,0)A9=(−1/2,1,0,1/2)A10=(−1/2,1,1/2,0)A11=(−1/2,1/2,1,0)A12=(−1/2,1/2,0,1)\begin{array}{llll} A_1=(1,1/2,1,-3/2) & A_2=(1/2,1,1,-3/2) & A_3=(1/2,1,-3/2,1) & A_4=(1,1/2,-3/2,1)\\ A_5=(1,-1/2,0,1/2) & A_6=(1,-1/2,1/2,0) & A_7=(1/2,-1/2,0,1) & A_8=(1/2,-1/2,1,0)\\ A_9=(-1/2,1,0,1/2) & A_{10}=(-1/2,1,1/2,0) & A_{11}=(-1/2,1/2,1,0) & A_{12}=(-1/2,1/2,0,1) \end{array}

and of their index-swapped counterparts B1,…,B12B_1,\dots,B_{12}, each BjB_j being obtained from AjA_j by permuting the indices 1 and 3 (for instance B2=(1,1,1/2,−3/2)B_2=(1,1,1/2,-3/2)). The theorem asserts that the Walsh biest maps Lp1×Lp2×Lp3L^{p_1}\times L^{p_2}\times L^{p_3} into Lp4′L^{p_4'} for every tuple in DD; the center point (1/2,1/2,1/2,−1/2)(1/2,1/2,1/2,-1/2) is the L2×L2×L2→L2/3L^2\times L^2\times L^2\to L^{2/3} bound. Nothing is asserted on the boundary ∂D\partial D: the vertices themselves, such as the L1×L2×L1→L2/5L^1\times L^2\times L^1\to L^{2/5} configuration corresponding to A1A_1, are precisely the endpoint questions.

Two extensions of the operator frame the open problem. First, for a parameter tuple β=(β1,β2,β3)\beta=(\beta_1,\beta_2,\beta_3) with nonzero entries one can reparametrize the ordering, defining (following the later literature)

Cβ(f1,f2,f3)(x)=∫ξ1/β1<ξ2/β2<ξ3/β3f^1(ξ1)f^2(ξ2)f^3(ξ3)e2πix(ξ1+ξ2+ξ3) dξ1dξ2dξ3,C^\beta(f_1,f_2,f_3)(x)=\int_{\xi_1/\beta_1<\xi_2/\beta_2<\xi_3/\beta_3}\hat f_1(\xi_1)\hat f_2(\xi_2)\hat f_3(\xi_3)e^{2\pi ix(\xi_1+\xi_2+\xi_3)}\,d\xi_1d\xi_2d\xi_3,

so that the biest is essentially C(1,1,1)C^{(1,1,1)}; the companion papers and the subsequent literature raise the question of which β\beta admit analogous estimates. The answer is governed by the degeneracy of the partial sums of β\beta: β\beta is called degenerate if some adjacent partial sum βi+βi+1\beta_i+\beta_{i+1} vanishes, semi-degenerate if some non-adjacent partial sum vanishes (for three inputs, the full sum β1+β2+β3\beta_1+\beta_2+\beta_3) while no adjacent one does, and non-degenerate if no partial sum ∑k=ijβk\sum_{k=i}^j\beta_k vanishes. Second, there is the higher-order ('triest' or 'multiest') hierarchy: the nn-linear operator TnT^n with symbol χ{ξ1<⋯<ξn}\chi_{\{\xi_1<\dots<\xi_n\}}, and its maximal Carleson-type variants

Mn(f1,…,fn)(x)=sup⁡N∣∫ξ1<⋯<ξn<N∏j=1nf^j(ξj) e2πix(ξ1+⋯+ξn) dξ1…dξn∣,M^n(f_1,\dots,f_n)(x)=\sup_N\Big|\int_{\xi_1<\dots<\xi_n<N}\prod_{j=1}^n\hat f_j(\xi_j)\,e^{2\pi ix(\xi_1+\dots+\xi_n)}\,d\xi_1\dots d\xi_n\Big|,

in which the upper cutoff of the ordering region is replaced by a supremum over all cutoffs - the ordering-analogue of the Carleson-Hunt theorem on almost-everywhere convergence of Fourier series. For a parameter tuple β=(β1,…,βn)\beta=(\beta_1,\dots,\beta_n) with nonzero entries, the reparametrized maximal variants MβnM^n_\beta are defined by the same supremum with the ordering constraint ξ1/β1<⋯<ξn/βn<N\xi_1/\beta_1<\dots<\xi_n/\beta_n<N, with β\beta non-degenerate in the sense above (no partial sum ∑k=ijβk\sum_{k=i}^{j}\beta_k vanishes), as in the multiest paper. The maximal operators are the ones needed for pointwise (not merely L2L^2) convergence of the eigenfunction expansions of general triangular systems. By contrast, the trilinear Hilbert transform - the analogue in which the ordering constraint is replaced by a condition of the form ξ1+2ξ2+3ξ3>0\xi_1+2\xi_2+3\xi_3>0 - is described in the source paper as just barely beyond the reach of known multilinear techniques and remains open.

Problem Statement

Determine the remaining mapping properties of the biest operator and its relatives - the problem left open by 'L^p estimates for the biest I. The Walsh case' (Muscalu, Tao and Thiele, 2004), which deferred three directions: the Fourier model, endpoint behavior on the boundary of DD, and the higher-order ('triest'/'multiest') analogues. One direction is now closed: the Fourier (non-dyadic) model of the biest itself is bounded on the same polytope DD (the 2004 sequel 'L^p estimates for the biest II. The Fourier case'), and its symbol class was later extended to arbitrary sufficiently regular product symbols (Jung, 2017). The classification question in part (2) below was posed as an open question by the subsequent literature (the 2003 Christ-Kiselev counterexample paper, the multiest paper arXiv:0712.2420, and Kesler's papers). The open problem consists of three linked parts.

(1) Endpoints of the biest. For the Walsh and the Fourier biest, determine the behavior exactly on the boundary ∂D\partial D of the open exponent polytope D=D′∩D′′D=D'\cap D'' of tuples α=(1/p1,1/p2,1/p3,1−1/p4′)\alpha=(1/p_1,1/p_2,1/p_3,1-1/p_4') in the hyperplane α1+α2+α3+α4=1\alpha_1+\alpha_2+\alpha_3+\alpha_4=1 - the polytope, constructed in the Background above from the twelve explicit vertices A1,…,A12A_1,\dots,A_{12} and their index-swapped counterparts B1,…,B12B_1,\dots,B_{12}, on which Theorem 1.8 of the source paper (and, for the Fourier model, Theorem 1.2 of the sequel) asserts the bounds under 1<p1,p2,p3≤∞1<p_1,p_2,p_3\le\infty and 0<p4′<∞0<p_4'<\infty: for which boundary tuples of DD in the closure of the admissible exponent class in which the sources pose their interpolation framework (Definition 2.1 of the source paper / Definition 3.1 of the sequel: αi<1\alpha_i<1, α1+α2+α3+α4=1\alpha_1+\alpha_2+\alpha_3+\alpha_4=1, at most one negative αi\alpha_i) does a strong-type or weak-type estimate Lp1×Lp2×Lp3→Lp4′L^{p_1}\times L^{p_2}\times L^{p_3}\to L^{p_4'} (or a refinement in Lorentz spaces Lp,rL^{p,r}, which refine the LpL^p scale by recording weak-type behavior) hold - including the vertex configurations such as α=(1,1/2,1,−3/2)\alpha=(1,1/2,1,-3/2), i.e. L1×L2×L1→L2/5L^1\times L^2\times L^1\to L^{2/5}, and the facets - and at which boundary points does boundedness fail? The boundary of DD also contains tuples with two negative αi\alpha_i, which lie outside the sources' framework (the estimate there is not defined) and are excluded from this question.

(2) The reparametrized family. For the operators CβC^\beta defined by the ordering constraint ξ1/β1<ξ2/β2<ξ3/β3\xi_1/\beta_1<\xi_2/\beta_2<\xi_3/\beta_3 with β\beta having nonzero entries, and for their higher-degree analogues, complete the delineation of the boundedness regime: characterize exactly the tuples β\beta, in every degree nn, for which Holder-type estimates hold, with the sharp range of exponents in each case. The known partial delineation - non-degenerate tuples satisfy a wide range of LpL^p estimates; degenerate tuples admit none; representative semi-degenerate tuples such as (1,1,−2)(1,1,-2) are bounded into LpL^p exactly for 1/2<p<∞1/2<p<\infty - is to be completed to all degeneracy patterns, including every configuration of vanishing partial sums in degrees n≥4n\ge4.

(3) Maximal ordering operators. Prove or disprove L2L^2 (and LpL^p) estimates for the maximal Carleson-type variants MnM^n and MβnM^n_\beta for n≥3n\ge3 - the supremum over frequency cutoffs of the ordering multipliers, with non-degenerate β\beta - whose boundedness was announced as forthcoming by the authors in the 2007 multiest paper and identified there as the exact reduction target for almost-everywhere convergence of the eigenfunction (WKB) expansions of general upper-triangular systems.

References

  1. Camil Muscalu, Terence Tao, and Christoph Thiele, 'L^p estimates for the biest I. The Walsh case', Mathematische Annalen 329(3):401-426, 2004. DOI: 10.1007/s00208-004-0518-1; arXiv:math/0102083; https://arxiv.org/abs/math/0102083
  2. Camil Muscalu, Terence Tao, and Christoph Thiele, 'L^p estimates for the biest II. The Fourier case', Mathematische Annalen 329(3):427-461, 2004. DOI: 10.1007/s00208-003-0508-8; arXiv:math/0102084; https://arxiv.org/abs/math/0102084
  3. Camil Muscalu, Terence Tao, and Christoph Thiele, 'A counterexample to a multilinear endpoint question of Christ and Kiselev', Mathematical Research Letters 10(2):237-246, 2003. DOI: 10.4310/MRL.2003.V10.N2.A10; arXiv:math/0108156; https://arxiv.org/abs/math/0108156
  4. Camil Muscalu, Terence Tao, and Christoph Thiele, 'Multi-linear multipliers associated to simplexes of arbitrary length', in Advances in Analysis, pp. 346-401, Princeton University Press, 2014. DOI: 10.1515/9781400848935-015; arXiv:0712.2420; https://arxiv.org/abs/0712.2420
  5. Camil Muscalu, Terence Tao, and Christoph Thiele, 'Multi-linear operators given by singular multipliers', Journal of the American Mathematical Society 15(2):469-496, 2002. DOI: 10.1090/s0894-0347-01-00379-4; https://doi.org/10.1090/s0894-0347-01-00379-4
  6. Robert Kesler, 'Mixed estimates for degenerate multi-linear operators associated to simplexes', Journal of Mathematical Analysis and Applications 424(1):344-360, 2015. DOI: 10.1016/j.jmaa.2014.11.005; arXiv:1311.2322; https://arxiv.org/abs/1311.2322
  7. Robert Kesler, 'Generic trilinear multipliers associated to degenerate simplexes', Collectanea Mathematica 70(1):153-186, 2019. DOI: 10.1007/s13348-018-0224-z; arXiv:1609.05946; https://arxiv.org/abs/1609.05946
  8. Robert Kesler, 'Unboundedness theorems for symbols adapted to large subspaces', Studia Mathematica 244(2):109-158, 2019. DOI: 10.4064/sm8750-7-2017; arXiv:1609.05954; https://arxiv.org/abs/1609.05954
  9. Robert Kesler, 'L^p estimates for semi-degenerate simplex multipliers', Revista Matematica Iberoamericana 36(1):99-158, 2020. DOI: 10.4171/rmi/1123; arXiv:1609.05964; https://arxiv.org/abs/1609.05964
  10. Joeun Jung, 'Iterated trilinear Fourier integrals with arbitrary symbols', Journal of Functional Analysis 273(10):3027-3060, 2017. DOI: 10.1016/j.jfa.2017.07.005; arXiv:1311.1574; https://arxiv.org/abs/1311.1574
  11. Michael Lacey and Christoph Thiele, 'L^p estimates on the bilinear Hilbert transform for 2<p<infty', Annals of Mathematics (2) 146(3):693-724, 1997. DOI: 10.2307/2952458; https://doi.org/10.2307/2952458
  12. Michael Lacey and Christoph Thiele, 'On Calderon's conjecture', Annals of Mathematics (2) 149(2):475-496, 1999. DOI: 10.2307/120971; https://doi.org/10.2307/120971
  13. Francesco Di Plinio and Christoph Thiele, 'Endpoint bounds for the bilinear Hilbert transform', Transactions of the American Mathematical Society 368(6):3931-3972, 2016. DOI: 10.1090/tran/6548; arXiv:1403.5978; https://arxiv.org/abs/1403.5978
  14. Ciprian Demeter and Francesco Di Plinio, 'Endpoint bounds for the quartile operator', Journal of Fourier Analysis and Applications 19(4):836-856, 2013. DOI: 10.1007/s00041-013-9275-4; arXiv:1206.3798; https://arxiv.org/abs/1206.3798
  15. Peter M. Luthy, 'On a biparameter maximal multilinear operator', Journal of Functional Analysis 268(5):1105-1152, 2015. DOI: 10.1016/j.jfa.2014.11.010; arXiv:1409.6763; https://arxiv.org/abs/1409.6763
  16. Camil Muscalu, 'Some remarks on the n-linear Hilbert transform for n>=4', Mathematical Research Letters 21(5):1155-1176, 2014. DOI: 10.4310/mrl.2014.v21.n5.a11; arXiv:1209.6391; https://arxiv.org/abs/1209.6391
  17. Cristina Benea and Itamar Oliveira, 'The oscillatory biest operator', preprint, 2024. arXiv:2412.18079; https://arxiv.org/abs/2412.18079
  18. Polona Durcik and Christoph Thiele, 'Singular Brascamp-Lieb: A Survey', Springer INdAM Series book chapter, pp. 321-349, 2021. DOI: 10.1007/978-3-030-72058-2_9; arXiv:1904.08844; https://arxiv.org/abs/1904.08844