Meissner Solution in Superconducting Systems

• Meissner solution is a configuration in superconductors where magnetic fields decay exponentially, characterized by the London penetration depth and vortex-free states.
• It underpins various phenomena from classical superconductivity to black-hole electrodynamics and photonic or spin analogues in insulating systems.
• The concept integrates boundary value approaches, variational principles, and thermodynamic analyses to explain flux expulsion and current screening under different conditions.

In the research literature, “Meissner solution” denotes several related but non-identical objects. In superconductivity it most commonly refers to the vortexless screening configuration solving the London or Ginzburg–Landau boundary-value problem, with magnetic field expelled from the bulk. In black-hole electrodynamics it refers either to horizon-flux expulsion in the extremal limit or, in astrophysical jet theory, to the split-monopole horizon topology proposed to evade that expulsion. More recent work extends the terminology to spin, holographic, nonreciprocal, and photonic settings in which a gauge field or a gauge-field analogue is screened by persistent currents or counterflow patterns (Yoshioka, 2012, Pan, 2018, Penna, 2014, Petrescu et al., 2013, Natsuume et al., 2022, Li et al., 3 Mar 2025).

1. Conventional superconducting meaning

In the standard superconductivity setting, the Meissner solution is the static London profile for a superconducting half-space. In London gauge, the supercurrent and vector potential satisfy

$$\mathbf{J}_s = -\frac{n_s e^2}{m}\,\mathbf{A}, \qquad \nabla \times \mathbf{J}_s = -\frac{n_s e^2}{m}\,\mathbf{B},$$

which combine with Maxwell’s equation to give

$$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$

For a semi-infinite superconductor occupying $x>0$ with an applied field parallel to the surface, the canonical Meissner solution is

$B_z(x) = B_0 e^{-x/\lambda}.$

This identifies $\lambda$ as the London penetration depth and the surface current layer as the locus of magnetic screening (Yoshioka, 2012).

A central distinction, emphasized repeatedly in the literature, is that the Meissner effect is not identical to perfect conductivity. In a classical perfect conductor, flux conservation and Lenz’s law would preserve whatever interior field existed when the conductivity became perfect, producing a history-dependent final state. By contrast, the Meissner state is an equilibrium state in which the interior field vanishes below the relevant critical field, independent of magnetic history (Yoshioka, 2012).

The energetic interpretation is disputed. Yoshioka argues that the Meissner profile cannot be derived from classical mechanics alone because expelling magnetic field from a volume requires positive work equal to the field energy removed, and because the stabilizing condensation energy is quantum mechanical (Yoshioka, 2012). By contrast, the review literature represented by Essén and Fiolhais argues that, once internal field energy and dissipationless current rearrangement are treated correctly, flux expulsion can be understood as approach to a magnetostatic energy minimum, with the London equation emerging from minimization of magnetic plus kinetic energy (Essen et al., 2011). A more recent thermodynamic treatment again reproduces the static London equation and the profile $B(x)=B_0 e^{-x/\lambda}$, but does so from the principle that superfluid flow minimizes energy at fixed entropy, together with a modified first London equation that allows field evolution during the transition (Attard, 9 Sep 2025).

2. Ginzburg–Landau, TDGL, and rigorous vortexless solutions

Within Ginzburg–Landau theory, “Meissner solution” acquires a precise variational meaning: it is a vortexless critical point or minimizer with positive order-parameter modulus. In the rigorous three-dimensional treatment of type-II superconductors, a Meissner state is described by a pair $(f,\mathbf A)$, where $f>0$ on $\Omega$ and $\mathbf A$ is defined on all of $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$0 with vanishing inner trace of the normal component on $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$1. The interior equations are

$$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$2

supplemented by the exterior condition $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$3 in $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$4, tangential continuity across $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$5, and $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$6 at infinity (Pan, 2018). That work shows that the Meissner solution is smooth in $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$7, that the exterior regularity of $$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$8 can be poor, and that the system naturally decomposes into an interior boundary-value problem and an exterior problem (Pan, 2018).

The large-$$\nabla^2 \mathbf{B} = \frac{1}{\lambda^2}\,\mathbf{B}, \qquad \lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}.$$9 limit yields a nonlinear curl system. As $x>0$0, the magnetic potential part of the Meissner solution converges in $x>0$1 to a solution of

$x>0$2

which is identified as the correct three-dimensional limit of the Meissner system (Pan, 2018).

In inhomogeneous extreme type-II superconductors with pinning term $x>0$3, the Meissner solution is the unique vortexless configuration that globally minimizes the Ginzburg–Landau energy below the first critical field $x>0$4, and a unique Meissner-type local minimizer persists for fields slightly below the superheating field $x>0$5 (Díaz-Vera et al., 2024). The analysis introduces the canonical vortexless pair $x>0$6, with $x>0$7 solving the pinned zero-field problem and

$x>0$8

The first critical field satisfies

$x>0$9

