Surface Maxwell Model Approaches

• Surface Maxwell Model is a family of frameworks that organizes electromagnetic fields via surface-localized data across gravitational, optical, and fluid systems.
• It highlights how surface data influences key phenomena such as modified gravitational inequalities, interface electromagnetic waves, and numerical boundary algorithms.
• The model integrates geometric, topological, and dynamic aspects from Einstein–Maxwell theories to discrete numerical implementations, enabling practical analysis of complex boundaries.

Searching arXiv for recent and relevant papers on "surface Maxwell model" and closely related formulations. “Surface Maxwell Model” is not a single canonical construction in the arXiv literature surveyed here. The term is used for several technically distinct frameworks in which Maxwell fields are organized by a distinguished surface: a compact two-surface in Einstein–Maxwell geometry, a material interface supporting surface electromagnetic waves, a cutoff surface in fluid/gravity duality, a wall in gas-surface interaction theory, a codimension-two conformal defect, or a moving boundary in electrodynamics and free-interface problems (Lee et al., 2020, Bliokh et al., 2019, Niu et al., 2011, Chen et al., 2022, Herzog et al., 2022, Wang, 2022). This suggests that the phrase is best understood as a family of surface-based Maxwell formulations rather than as a standardized model class.

1. Surface-organized Maxwell theories

Across the literature, the recurring structure is that the Maxwell sector does not enter only through bulk field equations or global charge data; it is resolved relative to a surface, and the surface data become dynamically or geometrically essential.

Usage	Surface object	Maxwell-sector role
Einstein–Maxwell strong-gravity surfaces	LTS, DTTS, AGPS	local energy density and pressure/tension correct areal bounds
Surface-wave optics	planar material interface	topological interface states; TE/TM surface Maxwell waves
Fluid/gravity and covariant phase space	cutoff surface, codimension-two boundary	forcing term, edge modes, $U(1)$ surface symmetry
Rarefied-gas theory	solid wall	diffuse/specular reflection with accommodation
Moving-medium electrodynamics	moving/deforming interface	flux transport through time-dependent surfaces

In the Einstein–Maxwell setting, the surface organizes gravitational inequalities. In optics and plasmonics, the surface is the locus of localized electromagnetic modes. In kinetic theory, the surface is a scattering boundary. In gauge/gravity and covariant phase space, the surface carries quasi-local stress tensors, edge modes, or symmetry generators (Lee et al., 2020, Lee et al., 2024, Bliokh et al., 2019, Nakata et al., 2022, Niu et al., 2011, Setare et al., 2018).

A common misconception is that “surface Maxwell model” must denote a surface electromagnetic constitutive law in the narrow engineering sense. The literature surveyed here indicates a broader usage: the phrase can denote geometric inequalities, topological interface theories, defect classifications, or numerical boundary algorithms, provided that the Maxwell sector is fundamentally organized by surface data rather than only by bulk variables. This is an interpretive synthesis of the cited usages.

2. Einstein–Maxwell geometry of strong-gravity surfaces

In general relativity, one important usage concerns compact two-surfaces that probe strong gravity. For a loosely trapped surface (LTS), a compact two-surface $S_0$ in a spacelike hypersurface $\Sigma$ satisfies

$k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$

where $k'$ is the derivative of the mean curvature along the outward spacelike normal. For a dynamically transversely trapping surface (DTTS), a closed orientable two-surface $\sigma_0\subset\Sigma$ admits a timelike hypersurface $S$ intersecting $\Sigma$ at $\sigma_0$ such that

$$\bar k=0,\qquad \max(\bar K_{ab}k^ak^b)=0,\qquad {}^{(3)}\bar{\pounds}_{\bar n}\bar k\le 0.$$

These notions were introduced to characterize strong-gravity regions without exact spherical symmetry (Lee et al., 2020).

The central Einstein–Maxwell result is that the Maxwell field changes the expected photon-sphere-type area bounds through local surface corrections. For LTSs,

$S_0$0

with

$S_0$1

For convex DTTSs under time symmetry,

$S_0$2

with

$S_0$3

The key point is that the electromagnetic field contributes not only through the total charge $S_0$4, but also through the energy density and the stress/pressure or tension of the field on the surface itself (Lee et al., 2020).

This local split is encoded in the Maxwell stress tensor

$S_0$5

together with

$S_0$6

and

$S_0$7

Because $S_0$8, the LTS bound is tighter than the naive Reissner–Nordström photon-sphere bound. By contrast, $S_0$9 can be positive or negative, so the DTTS bound can become stronger or weaker depending on the field orientation relative to the surface (Lee et al., 2020).

