Magnetic Steklov Operators

• Magnetic Steklov operators are spectral boundary operators that extend the classical Dirichlet-to-Neumann map by incorporating a magnetic potential in scalar, vector, or differential-form settings.
• They play a crucial role in spectral geometry, inverse problems, and electromagnetic PDEs, enabling explicit spectral calculations and variational characterizations via gauge invariance.
• The operators serve as key tools for analyzing eigenvalue asymptotics, isoperimetric inequalities, and eigenfunction localization, thereby advancing both theoretical and numerical investigations.

A magnetic Steklov operator is a spectral boundary operator that generalizes the classical Steklov (Dirichlet-to-Neumann) operator by incorporating a magnetic potential, acting on scalar, vector, or differential-form-valued fields, and arising in both scalar and electromagnetic PDEs on domains with boundary. In the presence of a magnetic field, the operator encodes how boundary values of a solution to an elliptic PDE determine the boundary values of its magnetic (co)normal derivatives. Magnetic Steklov operators appear naturally in the analysis of magnetic Laplacians, Schrödinger operators with magnetic/electric potentials, and Maxwell’s equations; they are also focal objects for spectral geometry, inverse problems, and boundary value analysis in both mathematics and mathematical physics.

1. Formal Definition and Variants

Let $(M, g)$ be a compact Riemannian manifold with smooth boundary $\partial M$, and let $A \in \Omega^1(M; \mathbb{R})$ denote a magnetic potential.

• For functions, the magnetic Laplacian is $L_A = (d + iA)^*(d + iA)$, acting on $u \in C^\infty(M)$, leading to the boundary value problem

$$\begin{cases} L_A u = 0 & \text{in } M, \ u|_{\partial M} = f \in C^\infty(\partial M). \end{cases}$$

The magnetic Steklov operator is then the Dirichlet-to-Neumann map:

$$\Lambda_A(f) = \left. (d + iA)u(\nu) \right|_{\partial M},$$

where $\nu$ is the outward unit normal to $\partial M$.

• For systems and differential forms, and in the electromagnetics (Maxwell) setting, the magnetic Steklov operator generalizes as the operator which, given tangential electric (or magnetic) fields on the boundary, returns the tangential magnetic (or electric) field, with the PDE incorporating terms involving the magnetic (vector) potential (Ferraresso et al., 2022, Halla, 2019, Chakradhar et al., 10 Nov 2025).
• On differential forms, the operator is defined via the magnetic Hodge Laplacian $\Delta_A = d_A \delta_A + \delta_A d_A$ ($\partial M$0) with appropriate absolute or relative boundary conditions, leading to a family of Steklov operators $\partial M$1 on $\partial M$2-forms (Chakradhar et al., 10 Nov 2025).

2. Spectral Theory and Analytic Structure

Magnetic Steklov operators are classical, self-adjoint pseudo-differential operators of order one on $\partial M$3, with a discrete real spectrum

$\partial M$4

when acting on functions ($\partial M$5); for vector or form-valued settings, negativity or sign symmetry may appear (Chakradhar et al., 2024, Ferraresso et al., 2022, Chakradhar et al., 10 Nov 2025). For Maxwell's equations in the self-adjoint case, the spectrum consists of three parts: essential spectrum at zero, an infinite negative sequence accumulating only at zero, and an infinite positive sequence accumulating only at infinity (Halla, 2019).

In all settings, there is an orthonormal basis of boundary eigenfunctions or eigenmodes, providing a natural Fourier basis for $\partial M$6 or the relevant Sobolev/trace spaces (Lamberti et al., 2020, Chakradhar et al., 10 Nov 2025). Explicit spectra are available in symmetric geometries (e.g., the disk, the ball, annuli) for certain choices of potentials, often involving Bessel or Laguerre functions, and their angular momentum decompositions reflect the interaction with the magnetic flux (Chakradhar et al., 2024, Chakradhar et al., 10 Nov 2025, Colbois et al., 2022).