so pinning modifies both the Meissner profile and the onset of vortices through $B_z(x) = B_0 e^{-x/\lambda}.$0 and $B_z(x) = B_0 e^{-x/\lambda}.$1 (Díaz-Vera et al., 2024).

A finite-size TDGL perspective sharpens the distinction between vortexless and bulk behavior. For homogeneous square samples in the Meissner state, the ideal bulk relation is

$B_z(x) = B_0 e^{-x/\lambda}.$2

in cgs units, but finite samples deviate from this line even without vortices. Using the initial Meissner slope

$B_z(x) = B_0 e^{-x/\lambda}.$3

the cited simulations define bulk-like behavior by $B_z(x) = B_0 e^{-x/\lambda}.$4. For $B_z(x) = B_0 e^{-x/\lambda}.$5, the inferred threshold lengths are extremely large, with $B_z(x) = B_0 e^{-x/\lambda}.$6 and $B_z(x) = B_0 e^{-x/\lambda}.$7, showing that a vortex-free Meissner solution need not yet realize bulk magnetization (Presotto et al., 2013).

3. Inhomogeneity, defects, and partial Meissner response

Once the penetration depth becomes spatially inhomogeneous, the Meissner solution is no longer described by a single homogeneous exponential. In the London model with $B_z(x) = B_0 e^{-x/\lambda}.$8, Kogan and Kirtley develop a first-order perturbation theory for the Meissner response of a half-space with planar or point defects. For a planar defect perpendicular to the surface,

$B_z(x) = B_0 e^{-x/\lambda}.$9

and the defect correction to the reflected field is expressed by an explicit Fourier-space kernel $\lambda$0, enabling forward modelling of scanning SQUID and MFM signals from local superfluid-density variations (Kogan et al., 2011). In that setting, “Meissner solution” is the local response of the inhomogeneous London equation rather than a single scalar penetration depth.

A complementary issue is partial expulsion during the superconducting transition of type-II niobium. A superconducting magnetic flux lens experiment on Nb sheet material quantifies the expelled fraction by collimating flux through a central aperture during cooldown. Across 669 transitions, the expulsion ratio $\lambda$1 increases approximately linearly with spatial temperature gradient $\lambda$2, with campaign-dependent slopes $\lambda$3–$\lambda$4, maximum $\lambda$5, and maximum $\lambda$6. The principal empirical conclusion is that expulsion improves with spatial temperature gradient and is independent on the cooling rate once $\lambda$7 is fixed (Ivanov et al., 2021). In this usage, the “Meissner solution” is an experimentally optimized, partially achieved state in which vortex trapping is minimized rather than mathematically absent.

4. Unconventional superconducting Meissner profiles

Several recent works generalize the Meissner solution by altering the screening kernel itself. In Nb thin films coated with a self-assembled monolayer of chiral molecules, low-energy muon spin rotation shows that the local field profile is no longer fitted by a single conventional slab solution. The fitted form is

$\lambda$8

with $\lambda$9, $B(x)=B_0 e^{-x/\lambda}$0, and $B(x)=B_0 e^{-x/\lambda}$1 for the chiral-molecule/Nb sample, compared with $B(x)=B_0 e^{-x/\lambda}$2 for bare Nb. The interpretation advanced there is that odd-frequency spin-triplet correlations generated by a spin-active chiral interface add a paramagnetic contribution to the conventional diamagnetic screening (Alpern et al., 2021).

In parity-mixed superconductors, Watanabe, Daido, and Yanase identify a nonreciprocal superfluid density $B(x)=B_0 e^{-x/\lambda}$3 and derive a sign-dependent penetration depth

$B(x)=B_0 e^{-x/\lambda}$4

The corresponding Meissner fields are

$B(x)=B_0 e^{-x/\lambda}$5

so magnetic screening becomes asymmetric under field reversal. For the UTe$B(x)=B_0 e^{-x/\lambda}$6 case study, the paper estimates $B(x)=B_0 e^{-x/\lambda}$7 without strong renormalization and $B(x)=B_0 e^{-x/\lambda}$8 for $B(x)=B_0 e^{-x/\lambda}$9, making the effect potentially measurable in heavy-fermion settings (Watanabe et al., 2021).

Holographic superconductors require a different modification: the boundary Maxwell field must be made dynamical. With the semiclassical boundary Maxwell equation imposed, the holographic Meissner effect reappears analytically, and for $(f,\mathbf A)$0 the penetration depth is

$(f,\mathbf A)$1

A distinctive result is that the bound current prevents the extreme Type I limit even as $(f,\mathbf A)$2, because $(f,\mathbf A)$3 remains finite in that limit (Natsuume et al., 2022).

A still more radical deformation is produced by a purely spatial Chern–Simons term in Landau–Ginzburg theory. There the magnetic penetration depth becomes complex, the field acquires periodic spatial oscillations rather than monotone decay, and a vortex solution with magnetic field inversion may emerge as parity breaking increases (Tao, 2016). This shifts the meaning of “Meissner solution” from pure expulsion to damped oscillatory screening.

5. Black-hole electrodynamics and the astrophysical “Meissner solution”