The same logic is extended by the theory of attractive gravity probe surfaces (AGPSs), which generalizes minimal-surface inequalities to five surface classes, including weak-gravity regions. In that framework, electric and magnetic charges, angular momentum, gravitational-wave energy, and matter contributions enter generalized Penrose-type areal inequalities, and one obtains a notion of surface extremality analogous to charged-black-hole extremality (Lee et al., 2024). In this usage, a “surface Maxwell model” is a geometry in which the electromagnetic sector constrains the allowed surface area and extremality threshold through both local field decomposition and global charges.

A numerical example in the Majumdar–Papapetrou two-black-hole spacetime shows that for common marginally DTTSs surrounding both black holes, $\Sigma$0 is very small, typically below $\Sigma$1, and $\Sigma$2. In that case the corrected bound and the naive photon-sphere area nearly coincide, but the correction is structurally present (Lee et al., 2020).

3. Surface Maxwell waves, topology, and defect triviality

In classical optics, “surface Maxwell model” can denote a topological reformulation of interface electromagnetic waves. For homogeneous, isotropic media with real-valued permittivity $\Sigma$3 and permeability $\Sigma$4, surface Maxwell waves are electromagnetic waves localized at a planar interface. Their origin is described through the spectrum of a medium-dependent helicity operator, which is generically non-Hermitian even in lossless media. The topological invariant is a $\Sigma$5 number, or equivalently a pair of $\Sigma$6 numbers,

$\Sigma$7

and bulk-boundary correspondence is expressed as

$\Sigma$8

In this framework, if only $\Sigma$9 changes sign one gets one TM mode, if only $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$0 changes sign one gets one TE mode, and if both $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$1 and $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$2 change sign there are two surface modes, one TE and one TM (Bliokh et al., 2019).

A complementary formulation identifies surface plasmon polaritons as degenerated electric zero modes protected by a hidden combined symmetry $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$3, where $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$4 is mirror reflection and $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$5 acts as

$k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$6

For a medium satisfying

$k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$7

the relevant topological indicator is defined from surface impedance,

$k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$8

The condition $k|_{S_0}>0,\qquad k'|_{S_0}\ge 0,$9 gives the surface-plasmon condition, and the bulk–edge argument is formulated directly in terms of the sign of the DC surface reactance. This theory also emphasizes that the zero modes exhibit vector-field rotation and are related to Keller–Dykhne duality (Nakata et al., 2022).

A different surface-theoretic result arises in Euclidean free Maxwell theory with a codimension-two defect. Reflection-positive, conformally invariant defects are severely constrained: only generalized free fields can appear in the defect operator product expansion of the bulk Maxwell field, and the resulting defect correlators satisfy Wick’s theorem. The allowed defect operators are a scalar with $k'$0, $k'$1, and defect vectors with $k'$2; more generally, the defect operator algebra closes only under generalized-free double-twist composites (Herzog et al., 2022). Thus, while interface Maxwell problems can exhibit nontrivial topological surface modes, conformal surface defects in free 4D Maxwell theory are “trivial” in the precise generalized-free-field sense.

4. Cutoff surfaces, edge modes, and surface effective theories

In fluid/gravity duality, a charged AdS black brane with a Dirichlet cutoff surface at $k'$3 yields an effective surface fluid. Under the non-relativistic long-wavelength scaling

$k'$4

the Maxwell equations at $k'$5 reduce to

$k'$6

so the Maxwell sector enforces incompressibility. The momentum constraint on the cutoff surface becomes the incompressible Navier–Stokes equation with forcing,

$k'$7

with

$k'$8

In Einstein–Maxwell theory,

$k'$9

independent of both the cutoff $\sigma_0\subset\Sigma$0 and the black brane charge, whereas in Gauss–Bonnet–Maxwell theory the ratio remains cutoff-independent but depends on the charge density and the Gauss–Bonnet coupling $\sigma_0\subset\Sigma$1 (Niu et al., 2011). In this usage, the Maxwell sector appears on the surface as an external body force rather than as an independent boundary gauge field.

A covariant-phase-space usage makes the Maxwell boundary degrees of freedom explicit. In Einstein–Maxwell theory, combined diffeomorphism and $\sigma_0\subset\Sigma$2 gauge transformations are written as $\sigma_0\subset\Sigma$3, and the ordinary symplectic potential is not invariant under these combined transformations. Introducing new boundary fields $\sigma_0\subset\Sigma$4 and $\sigma_0\subset\Sigma$5 extends phase space and produces edge modes. The resulting surface-preserving symmetry group is the semidirect sum of the two-dimensional diffeomorphism group on a spacelike codimension-two surface with $\sigma_0\subset\Sigma$6 and $\sigma_0\subset\Sigma$7, and the Casimir of $\sigma_0\subset\Sigma$8 is the area element (Setare et al., 2018).

In condensed-matter effective theory, the gapped surface of a three-dimensional topological insulator in an external electromagnetic field can be described by a Chern–Simons theory whose gauge connection takes values in the Maxwell algebra,

