Flexibility of Steklov eigenvalues via boundary homogenisation

Recently, D. Bucur and M. Nahon used boundary homogenisation to show the remarkable ﬂexibility of Steklov eigenvalues of planar domains. In the present paper we extend their result to higher dimensions and to arbitrary manifolds with boundary, even though in those cases the boundary does not generally exhibit any periodic structure. Our arguments use framework of variational eigenvalues and provides a diﬀerent proof of the original results. Furthermore, we present an application of this ﬂexibility to the optimisation of Steklov eigenvalues under perimeter constraint. It is proved that the best upper bound for normalised Steklov eigenvalues of surfaces of genus zero and any ﬁxed number of boundary components can always be saturated by planar domains. This is the case even though any actual maximiser (except for simply connected surfaces) is always far from being planar themselves. In particular, it yields sharp upper bound for the ﬁrst Steklov eigenvalue of doubly connected planar domains.


Introduction, main results and setting . Optimisation of Steklov eigenvalues
Let (M , g ) be a complete smooth Riemannian manifold, Ω ⊂ M be a domain with non-empty Lipschitz boundary and 0 ≡ β : ∂Ω → [0, ∞) be a non-negative function.We refer to such an Ω as a manifold with Lipschitz boundary; any abstract manifold with smooth boundary can be realised in this way.Consider the eigenvalue problem in Ω, ∂ ν u = σβu on ∂Ω.
( . ) Under some integrability conditions on β to be made explicit later (see Theorem .), the eigenvalues are discrete and form a sequence For every k, the naturally normalised eigenvalue is see [GKL , KM ] for a discussion around the naturality of that normalisation.The case β ≡ 1 is of particular interest and is referred to as the Steklov problem.The corresponding Steklov eigenvalues σ k (Ω, g , 1) are denoted simply as σ k (Ω, g ).For many known results and open questions about the Steklov problem, the reader can refer to the survey [GP ] and the references therein.In the present paper we are mainly concerned with the optimisation problem for normalised Steklov eigenvalues.
The first result of this type was obtained by Weinstock [Wei ] who proved that the round disk maximises the first normalised Steklov eigenvalue in the class of all bounded simply connected smooth planar domains.The optimisation problem for other topologies of domains in R 2 remains unsolved.At the same time, if one does not impose any assumptions on the topology of the planar domain, then the optimal upper bound for all normalised Steklov eigenvalues is for all k ∈ N σ k (Ω, g ) 8πk, was found in [GKL ].The main goal of the present paper is to apply the ideas of [BN ] to the optimisation problem for planar domains of fixed topology.Among other things, this allows us to determine the optimal upper bound for the first normalized Steklov eigenvalue in the class of planar domains with exactly 2 boundary components.
As a starting point, let us note that an examination of Weinstock's proof yields that the round disk continues to be the maximiser in the much larger class of all simply connected Riemannian surfaces.The main observation of the present paper is that the same holds for other topologies as well, namely, the optimal upper bound for normalized Steklov eigenvalues for planar domains of fixed topology does not increase after including arbitrary Riemannian surfaces of the same topological type.To give a precise statement, for any γ 0 and b 1, we let Ω γ,b be the compact connected surface with boundary of genus γ with b boundary components, and we define Theorem . .For every b 1 and k 0 one has where the supremum is taken over the set of all bounded Lipschitz domains in R 2 with b boundary components.
Remark . .It follows from our proof that in fact, Σ k (γ, b) is saturated by domains in a surface of constant curvature for every γ 0.
The quantities Σ k (γ, b) have received a lot of attention following the influential work of Fraser and Schoen [FS ], who established the connection between Σ k (γ, b) and free boundary minimal immersions of Ω g ,b into a Euclidean ball.In particular, they showed that for the annulus Ω 0,2 = A, Σ 1 (0, 2) is achieved by a metric g cc on the so-called critical catenoid in B 3 .Combining this result with Theorem .one obtains the following.
Corollary . .Let Ω ⊂ R 2 be a smooth bounded domain with 2 boundary components.Then one has The inequality is sharp, i.e. there exists a sequence of domains Ω n such that σ 1 (Ω n ) → σ 1 (A, g cc ).
Remark . .Theorem .and Corollary .can be extended to Lipschitz rather than smooth domains.In such a case, however, inequality ( . ) would stop being strict.In order to rule out the equality case one would need to show a regularity theorem for σ 1 -maximisers in the spirit of [KS , Theorem . ].
Many sequences of planar domains saturate bound ( .).For any bounded Ω ⊂ R 2 conformal to (A, g cc ) one can find a maximizing sequence Ω n such that Ω n → Ω in Hausdorff distance.Here is a concrete example of one of those maximising sequences, which follows from the proof of Theorem .and the geometry of g cc .Let t 1 be the unique solution of cosh t = t .Set Ω 0 = {z ∈ R 2 , r < |z| < R}, where log R r = 2t 1 .Then define Ω n ⊂ R 2 to be the (topological) annulus whose outer boundary is the same as Ω 0 , but whose inner boundary oscillates uniformly with period 2π/n, where the amplitude of the oscillations is chosen so that the length of the inner boundary component coincides with the length of the outer boundary.As n → ∞, the amplitude in this construction is of order O n −1 and the domains Ω n converge in the Hausdorff metric to Ω 0 while σ 1 (Ω n ) → σ 1 (A, g cc ) as n → ∞.

