This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Lower bounds for the first eigenvalue of the Laplacian with zero magnetic field in planar domains

Bruno Colbois and Alessandro Savo
Abstract

We study the Laplacian with zero magnetic field acting on complex functions of a planar domain Ω\Omega, with magnetic Neumann boundary conditions. If Ω\Omega is simply connected then the spectrum reduces to the spectrum of the usual Neumann Laplacian; therefore we focus on multiply connected domains bounded by convex curves and prove lower bounds for its ground state depending on the geometry and the topology of Ω\Omega. Besides the area, the perimeter and the diameter, the geometric invariants which play a crucial role in the estimates are the the fluxes of the potential one-form around the inner holes and the distance between the boundary components of the domain; more precisely, the ratio between its minimal and maximal width. Then, we give a lower bound for doubly connected domains which is sharp in terms of this ratio, and a general lower bound for domains with an arbitrary number of holes. When the inner holes shrink to points, we obtain as a corollary a lower bound for the first eigenvalue of the so-called Aharonov-Bohm operators with an arbitrary number of poles.


Classification AMS 20002000: 58J50, 35P15
Keywords: Magnetic Laplacian, spectrum, lowest eigenvalue, planar domains
Acknowledgments: Research partially supported by INDAM and GNSAGA of Italy

1 Introduction

1.1 Definitions and state of the art

Let Ω\Omega be a bounded, open, connected domain with smooth boundary Ω\partial\Omega in a Riemannian manifold (M,g)(M,g) and let AA be a smooth real one-form on Ω\Omega, (the potential one-form). Define a connection A\nabla^{A} on the space of complex-valued functions C(Ω,𝐂)C^{\infty}(\Omega,{\bf C}) as follows:

XAu=XuiA(X)u,\nabla^{A}_{X}u=\nabla_{X}u-iA(X)u,

for all vector fields XX on Ω\Omega, where \nabla is the Levi-Civita connection of MM. The magnetic Laplacian with potential AA is the operator acting on C(Ω,𝐂)C^{\infty}(\Omega,{\bf C}):

ΔA=(A)A.\Delta_{A}=(\nabla^{A})^{\star}\nabla^{A}.

In 𝐑n{\bf R}^{n} this gives explicitly, in the usual notation:

ΔA=(i+A)2,\Delta_{A}=(i\nabla+A^{\sharp})^{2},

where AA^{\sharp} is the dual vector field of AA, the vector potential. The two-form B=dAB=dA is the magnetic field; dually, in dimension 22, BB is the vector field B=curlAB={\rm curl}A^{\sharp}.

Scope of this paper is to discuss the spectrum of ΔA\Delta_{A} for planar domains. Hence in what follows we take Ω𝐑2\Omega\subset{\bf R}^{2}.

The spectrum of the magnetic Laplacian has been studied extensively for Dirichlet boundary conditions (u=0u=0 on Ω\partial\Omega), and we denote by λ1D(Ω,A)\lambda_{1}^{D}(\Omega,A) the first eigenvalue. First we remark that, thanks to the diamagnetic inequality, one always has:

λ1D(Ω,A)λ1D(Ω,0),\lambda_{1}^{D}(\Omega,A)\geq\lambda_{1}^{D}(\Omega,0),

and in particular λ1D(Ω,A)>0\lambda_{1}^{D}(\Omega,A)>0. For planar domains and constant magnetic field (that is, dA=BdA=B and |B|\lvert{B}\rvert constant), a Faber-Krahn inequality holds, in the sense that the first eigenvalue of a planar domain is minimized by that of the disk of the same area (see [5]). Estimates for sums of eigenvalues can be found in [9].

However in this paper we deal with magnetic Neumann boundary conditions, that is we impose NAu=0,\nabla^{A}_{N}u=0, on the boundary, where NN is the inner unit normal to Ω\partial\Omega. It is known that then ΔA\Delta_{A} admits a discrete spectrum

0λ1N(Ω,A)λ2N(Ω,A)0\leq\lambda_{1}^{N}(\Omega,A)\leq\lambda_{2}^{N}(\Omega,A)\leq\dots

diverging to ++\infty. The first eigenvalue has the following variational characterization:

λ1N(Ω,A)=minuH1(Ω){0}Ω|Au|2𝑑xΩ|u|2𝑑x.\lambda_{1}^{N}(\Omega,A)=\min_{u\in H_{1}(\Omega)\setminus\{0\}}\frac{\int_{\Omega}\lvert{\nabla^{A}u}\rvert^{2}dx}{\int_{\Omega}|u|^{2}dx}. (1)

For computing lower bounds the diamagnetic inequality is of no use; in fact it gives:

λ1N(Ω,A)λ1N(Ω,0)=0,\lambda_{1}^{N}(\Omega,A)\geq\lambda_{1}^{N}(\Omega,0)=0,

because λ1N(Ω,0)\lambda_{1}^{N}(\Omega,0) is simply the first eigenvalue of the usual Neumann Laplacian, which is zero (the associated eigenspace being spanned by the constant functions). There are fewer estimates in this regard; let us first discuss the case of a constant magnetic field |B|=B0>0\lvert{B}\rvert=B_{0}>0 on planar domains. The paper [4] gives a lower bound of λ1N(Ω,A)\lambda_{1}^{N}(\Omega,A) in terms of the inradius of Ω\Omega, λ1N(Ω,0)\lambda_{1}^{N}(\Omega,0) and of course B0B_{0}. Asymptotic expansions as |B|\lvert{B}\rvert\to\infty are obtained in [7]. We also mention the paper [6] which investigates the validity of a reverse Faber-Krahn inequality for constant magnetic field B0B_{0}, that is: is it true that λ1N(Ω,A)\lambda_{1}^{N}(\Omega,A) is always bounded above by that of a disk with equal volume ? It is proved there that this inequality is true when B0B_{0} is either sufficiently small or sufficiently large, but the general case is still open in the simply connected case.

\bullet\quadIn this paper we prove three lower bounds for the first eigenvalue of planar domains under Neumann conditions, when the magnetic field is identically zero. Since this will be the only boundary condition we consider, from now on we will simply write λ1(Ω,A)\lambda_{1}(\Omega,A) instead of λ1N(Ω,A)\lambda_{1}^{N}(\Omega,A).

Let us first clarify the circumstances under which the first eigenvalue might be positive even if the magnetic potential is a closed one-form on Ω\Omega. This is intimately related to a phenomenon in quantum mechanics predicted in 1959 and known as Aharonov-Bohm effect, which has also experimental evidence: a particle travelling a region in the plane might be affected by the magnetic field even if this is identically zero on its path. In fact what the particle ”feels” is not the magnetic field but, rather, the magnetic potential AA, provided that AA is closed but not exact, and that the flux of AA around the pole may assume non-integer values (see below for the precise condition).

Let us be more precise. From the definition we see that, if A=0A=0, the spectrum of ΔA\Delta_{A} coincides with the spectrum of the usual Laplacian under Neumann boundary conditions. The same is true when A=dfA=df is an exact one-form, by the well-known gauge invariance of the magnetic Laplacian. This fundamental property states that the spectrum of ΔA+df\Delta_{A+df} is the same as the spectrum of ΔA\Delta_{A}, for any fC(Ω)f\in C^{\infty}(\Omega), which follows from the identity:

ΔAeif=eifΔA+df\Delta_{A}e^{-if}=e^{-if}\Delta_{A+df}

showing that ΔA\Delta_{A} and ΔA+df\Delta_{A+df} are unitarily equivalent.

On the other hand, if the magnetic field B=dAB=dA is non-zero, then λ1(Ω,A)\lambda_{1}(\Omega,A) is strictly positive. One could then ask if λ1(Ω,A)\lambda_{1}(\Omega,A) has to vanish whenever the magnetic field is zero, that is, whenever AA is a closed one-form.

To that end, let cc be a closed curve in Ω\Omega (a loop). The quantity:

ΦcA=12πcA\Phi^{A}_{c}=\dfrac{1}{2\pi}\oint_{c}A

is called the flux of AA across cc (we assume that cc is travelled once, and we will not specify the orientation of the loop; this will not affect any of the statements, definitions or results which we will prove in this paper).

It turns out that

\bullet\quadλ1(Ω,A)=0\lambda_{1}(\Omega,A)=0 if and only if AA is closed and the cohomology class of AA is an integer, that is, the flux of AA around any loop is an integer.

This was first observed by Shigekawa [12] for closed manifolds, and then proved in [8] for manifolds with boundary. This remarkable feature of the magnetic Laplacian shows its deep relation with the topology of the underlying manifold Ω\Omega. In this paper we will focus precisely on the situation where the potential one form is closed, and we will then give two lower bounds for the first eigenvalue λ1(Ω,A)\lambda_{1}(\Omega,A).

Let us then recall a few previous results when the magnetic field is assumed to vanish. A lower bound for a general Riemannian cylinder (i.e. the surface 𝐒1×(0,L){\bf S}^{1}\times(0,L) endowed with a Riemannian metric) and zero magnetic field has been given in [3], and is somewhat the inspiration of this work: one of two main results here is in fact to improve such bound when Ω\Omega is a doubly connected planar domain.

Directly related to the Aharonov-Bohm effect, we mention the papers [1] and [10] which investigate the behavior of the spectrum of a domain with a pole Ω{a}\Omega\setminus\{a\} when the pole aa approaches the boundary, for Dirichlet boundary conditions. We remark here that the pole is a distinguished point a=(a1,a2)a=(a_{1},a_{2}) and the potential is the harmonic one-form:

Aa=12(x2a2(x1a1)2+(x2a2)2dx1+x1a1(x1a1)2+(x2a2)2dx2)A_{a}=\frac{1}{2}\Big{(}-\dfrac{x_{2}-a_{2}}{(x_{1}-a_{1})^{2}+(x_{2}-a_{2})^{2}}dx_{1}+\dfrac{x_{1}-a_{1}}{(x_{1}-a_{1})^{2}+(x_{2}-a_{2})^{2}}dx_{2}\Big{)}

which has flux 12\frac{1}{2} across any closed curve enclosing aa, giving rise to a magnetic field which is a Dirac distribution concentrated at the pole aa (therefore, the magnetic field indeed vanishes on Ω{a}\Omega\setminus\{a\}). The magnetic Laplacian ΔAa\Delta_{A_{a}} acting on Ω{a}\Omega\setminus\{a\} is often called an Aharonov-Bohm operator. One could think to a domain with a pole as a doubly connected domain for which the inner boundary curve shrinks to a point.

We will in fact give a lower bound for the first eigenvalue of Aharonov-Bohm operators with many poles, and Neumann boundary conditions (see Theorem 3).

The Aharonov-Bohm operators play an interesting role in the study of minimal partitions, see chapter 8 of [2].

For Neumann boundary conditions, we mention the paper [8], where the authors study the multiplicity and the nodal sets corresponding to the ground state for non-simply connected planar domains with harmonic potential. For doubly connected domains, it is shown that λ1(Ω,A)\lambda_{1}(\Omega,A) is maximal precisely when ΦA\Phi^{A} is congruent to 12\frac{1}{2} modulo integers (this fact is no longer true when there are more than two holes). The proof relies on a delicate argument involving the nodal line of a first eigenfunction and the conclusion does not follow from a specific comparison argument, or from an explicit lower bound.

The focus of this paper is on lower bounds for multiply connected planar domains and zero magnetic field defined by the closed potential form AA. By what we have just said, it is clear that estimating the first eigenvalue is a trivial problem if Ω\Omega is simply connected, because then any closed one-form is automatically exact, and therefore λ1(Ω,A)=0\lambda_{1}(\Omega,A)=0 by gauge invariance. Therefore, we restrict our study to domains with nn holes, with n1n\geq 1.

In this paper we will prove: an improved lower bound for doubly connected domains; a general lower bound for multiply connected domains with an arbitrary number of convex holes; a lower bound for a general convex domain with an arbitrary number of punctures. Let us describe these results in detail.

1.2 A lower bound for doubly connected domains

Let us start from doubly connected domains (n=1n=1) hence domains of type:

Ω=FG¯,\Omega=F\setminus\bar{G},