$\sigma_0\subset\Sigma$9

This construction incorporates both Lorentz and magnetic-translation symmetries of the surface states. The resulting boundary action contains the usual electromagnetic Chern–Simons term, geometric and torsional terms, and a non-minimal coupling

$S$0

interpreted as a relativistic version of the Wen–Zee term. In flat space, the theory reproduces the Hall response and supports a $S$1 chiral conformal field theory on defect lines of the gapped boundary (Palumbo, 2016).

5. Gas-surface interaction, moving boundaries, and free interfaces

In rarefied-gas dynamics, the Maxwell gas-surface interaction model is a wall boundary condition that mixes diffuse and specular reflection through an accommodation coefficient $S$2: $S$3 Here $S$4 gives full diffuse reflection and $S$5 gives pure specular reflection. The wall density $S$6 is fixed by the no-penetration condition, and off-grid specular reflection on unstructured velocity space is reconstructed by interpolation together with macro-conservation and micro-consistency corrections. This algorithm is designed for both DVM and UGKS and is validated on supersonic and hypersonic rarefied flows (Chen et al., 2022).

A related but distinct refinement is the Maxwell–Smoluchowski theory for Knudsen-regime transport in straight channels. There the tangential momentum accommodation coefficient $S$7 determines the self-diffusivity through

$S$8

The key result is that for sufficiently low roughness,

$S$9

where $\Sigma$0 is a flatness parameter, $\Sigma$1 is a shape parameter extracted from effective microgeometry, and $\Sigma$2 depends only on the channel cross-section. For a circular channel,

$\Sigma$3

For sphere-packing microgeometry,

$\Sigma$4

The classical accommodation coefficient is thereby geometrized (Chumley et al., 2023).

For charged media with time-dependent volume, shape, and boundary moving with an arbitrary slow velocity field $\Sigma$5, a lab-frame derivation from the integral laws yields generalized Maxwell equations

$\Sigma$6

$\Sigma$7

together with the continuity equation

$\Sigma$8

A further surface polarization term $\Sigma$9 is introduced through

$\sigma_0$0

to represent surface electrification and mechano-electrical coupling (Wang, 2022). In this usage, the model is surface-based because the derivation tracks fluxes through moving surfaces.

At a plasma–vacuum free boundary, the Maxwell sector can also be destabilizing. For the linearized ideal MHD–Maxwell free interface problem with surface tension, however, one obtains an a priori estimate in a conormal Sobolev space without assuming any stability conditions on the basic state. The stabilizing mechanism is the Laplace–Young term

$\sigma_0$1

which yields positive coercive control of $\sigma_0$2 at the boundary (Trakhinin, 2024).

6. Microscopic and discrete surface realizations

In nonlinear plasmonics, a Maxwell–hydrodynamic model for metallic metasurfaces provides a microscopic, time-domain realization of surface and bulk nonlinearities. The electromagnetic fields obey

$\sigma_0$3

while the electron fluid satisfies

$\sigma_0$4

$\sigma_0$5

Surface nonlinearity is not imposed through an explicit surface susceptibility; it emerges from charge accumulation and current discontinuity at metal boundaries, while bulk nonlinearity comes from convective acceleration, magnetic force, and charge-density coupling. The model is implemented in a Yee-grid FDTD framework with an ADE update and a two-step splitting of the current equation (Fang et al., 2017).

A fully different realization appears in the discrete exterior calculus formulation of semi-discrete Maxwell equations. In three dimensions, time remains continuous while space is discretized: $\sigma_0$6 The two-dimensional specialization is naturally interpreted as a Maxwell theory on a surface mesh. On a combinatorial two-dimensional torus, with $\sigma_0$7 and $\sigma_0$8, the equations reduce to a finite-dimensional first-order linear ODE system,

$\sigma_0$9

together with Gauss-law constraints. The construction preserves the coboundary identity $$\bar k=0,\qquad \max(\bar K_{ab}k^ak^b)=0,\qquad {}^{(3)}\bar{\pounds}_{\bar n}\bar k\le 0.$$0, a discrete Poynting theorem, gauge invariance, and a Laplacian/wave-equation structure, although the discrete Hodge star satisfies a shifted $$\bar k=0,\qquad \max(\bar K_{ab}k^ak^b)=0,\qquad {}^{(3)}\bar{\pounds}_{\bar n}\bar k\le 0.$$1 relation rather than the continuum identity (Sushch, 30 Oct 2025).

Taken together, these realizations show that the phrase “Surface Maxwell Model” spans several non-equivalent but structurally related programs. In every case, the decisive move is to treat a surface as the primary organizer of Maxwell data: as a geometric probe, interface, defect, wall, cutoff, moving boundary, or discrete manifold. The Maxwell sector then enters not merely through bulk field values, but through surface-localized geometry, topology, transport, symmetry, or numerical flux structure.