3. Key Operator-Theoretic and Geometric Properties

Gauge Invariance

The magnetic Steklov operator is invariant up to conjugacy under gauge transformations:

$\partial M$7

so only the magnetic field $\partial M$8 (the curvature of $\partial M$9) is physically and spectrally relevant in most regimes (Chakradhar et al., 2024, Ferreira et al., 17 Jan 2026). In the context of eigenvalue multiplicity and spectral determination, only gauge-equivalence classes of $A \in \Omega^1(M; \mathbb{R})$0 can be reconstructed from Steklov data.

Variational Principles

For smooth domains $A \in \Omega^1(M; \mathbb{R})$1,

$A \in \Omega^1(M; \mathbb{R})$2

characterizes the lowest Dirichlet-to-Neumann eigenvalue (with analogous forms for $A \in \Omega^1(M; \mathbb{R})$3-forms and systems) (Kachmar et al., 19 Feb 2026, Colbois et al., 2022, Chakradhar et al., 10 Nov 2025).

Asymptotics and Trace Formulas

The magnetic Steklov spectrum admits Weyl-type asymptotics, with leading behavior governed by boundary geometry and magnetic fluxes, and subleading/splitting terms encoding boundary normal jets of the magnetic and electric potentials (Cekić et al., 2024, Ferreira et al., 17 Jan 2026). On Anosov boundaries with simple length spectrum, precise wave trace expansions relate Steklov spectral singularities to closed geodesics and allow for full Taylor-jet (all normal derivatives) determination at the boundary. The subprincipal symbol of the operator encodes, to leading order, the restriction of $A \in \Omega^1(M; \mathbb{R})$4 and $A \in \Omega^1(M; \mathbb{R})$5 to the boundary (Ferreira et al., 17 Jan 2026).

4. Inverse Spectral Problems and Determination Results

Under non-degenerate geometric conditions, such as Anosov boundary with simple length spectrum, the full Steklov spectrum of $A \in \Omega^1(M; \mathbb{R})$6 uniquely determines:

• The boundary values and all boundary-normal derivatives (Taylor series) of the magnetic and electric potentials, up to the addition of a gauge exact form to $A \in \Omega^1(M; \mathbb{R})$7.
• Under analyticity, the full jets determine $A \in \Omega^1(M; \mathbb{R})$8 and $A \in \Omega^1(M; \mathbb{R})$9 uniquely near the boundary (Ferreira et al., 17 Jan 2026).

On compact surfaces with boundary, spectral asymptotics detect the set of boundary component lengths, magnetic parallel transport, and fluxes modulo obstructions from degenerate situations (coincidence of lengths or special fluxes), but the number of boundary components may not always be spectrally determined (Cekić et al., 2024).

For magnetic Steklov operators on differential forms, the possibility of reconstructing the magnetic potential up to gauge from the spectrum remains an open direction but is suggested by the structural analogy to scalar and vector-valued cases (Chakradhar et al., 10 Nov 2025).

5. Isoperimetric and Geometric Inequalities

Magnetic Steklov eigenvalues admit sharp isoperimetric inequalities extending the classical results of Szegő-Weinberger, Brock, Weinstock, and Fraser-Schoen:

• For planar domains with Aharonov–Bohm or uniform magnetic potential, the disk maximizes the first Steklov eigenvalue among domains of fixed area or perimeter, for magnetic fields of moderate strength (Colbois et al., 2022, Kachmar et al., 19 Feb 2026).
• Explicit bounds relate the first eigenvalue to magnetic Cheeger-type constants, combining geometric and topological information with the "frustration index" (a measure of magnetic holonomy) (Chakradhar et al., 2024).
• For annular or higher genus surfaces, sharp upper bounds for normalized eigenvalues are expressed in terms of conformal modulus, with the maximizers exhibiting connections to linear Weingarten, weighted-minimal, or free-boundary catenoid-type surfaces, and geometry depends on flux parameters (Provenzano et al., 2023).

