--- title: Jacobi-Steklov Spectrum Overview url: https://www.emergentmind.com/topics/jacobi-steklov-spectrum type: topic --- # Jacobi-Steklov Spectrum Overview The Jacobi–Steklov spectrum is the spectrum of a boundary spectral problem in which the spectral parameter appears in the boundary condition rather than in the interior equation. In its standard form, one considers a second-order elliptic operator \(J\), often a Jacobi operator such as \(J=\Delta+q(x)\), and a Steklov-type boundary condition such as \(\partial_\nu u+\beta(x)u=\sigma u\). The classical Steklov problem for the Laplacian, \[ \Delta u=0 \text{ in } \Omega,\qquad \partial_\nu u=\sigma u \text{ on } \partial\Omega, \] is the canonical simplest example. In smooth settings the associated boundary operator is a first-order elliptic self-adjoint pseudodifferential operator, so the spectrum is discrete, real, and unbounded, and it is naturally interpreted as the spectrum of a generalized Dirichlet-to-Neumann map [1411.6567, 2410.07462]. ## 1. Definition and operator-theoretic formulation For a smooth compact manifold with boundary \(M\), the classical Dirichlet-to-Neumann operator is defined by harmonic extension: if \(f\) is prescribed on \(\partial M\), one solves \[ \Delta(Hf)=0 \text{ in } M,\qquad Hf|_{\partial M}=f, \] and sets \[ \Lambda f=\partial_\nu(Hf). \] The Steklov eigenvalue problem is then \(\Lambda f=\sigma f\), equivalently \[ \Delta u=0 \text{ in } M,\qquad \partial_\nu u=\sigma u \text{ on } \partial M. \] Its spectrum is discrete and nonnegative, \[ 0=\sigma_0<\sigma_1\le \sigma_2\le \dots \nearrow \infty, \] under the regularity assumptions used in the standard Steklov theory [1411.6567]. The Jacobi–Steklov generalization replaces the interior Laplacian by a Jacobi-type operator. In the supplied formulations, one natural model is \[ J=\Delta+q(x), \] with boundary condition \[ \partial_\nu u+\beta(x)u=\sigma u. \] More generally, one may consider \[ J=-\operatorname{div}(A\nabla)+V, \] solve \(Ju=0\) with prescribed boundary values, and define the generalized Dirichlet-to-Neumann map \(\Lambda_J\) by \(f\mapsto B(H_Jf)\), where \(B\) is the boundary operator appearing in the Steklov-type condition. The resulting eigenvalues are the Jacobi–Steklov eigenvalues [1311.5533, 1411.6567]. Several standard generalizations already fit this scheme. A weighted Steklov problem is \[ \Delta f=0 \text{ in } M,\qquad \partial_n f=\sigma\,\delta f \text{ on } \Sigma, \] with boundary density \(\delta\), and a magnetic Steklov problem is \[ \Delta^\eta F=0,\qquad \nu\lrcorner d^\eta F=\sigma F \text{ on } \partial M, \] where \(\Delta^\eta\) is the magnetic Laplacian and \(T^\eta(f)=-\nu\lrcorner d^\eta \hat f\) is the magnetic Steklov operator [1103.2863, 2410.07462]. ## 2. Pseudodifferential structure and high-energy asymptotics For smooth boundary, the Dirichlet-to-Neumann operator is a first-order elliptic self-adjoint pseudodifferential operator whose principal symbol is \(|\xi|\) with respect to the induced boundary metric. This immediately places Steklov and Jacobi–Steklov problems in the standard microlocal framework for elliptic \(\Psi\)DOs and yields Weyl-type asymptotics for the counting function [1311.5533, 1411.6567]. In dimension \(n\), the classical Steklov counting function satisfies \[ \#\{\sigma_j<\sigma\} = \frac{\omega_{n-1}\,|\partial\Omega|}{(2\pi)^{n-1}}\sigma^{n-1}+O(\sigma^{n-2}), \] and equivalently \[ \sigma_j = 2\pi\left(\frac{j}{\omega_{n-1}|\partial\Omega|}\right)^{\frac{1}{n-1}}+O(1). \] On surfaces this simplifies drastically. If the boundary is connected and has length \(L\), then \[ \sigma_{2j}=\sigma_{2j-1}=\frac{2\pi}{L}j+O(j^{-\infty}), \] and for multiple boundary components the asymptotic model is a rearranged union of arithmetic progressions \(\frac{2\pi}{\ell(M_i)}\mathbb{N}\), each appearing twice [1311.5533]. For Jacobi–Steklov operators, the same principal-symbol argument implies that lower-order terms do not alter the order of growth. A standard consequence is that the leading Weyl term is governed by the boundary measure, while lower-order data enter only in subprincipal corrections. In the supplied formulations this is repeatedly expressed by saying that a generalized Dirichlet-to-Neumann operator for \(J\) remains an elliptic first-order operator with the same principal symbol structure; this suggests that the leading Weyl growth matches the classical Steklov law [1411.6567, 1311.5533]. A model beyond the smooth case is provided by spherical cylinders \(\Omega=B_R^m\times(0,L)\). For these domains, the Steklov counting function admits a two-term asymptotic expansion in which the second term contains both a curvature contribution from smooth boundary pieces and an edge contribution from the codimension-two strata where the lateral boundary meets the bases. The edge coefficient is represented by explicit angular and Beta-function integrals, for example \[ G_{n,1} = \frac{\Gamma(\tfrac12)^n}{2^n\Gamma(\tfrac n2)} \int_0^1 (1+x)^{-1/2}(1-x)^{\frac n2-1}\,dx, \] which is a universal right-angle edge constant in the second term [2506.20498]. ## 3. Spectral invariants of boundary geometry A central feature of Steklov-type spectra is that they are strongly boundary-sensitive. For smooth compact surfaces, the Steklov spectrum determines the number of boundary components and the lengths of those components. In particular, if \(\ell_{\max}\) is the length of the largest boundary component, then \[ \ell_{\max} = \frac{2\pi}{\limsup_{j\to\infty}(\sigma_{j+1}-\sigma_j)}. \] This is obtained from the refined asymptotic model \(\sigma_j=S(\rho)_j+O(j^{-\infty})\), where \(\rho_i=2\pi/\ell(M_i)\) and \(S(\rho)\) is the reordered union of the corresponding arithmetic progressions [1311.5533]. For compact Riemannian orbisurfaces, the boundary components come in two types: smooth circles and type II components, which are quotient circles with two reflector singularities. The Steklov spectrum determines the numbers of type I and type II boundary components and therefore detects the presence and number of orbifold singularities on the boundary. It also determines the boundary lengths modulo the explicit equivalence relation \[ L_1\sqcup L_1\sqcup 2\bar L_1 = L_2\sqcup L_2\sqcup 2\bar L_2, \] which is the precise ambiguity in the two-dimensional orbifold setting [1609.05142]. At the same time, the Steklov spectrum does not detect all singular geometric features. A flat disk can be Steklov isospectral to a cone: for the flat cone \(\mathcal C_k\) obtained as a quotient of a disk by a rotation of angle \(2\pi/k\), the Dirichlet-to-Neumann operator coincides with that of the flat disk \(D(2\pi)\). Consequently, the Steklov spectrum does not detect interior orbifold singularities and does not determine the orbifold Euler characteristic [1609.05142]. For convex polygons, the characteristic polynomial \(P_\Omega\) extracted from the Steklov spectrum encodes side lengths and angles through trigonometric frequencies and angle-dependent coefficients. Within this framework, almost all triangles are uniquely determined by their Steklov spectra within the class of all triangles; rectangles are uniquely determined among all triangles and convex quadrilaterals; and triangles and quadrilaterals are spectrally distinguished from simply connected smooth domains. In this sense, corners are spectrally visible in low polygonal complexity [2604.18977]. ## 4. Nonsmooth domains, edges, and singular spectral phenomena Nonsmooth boundary geometry changes the subleading structure of the spectrum. In piecewise smooth domains, the principal Weyl term remains determined by total boundary measure, but the next term collects contributions from smooth curvature and from edges. Spherical cylinders \(B\times I\) are a model example: the lateral boundary contributes curvature, the flat bases contribute differently, and the right-angle intersections produce explicit edge terms weighted by universal constants such as \(G_{n,1}\) [2506.20498]. Geometric singularities can also change the spectral type itself. In a bounded domain with a cusp of sharpness exponent \(m=2\), the continuous spectrum of the Steklov problem is \[ [\lambda_\dagger,\infty), \qquad \lambda_\dagger = \left(n-\frac32\right)^2\frac{|\omega|}{|\partial\omega|}. \] If the cusp is blunted at scale \(\varepsilon\), the spectrum on the truncated domains is discrete, but as \(\varepsilon\to0^+\) three kinds of behavior appear: stable eigenvalues, blinking eigenvalues, and gliding eigenvalues. Blinking eigenvalues reappear near any fixed \(\lambda^\flat>\lambda_\dagger\) along logarithmically spaced sequences \[ \varepsilon_k^\flat = \exp\Big\{-2\big(\tau_0(\lambda^\flat)\big)^{-1}\big((2k+1)\pi+\Theta(\lambda^\flat)\big)\Big\}, \] and this is part of the mechanism by which the discrete spectra of the blunted domains generate the continuous spectrum of the cuspidal limit [1807.00514]. Another singular phenomenon concerns perturbations of rough boundaries. Under Lipschitz graph convergence together with perimeter convergence, the Steklov resolvents converge compactly in the sense of Vainikko, so each eigenvalue converges and eigenfunctions converge strongly in \(H^1\) after identification by the connecting system. If the surface elements converge only weakly, the limit problem becomes a weighted Steklov problem \[ \Delta u=0,\qquad \partial_\nu u=\lambda\,y(x)\,u \text{ on } \partial\Omega, \] and if the boundary oscillates too strongly, all positive eigenvalues collapse to \(0\). In periodically corrugated domains \(g_\varepsilon(x')=\varepsilon^\alpha b(x'/\varepsilon)\), one obtains a trichotomy: stability for \(\alpha>1\), a weighted limit for \(\alpha=1\), and degeneration for \(0<\alpha<1\) [2103.04991]. ## 5. Bounds, comparison principles, and coarse control A large part of the modern theory concerns upper and lower bounds on normalized Steklov eigenvalues. For a bounded domain \(\Omega\subset N^{n+1}\) with boundary \(\Sigma\), the normalized eigenvalues are \[ \bar\sigma_k(\Omega)=\sigma_k(\Omega)|\Sigma|^{1/n}, \] and, in the weighted setting, \[ \bar\sigma_k(M,\delta)=\sigma_k(M,\delta)\,m(\Sigma,\delta)\,|\Sigma|^{1/n}. \] If the ambient metric is conformally equivalent to a complete metric with nonnegative Ricci curvature, then \[ \bar\sigma_k(\Omega,\delta) \le \frac{\gamma(n)}{I(\Omega)^{\frac{n-1}{n}}}\, k^{\frac{2}{n+1}}, \qquad I(\Omega)=\frac{|\Sigma|}{|\Omega|^{n/(n+1)}}. \] For bounded domains in \(\mathbb R^{n+1}\), \(\mathbb H^{n+1}\), or a hemisphere of \(\mathbb S^{n+1}\), this yields the uniform bound \[ \bar\sigma_k(\Omega)\le C_n\,k^{\frac{2}{n+1}}. \] On compact orientable surfaces of genus \(\gamma\), \[ \bar\sigma_k(M)\le C\left\lfloor\frac{\gamma+3}{2}\right\rfloor k. \] These inequalities are not Weyl-sharp in the exponent of \(k\), but they are uniform and geometric [1103.2863]. There are also direct comparison inequalities between Steklov eigenvalues and boundary Laplace eigenvalues. If \(\lambda_k(\Sigma)\) are the Laplace–Beltrami eigenvalues of the boundary and \[ \bar\lambda_k(\Sigma)=\lambda_k(\Sigma)|\Sigma|^{2/n}, \] then, under nonnegative Ricci curvature and a mild localization hypothesis, \[ \bar\lambda_k(\Sigma)\,\bar\sigma_l(\Omega) \le \kappa_n \left(\frac{|\Sigma|}{|\Omega|}\right)^{3/n} k^{2/n}l^{2/(n+1)}. \] This makes precise the principle that large boundary Laplace eigenvalues constrain normalized Steklov eigenvalues [1103.2863]. For manifolds with cylindrical boundary and bounded geometry, Steklov eigenvalues can also be compared to those of coarse discretizations. An \(\epsilon\)-discretization produces a graph with boundary \((\Gamma_M,V_\Sigma)\) roughly isometric to the manifold, and there exist constants \(a,b>0\), depending only on the geometry class and \(\epsilon\), such that \[ a<\frac{\sigma_2(M)}{\sigma_2(\Gamma_M,V_\Sigma)}