Steklov-Neumann Eigenproblem Analysis

• Steklov–Neumann eigenproblem is a boundary spectral problem where the spectral parameter appears in a Steklov condition on a designated boundary portion and a Neumann condition on the complement.
• It employs mixed formulations, operator-theoretic approaches, and variational methods to connect classical Neumann eigenvalues with boundary mass concentration limits.
• Applications include eigenvalue optimization, spectral comparisons in annuli, and numerical approximations that illustrate the interplay between boundary flux operators and mixed boundary conditions.

The Steklov–Neumann eigenproblem denotes a family of boundary spectral problems in which the spectral parameter appears in a boundary condition of Steklov type on only part of the boundary, while the complementary part carries a homogeneous Neumann condition. In the weighted formulation,

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma \rho\,u & \text{on }\partial\Omega, \end{cases}$$

the choice $\rho\equiv 1$ on a distinguished boundary portion and $\rho\equiv 0$ on the complement produces the mixed Steklov–Neumann, or sloshing, problem. A second, closely related line of work interprets Steklov eigenvalues as limits of Neumann eigenvalues when interior mass concentrates in a thin boundary layer, so that Steklov spectra arise as critical boundary-concentration endpoints of Neumann families (Girouard et al., 2014, Lamberti et al., 2014).

1. Mixed boundary formulations

The standard mixed formulation is obtained by decomposing the boundary as

$\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$

and imposing

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$

In the survey literature this is identified as a special case of the weighted Steklov problem, and specifically as the sloshing problem when $F$ is the free surface and $B$ the walls of the container (Girouard et al., 2014).

A particularly important geometric realization occurs on doubly connected domains. If

$\partial\Omega=\Gamma_1\cup\Gamma_2,$

with $\Gamma_1$ the inner boundary and $\Gamma_2$ the outer one, the mixed Steklov–Neumann problem takes the form

$\rho\equiv 1$0

Its weak formulation is

$\rho\equiv 1$1

and the spectrum is discrete:

$\rho\equiv 1$2

The first nontrivial eigenvalue is characterized by

$\rho\equiv 1$3

so orthogonality is imposed on the Steklov part of the boundary rather than in the bulk (Basak et al., 26 Mar 2026).

This mixed formulation should be distinguished from the pure Steklov problem on the full boundary, but it is not separate from weighted Steklov theory. In the framework

$\rho\equiv 1$4

mixed Steklov–Neumann conditions are encoded by a weight that vanishes on the Neumann portion. This equivalence is one of the main reasons the Steklov–Neumann problem is often treated inside weighted Steklov spectral geometry rather than as an isolated boundary value problem (Girouard et al., 2014).

2. Operator-theoretic and variational structure

For the classical Steklov problem on a compact Riemannian manifold with boundary $\rho\equiv 1$5,

$\rho\equiv 1$6

the spectrum coincides with that of the Dirichlet-to-Neumann operator

$\rho\equiv 1$7

where $\rho\equiv 1$8 is the harmonic extension of $\rho\equiv 1$9. In the smooth case, $\rho\equiv 0$0 is a first-order elliptic pseudodifferential operator with the same principal symbol as $\rho\equiv 0$1. This viewpoint extends directly to mixed problems by restricting the spectral boundary condition to a distinguished subset of the boundary (Girouard et al., 2014).

The spectrum is discrete whenever the trace operator

$\rho\equiv 0$2

is compact; the survey states that this holds, for instance, if $\rho\equiv 0$3 has Lipschitz boundary. For the pure Steklov problem the min–max principle is

$\rho\equiv 0$4

with orthogonality to constants on the boundary. In the mixed setting on doubly connected domains, the corresponding denominator is supported only on the Steklov boundary $\rho\equiv 0$5 (Girouard et al., 2014, Basak et al., 26 Mar 2026).

A concrete model reduction occurs for the square. By diagonal symmetries, the Steklov problem on the square decomposes into four mixed Steklov problems on a right isosceles triangle; in each reduced problem the Steklov condition is imposed on the hypotenuse, and Dirichlet or Neumann conditions are imposed on the legs according to parity under reflection. One of these reductions is explicitly a sloshing, hence mixed Steklov–Neumann, problem. In the even–even class, the restrictions to the hypotenuse are eigenfunctions of the free beam equation