. Flexibility of the Steklov spectrum
Theorem .can be proved by using as a main tool the material already contained in [BN ].Despite that, we take this opportunity to give an alternative proof using the framework of measure eigenvalues developed in [GKL ].This allows and to extend the results of [BN ] to higher dimension and in a geometric context.
We first make the observation that every compact connected smooth manifold with boundary can be realised as a bounded smooth domain in a complete Riemannian manifold (M , g ).Through this equivalence, we define manifolds with Lipschitz boundary as bounded Lipschitz domains in a complete Riemannian manifold.The weighted Steklov problem ( . ) can be defined for those manifolds as well, the normal derivative being only well-defined almost everywhere.
We prove the following flexibility result for Steklov eigenvalues, which was first observed in [BN ] for planar domains.

Theorem . . Let Ω be a compact connected Riemannian manifold with Lipschitz boundary and let
Then, there exists a family of domains Ω ε ⊂ Ω with Lipschitz boundary such that .As ε → 0, ∂Ω ε → ∂Ω in the Hausdorff distance.
As with [BN ], the proof is based on homogenisation of the boundary.However, when d 3 the boundary may no longer carry a periodic structure which means that classical homogenisation constructions do not work in that setting.Instead, we adapt the geometric homogenisation ideas from [GL ], which do not require any periodic structure.Furthermore, we interpret the statement ( ) of Theorem . in the formalism of variational eigenvalues, which in turn allows us to apply the general convergence results presented in [GKL ].In particular, this approach results in a more streamlined proof compared to [BN ].Let us note that boundary homogenisation of the Steklov problem in dimension d 3 was studied by Ferrero-Lamberti in [FL ], however as with most boundary homogenisation results it required domains in Euclidean space to be of product type; we make no such geometric assumptions.
Theorem . is a consequence of Theorem ., Koebe uniformization theorem and conformal invariance of Steklov eigenvalues for d = 2.For d 3, Steklov eigenvalues are no longer conformally invariant, but one still has the following corollary of Theorem . .Corollary . .Let (M , g ) be a closed Riemannian manifold.Then for any k 0 one has where Ω varies over all smooth domains Ω ⊂ M .

Flexibility of the Steklov spectrum
Informally, this corollary states that the introduction of density does not change the optimal upper bound for the normalized Steklov eigenvalues.At the same time, the problem with density is more natural from the geometric viewpoint [KM ].

. Plan of the paper
In Section , we prove Theorem .and its Corollary .using conformal changes of variable and assuming Theorem . .Then, in Section we prove Theorem . .This is done by first assuming that Ω and β are smooth, using a geometric homogenisation procedure on the boundary.Then, we relax the smoothness assumption and in turn approximate eigenvalues for singular densities, then domains with Lipschitz boundary, in the end extracting a diagonal subsequence from these procedures.

. Notation
We make extensive use throughout the paper of Landau's asymptotic notation.We write The limit in that last bullet point will be either as a parameter tends to 0 or ∞ and will be clear from context.The use of a subscript, for instance f 1 ≪ Ω f 2 means that the constant C or the quantities involved in the definition of the limit may depend on the subscript.We make use of a generalisation of L p spaces, called Orlicz spaces.Given Φ an increasing, nonegative convex function on [0, ∞), Φ(L)(Ω) is the space In addition to Φ(x) = x p (which corresponds to L p spaces), we also will refer to the case Φ(x) = e x , denoted exp L, Φ(x) = x log(1 + x), denoted L log L which is dual to expL, and Φ(x) = x 2 log(1 + x) −1/2 denoted L 2 (log L) −1/2 .For a reference on Orlicz space, see [BS ].

Conformal changes of the metric
In this section we prove Theorem .and its corollary assuming Theorem . .We start by introducing the notion of variational eigenvalues and look at how they behave under a conformal change of variables.