6. Model Problems and Explicit Solution Regimes

In planar or Euclidean geometry, explicit magnetic Steklov spectra are computable for disks, balls, or annuli, with the eigenvalues given by closed formulas involving magnetic flux parameters, angular momentum quantum numbers, and special functions (Bessel, Laguerre):

• On the Euclidean disk with Aharonov–Bohm potential of flux $L_A = (d + iA)^*(d + iA)$0 at the center, $L_A = (d + iA)^*(d + iA)$1, shifting the angular Fourier modes (Colbois et al., 2022).
• For the ball in 2D or 4D, with Killing field potential, Laguerre polynomial structure appears in the eigenvalues and multiplicities (Chakradhar et al., 2024, Chakradhar et al., 10 Nov 2025).
• In Maxwell or electromagnetic settings, vector spherical harmonics permit explicit computation in the unit ball, separating TE and TM modes (Ferraresso et al., 2022).

Strong and weak magnetic field limits yield contrasting asymptotics for the lowest eigenvalues in exterior domains, with flux dependence entering in distinct orders of the expansion, reflecting semiclassical and tunneling effects (Helffer et al., 25 Aug 2025).

7. Exponential Localization and Eigenfunction Behavior

For strong magnetic fields, ground-state Steklov eigenfunctions exhibit exponential concentration near boundary points where the magnetic field vanishes to maximal order:

$L_A = (d + iA)^*(d + iA)$2

with $L_A = (d + iA)^*(d + iA)$3 an Agmon-type distance to the minimal magnetic well on the boundary. This boundary-layer phenomenon generalizes previous semiclassical results and quantifies the localization in terms of field vanishing order (Shen, 18 Nov 2025).

8. Extensions: Forms, Systems, and Electromagnetic Settings

• Magnetic Steklov operators extend to arbitrary differential forms, with boundary conditions, variational formulations, and spectra exhibiting rich structure, including explicit violation of the diamagnetic inequality known from scalar Laplacians (Chakradhar et al., 10 Nov 2025).
• For Maxwell equations, Steklov-type problems (sometimes called electromagnetic or magnetic Steklov eigenproblems) act on tangential field traces, with operator pencils in $L_A = (d + iA)^*(d + iA)$4 and sophisticated spectral structure including block operator representations and Fredholm theory (Ferraresso et al., 2022, Halla, 2019, Halla, 2019, Lamberti et al., 2020).
• Modified versions (using tangential projections or Helmholtz decompositions) ensure desirable compactness and variational properties for numerical and inverse applications (Halla, 2019).

Table: Selected Model and Theoretical Results

Setting	Main Explicit Spectral Result	Reference
Planar disk with AB flux $L_A = (d + iA)^*(d + iA)$5	$L_A = (d + iA)^*(d + iA)$6	(Colbois et al., 2022)
Manifold with Anosov boundary	Steklov spectrum determines full boundary Taylor jets of $L_A = (d + iA)^*(d + iA)$7	(Ferreira et al., 17 Jan 2026)
Surface, $L_A = (d + iA)^*(d + iA)$8 boundary components	Two arithmetic progressions per component (flux shift)	(Cekić et al., 2024)
$L_A = (d + iA)^*(d + iA)$9-forms on ball, potential $u \in C^\infty(M)$0	Laguerre polynomial formulas for Steklov eigenvalues	(Chakradhar et al., 10 Nov 2025)
Maxwell (unit ball, vector fields)	TE/TM modes via spherical Bessel functions	(Ferraresso et al., 2022)

These results collectively establish the central role of magnetic Steklov operators in boundary spectral geometry and mathematical physics, unifying abstract analytic properties, explicit spectral calculations, inverse boundary problems, and geometric extremal theory. The interplay between gauge, global topology, boundary geometry, and magnetic structure produces a spectrum of phenomena not present in the classical, non-magnetic setting.