$\rho\equiv 0$6

with boundary conditions

$\rho\equiv 0$7

showing that mixed Steklov–Neumann problems can sometimes be reduced to one-dimensional self-adjoint spectral problems (Girouard et al., 2014).

3. Neumann-to-Steklov limits

A distinct but mathematically adjacent meaning of the Steklov–Neumann theme arises from boundary mass concentration. Let

$\rho\equiv 0$8

and define the concentrating density

$\rho\equiv 0$9

Then one studies the Neumann problem

$\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$0

For each fixed $\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$1, the corresponding Rayleigh quotient converges to the Steklov quotient, and for bounded $\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$2 domains one has

$\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$3

Thus the Steklov spectrum is the spectral limit of Neumann problems with fixed total mass concentrating near the boundary (Lamberti et al., 2014, Lamberti et al., 2016).

On the unit ball this limit is explicit. For constant boundary density $\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$4, the Steklov eigenvalues are

$\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$5

and the eigenfunctions are homogeneous harmonic polynomials of degree $\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$6. Separation of variables and Bessel-function analysis yield the first-order expansion

$\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$7

hence

$\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$8

Therefore every positive branch is strictly increasing for sufficiently small $\partial\Omega=F\cup B,\qquad F\cap B=\varnothing,$9, and the limiting Steklov eigenvalue is a local minimizer of the corresponding Neumann branch (Lamberti et al., 2016).

This asymptotic picture explains the description of Steklov eigenvalues as “critical Neumann eigenvalues.” They are not merely formal boundary analogues of Neumann spectra; they are realized as singular endpoint values of Neumann problems, and on the ball they sit at the bottom of nearby concentrating-mass branches (Lamberti et al., 2014).

4. Annuli, holes, and symmetry

The mixed Steklov–Neumann problem is especially tractable on concentric annuli. For

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$0

the eigenvalues are

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$1

with eigenfunctions

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$2

where $$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$3 are spherical harmonics of degree $$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$4. The first nonzero eigenvalue corresponds to $$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$5, has multiplicity $$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$6, and the associated eigenfunctions can be written in Cartesian form (Basak et al., 26 Mar 2026).

Among eccentric annuli of the form

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$7

with fixed radii and $$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$8, the first nonzero mixed eigenvalue satisfies

$$\begin{cases} \Delta u=0 & \text{in }\Omega,\ \partial_\nu u=\sigma u & \text{on }F,\ \partial_\nu u=0 & \text{on }B. \end{cases}$$9

with equality only for the concentric configuration. This is a sharp extremal result for the first mixed Steklov–Neumann eigenvalue in that class (Basak et al., 26 Mar 2026).

When the inner spherical hole shrinks,

$F$0

one has

$F$1

where $F$2 is the first nontrivial Steklov eigenvalue of the unperforated outer domain. After harmonic extension into the hole, suitably normalized eigenfunctions converge strongly in $F$3. The corresponding first mixed eigenfunction has exactly two nodal domains (Basak et al., 26 Mar 2026).

Higher mixed eigenvalues on domains with holes also admit annular comparison under symmetry assumptions. If

$F$4

and $F$5 is connected, smooth, centered at the origin, and symmetric of order $F$6, then for the concentric annulus

$F$7

of the same volume,

$F$8

The proof uses the $F$9 annular eigenfunctions together with cancellation of mixed moments enforced by order-$B$0 symmetry. The same work gives counterexamples with only order-$B$1 symmetry, so the symmetry assumption is structurally essential rather than cosmetic (Basak et al., 2024).

5. Perturbation, optimization, and criticality

One optimization problem varies the boundary partition while keeping the geometry fixed. Starting from a decomposition

$B$2

one inserts a small interval $B$3 of length $B$4 and replaces it by Neumann boundary:

$B$5

If $B$6 is an eigenvalue of multiplicity $B$7 and at least one eigenfunction in the eigenspace does not vanish at the center $B$8 of the inserted interval, then

$B$9

The associated Green’s function satisfies

$\partial\Omega=\Gamma_1\cup\Gamma_2,$0