. Function spaces and variational eigenvalues
We study the weighted Steklov problem .through the formalism developed in [GKL ], see also [Kok , KLP ].For any domain with Lipschitz boundary Ω ⊂ M and any Radon measure µ supported on Ω, we define the Sobolev spaces W 1,p (Ω, µ) as the closure of C ∞ (Ω) under the norm we write W 1,p (Υ) := W 1,p (Υ, dv g ) for the usual Sobolev space.
We say that a measure µ is admissible if the trace operator T µ : W 1,2 (Ω) → L 2 (Ω, µ) is compact.We note that under such conditions W 1,2 (Ω, µ) is isomorphic to W 1,2 (Ω), see [GKL , Theorems . and . ] For an admissible measure µ and f ∈ C ∞ (Ω) we define the Rayleigh quotients From this Rayleigh quotient we define the variational eigenvalues where the infimum is taken over all (k +1)-dimensional subspaces F k+1 ⊂ C ∞ (Ω) that remain (k +1)dimensional in L 2 (Ω, µ).Admissibility of µ ensures that the variational eigenvalues are discrete and form a sequence (see [GKL , Proposition . ]) The main example of variational eigenvalues employed in the present paper is the following.

. Conformal optimisation
We are now ready to prove the optimisation theorems for d = 2 under the assumption of Theorem . .
Proof of Corollary . .The inequality ( . ) and its sharpness follows immediately from Theorem .and [FS , Theorem . ].It remains to show that the equality can not be achieved by a smooth domain Ω.Suppose that it does, then by [FS , Theorem . ] there exists ω ∈ C ∞ (A) such that ω = 0 on ∂A and (Ω, g 0 ) is isometric to (A, e −2ω g cc ), where g cc is a metric on a free boundary minimal annulus in B 3 .Then the formula for Gauss curvature in a conformal metric implies that ω is solution to the following problem ( . ) Let κ and κ cc be the geodesic curvture of Ω and critical catenoid respectively.Recall that the isometry group of the critical catenoid acts transitively on its boundary.Thus, κ cc is constant.Similarly, since the solution to ( . ) is unique, the function ∂ ν ω is also constant on ∂A.Then one has κ = κ cc −∂ ν ω is also constant.The only curves of constant geodesic curvature κ on R 2 are circles of radius κ −1 .Hence ∂Ω consists of two circles of the same radius, which is impossible.

Flexibility of the spectrum
In this section we prove Theorem ., first under the assumptions that ∂Ω is smooth and β > 0 is a smooth density, then under the weaker assumption that ∂Ω is Lipschitz and β is in an appropriate integrability class.The section is organised as follows.First, we describe the boundary homogenisation construction yielding the appropriate domains Ω ε .Then, we briefly recall abstract tools defined in [GKL ] to study eigenvalue continuity results, and we use them in order to obtain continuity of the Steklov eigenvalues of Ω ε to weighted Steklov eigenvalue on Ω.Finally, we extend the results to the rough case.

. Boundary homogenisation
This construction combines elements found in [GL , Section ] (for the geometric distribution of the perturbations) and in [BN ] (for the type of perturbation).A distinction from the construction in [BN ] is that the approximation is done "from the inside", allowing us to perform the construction intrinsically in the geometric setting.In this subsection we assume that Ω has smooth boundary and 0 < β ∈ C ∞ (∂Ω).This assumption will be relaxed later in Section . .Invariance of normalised eigenvalues under scaling of the density allows us to furthermore assume that β 1.
Let h be the induced metric on ∂Ω, and assume that ε > 0 is small enough that h is uniformly almost Euclidean in balls of radius 3ε.In other words assume that in geodesic polar coordinates around any z ∈ ∂Ω, h reads where g S d−2 is the round metric on the d − 2-dimensional sphere, and r is a symmetric 2-tensor such that For every ε > 0, let S ε be a maximal ε-separated subset of ∂Ω and let V ε be the Voronoȋ tesselation associated with S ε , i.e.V ε := V ε z : z ∈ S ε , where and the distance is computed with respect to the metric h.We construct a sequence of domains Ω ε ⊂ Ω in the following way.For every z ∈ S ε and θ ∈ S d−2 , let ρ θ,z be the distance from z to ∂V ε z along the geodesic starting with direction θ.Then, define w ε Then, w ε z is piecewise smooth, vanishes on ∂V ε z and satisfies the estimates For any smooth nonnegative function α : ∂Ω → R, we have that In a neighbourhood of size 2ε α ∞ of the boundary ∂Ω, write Fermi coordinates as x = (y, t ), where t is the distance along the unit speed geodesic normal to the boundary at y. Define and Finally, we define Ω ε as which has boundary We note that the family Ω ε has equi-Lipschitz boundary, with the constant depending only on g , ∂Ω, and α.Furthermore, Finally, for almost every y ∈ V ε z , if x = (y, t ) ∈ Z ε z then the area element of ∂Ω ε at x is given by ( . ) We choose which we can do since we assumed β 1.

