Thin triangles

The purpose of this page is to establish

Theorem: the hot spots conjecture is true for acute-angled triangles ABC when [math]\displaystyle{ \angle BAC }[/math] is sufficiently small.

Let us write [math]\displaystyle{ \varepsilon := \angle BAC }[/math], which we view as being small. We may normalize A = (0,0) and B = (1,0). Since the other two angles [math]\displaystyle{ \angle ABC, ACB }[/math] are [math]\displaystyle{ \pi/2 - O(\varepsilon) }[/math], we may also normalize [math]\displaystyle{ C = (1 + O(\varepsilon^2), \varepsilon + O(\varepsilon^2)) }[/math].

Eigenvalue bound

For the sector of radius 1 and aperture [math]\displaystyle{ \varepsilon }[/math], we know (as discussed at the hot spots conjecture) that the second eigenvalue is [math]\displaystyle{ j_1^{-2} }[/math], where [math]\displaystyle{ j_1=3.8317\ldots }[/math] is the first solution to [math]\displaystyle{ J'_0(j_1)=0 }[/math] (or equivalently [math]\displaystyle{ J_1(j_1)=0 }[/math]). We claim that for the triangle ABC, the second eigenvalue is [math]\displaystyle{ j_1^{-2}+O(\varepsilon) }[/math].

First, the upper bound. We can take the second eigenfunction [math]\displaystyle{ J_0(\sqrt{\lambda}r) }[/math] for an inscribed sector to ABC of radius [math]\displaystyle{ 1-\varepsilon }[/math] and aperture [math]\displaystyle{ \varepsilon }[/math], with [math]\displaystyle{ \lambda = (1-\varepsilon)^{-2} j_1^{-2} }[/math], and extend it smoothly to a function on ABC obeying the Neumann boundary conditions. It will not quite mean zero, but we can subtract a constant of size [math]\displaystyle{ O(\varepsilon) }[/math] to make it of mean zero. The Rayleigh quotient for this object can be shown to be [math]\displaystyle{ j_1^{-2} + O(\varepsilon) }[/math], providing the upper bound.

The lower bound can be proven by similar methods once we get good enough C^2 type bounds on eigenfunctions on ABC, but for now let us just establish the weaker bound [math]\displaystyle{ \lambda \gg 1 }[/math]. It suffices to establish the Poincare inequality

[math]\displaystyle{ \int_{ABC} |u - c|^2 \ll \int_{ABC} |\nabla u|^2 }[/math]

for all u on ABC (not necessarily obeying the Neumann boundary condition) and some constant c depending on u. But one can rescale ABC to, say, the unit equilateral triangle in which case the bound is classical, and note that the constants for the rescaling are favorable.

Regularity of eigenfunctions

Let u be the second eigenfunction on ABC, normalised so that

[math]\displaystyle{ \int_{ABC} |u|^2 = |ABC|. }[/math]

By repeated integration by parts (as discussed at the hot spots conjecture) we have

[math]\displaystyle{ \int_{ABC} |\nabla^j u|^2 = \lambda^j |ABC| }[/math] (*)

for j=0,1,2,3. In particular, [math]\displaystyle{ \nabla^j u }[/math] is O(1) on the average on ABC for j=0,1,2,3. In the middle third of the triangle, we can use elliptic regularity (after reflecting the triangle across AB and AC many times so that a disk of radius ~1 sits inside the domain) to in fact conclude that all derivatives are O(1) in this region.

By working in polar coordinates around A, we may expand

[math]\displaystyle{ u(r,\theta) = \sum_{k=0}^\infty c_k J_{\pi k/\varepsilon}(\sqrt{\lambda} r) \cos(\pi k / \varepsilon) }[/math]

for some coefficients [math]\displaystyle{ c_k }[/math]. Because all derivatives of u are O(1) in the middle third of the triangle, we see that [math]\displaystyle{ c_k J_{\pi k/\varepsilon}(\sqrt{\lambda} r) }[/math] is rapidly decreasing in k in this middle third, which from the asymptotics of Bessel functions gives excellent regularity on the left third of the triangle; in particular it is not difficult to get uniform bounds on the C^3 norm in this third.

The situation is more delicate on the right third of the triangle. We begin in the region of size [math]\displaystyle{ O(\varepsilon) }[/math] around the edge BC. From (*) we know that [math]\displaystyle{ \nabla^3 u = O(\varepsilon^{-1/2}) }[/math] on average in this region, while from the Poincare inequality (or Sobolev inequality) one can show that [math]\displaystyle{ \nabla^2 u = O(1) }[/math] on average in this region. By using Bessel function expansions around B and C one can then show that [math]\displaystyle{ \nabla^2 u = O(1) }[/math] pointwise here; in the rest of the triangle one can use elliptic regularity to also get [math]\displaystyle{ \nabla^2 u = O(1) }[/math] first on average in regions of diameter [math]\displaystyle{ \varepsilon }[/math], and then pointwise by elliptic regularity (and reflection). Thus u is bounded uniformly in C^2.

This is already enough regularity on u to extend u to a circumscribing sector and show that [math]\displaystyle{ \lambda = j_0^{-2} + O(\varepsilon) }[/math]. Now we can relate u more carefully with a Bessel function. We again work in polar coordinates around the origin A. If we write

[math]\displaystyle{ u_0(r) := \frac{1}{\varepsilon} \int_0^\varepsilon u(r,\theta)\ d\theta }[/math]

for the averaged radial component of u, defined for [math]\displaystyle{ 0 \leq r \leq 1-\varepsilon }[/math] (say), then from integration by parts we see that

[math]\displaystyle{ u''_0(r) + \frac{1}{r} u'_0(r) + \lambda u_0(r) = 0, }[/math]

thus u_0 is a multiple [math]\displaystyle{ c_0 J_0(\sqrt{\lambda} r) }[/math] of the Bessel function [math]\displaystyle{ J_0(\sqrt{\lambda} r) }[/math] (recall that u_0 has to be continuous at the origin). From the uniform C^2 bounds, we have [math]\displaystyle{ u(r,\theta) = u_0(r) + O(\varepsilon) }[/math] for all [math]\displaystyle{ 0 \leq \theta \leq \varepsilon }[/math] and [math]\displaystyle{ 0 \leq r \leq 1-\varepsilon }[/math]. Given the normalisation of u, this implies that

[math]\displaystyle{ c_0^2 \int_0^1 J_0(\sqrt{\lambda} r)^2\ \varepsilon r dr = \frac{1}{2} \varepsilon^2 + O(\varepsilon^3) }[/math]

which implies that [math]\displaystyle{ c_0 = \pm c_* + O(\varepsilon) }[/math] where [math]\displaystyle{ c_* }[/math] is the positive absolute constant

[math]\displaystyle{ c_* := (2 \int_0^1 J_0(r/j_1)^2\ r dr)^{-1/2}. }[/math]

Without loss of generality we may take the positive sign, thus [math]\displaystyle{ c_0 = c_* + O(\varepsilon) }[/math], and so we have the asymptotic

[math]\displaystyle{ u(r,\theta) = c_0 J_0(r/j_1) + O(\varepsilon) }[/math] (**)

throughout the triangle. This is already enough to show that the maximum can only occur within [math]\displaystyle{ O(\varepsilon^{1/2}) }[/math] of A, and the minimum can only occur within [math]\displaystyle{ O(\varepsilon^{1/2}) }[/math] of B or C.

Now let us show that the maximum can only occur at A. If for contradiction there is a nearby point P to A which also attains the maximum, then [math]\displaystyle{ \partial_r u = 0 }[/math] at P, but also [math]\displaystyle{ \partial_r u=0 }[/math] at A. Hence, by Rolle's theorem, one has [math]\displaystyle{ \partial_{rr} u = 0 }[/math] at some intermediate point Q on the interval AP. But from the Bessel expansion we see that [math]\displaystyle{ \partial_{rr} u }[/math] is bounded away from zero at A, and so from the uniform C^3 bounds we obtain a contradiction if [math]\displaystyle{ \varepsilon }[/math] is small enough.

Now we assert that the minimum can only occur on BC. Let n be the unit normal to BC, thus [math]\displaystyle{ \partial_n u = 0 }[/math] on BC. If for contradiction there is a nearby point Q to BC which also attains the minimum, then [math]\displaystyle{ \partial_n u = 0 }[/math] at P, so by Rolle's theorem we have [math]\displaystyle{ \partial_{nn} u = 0 }[/math] for some intermediate point Q on the line segment from P to BC in the direction n. As n differs from the radial direction by [math]\displaystyle{ O(\varepsilon) }[/math], we conclude from the uniform C^2 bounds that [math]\displaystyle{ \partial_{rr} u = O(\varepsilon) }[/math] at Q.

We will insert this into the equation

[math]\displaystyle{ \partial_{rr} u + \frac{1}{r} \partial_r u + \frac{1}{r^2} \partial_{\theta\theta} u + \lambda u = 0 }[/math]

which is the eigenfunction equation in polar coordinates. As Q is near to BC, we see from (**) that u is comparable to 1, so [math]\displaystyle{ \lambda u }[/math] is comparable to 1 also. From the Neumann boundary condition we have [math]\displaystyle{ \partial_r u = O(\varepsilon) }[/math] on BC, and so since Q is within [math]\displaystyle{ O(\varepsilon^{1/2}) }[/math] of BC we conclude from the uniform C^2 bounds that [math]\displaystyle{ \partial u = O(\varepsilon^{1/2}) }[/math] at Q. Putting all this together, we see that

[math]\displaystyle{ |\partial_{\theta\theta} u| \sim 1 }[/math]

at Q.

To be continued...