which yields an algorithm for placing and enlarging Neumann pieces so that a prescribed parameter becomes close to a Steklov–Neumann eigenvalue and the Green’s function becomes large (Ammari et al., 2019).

A different optimization theory varies the boundary mass density in the weighted Steklov problem

$\partial\Omega=\Gamma_1\cup\Gamma_2,$1

For a finite cluster $\partial\Omega=\Gamma_1\cup\Gamma_2,$2 of eigenvalues, the elementary symmetric functions $\partial\Omega=\Gamma_1\cup\Gamma_2,$3 are real-analytic on the noncollision set $\partial\Omega=\Gamma_1\cup\Gamma_2,$4, and their Fréchet derivatives are

$\partial\Omega=\Gamma_1\cup\Gamma_2,$5

Under fixed total boundary mass

$\partial\Omega=\Gamma_1\cup\Gamma_2,$6

critical densities are characterized by

$\partial\Omega=\Gamma_1\cup\Gamma_2,$7

On the unit ball, the constant density

$\partial\Omega=\Gamma_1\cup\Gamma_2,$8

is critical for all $\partial\Omega=\Gamma_1\cup\Gamma_2,$9 (Lamberti et al., 2014).

This last fact is one of the main structural differences between Steklov and standard Dirichlet–Neumann spectral optimization. The paper explicitly contrasts the Steklov case with a “maximum principle” for many other elliptic spectral optimization problems, where analogous critical density conditions are typically impossible. In the Steklov setting the condition lives on the boundary, and in highly symmetric domains such as the ball it can be satisfied. This is another sense in which Steklov eigenvalues behave as critical Neumann eigenvalues rather than as ordinary boundary counterparts of Dirichlet or Neumann spectra (Lamberti et al., 2014).

6. Computation and broader generalizations

For the classical Steklov operator on polygonal planar domains, conforming finite elements can produce certified lower bounds. The key abstract estimate is

$\Gamma_1$0

where $\Gamma_1$1 is the conforming finite-element eigenvalue and $\Gamma_1$2 is a computable projection-error constant obtained from an auxiliary nonhomogeneous Neumann problem and a hypercircle construction. In the local analysis of the trace constant, the paper also introduces a mixed Steklov–Neumann problem on a triangle, showing that even certification techniques for pure Steklov spectra naturally generate mixed Steklov–Neumann subproblems (Nakano et al., 2020).

Boundary-only methods are also available. For smooth simply connected planar domains, the Dirichlet-to-Neumann operator can be written in terms of a generalized conjugation operator:

$\Gamma_1$3

After Nyström discretization of the boundary integral equation and Fourier differentiation, one obtains a dense algebraic eigenvalue problem

$\Gamma_1$4

or, in the notation of the paper,

$\Gamma_1$5

This formulation treats interior and exterior Steklov problems in a unified way and reconstructs eigenfunctions by harmonic extension from their boundary traces (Swan et al., 7 Apr 2026).

The Neumann-to-Steklov principle extends far beyond the linear Laplacian. For admissible possibly irregular domains obtained from the unit ball by trace-compatible Sobolev homeomorphisms, the first nontrivial weighted $\Gamma_1$6-Neumann eigenvalue with concentrating bulk weight $\Gamma_1$7 satisfies

$\Gamma_1$8

and normalized minimizers converge, up to subsequences, strongly in $\Gamma_1$9 to weighted Steklov minimizers. Equivalently, the best constants in weighted Poincaré inequalities converge to the best constants in weighted trace inequalities (Menovschikov, 10 May 2026).

Beyond the Laplace setting, related Steklov-type boundary spectra appear for the modified Helmholtz equation, where polygonal corners numerically produce asymptotics

$\Gamma_2$0

and for Maxwell’s equations, where the boundary relation

$\Gamma_2$1

defines a compact NtD-based Steklov-type theory for tangential fields. Discrete analogues also exist on finite subgraphs of $\Gamma_2$2, where the Dirichlet-to-Neumann operator defines a lattice Steklov spectrum (Chaigneau et al., 2023, Lamberti et al., 2020, Han et al., 2019). These developments suggest that the Steklov–Neumann paradigm is best understood not as a single boundary value problem, but as a broad spectral mechanism linking boundary flux operators, mixed boundary partitions, and Neumann-to-boundary concentration limits.