. Continuity of eigenvalues -the smooth setting
We start by introducing conditions under which which variational eigenvalues are continuous with respect to the measures used to define them.For n ∈ N, let Ω n ⊂ Ω.Let µ n , µ be Radon measures supported respectively on Ω n , Ω, we introduce the following three conditions: (M ) µ n * − µ as measures on Ω and Vol g (Ω \ Ω n ) → 0; (M ) the measures µ, µ n are admissible for all n; (M ) there is an equibounded family of extension maps The following proposition appears as [GKL , Proposition . ].
Proposition . .Suppose that Ω n ⊂ Ω is a sequence of domains and µ, µ n are Radon measures on respectively Remark . .The space W 1,2,−1/2 (Ω) is the space of all functions in L 2 (log L) −1/2 such that their distributional gradient also belongs in that space.It is a space which is contained W 1,p (Ω) for all 1 p < 2 so that the convergence in the previous theorem can be verified in the dual of any of those spaces.
Condition (M ) follows from the trace inequality and the fact that β 1. Condition (M ) follows from the fact that for all ε, Ω ε are Lipschitz domains whose Lipschitz constant is controlled by C Ω sup x∈Ω β(x), and C Ω depends on Ω through the metric.

Mikhail Karpukhin and Jean Lagacé
We claim that the operator Changing the order of integration and integrating over y completes the proof of the claim.Thus, one has By monotonicity, we have that This completes the proof that ( . ) holds, which was enough for our purposes, and the proof of Theorem .under smoothness assumptions is complete.
. Continuity of eigenvalues -the singular setting

. . Singular densities
We first give a condition on β so that βH d−1 ⌊ ∂Ω is an admissible measure.
Lemma . .Suppose that d 3 (respectively d = 2) and that 0 Proof.The case d = 2 is proven in [KLP , Proposition . ], by factoring T β through the bounded trace W 1,2 (Ω) → expL 2 (∂Ω) and appropriate multiplication operators, so that T β is seen to be a norm limit of compact operators.The case d 3 is dealt with in the same way, using instead the bounded trace W 1,2 (Ω) → L 2(d−1) d−2 (∂Ω) given by Gagliardo's trace theorem [Gag ].

. . Lipschitz boundary
In order to study convergence of eigenvalues of domains with Lipschitz boundary, we need [BGT , Theorem . ].Note that this result is proven in the Euclidean setting, but its proof extends to the Riemannian setting directly, see [GL , Lemma . ] for an adaptation to the Riemannian setting of the only part of the proof which is not completely local.
Proof of Theorem .for manifolds with Lipschitz boundary.We first prove that we can any compact manifold with Lipschitz boundary with a sequence of domains with smooth boundary in such a way that the Steklov eigenvalues are stable.
Let Ω be a compact manifold with Lipschitz boundary.Following [Ver , Theorem A. ], (see [MT , Appendix A] for a discussion of the adaptation to the Riemannian case), there exists a sequence of smooth domains Ω n ⊂ Ω converging to Ω such that the boundaries of Ω, Ω n may be respectively parametrised by a finite number of equi-Lipschitz maps γ j , γ j ,n such that γ j ,n → γ j uniformly.This implies in particular that if T n : BV(Ω n ) → L 1 (∂Ω n ) is the trace operator, their norms remains uniformly bounded since it can be estimated in terms of the Lipschitz constants of ∂Ω n and the volume of Ω [AG ].In particular, it follows from Proposition .that for all k ∈ N, σ k (Ω n , g ) → σ k (Ω, g ).
It also follows from [Ver , Theorem A. ] that there are bi-Lipschitz homeomorphisms Φ n from ∂Ω n to ∂Ω, whose bi-Lipschitz character is preserved uniformly in n.In particular, Φ n induces an isomorphism Φ * n : L p (∂Ω) → L p (∂Ω) for every p ∈ [1, ∞], whose norms are uniformly bounded in n.Therefore, extracting a diagonal subsequence from -first finding a sequence of domains Ω n with smooth boundary converging to Ω; -then approximating the weight β • Φ n ∈ L d−1 (∂Ω n ) by smooth weights β m,n ; -finally finding a sequence of domains with Lipschitz boundary Ω ε m,n so that H d−1 ⌊ ∂Ω ε m,n converges to β m,n H d−1 ⌊ ∂Ω n ; provides us with the required sequence of domains proving our claim.