Rauch's hot spots conjecture predicts that on a bounded insulated domain the eigenfunction of the first nonzero Neumann eigenvalue attains its hottest and coldest points on the boundary; the strong form requires the extrema to be attained only on the boundary. Kennedy and Rohleder proved the strong form for multidimensional lip domains, and for the distinguished eigenfunctions on domains symmetric in all coordinate hyperplanes whose intersection with an orthant is a lip domain, using a purely variational, vector-valued Laplacian whose spectrum encodes the Neumann and mixed Dirichlet-Neumann spectra. This record poses the problem left open by that work: establish the strong hot spots property for explicit classes of bounded piecewise-smooth domains in Rd\mathbb{R}^d, d≥3d\ge3, strictly beyond those two classes — with the flagship instance being bounded convex domains in R3\mathbb{R}^3 and other low dimensions — or refute it for a natural candidate class by an explicit counterexample. Related open questions concern the spectral theory of the framework itself: the characterization of the divergence-free sector A2A_2 in dimensions d≥4d\ge4, the comparison τ1(C)>μ2\tau_1(C)>\mu_2 between the first curl-curl and first nonzero Neumann eigenvalues beyond strictly convex C1,1C^{1,1} and lip domains, and the minimal boundary regularity; progress on these is partial progress on the framework side, not an answer to the main question. A precise nonempty open core remains: convex counterexamples exist only for all sufficiently large dimensions (de Dios Pont, 2024), and the planar convex case is still open.

Background

Consider the heat equation on a bounded domain Ω⊂Rd\Omega\subset\mathbb{R}^d with a perfectly insulated boundary, i.e. the homogeneous Neumann condition ∂νu=0\partial_\nu u=0 on ∂Ω\partial\Omega, where ν\nu is the outer unit normal. The Neumann Laplacian −ΔN-\Delta_N has eigenvalues 0=μ1<μ2≤μ3≤…0=\mu_1<\mu_2\le\mu_3\le\dots; after the constants, the first nontrivial eigenfunction ψ2\psi_2 (with eigenvalue μ2\mu_2) governs the long-time behavior of solutions. The physical intuition that the hottest and coldest spots of an insulated body migrate to its surface led Rauch in the 1970s to the hot spots conjecture: ψ2\psi_2 attains its maximum and minimum on ∂Ω\partial\Omega. The strong form asks that these extrema be attained only on ∂Ω\partial\Omega. The conjecture is false for completely general domains: Burdzy and Werner (1999) constructed planar domains with holes whose second Neumann eigenfunctions have interior extrema, so research targets structured classes. In the plane the conjecture is known for long and thin domains (R. Bañuelos and K. Burdzy; R. Atar and K. Burdzy), lip domains (Atar–Burdzy: domains enclosed between two Lipschitz graphs with Lipschitz constant at most one, proved via couplings of reflected Brownian motion, i.e. the probabilistic representation of the Neumann heat flow), domains with two axes of symmetry (Jerison–Nadirashvili), convex domains with one axis of symmetry (Pascu), and Euclidean triangles (Judge–Mondal); Steinerberger showed that in convex planar domains hot spots lie within one inradius of a 'tip' of the domain. The case of general convex planar domains remains open. In dimensions d≥3d\ge3 positive results are far sparser: besides a classical observation for cylindrical domains, the main ones are the lip-type domains of Yang (via reflected Brownian motion) and the deformation-based results of Chen–Li–Wang for certain symmetric domains. Kennedy and Rohleder introduced a purely variational approach: a self-adjoint vector-valued Laplacian AA acting on vector fields on Ω\Omega, whose quadratic form contains a boundary term involving the shape operator of ∂Ω\partial\Omega (the matrix of normal derivatives of ν\nu, i.e. the curvature of the boundary), and whose spectrum contains the nonzero Neumann and mixed Dirichlet–Neumann eigenvalues; the eigenfunction for the lowest eigenvalue is essentially ∇ψ2\nabla\psi_2, so sign properties of this vector field (coordinate-wise monotonicity of ψ2\psi_2) yield the strong hot spots conclusion. The underlying function space splits into a gradient part and a divergence-free part, and AA correspondingly splits; the divergence-free sector A2A_2 is well understood only in low dimensions: for simply connected planar Ω\Omega the spectrum of A2A_2 coincides with the Dirichlet spectrum, while in R3\mathbb{R}^3, A2A_2 equals the curl-curl operator Cu=curl⁡curl⁡uCu=\operatorname{curl}\operatorname{curl}u (the Maxwell eigenvalue operator with the boundary condition u×ν=0u\times\nu=0 on tangential fields with div⁡u=0\operatorname{div}u=0). Kennedy and Rohleder proved the strong form for bounded lip domains in Rd\mathbb{R}^d (after rotation, the outer normal has exactly two nonzero components of opposite sign, or equals ±ej\pm e_j; together with piecewise C∞C^\infty faces and a uniform exterior ball condition, meaning a ball of fixed radius rolls freely outside Ω\Omega), and a companion result for domains symmetric with respect to every coordinate hyperplane whose intersection with some orthant is a lip domain, for the distinguished eigenfunctions that are antisymmetric in one and symmetric in the remaining hyperplanes. Their method also yields τ1(C)>μ2\tau_1(C)>\mu_2 for lip domains in R3\mathbb{R}^3, where τ1(C)\tau_1(C) is the first curl-curl eigenvalue; Zhang had proved the same comparison for strictly convex C1,1C^{1,1} domains. Very recently de Dios Pont constructed smooth, centrally symmetric convex domains in Rd\mathbb{R}^d, for all sufficiently large dd, whose second Neumann eigenfunctions attain interior extrema, refuting the conjecture for convex sets in high dimensions; the higher-dimensional question therefore survives only for structured classes and for low dimensions such as d=3d=3.