where FF and GG are open and smooth. We assume FF and GG convex. Let ΦA\Phi^{A} be the flux of the closed potential AA around the inner boundary curve G\partial G: by Shigekawa result, the lower bound is simply zero when ΦA\Phi^{A} is an integer. Then, to hope for a positive lower bound, we need to measure how much ΦA\Phi^{A} is far from being an integer, and the natural invariant will then be:

d(ΦA,𝐙)=min{|ΦAk|:k𝐙}.d(\Phi^{A},{\bf Z})=\min\{\lvert{\Phi^{A}-k}\rvert:k\in{\bf Z}\}.

The second important ingredient for our lower bounds is the ratio βB\frac{\beta}{B} between the minimal width and the maximal width of Ω\Omega. To be more precise, let us say that the line segment σΩ\sigma\subset\Omega is an orthogonal ray if it hits the inner boundary G\partial G orthogonally. By definition, the minimal width β\beta (resp. maximal width BB) of Ω\Omega is the minimal (resp. maximal ) length of an orthogonal ray contained in Ω\Omega:

Refer to caption
Figure 1: The minimal width β\beta and the maximal width BB of an annulus Ω=FG¯\Omega=F\setminus\bar{G}

Note that the ratio βB\frac{\beta}{B} is invariant under homotheties, and reaches its largest value 11 whenever the boundary components are parallel curves.

In Theorem 2 of [3] we prove the lower bound:

λ1(Ω,A)4π2|F|2β(Ω)2B(Ω)2d(ΦA,𝐙)2.\lambda_{1}(\Omega,A)\geq\dfrac{4\pi^{2}}{\lvert{\partial F}\rvert^{2}}\dfrac{\beta(\Omega)^{2}}{B(\Omega)^{2}}d(\Phi^{A},{\bf Z})^{2}. (2)

We insist on the fact that if βB\frac{\beta}{B} is bounded below away from zero we get a positive, uniform lower bound even if β\beta tends to zero. Think for example to a concentric annulus Ω\Omega of radii 11 and 1+β1+\beta; then βB=1\frac{\beta}{B}=1 and as β0\beta\to 0 the lower bound will approach 4π2d(ΦA,𝐙)24\pi^{2}d(\Phi^{A},{\bf Z})^{2}, a positive number, which is just the first eigenvalue of the unit circle.

This means that (for fixed perimeter) in order to get λ1\lambda_{1} small, the ratio βB\frac{\beta}{B} (and not just β\beta) has to be small.

Sharpness in terms of βB\frac{\beta}{B}. In [3] we showed that if βB\frac{\beta}{B} is small then the first eigenvalue could indeed be small. We then looked for an example which could show that the dependance on β2B2\frac{\beta^{2}}{B^{2}} is sharp, and we could not find it. Rather, in Examples 14 and 15 in [3], we constructed examples of domains such that BB is bounded below, say by 11, |F|\lvert{\partial F}\rvert is bounded above, β\beta goes to zero and λ1(Ω,A)\lambda_{1}(\Omega,A) goes to zero proportionally to β\beta, for any non-integral flux. Therefore, if one could replace β2B2\frac{\beta^{2}}{B^{2}} by the linear factor βB\frac{\beta}{B} in (2), one would obtain a sharp inequality (with respect to βB\frac{\beta}{B}). See Figure 2 below for the example which shows sharpness.

This is in fact possible, and the theorem which follows should be regarded as the first main theorem of this paper.

Theorem 1.

Let Ω=FG¯\Omega=F\setminus\bar{G} be an annulus in the plane, with FF and GG convex with piecewise-smooth boundary. Let AA be a closed one-form with flux ΦA\Phi^{A} around the inner hole GG. Then:

a)   One has the lower bound:

λ1(Ω,A)π28|F|2|F|2D(F)4β(Ω)B(Ω)d(ΦA,𝐙)2.\lambda_{1}(\Omega,A)\geq\dfrac{\pi^{2}}{8}\cdot\dfrac{\lvert{F}\rvert^{2}}{\lvert{\partial F}\rvert^{2}D(F)^{4}}\cdot\frac{\beta(\Omega)}{B(\Omega)}\cdot d(\Phi^{A},{\bf Z})^{2}.

where β(Ω)\beta(\Omega) and B(Ω)B(\Omega) are, respectively, the minimum and maximum width of Ω\Omega, and D(F)D(F) is the diameter of FF.

b)   If the outer boundary F\partial F is smooth, and if β(Ω)\beta(\Omega) is less than the injectivity radius of the normal exponential map of F\partial F, then we have the simpler lower bound:

λ1(Ω,A)π2|F|2β(Ω)B(Ω)d(ΦA,𝐙)2,\lambda_{1}(\Omega,A)\geq\dfrac{\pi^{2}}{\lvert{\partial F}\rvert^{2}}\frac{\beta(\Omega)}{B(\Omega)}d(\Phi^{A},{\bf Z})^{2}, (3)

Note that, modulo a factor of 44, b) is formally identical to (2) with β/B\beta/B replacing β2/B2\beta^{2}/B^{2}.

We observe that there is no positive constant cc such that

β(Ω)B(Ω)c|F|2D(F)4\frac{\beta(\Omega)}{B(\Omega)}\geq c\dfrac{\lvert{F}\rvert^{2}}{D(F)^{4}}

for all doubly convex annuli in the plane (otherwise, the lower bound would be independent on the inner hole, and this is impossible). This means that Theorem 1 is not a trivial consequence of (2).

In fact, the proof of Theorem 1 uses a suitable partition of Ω\Omega into overlapping annuli for which βB\frac{\beta}{B} is, so to speak, as small as possible (see Section 2 below, and in particular Figure 3 for an example). Recall the δ\delta-interior ball condition:

given xFx\in\partial F, there is a ball of radius δ\delta tangent to F\partial F at xx and entirely contained in FF.

Here and for further applications, we say that the injectivity radius of F\partial F is Inj(F){\rm Inj}(\partial F) if FF satisfies the δ\delta-interior ball condition for any δInj(F)\delta\leq{\rm Inj}(\partial F). If F\partial F is smooth, its injectivity radius is positive.

Finally we picture below the family of domains Ωϵ\Omega_{\epsilon} which realize sharpness. Ωϵ\Omega_{\epsilon} is the difference between two rectangles with parallel sides, with boundaries being ϵ\epsilon units apart. Hence β(Ωϵ)=ϵ\beta(\Omega_{\epsilon})=\epsilon and B(Ωϵ)B(\Omega_{\epsilon}) is uniformly bounded above by 5\sqrt{5}.

Refer to caption
Figure 2: The domain Ωϵ\Omega_{\epsilon} has minimal width ϵ\epsilon and lowest eigenvalue going to zero proportionally with ϵ\epsilon.

We show in Section 2.8 that

π23605d(ΦA,𝐙)2λ1(Ωϵ,A)ϵ110,\frac{\pi^{2}}{360\sqrt{5}}d(\Phi^{A},{\bf Z})^{2}\leq\dfrac{\lambda_{1}(\Omega_{\epsilon},A)}{\epsilon}\leq\frac{1}{10},

so that λ1(Ω,A)\lambda_{1}(\Omega,A) goes to zero proportionally to ϵβB\epsilon\sim\frac{\beta}{B}.

1.3 A general lower bound for multiply connected domains

Now let Ω\Omega be an nn-holed planar domain, which we write as follows:

Ω=F(G¯1G¯n)\Omega=F\setminus(\bar{G}_{1}\cup\dots\cup\bar{G}_{n}) (4)

where the inner holes G1,,GnG_{1},\dots,G_{n} are smooth, open and disjoint. We furthermore assume that F,G1,,GnF,G_{1},\dots,G_{n} are convex. Note that:

Ω=FG1Gn.\partial\Omega=\partial F\cup\partial G_{1}\cup\dots\cup\partial G_{n}.

We will call G1Gn\partial G_{1}\cup\dots\cup\partial G_{n} the inner boundary of Ω\Omega. The minimal and maximal widths of Ω\Omega are defined as in the case n=1n=1, namely β\beta is the minimal length of a line segment contained in Ω\Omega and hitting the inner boundary orthogonally, and the maximal length of such line segments is by definition the maximal width BB.

It is clear that we could replace B(Ω)B(\Omega) by the diameter of FF, and β(Ω)\beta(\Omega) by the invariant:

β~(Ω)=min{d(Gj,Gk),d(Gh,F):jk,h=1,,n}.\tilde{\beta}(\Omega)=\min\{d(\partial G_{j},\partial G_{k}),d(\partial G_{h},\partial F):j\neq k,h=1,\dots,n\}.

In this section we give a lower bound of λ1(Ω,A)\lambda_{1}(\Omega,A) when Ω\Omega has an arbitrary number of convex holes.

Here is the estimate.

Theorem 2.

Let Ω=F(G¯1G¯n)\Omega=F\setminus(\bar{G}_{1}\cup\dots\cup\bar{G}_{n}) be an n-holed planar domain, where F,G1,,GnF,G_{1},\dots,G_{n} are smooth, open and convex. Let AA be a closed potential having flux Φj\Phi_{j} around the jj-th inner boundary curve Gj\partial G_{j}, for j=1,,nj=1,\dots,n, and let γ=minj=1,,nd(Φj,𝐙).\gamma=\min_{j=1,\dots,n}d(\Phi_{j},{\bf Z}). Then we have:

λ1(Ω,A)π22(|F|+2πB(Ω))2β(Ω)4B(Ω)4γ2.\lambda_{1}(\Omega,A)\geq\dfrac{\pi^{2}}{2\Big{(}\lvert{\partial F}\rvert+2\pi B(\Omega)\Big{)}^{2}}\dfrac{\beta(\Omega)^{4}}{B(\Omega)^{4}}\cdot\gamma^{2}. (5)

where β(Ω)\beta(\Omega) and B(Ω)B(\Omega) are, respectively, the minimal and maximal width of Ω\Omega.

The proof uses a suitable decomposition of Ω\Omega into a finite union of annuli, and a lower bound proved in [3] for annuli whose outer boundary is star-shaped with respect to the inner boundary curve. A stronger estimate is proved when the inner holes are disks of the same radius (see Theorem 13).

1.4 A lower bound for Aharonov-Bohm operators with many poles

The power β4B4\frac{\beta^{4}}{B^{4}} in the previous estimate is probably not sharp; it appears to be there for technical reasons. By shrinking the inner boundary curves to points we obtain an estimate in terms of β2B2\frac{\beta^{2}}{B^{2}}, which has an interesting interpretation in terms of Aharonov-Bohm operators with many poles.

Precisely, we fix a convex domain Ω\Omega and choose nn points inside it, say 𝒫={p1,,pn}\mathcal{P}=\{p_{1},\dots,p_{n}\}. Consider the punctured domain Ω𝒫\Omega\setminus\mathcal{P}. Given a closed one-form AA, we define:

λ1(Ω𝒫,A)=lim infδ0λ1(Ω𝒫(δ),A)\lambda_{1}(\Omega\setminus\mathcal{P},A)=\liminf_{\delta\to 0}\lambda_{1}(\Omega\setminus\mathcal{P}(\delta),A)

where 𝒫(δ)\mathcal{P}(\delta) is the δ\delta-neighborhood of 𝒫\mathcal{P} (it obviously consists of a finite set of disks of radius δ\delta). It is not our scope in this paper to investigate the convergence in terms of δ\delta; however, what we are looking at could be interpreted as the first eigenvalue of a Aharonov-Bohm operator with poles p1,,pnp_{1},\dots,p_{n} and Neumann boundary conditions. The proof of the theorem in the previous section simplifies, to give a general lower bound in terms of the distance between the poles, and the distance of each pole to the boundary. To that end, define:

{β(𝒫)=min{d(pj,pk),d(pm,Ω):pjpk,pm𝒫}B(𝒫)=max{d(pj,pk),d(pm,Ω):pjpk,pm𝒫}\left\{\begin{aligned} &\beta(\mathcal{P})=\min\{d(p_{j},p_{k}),d(p_{m},\partial\Omega):p_{j}\neq p_{k},p_{m}\in\mathcal{P}\}\\ &B(\mathcal{P})=\max\{d(p_{j},p_{k}),d(p_{m},\partial\Omega):p_{j}\neq p_{k},p_{m}\in\mathcal{P}\}\end{aligned}\right.

Of course B(𝒫)B({\mathcal{P}}) could be conveniently bounded above by the diameter of Ω\Omega. Let AA be as usual a closed one-form having flux Φj\Phi_{j} around the pole pjp_{j}. Then we have:

Theorem 3.

Let Ω\Omega be a convex domain and 𝒫={p1,,pn}\mathcal{P}=\{p_{1},\dots,p_{n}\} a finite set of poles. For the punctured domain Ω𝒫\Omega\setminus\mathcal{P} we have the bound:

λ1(Ω𝒫,A)4π2|Ω|2β(𝒫)2B(𝒫)2γ2.\lambda_{1}(\Omega\setminus\mathcal{P},A)\geq\dfrac{4\pi^{2}}{\lvert{\partial\Omega}\rvert^{2}}\dfrac{\beta(\mathcal{P})^{2}}{B({\mathcal{P}})^{2}}\gamma^{2}.

where γ=minj=1,,nd(Φj,𝐙)\gamma=\min_{j=1,\dots,n}d(\Phi_{j},{\bf Z}), and Φj\Phi_{j} is the flux of the closed potential AA around pjp_{j}.

\bullet\quadIn a forthcoming paper, we will give upper bounds for the Laplacian with zero magnetic field on multiply connected planar domains, which are closely related to the topology (number of holes) of the domain.

The rest of the paper is devoted to the proof of Theorems 1,2 and 3.

2 Proof of Theorem 1

The proof depends on a suitable way to partition our domain Ω\Omega. We first remark the simple fact that the first eigenvalue of a domain is controlled from below by the smallest first eigenvalue of the subdomains of a partition of Ω\Omega (Proposition 4). Then, we need to extend inequality (2) to piecewise-smooth boundaries, see Section 2.2. In Section 2.3 we state our main geometric facts, Lemma 6 and Lemma 7, and then we prove Theorem 1 (see Section 2.4). Finally, in Section 2.5, we define the partition and we prove Lemma 6 and Lemma 7.

2.1 A simple lemma

We say that the family of open subdomains {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} is a partition of Ω\Omega, if Ω¯=Ω¯1Ω¯n\bar{\Omega}=\bar{\Omega}_{1}\cup\dots\cup\bar{\Omega}_{n}. Thus, the members of the partition might overlap and some of the intersections ΩjΩj\Omega_{j}\cap\Omega_{j} could have positive measure. If furthemore ΩjΩk\Omega_{j}\cap\Omega_{k} is empty for all jkj\neq k then we say that the partition is disjoint. We observe the following standard fact whose proof is easy:

Proposition 4.

Let {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} be a partition of the domain Ω\Omega. Let AA be any closed potential. Then, there is an index k=1,,nk=1,\dots,n such that

λ1(Ω,A)1nλ1(Ωk,A).\lambda_{1}(\Omega,A)\geq\dfrac{1}{n}\lambda_{1}(\Omega_{k},A). (6)

If the partition is disjoint, then:

λ1(Ω,A)minj=1,,nλ1(Ωj,A).\lambda_{1}(\Omega,A)\geq\min_{j=1,\dots,n}\lambda_{1}(\Omega_{j},A). (7)
Proof.

We start proving (6). Let uu be an eigenfunction associated to λ1(Ω,A)\lambda_{1}(\Omega,A). We use it as test-function for λ1(Ωj,A)\lambda_{1}(\Omega_{j},A) and obtain, for all jj:

λ1(Ωj,A)Ωj|u|2Ωj|Au|2.\lambda_{1}(\Omega_{j},A)\int_{\Omega_{j}}\lvert{u}\rvert^{2}\leq\int_{\Omega_{j}}\lvert{\nabla^{A}u}\rvert^{2}. (8)

Now

Ω|u|2j=1nΩj|u|2nΩk|u|2\int_{\Omega}\lvert{u}\rvert^{2}\leq\sum_{j=1}^{n}\int_{\Omega_{j}}\lvert{u}\rvert^{2}\leq n\int_{\Omega_{k}}\lvert{u}\rvert^{2}

where the index kk is chosen so that Ωk|u|2\int_{\Omega_{k}}\lvert{u}\rvert^{2} is maximum among all j=1,,nj=1,\dots,n. Then:

λ1(Ωk,A)Ω|u|2\displaystyle\lambda_{1}(\Omega_{k},A)\int_{\Omega}\lvert{u}\rvert^{2} nλ1(Ωk,A)Ωk|u|2\displaystyle\leq n\lambda_{1}(\Omega_{k},A)\int_{\Omega_{k}}\lvert{u}\rvert^{2}
=nΩk|Au|2\displaystyle=n\int_{\Omega_{k}}\lvert{\nabla^{A}u}\rvert^{2}
nΩ|Au|2\displaystyle\leq n\int_{\Omega}\lvert{\nabla^{A}u}\rvert^{2}
=nλ1(Ω,A)Ω|u|2\displaystyle=n\lambda_{1}(\Omega,A)\int_{\Omega}\lvert{u}\rvert^{2}

That is: λ1(Ωk,A)nλ1(Ω,A),\lambda_{1}(\Omega_{k},A)\leq n\lambda_{1}(\Omega,A), which is the assertion.

For the proof of (7), let λmin=minj=1,,nλ1(Ωj,A)\lambda_{\rm min}=\min_{j=1,\dots,n}\lambda_{1}(\Omega_{j},A). From (8) we have, for all jj:

Ωj|Au|2λ1(Ωj,A)Ωj|u|2λminΩj|u|2\int_{\Omega_{j}}\lvert{\nabla^{A}u}\rvert^{2}\geq\lambda_{1}(\Omega_{j},A)\int_{\Omega_{j}}\lvert{u}\rvert^{2}\geq\lambda_{\rm min}\int_{\Omega_{j}}\lvert{u}\rvert^{2}

We now sum over j=1,,nj=1,\dots,n and obtain Ω|Au|2λminΩ|u|2.\int_{\Omega}\lvert{\nabla^{A}u}\rvert^{2}\geq\lambda_{\rm min}\int_{\Omega}\lvert{u}\rvert^{2}. As uu is a first eigenfunction the left-hand side is precisely λ1(Ω,A)Ω|u|2\lambda_{1}(\Omega,A)\int_{\Omega}\lvert{u}\rvert^{2}, and the inequality follows.

2.2 Convex annuli with piecewise-smooth boundary

From now on Ω\Omega will be an annulus in the plane with boundary components Γint,Γext\Gamma_{\rm int},\Gamma_{\rm ext} which we assume convex and piecewise-smooth. We will write Ω=FG¯\Omega=F\setminus\bar{G} where FF and GG are open, convex, with piecewise smooth boundary. In that case Γint=G\Gamma_{\rm int}=\partial G and Γext=F\Gamma_{\rm ext}=\partial F.

Let pp be a point of G\partial G where G\partial G is not smooth (pp will then be called a vertex). The normal cone of GG at pp is the set

NG(p)={x𝐑2:x,yp0,for allyG}.N_{G}(p)=\{x\in{\bf R}^{2}:\langle{x},{y-p}\rangle\leq 0,\quad\text{for all}\quad y\in G\}.

Then NG(p)N_{G}(p) is the closed exterior wedge bounded by the normal lines to the two smooth curves concurring at pp, its boundary is the broken line depicted in the figure below. Call αp\alpha_{p} its angle at pp.

\bullet\quadWe remark the obvious fact that 0<αp<π0<\alpha_{p}<\pi.

We now define the minimum and maximum width in the piecewise-smooth case. These are defined in (9) and depicted in the Figure 3 below.

Refer to caption
Figure 3: A vertex pp of G\partial G and its normal cone NG(p)N_{G}(p).

For a unit vector vv applied in pp and pointing inside NG(p)N_{G}(p) we let γp,v(t)=p+tv\gamma_{p,v}(t)=p+tv denote the ray exiting pp in the direction vv, and let Q(p,v)Q(p,v) be the intersection of γp,v\gamma_{p,v} with Γext=F\Gamma_{\rm ext}=\partial F. We define:

{β(p)=infvNG(p)d(p,Q(p,v))B(p)=supvNG(p)d(p,Q(p,v))\left\{\begin{aligned} &\beta(p)=\inf_{v\in N_{G}(p)}d(p,Q(p,v))\\ &B(p)=\sup_{v\in N_{G}(p)}d(p,Q(p,v))\end{aligned}\right. (9)

We notice that at a smooth point qq the cone at qq degenerates to the normal segment at qq. Hence at a smooth point qq one has

β(q)=B(q).\beta(q)=B(q).

We now define

{β(Ω)=inf{β(p):pG}B(Ω)=sup{β(p):pG}\left\{\begin{aligned} &\beta(\Omega)=\inf\{\beta(p):p\in\partial G\}\\ &B(\Omega)=\sup\{\beta(p):p\in\partial G\}\end{aligned}\right. (10)

β(Ω)\beta(\Omega) and B(Ω)B(\Omega) will be called the minimum width and, respectively, the maximum width of Ω\Omega. We remark that when the two boundary components are smooth and parallel then β=B\beta=B and the ratio βB\frac{\beta}{B} assumes its largest possible value, which is 11.

As a first step in the proof of Theorem 1, we extend the inequality (2) to the piecewise-smooth case.

Theorem 5.

Let Ω=FG¯\Omega=F\setminus\bar{G} be an annulus in the plane whose boundary components are convex and piecewise smooth. Let β=β(Ω)\beta=\beta(\Omega) and B=B(Ω)B=B(\Omega) be the invariants defined in (10). Then for any closed potential having flux ΦA\Phi^{A} around the inner boundary curve one has the lower bound:

λ1(Ω,A)4π2|F|2β(Ω)2B(Ω)2d(ΦA,𝐙)2,\lambda_{1}(\Omega,A)\geq\dfrac{4\pi^{2}}{\lvert{\partial F}\rvert^{2}}\dfrac{\beta(\Omega)^{2}}{B(\Omega)^{2}}d(\Phi^{A},{\bf Z})^{2},

where |F|\lvert{\partial F}\rvert is the length of the outer boundary.

Proof.

First, Ω\Omega admits an exhaustion by convex annuli with C1C^{1}-boundary, say {Ωϵ:ϵ>0}\{\Omega_{\epsilon}:\epsilon>0\}. By that we mean:

a) Ωϵ=FϵG¯ϵ\Omega_{\epsilon}=F_{\epsilon}\setminus\bar{G}_{\epsilon} where FϵF_{\epsilon} and GϵG_{\epsilon} are convex and have C1C^{1}-smooth boundary;

b) FϵFF_{\epsilon}\subseteq F and GϵGG_{\epsilon}\supseteq G so that ΩϵΩ\Omega_{\epsilon}\subseteq\Omega;

c) Ω=ϵ>0Ωϵ\Omega=\cup_{\epsilon>0}\Omega_{\epsilon} and in particular limϵ0|ΩΩϵ|=0\lim_{\epsilon\to 0}\lvert{\Omega\setminus\Omega_{\epsilon}}\rvert=0.

To construct FϵF_{\epsilon} we round off corners at distance ϵ\epsilon to each of the vertices of F\partial F; to construct GϵG_{\epsilon} we just take the convex domain bounded by the ϵ\epsilon-neighborhood of GG.

Let uu be an eigenfunction associated to λ1(Ω,A)\lambda_{1}(\Omega,A); by restriction we obtain a test-function for Ωϵ\Omega_{\epsilon}, hence by the min-max principle:

Ωϵ|Au|2Ωϵ|u|2λ1(Ωϵ,A).\dfrac{\int_{\Omega_{\epsilon}}\lvert{\nabla^{A}u}\rvert^{2}}{\int_{\Omega_{\epsilon}}\lvert{u}\rvert^{2}}\geq\lambda_{1}(\Omega_{\epsilon},A).

Let (Ω)\mathcal{L}(\Omega) be the functional:

(Ω)=4π2|F|2β(Ω)2B(Ω)2d(ΦA,𝐙)2\mathcal{L}(\Omega)=\dfrac{4\pi^{2}}{\lvert{\partial F}\rvert^{2}}\dfrac{\beta(\Omega)^{2}}{B(\Omega)^{2}}d(\Phi^{A},{\bf Z})^{2}

We can apply (2) and obtain λ1(Ωϵ)(Ωϵ)\lambda_{1}(\Omega_{\epsilon})\geq\mathcal{L}(\Omega_{\epsilon}) because Ωϵ\Omega_{\epsilon} has smooth boundary; then, for all ϵ>0\epsilon>0:

Ωϵ|Au|2Ωϵ|u|2(Ωϵ).\dfrac{\int_{\Omega_{\epsilon}}\lvert{\nabla^{A}u}\rvert^{2}}{\int_{\Omega_{\epsilon}}\lvert{u}\rvert^{2}}\geq\mathcal{L}(\Omega_{\epsilon}).

We now pass to the limit as ϵ0\epsilon\to 0 on both sides; as (Ωϵ)(Ω)\mathcal{L}(\Omega_{\epsilon})\to\mathcal{L}(\Omega) (as we can easily see from the definitions in (10)), we obtain the assertion: λ1(Ω)(Ω).\lambda_{1}(\Omega)\geq\mathcal{L}(\Omega).

2.3 Preparatory results

In this section we state the two main technical lemmas; the partition of the annulus Ω\Omega will be defined in Section 2.5.

So let Ω=FG¯\Omega=F\setminus\bar{G} be an annulus as above and let β(0,β(Ω)]\beta\in(0,\beta(\Omega)]. We consider the distance functions:

ρ1,ρ2:F[0,),\rho_{1},\rho_{2}:F\to[0,\infty),

where ρ1(x)=d(x,G)\rho_{1}(x)=d(x,G) and ρ2(x)=d(x,F)\rho_{2}(x)=d(x,\partial F). Fix a parameter β>0\beta>0. As GG is convex, with piecewise-smooth boundary, it is well-known that the equidistants {ρ1=β}\{\rho_{1}=\beta\} are C1C^{1}-smooth curves. We say that the parameter β\beta is regular if the equidistant {ρ2=β}\{\rho_{2}=\beta\} is a piecewise-smooth curve. Following Appendix 2 in [11], we know that the set of regular parameters has full measure in (0,β(Ω)](0,\beta(\Omega)]; as a consequence, there exists a sequence of regular parameters {βj}β(Ω)\{\beta_{j}\}\to\beta(\Omega) as jj\to\infty.

\bullet\quadBy using an obvious limiting procedure, from now on we take

β=β(Ω)\beta=\beta(\Omega)

and can assume that it is a regular parameter, so that ρ2=β\rho_{2}=\beta is a piecewise-smooth curve.

Lemma 6.

Let Ω=FG¯\Omega=F\setminus\bar{G} be an annulus in the plane with FF and GG convex with piecewise-smooth boundary, and let β=β(Ω)\beta=\beta(\Omega) (which, by assumption, is a regular parameter). Then Ω\Omega admits a partition {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} into (overlapping) subdomains Ωk\Omega_{k} with the following properties.

a)   Ωk\Omega_{k} is an annulus bounded by two convex piecewise-smooth curves, that is, Ω=FkG¯k\Omega=F_{k}\setminus\bar{G}_{k} with FkF_{k} and GkG_{k} convex, and GkG_{k} contains GG (see figure in Section 2.5 below).

b)   The number nn of annuli in the partition can be taken so that:

n2B(Ω)β.n\leq\dfrac{2B(\Omega)}{\beta}.

We estimate the ratio βB\frac{\beta}{B} of each piece as follows.

Lemma 7.

Let {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} be the partition in the previous lemma. For all k=1,,nk=1,\dots,n one has the following facts.

a)   |Fk||F|\lvert{\partial F_{k}}\rvert\leq\lvert{\partial F}\rvert and β(Ωk)=β\beta(\Omega_{k})=\beta.