In Kerr spacetimes, the “black hole Meissner effect” is a geometric phenomenon: as $(f,\mathbf A)$4, smooth stationary axisymmetric electromagnetic fields such as Wald’s vacuum solution lose horizon-threading flux. Penna’s analysis argues that this expulsion is not universal. The key observation is that the proper radial distance across the extremal throat diverges, so a stationary axisymmetric field with nonzero transverse component at the horizon is generically expelled by the no-monopoles constraint. The exceptional case is a field that becomes entirely radial at the horizon. Penna proposes the split-monopole as the corresponding “Meissner solution,” with flux function

$(f,\mathbf A)$5

in the northern hemisphere. Because $(f,\mathbf A)$6 and $(f,\mathbf A)$7 remains finite on the horizon, the configuration evades the cylinder argument, satisfies Znajek regularity, and preserves Blandford–Znajek power

$(f,\mathbf A)$8

up to near-extremal spin (Penna, 2014).

The same literature emphasizes that ordinary Blandford–Znajek jets are impossible if the Meissner effect operates and expels the field, because horizon-threading flux is required both geometrically and in the membrane paradigm (Penna, 2014). GRMHD simulations reportedly show no jet quenching up to $(f,\mathbf A)$9 and spontaneously develop a large split-monopole component near the horizon, consistent with that proposal (Penna, 2014).

A complementary force-free analysis by Pan and Yu recasts the steady axisymmetric Kerr magnetosphere equation into a Poisson-like form

$f>0$0

and proves that a Kerr black-hole force-free magnetosphere does not possess the Meissner effect when

$f>0$1

Under those monotonicity conditions, magnetic flux threading the inner light surface must also thread the horizon, so flux expulsion is ruled out in the force-free setting (Gong et al., 2016).

In strong magnetic fields, the extremal magnetized Kerr–Newman near-horizon geometry becomes that of an unmagnetized extremal Kerr–Newman solution with effective parameters $f>0$2. The horizon magnetic flux through the upper hemisphere is

$f>0$3

equivalently

$f>0$4

and vanishes iff $f>0$5. In that strong-field near-horizon setting, complete flux expulsion is therefore restricted to special subfamilies rather than generic extremal magnetized solutions (Hejda et al., 2016).

6. Insulating and photonic analogues

The term has also migrated into systems where electric charge transport is absent or secondary. In the two-species bosonic lattice model of a bosonic Mott insulator, the total-density sector is gapped while the relative-phase sector remains coherent because of interspecies conversion $f>0$6. The effective spin gauge field is $f>0$7, and in the Mott phase the pinned relative phase yields

$f>0$8

In the Meissner state, $f>0$9, so

$\Omega$0

with spin penetration length

$\Omega$1

Here the “Meissner solution” is screening in the spin channel inside a charge-insulating phase, made possible by spin-charge separation and a pseudospin superfluid (Petrescu et al., 2013).

An explicitly photonic version appears in the quantum Rabi zigzag chain with staggered flux. In the Meissner superradiant phase, persistent chiral edge currents flow on the even and odd sublattices with opposite sign, satisfying

$\Omega$2

The cancellation is mediated by coherent cavity displacements rather than by Maxwell back-action, and the authors characterize it as directly analogous to Meissner surface currents. In this setting the “Meissner solution” is a many-body current pattern with global current cancellation, not literal magnetic-field expulsion (Li et al., 3 Mar 2025).

7. Conceptual issues, terminology, and recurring structure

A recurrent source of confusion is that the same label spans distinct ontological levels. In the superconducting literature, the Meissner solution may mean the London exponential profile, a vortexless Ginzburg–Landau minimizer, a finite-size TDGL screening branch, or an experimentally partial expulsion state shaped by vortex pinning and cooldown gradients (Yoshioka, 2012, Presotto et al., 2013, Pan, 2018, Ivanov et al., 2021). In black-hole physics, by contrast, the effect is geometric rather than material: it depends on near-horizon throat structure, stationarity, axisymmetry, and force-free regularity, not on condensation or diamagnetism (Penna, 2014, Gong et al., 2016).

Another recurring dispute concerns mechanism. One line of work treats the Meissner profile as classically derivable from correct magnetic and kinetic energy accounting, while another insists that equilibrium flux expulsion requires quantum condensation energy and therefore cannot be understood classically (Essen et al., 2011, Yoshioka, 2012). The literature represented here does not resolve that disagreement; it shows instead that “Meissner solution” is a term whose content depends on the modelling level at which one asks the question.

Across these settings, a common mathematical pattern nevertheless recurs. The solution is typically a distinguished branch selected by one of three criteria: a boundary-value problem for a screening equation, a variational principle favoring a vortexless or current-carrying state, or a regularity/topology constraint that preserves flux threading while excluding singular behavior. This suggests that the broadest encyclopedic meaning of “Meissner solution” is not a single formula, but a class of screening or counterflow configurations that realize magnetic exclusion—or a precise analogue of it—within the constitutive rules of the system under study.