Problem Statement

Let Ω⊂Rd\Omega\subset\mathbb{R}^d, d≥3d\ge3, be a bounded, connected, piecewise smooth Lipschitz domain (boundary a finite union of C∞C^\infty faces, satisfying a uniform exterior ball condition — a ball of some fixed positive radius can be placed outside Ω\Omega touching each boundary point — as in the source's Hypothesis 1.2). Say that Ω\Omega has the strong hot spots property if every eigenfunction ψ2\psi_2 of the Neumann Laplacian for the first nonzero eigenvalue μ2\mu_2 attains its maximum and its minimum only on ∂Ω\partial\Omega. The common objective is to determine how far the strong hot spots conjecture — and the variational vector-operator framework that currently proves it — extends among higher-dimensional domains, beyond the classes settled by Kennedy and Rohleder: multidimensional lip domains (settled in full), and domains symmetric in every coordinate hyperplane whose intersection with some orthant is a lip domain (settled for the eigenfunctions that are antisymmetric in one coordinate hyperplane and symmetric in the remaining ones). Main question. (1) Prove: establish the strong hot spots property for an explicit, natural class of such domains that is not contained in the union of the two settled classes and is genuinely open — the property must not already follow from established results for previously settled families (lip domains, the settled orthant-symmetric class, planar domains with two axes of symmetry, Euclidean triangles, long-and-thin domains), nor from cylindrical products of such families; the flagship open instance is the class of bounded convex domains in R3\mathbb{R}^3 that satisfy the domain regularity above, and establishing the property for this flagship class alone resolves this part; more generally, deciding for each fixed dimension dd below the unknown threshold beyond which de Dios Pont's counterexamples apply whether all bounded convex domains in Rd\mathbb{R}^d that satisfy the domain regularity above have the property is a strictly stronger statement whose resolution also settles this part. (2) Refute: exhibit an explicit domain belonging to a natural candidate class not covered by known counterexamples, together with an eigenfunction for μ2\mu_2 that attains a strict extremum at an interior point. Related open questions on the spectral side of the variational framework are: a characterization of the divergence-free sector A2A_2 of the vector-valued operator AA in dimensions d≥4d\ge4 — in d=3d=3 it is the curl-curl operator CC — and of its spectral relation to the Dirichlet and Neumann Laplacians; an extension of the comparison τ1(C)>μ2\tau_1(C)>\mu_2 from strictly convex C1,1C^{1,1} domains and lip domains in R3\mathbb{R}^3 to broader geometric classes and to higher dimensions, with sharp conditions for strictness; and the minimal regularity of the boundary faces — currently piecewise C∞C^\infty, weakenable to C3C^3 — under which the conclusions remain valid. These are independent questions: progress on them constitutes partial progress on the framework side, and none of them is a prerequisite for, or a substitute for, an answer to the main question — indeed the source's all-dimensional variational proof (including for lip domains in d≥4d\ge4) uses no further properties of A2A_2.

References

  1. James B. Kennedy and Jonathan Rohleder, On the hot spots conjecture in higher dimensions, Transactions of the American Mathematical Society, Series B 13 (2026), no. 4, 108-131, DOI 10.1090/btran/245, https://doi.org/10.1090/btran/245; arXiv:2410.00816, https://arxiv.org/abs/2410.00816
  2. Jaume de Dios Pont, Convex sets can have interior hot spots, preprint (2024), arXiv:2412.06344, https://arxiv.org/abs/2412.06344
  3. Jaume de Dios Pont, Alexander W. Hsu, Mitchell A. Taylor, Sharp bounds on the failure of the hot spots conjecture, preprint (2025), arXiv:2508.16321, https://arxiv.org/abs/2508.16321
  4. Magnus Fries, Magnus Goffeng, Germán Miranda, Reproving Friedlander's inequality with the de Rham complex, Journal of Spectral Theory 15 (2025), no. 4, 1593-1613, DOI 10.4171/JST/582, https://doi.org/10.4171/JST/582; arXiv:2412.03369, https://arxiv.org/abs/2412.03369
  5. Jonathan Rohleder, A new approach to the hot spots conjecture, preprint (2021, revised 2023), arXiv:2106.05224, https://arxiv.org/abs/2106.05224
  6. Jonathan Rohleder, Curl curl versus Dirichlet Laplacian eigenvalues, Bulletin of the London Mathematical Society 57 (2025), no. 9, 2738-2747, DOI 10.1112/blms.70121, https://doi.org/10.1112/blms.70121; arXiv:2404.11301, https://arxiv.org/abs/2404.11301
  7. Zhibing Zhang, Comparison results for eigenvalues of curl curl operator and Stokes operator, Zeitschrift für angewandte Mathematik und Physik 69 (2018), no. 4, article 104, DOI 10.1007/s00033-018-0997-7, https://doi.org/10.1007/s00033-018-0997-7; arXiv:1807.08154, https://arxiv.org/abs/1807.08154
  8. Rodrigo Bañuelos and Krzysztof Burdzy, On the "hot spots" conjecture of J. Rauch, Journal of Functional Analysis 164 (1999), no. 1, 1-33, DOI 10.1006/jfan.1999.3397, https://doi.org/10.1006/jfan.1999.3397
  9. Krzysztof Burdzy and Wendelin Werner, A counterexample to the "Hot Spots" conjecture, Annals of Mathematics (2) 149 (1999), no. 1, 309-317, DOI 10.2307/121027, https://doi.org/10.2307/121027
  10. David Jerison and Nikolai Nadirashvili, The "hot spots" conjecture for domains with two axes of symmetry, Journal of the American Mathematical Society 13 (2000), no. 4, 741-772, DOI 10.1090/s0894-0347-00-00346-5, https://doi.org/10.1090/s0894-0347-00-00346-5
  11. Rami Atar and Krzysztof Burdzy, On Neumann eigenfunctions in lip domains, Journal of the American Mathematical Society 17 (2004), no. 2, 243-265, DOI 10.1090/s0894-0347-04-00453-9, https://doi.org/10.1090/s0894-0347-04-00453-9
  12. Chris Judge and Sugata Mondal, Euclidean triangles have no hot spots, Annals of Mathematics (2) 191 (2020), no. 1, 167-211, DOI 10.4007/annals.2020.191.1.3, https://doi.org/10.4007/annals.2020.191.1.3
  13. Mihai Pascu, Scaling coupling of reflecting Brownian motions and the hot spots problem, Transactions of the American Mathematical Society 354 (2002), no. 11, 4681-4702, DOI 10.1090/s0002-9947-02-03020-9, https://doi.org/10.1090/s0002-9947-02-03020-9
  14. Stefan Steinerberger, Hot spots in convex domains are in the tips (up to an inradius), Communications in Partial Differential Equations 45 (2020), no. 6, 641-654, DOI 10.1080/03605302.2020.1750427, https://doi.org/10.1080/03605302.2020.1750427
  15. Hongbin Chen, Yi Li and Lihe Wang, Monotone properties of the eigenfunction of Neumann problems, Journal de Mathématiques Pures et Appliquées 130 (2019), 112-129, DOI 10.1016/j.matpur.2019.01.013, https://doi.org/10.1016/j.matpur.2019.01.013
  16. Peng-fei Yang, The hot spots conjecture on a class of domains in R^n with n >= 3, Acta Mathematica Sinica, English Series 27 (2011), no. 4, 639-646, DOI 10.1007/s10255-011-0112-8, https://doi.org/10.1007/s10255-011-0112-8
  17. Bobo Hua and Jin Sun, The hot spots conjecture on Gaussian spaces, preprint (2024, revised 2026), arXiv:2412.20663, https://arxiv.org/abs/2412.20663
  18. Lawford Hatcher, Hot spots on cones and warped product manifolds, preprint (2025, revised 2026), arXiv:2508.18054, https://arxiv.org/abs/2508.18054
  19. Lawford Hatcher, Hot spots in domains of constant curvature, preprint (2025), arXiv:2508.13353, https://arxiv.org/abs/2508.13353
  20. Lawford Hatcher, Hot spots in convex hyperbolic planar domains with small eigenvalues, preprint (2026), arXiv:2605.21621, https://arxiv.org/abs/2605.21621
  21. Haiyun Deng, Xuyong Jiang, Xiaoping Yang, Critical points of second Neumann eigenfunctions on some convex domains in two-dimensional space forms, preprint (2026), arXiv:2607.17882, https://arxiv.org/abs/2607.17882
  22. Jonathan Rohleder, A note on hot-spots free subregions of convex domains, preprint (2025, revised 2026), arXiv:2506.22184, https://arxiv.org/abs/2506.22184