b)   The following estimate holds:

β(Ωk)B(Ωk)14|F|D(F)2,\dfrac{\beta(\Omega_{k})}{B(\Omega_{k})}\geq\dfrac{1}{4}\dfrac{\lvert{F}\rvert}{D(F)^{2}}, (11)

where D(F)D(F) is the diameter of FF.

c)   If β(Ω)\beta(\Omega) is less than the injectivity radius of F\partial F, then the following simpler lower bound holds for all kk:

β(Ωk)B(Ωk)12.\dfrac{\beta(\Omega_{k})}{B(\Omega_{k})}\geq\dfrac{1}{\sqrt{2}}.

The proof of Lemma 6 and Lemma 7 involves rather simple geometric constructions, but there are some delicate points to take care of, and will be done in Section 2.5. In fact, these two lemmas make it possible to write Ω\Omega as a union of subset Ωk\Omega_{k} such that the ratio β(Ωk)B(Ωk)\frac{\beta(\Omega_{k})}{B(\Omega_{k})} is bounded below, which make the proof of Theorem 1 quite easy, as follows.

2.4 Proof of Theorem 1

We use the partition {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} of Lemma 6. Let AA be a closed potential having flux ΦA\Phi^{A} around the inner boundary curve G\partial G; then, AA has the same flux around the inner component of Ωk\Omega_{k}, by Lemma 6a. By Proposition 4a there exists k{1,,n}k\in\{1,\dots,n\} such that

λ1(Ω,A)1nλ1(Ωk,A).\lambda_{1}(\Omega,A)\geq\dfrac{1}{n}\lambda_{1}(\Omega_{k},A). (12)

By Theorem 5 applied to Ω=Ωk\Omega=\Omega_{k} we see:

λ1(Ωk,A)4π2|Fk|2β(Ωk)2B(Ωk)2d(ΦA,𝐙)2.\lambda_{1}(\Omega_{k},A)\geq\dfrac{4\pi^{2}}{\lvert{\partial F_{k}}\rvert^{2}}\dfrac{\beta(\Omega_{k})^{2}}{B(\Omega_{k})^{2}}d(\Phi^{A},{\bf Z})^{2}.

By b) of Lemma 7 we see:

β(Ωk)2B(Ωk)2116|F|2D(F)4.\dfrac{\beta(\Omega_{k})^{2}}{B(\Omega_{k})^{2}}\geq\dfrac{1}{16}\dfrac{\lvert{F}\rvert^{2}}{D(F)^{4}}. (13)

This, together with the inequality |Fk||F|\lvert{\partial F_{k}}\rvert\leq\lvert{\partial F}\rvert, gives:

λ1(Ωk,A)π24|F|2|F|2D(F)4d(ΦA,𝐙)2.\lambda_{1}(\Omega_{k},A)\geq\dfrac{\pi^{2}}{4}\cdot\dfrac{\lvert{F}\rvert^{2}}{\lvert{\partial F}\rvert^{2}D(F)^{4}}\cdot d(\Phi^{A},{\bf Z})^{2}.

We insert this inequality in (12) and use the inequality 1nβ2B(Ω)\frac{1}{n}\geq\frac{\beta}{2B(\Omega)} (see Lemma 6b) to conclude that

λ1(Ω,A)\displaystyle\lambda_{1}(\Omega,A) β2B(Ω)λ1(Ωk,A)\displaystyle\geq\frac{\beta}{2B(\Omega)}\lambda_{1}(\Omega_{k},A)
π28|F|2|F|2D(F)4βB(Ω)d(ΦA,𝐙)2.\displaystyle\geq\dfrac{\pi^{2}}{8}\cdot\dfrac{\lvert{F}\rvert^{2}}{\lvert{\partial F}\rvert^{2}D(F)^{4}}\cdot\dfrac{\beta}{B(\Omega)}\cdot d(\Phi^{A},{\bf Z})^{2}.

This proves part a) of Theorem 1.

If β(Ω)\beta(\Omega) is less than the injectivity radius of F\partial F we proceed as before, using the lower bound β(Ωk)B(Ωk)12\frac{\beta(\Omega_{k})}{B(\Omega_{k})}\geq\frac{1}{\sqrt{2}} proved in Lemma 7c. We arrive easily at the inequality:

λ1(Ω,A)π2|F|2β(Ω)B(Ω)d(ΦA,𝐙)2,\lambda_{1}(\Omega,A)\geq\dfrac{\pi^{2}}{\lvert{\partial F}\rvert^{2}}\frac{\beta(\Omega)}{B(\Omega)}d(\Phi^{A},{\bf Z})^{2},

which is Theorem 1b).

2.5 The partition of Ω\Omega and the proof of Lemma 6

We start by showing the partition on a particular example, see Figure 4 below. The initial domain is a triangle FF minus a disk GG and β=β(Ω)\beta=\beta(\Omega). We draw the first three pieces Ω1,Ω2,Ω3\Omega_{1},\Omega_{2},\Omega_{3} and then the last one, which is Ω6\Omega_{6} and which coincides with the β\beta-neighborhood of the exterior boundary F\partial F (this is always the case). Note that the pieces overlap, hence the partition is not disjoint.

Refer to caption
Figure 4: The pieces Ω1,Ω2,Ω3\Omega_{1},\Omega_{2},\Omega_{3} and the last piece Ω6={ρ2<β}\Omega_{6}=\{\rho_{2}<\beta\} when the initial domain is the triangle minus the small disk and β=β(Ω)\beta=\beta(\Omega)

.

We now proceed to construct the partition in general. Let then Ω\Omega be convex annulus Ω=FG¯\Omega=F\setminus\bar{G} as above, and consider the distance functions:

ρ1,ρ2:F[0,),\rho_{1},\rho_{2}:F\to[0,\infty),

where ρ1(x)=d(x,G)\rho_{1}(x)=d(x,G) and ρ2(x)=d(x,F)=d(x,Fc)\rho_{2}(x)=d(x,\partial F)=d(x,F^{c}).

\bullet\quadAt step 1, we let F1={ρ1<β}F_{1}=\{\rho_{1}<\beta\}, G1=GG_{1}=G and Ω1=F1G¯1\Omega_{1}=F_{1}\setminus\bar{G}_{1}. That is, Ω1\Omega_{1} is simply the subset of FF at distance less than β\beta to GG.

\bullet\quadAt step 2, we let F2={ρ1<2β}F_{2}=\{\rho_{1}<2\beta\} and G2={ρ1<β}{ρ2>β}G_{2}=\{\rho_{1}<\beta\}\cap\{\rho_{2}>\beta\} and define Ω2=F2G¯2\Omega_{2}=F_{2}\setminus\bar{G}_{2}.

\bullet\quadAt the arbitrary step kk, we let Fk={ρ1<kβ}F_{k}=\{\rho_{1}<k\beta\} and Gk={ρ1<(k1)β}{ρ2>β}G_{k}=\{\rho_{1}<(k-1)\beta\}\cap\{\rho_{2}>\beta\}, and define:

Ωk=FkG¯k.\Omega_{k}=F_{k}\setminus\bar{G}_{k}.

Observe that for any choice of positive numbers a,ba,b the sets {ρ1<a}\{\rho_{1}<a\} and {ρ2<b}\{\rho_{2}<b\} are convex. Therefore, both FkF_{k} and GkG_{k} are convex; moreover GkFkG_{k}\subset F_{k} and then Ωk\Omega_{k} is an annulus. This proves part a) of Lemma 6.

For b) we first prove the following fact:

Fact. Let n[B(Ω)β]n\doteq\Big{[}\dfrac{B(\Omega)}{\beta}\Big{]} (the smallest integer greater than or equal to B(Ω)β\frac{B(\Omega)}{\beta}). Then Fn=FF_{n}=F and Gn={ρ2>β}G_{n}=\{\rho_{2}>\beta\}. In particular, Ωn={ρ2<β}\Omega_{n}=\{\rho_{2}<\beta\} and then, starting from nn, the sequence Ωn\Omega_{n} stabilizes: Ωn=Ωn+1=\Omega_{n}=\Omega_{n+1}=\dots.

For the proof we first observe that, from the definition of B(Ω)B(\Omega), we have F{ρ1<B(Ω)}F\subseteq\{\rho_{1}<B(\Omega)\}; then, if we fix nB(Ω)βn\geq\frac{B(\Omega)}{\beta} we have by definition FFnF\subseteq F_{n} hence F=FnF=F_{n}. To show that Gn={ρ2>β}G_{n}=\{\rho_{2}>\beta\} it is enough to show:

{ρ2>β}{ρ1<(n1)β}.\{\rho_{2}>\beta\}\subseteq\{\rho_{1}<(n-1)\beta\}. (14)

In fact, if not, there would be a point xFx\in F such that d(x,F)=ρ2(x)>βd(x,\partial F)=\rho_{2}(x)>\beta and ρ1(x)(n1)β\rho_{1}(x)\geq(n-1)\beta. Let yGy\in\partial G be a point at minimum distance to xx, and prolong the segment from yy to xx till it hits Γext=F\Gamma_{\rm ext}=\partial F at the point zz. It is clear that then d(y,z)=d(y,x)+d(x,z)>nβd(y,z)=d(y,x)+d(x,z)>n\beta. By definition of B(Ω)B(\Omega) we then have:

B(Ω)d(y,z)>nβ,B(\Omega)\geq d(y,z)>n\beta,

which contradicts the definition of nn. Hence (14) holds.

We now prove part b) of Lemma 6. Observe that Ω¯=k=1nΩ¯k\bar{\Omega}=\cup_{k=1}^{n}\bar{\Omega}_{k} and (by the definition of nn) B(Ω)βn1\frac{B(\Omega)}{\beta}\geq n-1. Since n1n2n-1\geq\frac{n}{2} for all n2n\geq 2 we see that, for all nn:

B(Ω)βn2,\frac{B(\Omega)}{\beta}\geq\frac{n}{2},

which gives the assertion.

2.6 Proof of Lemma 7

We now study the typical piece Ωk=FkG¯k\Omega_{k}=F_{k}\setminus\bar{G}_{k} in the partition. Observe that

Fk=1Fk2Fk,Gk=1Gk2Gk,\partial F_{k}=\partial_{1}F_{k}\cup\partial_{2}F_{k},\quad\partial G_{k}=\partial_{1}G_{k}\cup\partial_{2}G_{k},

where

{1Fk={ρ1=kβ}F¯2Fk={ρ1kβ}F{1Gk={ρ1=(k1)β}{ρ2β}2Gk={ρ1(k1)β}{ρ2=β}\left\{\begin{aligned} &\partial_{1}F_{k}=\{\rho_{1}=k\beta\}\cap\bar{F}\\ &\partial_{2}F_{k}=\{\rho_{1}\leq k\beta\}\cap\partial F\end{aligned}\right.\quad\left\{\begin{aligned} &\partial_{1}G_{k}=\{\rho_{1}=(k-1)\beta\}\cap\{\rho_{2}\geq\beta\}\\ &\partial_{2}G_{k}=\{\rho_{1}\leq(k-1)\beta\}\cap\{\rho_{2}=\beta\}\end{aligned}\right.

(some of these boundary pieces may be empty). As the equidistants {ρ1=r}\{\rho_{1}=r\} are C1C^{1}-smooth, and the equidistant {ρ2=β}\{\rho_{2}=\beta\} is piecewise smooth, we see that Fk\partial F_{k} and Gk\partial G_{k} are both piecewise smooth, hence

\bullet\quadΩk\partial\Omega_{k} is piecewise smooth.

The inner boundary is Gk\partial G_{k}, it is piecewise smooth with vertices in the set

S={ρ1=(k1)β}{ρ2=β}.S=\{\rho_{1}=(k-1)\beta\}\cap\{\rho_{2}=\beta\}.

Now we have to estimate the ratio β(Ωk)/B(Ωk)\beta(\Omega_{k})/B(\Omega_{k}) for a fixed k=1,,nk=1,\dots,n. Recall the function

β(x):Gk𝐑\beta(x):\partial G_{k}\to{\bf R}

defined in (9). First notice that the regular parts

1,regFk={ρ1=kβ}Fand1,regGk={ρ1=(k1)β}{ρ2>β}\partial_{1,{\rm reg}}F_{k}=\{\rho_{1}=k\beta\}\cap F\quad\text{and}\quad\partial_{1,\rm reg}G_{k}=\{\rho_{1}=(k-1)\beta\}\cap\{\rho_{2}>\beta\}

are parallel, at distance β\beta to each other. Hence

β(x)B(x)=1\frac{\beta(x)}{B(x)}=1 (15)

at all points x1,regGkx\in\partial_{1,\rm reg}G_{k}. Similarly, the regular sets

2,regFk={ρ1<kβ}Fand2,regGk={ρ1<(k1)β}{ρ2=β}\partial_{2,{\rm reg}}F_{k}=\{\rho_{1}<k\beta\}\cap\partial F\quad\text{and}\quad\partial_{2,\rm reg}G_{k}=\{\rho_{1}<(k-1)\beta\}\cap\{\rho_{2}=\beta\}

are parallel at distance β\beta and β(x)B(x)=1\frac{\beta(x)}{B(x)}=1 on 2,regGk\partial_{2,\rm reg}G_{k}.

Therefore, it only remains to control the ratio β(x)B(x)\frac{\beta(x)}{B(x)} at the vertices of Gk\partial G_{k}, which are finite, say p1,,pmp_{1},\dots,p_{m}.

Each break point pjp_{j} gives rise to a corresponding wedge W(pj)NG(pj)ΩW(p_{j})\doteq N_{G}(p_{j})\cap\Omega. In Figure 5 we enlarge the domain Ω2\Omega_{2} relative to the partition of Figure 4 and we show its set of wedges. In other words, every annulus Ωk\Omega_{k} is made up of strips of constant width β\beta and wedges, and we need to control β(x)B(x)\frac{\beta(x)}{B(x)} only at the wedges.

Refer to caption
Figure 5: On the left: the piece Ω2\Omega_{2} and its wedges. On the right: the last piece Ω6={ρ2<β}\Omega_{6}=\{\rho_{2}<\beta\}; the ratio βB\frac{\beta}{B} is small at the upper wedge, because the angle of the wedge at its break point is near π\pi.

Typically, the ratio βB\frac{\beta}{B} is small when there are small angles; nevertheless, this ratio is controlled from below by the diameter and the volume of the outer domain, as we will see in the next section.

As FkF_{k} is a convex subset of FF, we see that |Fk||F|\lvert{\partial F_{k}}\rvert\leq\lvert{\partial F}\rvert. Now it is clear from the construction that β(x)β\beta(x)\geq\beta for all xFkx\in\partial F_{k}; moreover, the inequality is attained. Therefore

β(Ωk)=β\beta(\Omega_{k})=\beta

for all kk. This proves part a) of Lemma 7.

2.7 End of proof of Lemma 7

The estimate βB\frac{\beta}{B} on the wedges of the generic piece Ωk\Omega_{k} will be a consequence of Lemma 8 below.

We recall that the cut-locus of F\partial F is the closure of the set of all points which can be joined to F\partial F by at least two minimizing segments; moreover, the injectivity radius of F\partial F is the minimum distance of F\partial F to the cut-locus. If F\partial F is smooth, its injectivity radius is positive. Finally the distance function d(,F)d(\cdot,\partial F) is smooth outside the cut-locus.

Then, we fix a piece Ωk\Omega_{k} and recall that β(Ωk)=β\beta(\Omega_{k})=\beta. For simplicity we write B(Ωk)=BB(\Omega_{k})=B and recall that, by definition, we have BβB\geq\beta.

Let {p1,,pn}\{p_{1},\dots,p_{n}\} be the vertices of Gk\partial G_{k}. For pp in this set, write p=γ1γ2p=\gamma_{1}\cap\gamma_{2}, where γ1,γ2\gamma_{1},\gamma_{2} are the arcs concurring at pp. Note that either γj\gamma_{j} is an equidistant to G\partial G, that is, is a subset of ρ1=(k1)β\rho_{1}=(k-1)\beta, and in that case we say that γj\gamma_{j} is parallel to G\partial G, or γj\gamma_{j} is an equidistant to F\partial F, that is, is a subset of ρ2=β\rho_{2}=\beta, and in that case we say that γj\gamma_{j} is parallel to F\partial F. There are two possibilities:

Type 1. The vertex p=γ1γ2p=\gamma_{1}\cap\gamma_{2} where γ1\gamma_{1} is parallel to G\partial G and γ2\gamma_{2} is parallel to F\partial F;

Type 2. γ1\gamma_{1} and γ2\gamma_{2} are both parallel to F\partial F.

Note that the second type corresponds to the situation where the vertex pp is a point of the cut-locus of F\partial F. The situation where γ1\gamma_{1} and γ2\gamma_{2} are both parallel to G\partial G does not occur, because then pp would belong to the cut locus of G\partial G; however the cut-locus of a convex domain is always contained in the interior of the domain; as pp is outside GG this is impossible.

\bullet\quadFor the partition in the example and its piece Ω2\Omega_{2} (see Figure 5), the vertices p1p_{1} and p3p_{3} are of type 1, while the vertex p2p_{2} is of type 2.

Lemma 8.

a)   If pp is of type 1 then the interior angle of GkG_{k} at pp is larger than or equal to π2\frac{\pi}{2}, hence the angle of the wedge W(p)W(p) at pp is at most π2\frac{\pi}{2}. Consequently,

βB(p)12.\frac{\beta}{B(p)}\geq\frac{1}{\sqrt{2}}.

b)   If pp is of type 2, then pp is in the cut-locus of F\partial F and one has:

βB(p)14|F|D(F)2.\frac{\beta}{B(p)}\geq\dfrac{1}{4}\dfrac{\lvert{F}\rvert}{D(F)^{2}}.

c)   If β=β(Ω)\beta=\beta(\Omega) is less than the injectivity radius of F\partial F then the estimate in a) will hold at all vertices of the decomposition.

Since the lower bound in b) is always weaker than that in a), we have b) at all vertices of Gk\partial G_{k}. It is clear that Lemma 8 completes the proof of Lemma 7.

Proof of Lemma 8a) If pp is of type 1 then pp is not on the cut-locus of F\partial F, hence ρ2\nabla\rho_{2} exists and is a well-defined unit vector in a neighborhood of pp. Note that ρ2(p)\nabla\rho_{2}(p) points in the direction where the distance to F\partial F increases (obviously an analogous observation holds for ρ1\nabla\rho_{1}). Now observe that the angle of the wedge W(p)W(p) is the angle between the vectors ρ1\nabla\rho_{1} and ρ2-\nabla\rho_{2} (see Figure 6).

Refer to caption
Figure 6: Proof of Lemma 8a

Hence, it is enough to show that the quantity

ρ1(p),ρ2(p)=c(p)\langle{\nabla\rho_{1}(p)},{\nabla\rho_{2}(p)}\rangle=c(p)

is non-positive. Assume on the contrary that c(p)>0c(p)>0. We let α(t)\alpha(t) denote the segment which minimizes distance from pp to G\partial G (parametrized by arc-length); hence α(t)=ρ1(α(t))\alpha^{\prime}(t)=-\nabla\rho_{1}(\alpha(t)). We let f(t)f(t) be the function which measures distance from α(t)\alpha(t) to F\partial F, so that:

f(t)=ρ2(α(t)).f(t)=\rho_{2}(\alpha(t)).

Now f(t)=ρ2(α(t)),α(t)=ρ2(α(t)),ρ1(α(t)f^{\prime}(t)=\langle{\nabla\rho_{2}(\alpha(t))},{\alpha^{\prime}(t)}\rangle=-\langle{\nabla\rho_{2}(\alpha(t))},{\nabla\rho_{1}(\alpha(t)}\rangle. In particular,

f(0)=c(p)<0.f^{\prime}(0)=-c(p)<0.

As f(0)=βf(0)=\beta, this means that for small tt one has ρ2(α(t))<β\rho_{2}(\alpha(t))<\beta, but this impossible because α(t)Gk\alpha(t)\in G_{k}, and all points of GkG_{k} are, by definition, at distance at least β\beta to F\partial F.

Hence cp0c_{p}\leq 0 and the angle of the wedge at pp is at most π2\frac{\pi}{2}. Now, the wedge W(p)W(p) is contained in the polygon with vertices p,q,s,rp,q,s,r as in the picture, hence B(p)d(p,s)2βB(p)\leq d(p,s)\leq\sqrt{2}\beta because the angle at pp is at most π2\frac{\pi}{2}, the angles at rr and ss are π2\frac{\pi}{2}, and d(p,r)=d(p,q)=βd(p,r)=d(p,q)=\beta.

Proof of b). Let pp be a vertex of type 2: then, the two arcs concurring at pp are parallel to F\partial F, and pp belongs to the cut locus of F\partial F. The boundary of the wedge W(p)W(p) is made of two distinct segment of the same length β\beta minimizing distance to F\partial F. Then, W(p)W(p) is contained in a wedge of the last member of the partition, that is, {ρ2<β}\{\rho_{2}<\beta\}. As B(p)B(p) depends only on W(p)W(p), we could as well estimate the ratio βB(p)\frac{\beta}{B(p)} by estimating the corresponding ratio for the wedges of {ρ2<β}\{\rho_{2}<\beta\}, which will allow to express βB(p)\frac{\beta}{B(p)} in terms of the geometry of {ρ2<β}\{\rho_{2}<\beta\}, hence, in terms of the geometry of FF.

The relevant picture is shown below (see Figure 7), in which we evidence such an edge WW (dark shadowed in the picture): it has its vertex in the point pp of ρ2=β\rho_{2}=\beta; we let qWq\in W be a point such that d(p,q)=Bd(p,q)=B. We omit to draw the inner boundary as it will play no role in the proof.

Let TT be the triangle with dotted boundary, with a vertex in qq and such that FWTF\setminus W\subseteq T. As ϕ+γ\phi+\gamma is the exterior angle at a vertex of the piecewise-smooth curve ρ2=β\rho_{2}=\beta, we see that ϕ+γπ\phi+\gamma\leq\pi.

Refer to caption
Refer to caption
Figure 7: Estimate of β/B\beta/B on a typical wedge

Each of the angles ψ\psi and α\alpha is less than π2\frac{\pi}{2}. Consider the circle with center pp and radius β\beta. If qq is inside this circle then B<βB<\beta which is impossible. Then qq is outside the circle; α\alpha and ψ\psi are, each, less then the corresponding angles at the vertex qq^{\prime} obtained by projecting qq on the circle. It is clear that each of these two angles is less than π/2\pi/2.

Let ww be the angle at the vertex ss. Then

βsinψ=BsinwB,\dfrac{\beta}{\sin\psi}=\dfrac{B}{\sin w}\geq B, (16)

and similarly βsinαB.\frac{\beta}{\sin\alpha}\geq B. If ψπ4\psi\geq\frac{\pi}{4} or απ4\alpha\geq\frac{\pi}{4} then βB12\frac{\beta}{B}\geq\frac{1}{\sqrt{2}} and we are finished because 12>|F|4D(F)2\frac{1}{\sqrt{2}}>\frac{\lvert{F}\rvert}{4D(F)^{2}}.

\bullet\quadHence we can assume from now on αψπ4\alpha\leq\psi\leq\frac{\pi}{4}.

Lemma 9.

In the above notation we have |T|24|F|.\lvert{T}\rvert\geq\dfrac{\sqrt{2}}{4}\lvert{F}\rvert.

Proof.

We first remark that π4ϕ3π4\frac{\pi}{4}\leq\phi\leq\frac{3\pi}{4}. In fact, we have ψπ4\psi\leq\frac{\pi}{4} and wπ2w\leq\frac{\pi}{2} so that ϕπ4.\phi\geq\frac{\pi}{4}. On the other hand, the same argument applies to γ\gamma, that is, γπ4\gamma\geq\frac{\pi}{4}. Therefore, as ϕ+γπ\phi+\gamma\leq\pi we conclude ϕ3π4.\phi\leq\frac{3\pi}{4}. The same bounds are satisfied by γ\gamma.

As WW is contained in the union of two parallelograms with sides β\beta and BB we see:

|W|2Bβ.\lvert{W}\rvert\leq 2B\beta. (17)

We set X1=FWX_{1}=F\setminus W; we also let P2P_{2} be the convex polygon with vertices p,q,r,sp,q,r,s and P1=TP2P_{1}=T\setminus P_{2} (see Figure 8).

Refer to caption
Figure 8: The polygons P1P_{1} and P2P_{2}.

We have disjoint unions:

{F=X1WT=P1P2\left\{\begin{aligned} &F=X_{1}\cup W\\ &T=P_{1}\cup P_{2}\end{aligned}\right.

We will show that

{|P1||X1|24|X1||P2|24|W|\left\{\begin{aligned} &\lvert{P_{1}}\rvert\geq\lvert{X_{1}}\rvert\geq\dfrac{\sqrt{2}}{4}\lvert{X_{1}}\rvert\\ &\lvert{P_{2}}\rvert\geq\dfrac{\sqrt{2}}{4}\lvert{W}\rvert\end{aligned}\right.

and the assertion will follow by summing up the two inequalities. Now the first inequality is obvious, because X1P1X_{1}\subseteq P_{1}. By the bounds for ϕ\phi and γ\gamma we see that sinϕ\sin\phi and sinγ\sin\gamma are both, at least, 12\frac{1}{\sqrt{2}}. Then:

|P2|\displaystyle\lvert{P_{2}}\rvert =12Bβsinϕ+12Bβsinγ\displaystyle=\dfrac{1}{2}B\beta\sin\phi+\dfrac{1}{2}B\beta\sin\gamma
22Bβ\displaystyle\geq\dfrac{\sqrt{2}}{2}B\beta

Combining the two estimates we see that |P2|24|W|\lvert{P_{2}}\rvert\geq\dfrac{\sqrt{2}}{4}\lvert{W}\rvert as asserted. ∎

End of proof of Lemma 7. Refer to Figure 8. We can assume that αψπ4\alpha\leq\psi\leq\frac{\pi}{4}. We let δ\delta be the length of the segment joining qq and uu (which meets the side opposite to qq orthogonally, by definition), so that:

|T|=12δ2(tanα+tanψ)δ2tanψ.\lvert{T}\rvert=\frac{1}{2}\delta^{2}(\tan\alpha+\tan\psi)\leq\delta^{2}\tan\psi.

The assumptions give tanψ2sinψ\tan\psi\leq\sqrt{2}\sin\psi, so that |T|2δ2sinψ.\lvert{T}\rvert\leq\sqrt{2}\delta^{2}\sin\psi. Using the lower bound for |T|24|F|\lvert{T}\rvert\geq\frac{\sqrt{2}}{4}\lvert{F}\rvert proved before, we have sinψ|F|4δ2\sin\psi\geq\dfrac{\lvert{F}\rvert}{4\delta^{2}} and then, from (16):

βBsinψ|F|4δ2|F|4D(F)2,\frac{\beta}{B}\geq\sin\psi\geq\dfrac{\lvert{F}\rvert}{4\delta^{2}}\geq\dfrac{\lvert{F}\rvert}{4D(F)^{2}},

the last inequality holding because evidently δD(F)\delta\leq D(F). This proves Lemma 8b and, with it, Lemma 7 is completely proved.

2.8 Example showing sharpness

This example is taken from [3], we repeat it below for the sake of clarity. Its scope is to show that the inequality of Theorem 1 is sharp in βB\frac{\beta}{B}.

We take FF to be the rectangle [4,4]×[0,4][-4,4]\times[0,4], Gϵ=[3,3]×[ϵ,2]G_{\epsilon}=[-3,3]\times[\epsilon,2] and consider the doubly connected domain:

Ωϵ=FG¯ϵ.\Omega_{\epsilon}=F\setminus\bar{G}_{\epsilon}.

We refer to the picture in the Introduction. We let AA be any closed 11-form. As a direct consequence of the gauge invariance of the magnetic Laplacian, it is proved in [3] that, for any planar domain Ω\Omega one has:

λ1(Ω,A)ν1(D),\lambda_{1}(\Omega,A)\leq\nu_{1}(D), (18)

where DD is any closed, simply connected subdomain of Ω¯\bar{\Omega}, and where ν1(D)\nu_{1}(D) denotes the first eigenvalue of the usual Laplacian with Neumann boundary conditions on DΩ\partial D\cap\partial\Omega and with Dirichlet boundary conditions on DΩ\partial D\cap\Omega.

Given our choice of Ωϵ\Omega_{\epsilon}, we remove from it the rectangle (1,1)×(0,ϵ)(-1,1)\times(0,\epsilon) to get the simply connected subdomain called DϵD_{\epsilon}. We estimate ν1(Dϵ)\nu_{1}(D_{\epsilon}) by taking the test-function as follows:

ϕ(x,y)={1on the complement of [2,1]×[0,ϵ][1,2]×[0,ϵ] x1on [1,2]×[0,ϵ]1xon [2,1]×[0,ϵ]\phi(x,y)=\left\{\begin{aligned} &1\quad\text{on the complement of $[-2,-1]\times[0,\epsilon]\cup[1,2]\times[0,\epsilon]$ }\\ &x-1\quad\text{on $[1,2]\times[0,\epsilon]$}\\ &1-x\quad\text{on $[-2,-1]\times[0,\epsilon]$}\end{aligned}\right.

It is readily seen that Dϵ|ϕ|2=2ϵ\int_{D_{\epsilon}}\lvert{\nabla\phi}\rvert^{2}=2\epsilon, while Dϵϕ2C>0\int_{D_{\epsilon}}\phi^{2}\geq C>0 (note that C>20C>20). Therefore:

ν1(Dϵ)ϵ10.\nu_{1}(D_{\epsilon})\leq\dfrac{\epsilon}{10}.

Given (18) we conclude that:

λ1(Ωϵ,A)β(Ωϵ)10=ϵ10\lambda_{1}(\Omega_{\epsilon},A)\leq\dfrac{\beta(\Omega_{\epsilon})}{10}=\dfrac{\epsilon}{10} (19)

Applying our lower bound in Theorem 1 to this case, we have

β=ϵ;B=5;|F|=32;|F|=24;D(F)=45,\beta=\epsilon;\ B=\sqrt{5};\ |F|=32;\ |\partial F|=24;\ D(F)=4\sqrt{5},

and we obtain

λ1(Ωϵ,A)π23605d(ΦA,𝐙)2ϵ.\lambda_{1}(\Omega_{\epsilon},A)\geq\frac{\pi^{2}}{360\sqrt{5}}d(\Phi^{A},{\bf Z})^{2}\,\epsilon. (20)

We now observe that the minimum width β(Ωϵ)=ϵ\beta(\Omega_{\epsilon})=\epsilon, by construction, and that B(Ωϵ)B(\Omega_{\epsilon}) is bounded above by 44. Taking into account (19) and (20) we see that λ1(Ωϵ)\lambda_{1}(\Omega_{\epsilon}) goes to zero proportionally to ϵβ(Ωϵ)B(Ωϵ)\epsilon\sim\frac{\beta(\Omega_{\epsilon})}{B(\Omega_{\epsilon})}.

3 Proof of Theorem 2

Let Ω\Omega be an nn-holed domain, which we write: Ω=F(G¯1G¯n),\Omega=F\setminus(\bar{G}_{1}\cup\dots\cup\bar{G}_{n}), with F,G1,,GnF,G_{1},\dots,G_{n} smooth, open and convex.

From now on we denote Γj=Gj\Gamma_{j}=\partial G_{j}. The idea is to use a suitable partition of Ω\Omega by annuli Ωj\Omega_{j} whose boundary is either a piece of F\partial F or is an equidistant curve from two interior boundary curves; each Ωj\Omega_{j} is an annulus Ωj=FjGj\Omega_{j}=F_{j}\setminus G_{j} with piecewise-smooth exterior boundary Fj\partial F_{j} which is star-shaped with respect to Γj=Gj\Gamma_{j}=\partial G_{j}. We can then apply a theorem in [3] and obtain the uniform bound, valid for all jj:

λ1(Ωj,A)2π29(|F|+2πB(Ω))2β(Ω)4B(Ω)4γ2.\lambda_{1}(\Omega_{j},A)\geq\dfrac{2\pi^{2}}{9\Big{(}\lvert{\partial F}\rvert+2\pi B(\Omega)\Big{)}^{2}}\dfrac{\beta(\Omega)^{4}}{B(\Omega)^{4}}\cdot\gamma^{2}. (21)

As the bound holds for all subdomains of a disjoint partition it holds a fortiori for Ω\Omega, thanks to Proposition 4.

3.1 The partition of Ω\Omega

We start by giving in Figure 9 below the picture of the partition {Ω1,Ω2,Ω3}\{\Omega_{1},\Omega_{2},\Omega_{3}\} when Ω\Omega has three holes. The inner boundary of each piece Ωj\Omega_{j} is made of equidistant sets from two suitable holes.

Refer to caption
Figure 9: The partition {Ω1,Ω2,Ω3}\{\Omega_{1},\Omega_{2},\Omega_{3}\} of a domain Ω\Omega with three holes. The curves Γjk\Gamma_{jk} are equidistant sets.

Here is the construction. For each j=1,,nj=1,\dots,n we consider the non-empty open set:

Fj={xF:d(x,Gj)<d(x,Gk)for allkj}.F_{j}=\{x\in F:d(x,G_{j})<d(x,G_{k})\quad\text{for all}\quad k\neq j\}.

If we set:

Hjk={x𝐑2:d(x,Gj)<d(x,Gk)}H_{jk}=\{x\in{\bf R}^{2}:d(x,G_{j})<d(x,G_{k})\} (22)

we see that we can write

Fj=kj(HjkF).F_{j}=\cap_{k\neq j}(H_{jk}\cap F).

It is clear that F¯=j=1nF¯j\bar{F}=\cup_{j=1}^{n}\bar{F}_{j}. We remark that Hjk\partial H_{jk} is the equidistant set from GjG_{j} and GkG_{k}:

Hjk={x𝐑2:d(x,Gj)=d(x,Gk)}.\partial H_{jk}=\{x\in{\bf R}^{2}:d(x,G_{j})=d(x,G_{k})\}. (23)

We have the following general fact.

Lemma 10.

Let G1G_{1} and G2G_{2} be disjoint smooth convex domains. Then the equidistant set H12\partial H_{12} as above is a smooth curve.

Proof.

Let A=𝐑2(G1G2)A={\bf R}^{2}\setminus(G_{1}\cup G_{2}) and let ρj\rho_{j} be the distance function to GjG_{j}, j=1,2j=1,2. The convexity of GjG_{j} implies that ρj\rho_{j} is smooth on the complement of GjG_{j}, so that ρj\rho_{j} is smooth on AA. Let f=ρ1ρ2f=\rho_{1}-\rho_{2}, so that H12\partial H_{12} is the zero set of ff. One has f=ρ1ρ2\nabla f=\nabla\rho_{1}-\nabla\rho_{2} and it is enough to show that f\nabla f has no critical points on {f=0}\{f=0\}. Let pp be a point in H12=H21\partial H_{12}=\partial H_{21} and let γ1\gamma_{1} be the line segment which minimizes the distance from pp to G1G_{1}. One has: γ1(t)=ptρ1(p)\gamma_{1}(t)=p-t\nabla\rho_{1}(p). The corresponding minimizing segment from pp to G2G_{2} is then γ2(t)=ptρ2(p)\gamma_{2}(t)=p-t\nabla\rho_{2}(p). If ρ1(p)=ρ2(p)\nabla\rho_{1}(p)=\nabla\rho_{2}(p) then γ1(t)=γ2(t)\gamma_{1}(t)=\gamma_{2}(t) and, as d(p,G1)=d(p,G2)d(p,G_{1})=d(p,G_{2}), the two minimizing segments would have the same foot qq, which would then belong to both G1G_{1} and G2G_{2}: but this impossible because G1G_{1} and G2G_{2} are disjoint.

Hence, on {f=0}\{f=0\} one has f0\nabla f\neq 0 which proves smoothness. ∎

As G1,,GkG_{1},\dots,G_{k} are disjoint we see that GjFjG_{j}\subset F_{j} and then we can introduce the annulus

ΩjFjG¯j,\Omega_{j}\doteq F_{j}\setminus\bar{G}_{j},

that is:

Ωj={xΩ:d(x,Gj)<d(x,Gk)for allkj}.\Omega_{j}=\{x\in\Omega:d(x,G_{j})<d(x,G_{k})\quad\text{for all}\quad k\neq j\}.

The family {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} gives rise to a disjoint partition of Ω\Omega, as the next lemma shows.

Lemma 11.

The following properties hold:

a)   Ω¯=j=1nΩ¯j\bar{\Omega}=\cup_{j=1}^{n}\bar{\Omega}_{j}.

b)   For jkj\neq k one has that ΩjΩk=\Omega_{j}\cap\Omega_{k}=\emptyset and Ω¯jΩ¯k\bar{\Omega}_{j}\cap\bar{\Omega}_{k} is a smooth curve (eventually empty).

c)   Ωj\Omega_{j} is an annulus with smooth inner boundary GjG_{j} and piecewise smooth outer boundary Fj\partial F_{j}. Moreover:

Fj=(kjΓjk)(FF¯j),\partial F_{j}=\Big{(}\cup_{k\neq j}\Gamma_{jk}\Big{)}\cup(\partial F\cap\bar{F}_{j}),

where Γjk=F¯jF¯k=Ω¯jΩ¯k\Gamma_{jk}=\bar{F}_{j}\cap\bar{F}_{k}=\bar{\Omega}_{j}\cap\bar{\Omega}_{k} is contained in the equidistant curve Hjk\partial H_{jk} from GjG_{j} and GkG_{k}.

Note that actually Ω¯jΩ¯k=F¯jF¯k\bar{\Omega}_{j}\cap\bar{\Omega}_{k}=\bar{F}_{j}\cap\bar{F}_{k}. The proof of the lemma is clear from the definitions.

3.2 Estimate of λ1(Ωj,A)\lambda_{1}(\Omega_{j},A)

As the partition of Lemma 11 is disjoint, from Proposition 4 we have:

λ1(Ω,A)=minj=1,,nλ1(Ωj,A).\lambda_{1}(\Omega,A)=\min_{j=1,\dots,n}\lambda_{1}(\Omega_{j},A).

Therefore, in this section, we estimate the first eigenvalue of the generic member of the partition. To that end, recall a relevant theorem from [3]. Let Ω1=F1G1\Omega_{1}=F_{1}\setminus G_{1} be an annulus with inner boundary curve Γ1=G1\Gamma_{1}=\partial G_{1}, where G1G_{1} is smooth and convex. For xΓ1x\in\Gamma_{1} and t0t\geq 0, consider the segment γ(t)=x+tN(x)\gamma(t)=x+tN(x) where N(x)N(x) is the unit normal to Γ1\Gamma_{1} oriented outside G1G_{1}. Let Q(x)Q(x) be the first intersection of γx(t)\gamma_{x}(t) with the outer boundary curve F1\partial F_{1}, and let θx\theta_{x} be the angle between γx\gamma^{\prime}_{x} and the outer normal ν\nu to F1F_{1} at Q(x)Q(x). We set:

m(Ω1)minxΓ1cosθx.m(\Omega_{1})\doteq\min_{x\in\Gamma_{1}}\cos\theta_{x}.

We recall that Ω1\Omega_{1} is said to be starlike with respect to Γ1\Gamma_{1} if, for any yF1y\in F_{1}, the segment minimizing distance from yy to Γ1\Gamma_{1} is entirely contained in Ω1\Omega_{1}.

We also set:

{β(Ω1)min{d(x,Q(x):xΓ1}B(Ω1)max{d(x,Q(x):xΓ1}\left\{\begin{aligned} &\beta(\Omega_{1})\doteq\min\{d(x,Q(x):x\in\Gamma_{1}\}\\ &B(\Omega_{1})\doteq\max\{d(x,Q(x):x\in\Gamma_{1}\}\end{aligned}\right.

which are called, respectively, the minimum and maximum width of Ω1\Omega_{1}.

The estimate in Theorem 2 of [3] says that:

λ1(Ω1,A)4π2|F1|2β(Ω1)m(Ω1)B(Ω1)d(ΦA,𝐙)2.\lambda_{1}(\Omega_{1},A)\geq\dfrac{4\pi^{2}}{\lvert{\partial F_{1}}\rvert^{2}}\dfrac{\beta(\Omega_{1})m(\Omega_{1})}{B(\Omega_{1})}d(\Phi^{A},{\bf Z})^{2}. (24)

We will apply (24) to each annulus Ωj\Omega_{j} in the above partition of Ω\Omega. We start from:

Lemma 12.

Let Ω=F(G1Gn)\Omega=F\setminus(G_{1}\cup\dots\cup G_{n}) and let Ωj=FjGj\Omega_{j}=F_{j}\setminus G_{j} be a piece in the partition defined above. Then:

a) Ωj\Omega_{j} is an annulus which is starlike with respect to Gj\partial G_{j}, and moreover:

m(Ωj)β(Ωj)2B(Ωj)m(\Omega_{j})\geq\dfrac{\beta(\Omega_{j})}{2B(\Omega_{j})}

b) One has the estimate:

|Fj|2B(Ωj)β(Ωj)(|Gj|+2πB(Ωj))\lvert{\partial F_{j}}\rvert\leq\dfrac{2B(\Omega_{j})}{\beta(\Omega_{j})}\Big{(}\lvert{\partial G_{j}}\rvert+2\pi B(\Omega_{j})\Big{)}

Lemma 12 allows to prove Theorem 2 as follows. We apply (24) to Ωj\Omega_{j} and get:

λ1(Ωj,A)π22(|Gj|+2πB(Ω))2β(Ωj)4B(Ωj)4γ2.\lambda_{1}(\Omega_{j},A)\geq\dfrac{\pi^{2}}{2\Big{(}\lvert{\partial G_{j}}\rvert+2\pi B(\Omega)\Big{)}^{2}}\dfrac{\beta(\Omega_{j})^{4}}{B(\Omega_{j})^{4}}\cdot\gamma^{2}.

To make the lower bound independent on jj, it is enough to observe that β(Ωj)β(Ω),B(Ωj)B(Ω)\beta(\Omega_{j})\geq\beta(\Omega),B(\Omega_{j})\leq B(\Omega) and |Gj||F|\lvert{\partial G_{j}}\rvert\leq\lvert{\partial F}\rvert. Then we get:

λ1(Ωj,A)π22(|F|+2πB(Ω))2β(Ω)4B(Ω)4γ2.\lambda_{1}(\Omega_{j},A)\geq\dfrac{\pi^{2}}{2\Big{(}\lvert{\partial F}\rvert+2\pi B(\Omega)\Big{)}^{2}}\dfrac{\beta(\Omega)^{4}}{B(\Omega)^{4}}\cdot\gamma^{2}.

which is the final step of the proof.

Then, it remains to prove Lemma 12.

3.3 Proof of Lemma 12a

It is enough to prove it for j=1j=1. We first prove that Ω1=F1G¯1\Omega_{1}=F_{1}\setminus\bar{G}_{1} is star shaped with respect to Γ1=G1\Gamma_{1}=\partial G_{1}. Let yF1y\in\partial F_{1} and let σ\sigma be the segment starting at yy and minimizing distance to Γ1\Gamma_{1}: let xΓ1x\in\Gamma_{1} be the foot of σ\sigma.

\bullet\quadNote that, as yF1y\in\partial F_{1}, we must have d(y,Γ1)d(y,Γk)d(y,\Gamma_{1})\leq d(y,\Gamma_{k}) for all k1k\neq 1.

We prove that σ\sigma is entirely contained in Ω1\Omega_{1}. Assume by contradiction that there is zσz\in\sigma such that zΩ1z\notin\Omega_{1}. Then zΩhz\in\Omega_{h} for some h1h\neq 1, and there exists qΓhq\in\Gamma_{h} with d(z,q)<d(z,x)d(z,q)<d(z,x). But then:

d(y,x)\displaystyle d(y,x) =d(y,z)+d(z,x)\displaystyle=d(y,z)+d(z,x)
>d(y,z)+d(z,q)\displaystyle>d(y,z)+d(z,q)
d(y,q)\displaystyle\geq d(y,q)

that is, d(y,x)>d(y,q)d(y,x)>d(y,q) and this means that d(y,Γ1)>d(y,Γh)d(y,\Gamma_{1})>d(y,\Gamma_{h}), which contradicts the assumption. Hence Ω1\Omega_{1} is star shaped.

We now estimate cosθx\cos\theta_{x}, and for convenience we refer to the picture below, Figure 10.

Refer to caption
Figure 10: The curve is the equidistant Γ12\Gamma_{12} from Γ1\Gamma_{1} and Γ2\Gamma_{2}. The tangent to Γ12\Gamma_{12} at p=Q(x)p=Q(x) is the line through pp orthogonal to ν\nu. It cuts the angle between ρ1\nabla\rho_{1} and ρ2\nabla\rho_{2} in half.

Let xΓ1x\in\Gamma_{1} and draw the segment γx(t)=x+tNx\gamma_{x}(t)=x+tN_{x} where NxN_{x} is the unit normal vector to Γ1\Gamma_{1} pointing outside G1G_{1}. It hits F1\partial F_{1} at the point p=Q(x)p=Q(x). If pFp\in\partial F we proceed as in [3] (because F\partial F is convex) and get

cosθxβ(Ω1)B(Ω1).\cos\theta_{x}\geq\dfrac{\beta(\Omega_{1})}{B(\Omega_{1})}.

If pFp\notin\partial F (as in the picture) then pΓ1kp\in\Gamma_{1k} for some k1k\neq 1; we can assume that k=2k=2. Let ww be the point in Γ2\Gamma_{2} such that d(p,x)=d(p,w)=dd(p,x)=d(p,w)=d; observe that cosθx=ρ1,ν\cos\theta_{x}=\langle{\nabla\rho_{1}},{\nu}\rangle where ν\nu is the normal to Γ12\Gamma_{12} at pp pointing away from Γ1\Gamma_{1}. Observe that ν\nu is the unit vector in the direction of ρ1ρ2\nabla\rho_{1}-\nabla\rho_{2}, and that ρ1+ρ2\nabla\rho_{1}+\nabla\rho_{2} is tangent to Γ12\Gamma_{12} at pp. If 2ϕ2\phi is the angle between ρ1\nabla\rho_{1} and ρ2\nabla\rho_{2} then we see that 2ϕ+2θx=π2\phi+2\theta_{x}=\pi, that is

θx=π2ϕ.\theta_{x}=\frac{\pi}{2}-\phi.

Consider the triangle with vertices x,w,px,w,p; it is isosceles on the basis xwxw, (whose length is denoted 2a2a); its height is part of the tangent line to the equidistant at pp. One sees that

dcosθx=dsinϕ=ad\cos\theta_{x}=d\sin\phi=a

hence

cosθx=ad.\cos\theta_{x}=\frac{a}{d}.

Now 2a=d(x,w)β(Ω1)2a=d(x,w)\geq\beta(\Omega_{1}) by definition of β(Ω1)\beta(\Omega_{1}); as the segment joining xx and pp is entirely contained in Ω1\Omega_{1} we see that dB(Ω1)d\leq B(\Omega_{1}). Hence

cosθxβ(Ω1)2B(Ω1)\cos\theta_{x}\geq\frac{\beta(\Omega_{1})}{2B(\Omega_{1})}

as asserted.

3.4 Proof of Lemma 12b

Recall that the typical piece of the decomposition is Ωj=FjGj\Omega_{j}=F_{j}\setminus G_{j}. We need to estimate |Fj|\lvert{\partial F_{j}}\rvert; this is a bit more difficult now because FjF_{j} is no longer convex (there are circumstances under which each FjF_{j} is actually convex - for example, when all holes are disks of the same radius - and we will discuss this case in the next section, to obtain a simpler final estimate).

Set j=1j=1 for concreteness. We apply Green formula to the function ρ1(x)=d(x,G1)\rho_{1}(x)=d(x,\partial G_{1}). Note that Δρ1(x)\Delta\rho_{1}(x) is the curvature at xx of the equidistant to G1\partial G_{1} through xx; as G1\partial G_{1} is convex one has Δρ10\Delta\rho_{1}\leq 0 on the complement of G1G_{1}. By Green formula:

Ω1Δρ1=Ω1ρ1,N,\int_{\Omega_{1}}\Delta\rho_{1}=\int_{\partial\Omega_{1}}\langle{\nabla\rho_{1}},{N}\rangle,

where NN is the inner unit normal. We let D(G1,B)D(G_{1},B) denote the BB-neighborhood of G1G_{1}, so that F1D(G1,B)F_{1}\subseteq D(G_{1},B) by the definition of BB. Since Δρ10\Delta\rho_{1}\leq 0:

Ω1Δρ1D(G1,B)Δρ1.\int_{\Omega_{1}}\Delta\rho_{1}\geq\int_{D(G_{1},B)}\Delta\rho_{1}.

By co-area formula:

D(G1,B)Δρ1=0Bρ1=rΔρ1𝑑r=2πB\int_{D(G_{1},B)}\Delta\rho_{1}=\int_{0}^{B}\int_{\rho_{1}=r}\Delta\rho_{1}\,dr=-2\pi B

because ρ1=rΔρ1=2π\int_{\rho_{1}=r}\Delta\rho_{1}=-2\pi for all rr (we are integrating the opposite of the curvature of a closed curve, and we always obtain 2π-2\pi). Therefore:

Ω1Δρ12πB.\int_{\Omega_{1}}\Delta\rho_{1}\geq-2\pi B. (25)

On the other hand Ω1=G1F1\partial\Omega_{1}=\partial G_{1}\cup\partial F_{1}. Hence:

Ω1ρ1,N=G1ρ1,N+F1ρ1,N.\int_{\partial\Omega_{1}}\langle{\nabla\rho_{1}},{N}\rangle=\int_{\partial G_{1}}\langle{\nabla\rho_{1}},{N}\rangle+\int_{\partial F_{1}}\langle{\nabla\rho_{1}},{N}\rangle.

The first piece is |G1|\lvert{\partial G_{1}}\rvert. On the outer boundary F1\partial F_{1} we see that:

ρ1,N=cosθxβ2B,\langle{\nabla\rho_{1}},{N}\rangle=-\cos\theta_{x}\leq-\dfrac{\beta}{2B},

where θx\theta_{x} is as in the proof of part a), and the inequality then follows from part a). Then:

Ω1ρ1,N|G1|β2B|F1|,\int_{\Omega_{1}}\langle{\nabla\rho_{1}},{N}\rangle\leq\lvert{\partial G_{1}}\rvert-\frac{\beta}{2B}\lvert{\partial F_{1}}\rvert,

and given (25) we obtain 2πB|G1|β2B|F1|-2\pi B\leq\lvert{\partial G_{1}}\rvert-\frac{\beta}{2B}\lvert{\partial F_{1}}\rvert, that is:

|F1|2Bβ(|G1|+2πB),\lvert{\partial F_{1}}\rvert\leq\frac{2B}{\beta}(\lvert{\partial G_{1}}\rvert+2\pi B),

which gives the assertion.

4 Proof of Theorem 3

About the partition {Ω1,,Ωn}\{\Omega_{1},\dots,\Omega_{n}\} of the previous section for domains with nn holes, we remark that if the inner holes G1,,GnG_{1},\dots,G_{n} are disks of the same radius rr, then the equidistant set between any pair of them is simply a straight line, and therefore each Γjk=Ω¯jΩ¯k\Gamma_{jk}=\bar{\Omega}_{j}\cap\bar{\Omega}_{k} is a line segment; moreover the subdomains F1,,FnF_{1},\dots,F_{n} are all convex: see Figure 11 which illustrates the partition when all holes shrink to a point.

Refer to caption
Figure 11: The partition {Ω1,Ω2,Ω3}\{\Omega_{1},\Omega_{2},\Omega_{3}\} for a domain Ω\Omega punctured at the points p1,p2,p3p_{1},p_{2},p_{3}.

We can directly apply (2) to each Ωj\Omega_{j} and obtain:

λ1(Ωj,A)4π2|Fj|2β(Ωj)2B(Ωj)2d(Φj,𝐙)2\lambda_{1}(\Omega_{j},A)\geq\dfrac{4\pi^{2}}{\lvert{\partial F_{j}}\rvert^{2}}\dfrac{\beta(\Omega_{j})^{2}}{B(\Omega_{j})^{2}}\cdot d(\Phi_{j},{\bf Z})^{2}

As FjF_{j} is convex, we have |Fj||F|\lvert{\partial F_{j}}\rvert\subseteq\lvert{\partial F}\rvert and therefore we arrive at the following estimate.

Theorem 13.

Let Ω=F(G1Gn)\Omega=F\setminus(G_{1}\cup\dots\cup G_{n}) with FF convex and G1,,GnG_{1},\dots,G_{n} being disjoint disks of center, respectively, p1,,pnp_{1},\dots,p_{n} and common radius r>0r>0. Then:

λ1(Ω,A)4π2|F|2β(Ω)2B(Ω)2γ2\lambda_{1}(\Omega,A)\geq\dfrac{4\pi^{2}}{\lvert{\partial F}\rvert^{2}}\dfrac{\beta(\Omega)^{2}}{B(\Omega)^{2}}\cdot\gamma^{2}

with γ=minj=1,,nd(Φj,𝐙)\gamma=\min_{j=1,\dots,n}d(\Phi_{j},{\bf Z}).

We remark that if we let r0r\to 0 in Theorem 13 we get the lower bound for the punctured domain Ω{p1,,pn}\Omega\setminus\{p_{1},\dots,p_{n}\} as in Theorem 3.

References

  • [1] Laura Abatangelo, Veronica Felli, Benedetta Noris, and Manon Nys. Sharp boundary behavior of eigenvalues for Aharonov-Bohm operators with varying poles. J. Funct. Anal., 273(7):2428–2487, 2017.
  • [2] Virginie Bonnaillie-Noël and Bernard Helffer. Nodal and spectral minimal partitions—the state of the art in 2016. In Shape optimization and spectral theory, pages 353–397. De Gruyter Open, Warsaw, 2017.
  • [3] Bruno Colbois and Alessandro Savo. Lower bounds for the first eigenvalue of the magnetic Laplacian. J. Funct. Anal., 274(10):2818–2845, 2018.
  • [4] Tomas Ekholm, Hynek Kovarik, and Fabian Portmann. Estimates for the lowest eigenvalue of magnetic Laplacians. J. Math. Anal. Appl., 439(1):330–346, 2016.
  • [5] L. Erdös. Rayleigh-type isoperimetric inequality with a homogeneous magnetic field. Calc. Var. Partial Differential Equations, 4:283–292, 1996.
  • [6] S. Fournais and B. Helffer. Inequalities for the lowest magnetic Neumann eigenvalue. Lett. Math. Phys., 109(7):1683–1700, 2019.
  • [7] Soeren Fournais and Bernard Helffer. Accurate eigenvalue asymptotics for the magnetic Neumann Laplacian. Ann. Inst. Fourier (Grenoble), 56(1):1–67, 2006.
  • [8] B. Helffer, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and M. P. Owen. Nodal sets for groundstates of Schrödinger operators with zero magnetic field in non-simply connected domains. Comm. Math. Phys., 202:629–649, 1999.
  • [9] R. S. Laugesen and B. A. Siudeja. Magnetic spectral bounds on starlike plane domains. ESAIM Control Optim. Calc. Var., 21(3):670–689, 2015.
  • [10] Benedetta Noris, Manon Nys, and Susanna Terracini. On the Aharonov-Bohm operators with varying poles: the boundary behavior of eigenvalues. Comm. Math. Phys., 339(3):1101–1146, 2015.
  • [11] Alessandro Savo. Lower bounds for the nodal length of eigenfunctions of the Laplacian. Ann. Global Anal. Geom., 19(2):133–151, 2001.
  • [12] Ichiro Shigekawa. Eigenvalue problems for the Schrödinger operator with the magnetic field on a compact Riemannian manifold. J. Funct. Anal., 75(1):92–127, 1987.

Bruno Colbois
Université de Neuchâtel, Institut de Mathématiques
Rue Emile Argand 11
CH-2000, Neuchâtel, Suisse

Alessandro Savo
Dipartimento SBAI, Sezione di Matematica, Sapienza Università di Roma
Via Antonio Scarpa 16
00161 Roma, Italy