Solenoidal Ohm's Law

• Solenoidal Ohm's law is a divergence-free constitutive formulation that attaches Ohm’s law to incompressible current sectors across various electrodynamic systems.
• It underpins advanced models like the incompressible Navier–Stokes–Maxwell framework by integrating pressure corrections to enforce current conservation and ensure energy identity balance.
• The formulation guides applications in nanoplasmonics and steady conduction, where enforcing a divergence-free current improves numerical stability and resolves field singularities.

Searching arXiv for papers directly relevant to “Solenoidal Ohm’s Law” and adjacent formulations. First, I’ll look for papers explicitly using the phrase or closely related formulations. Solenoidal Ohm’s law denotes a family of constitutive formulations in which Ohm’s law is attached to a divergence-free current, or to the divergence-free sector of a broader electrodynamic response. In the literature, the phrase is not used uniformly. It appears explicitly in incompressible Navier–Stokes–Maxwell theory as an Ohm law for a current $j$ satisfying $\nabla\cdot j=0$, while in other settings it is an interpretive label for either a transverse current sector that remains exactly Ohmic or a limiting current rendered divergence-free by continuity and Maxwell constraints (Guo et al., 16 Jul 2025, Mortensen et al., 2012, Lei et al., 2023).

1. Terminology and defining structure

In its most explicit contemporary form, solenoidal Ohm’s law is the incompressible current relation

$$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$

or equivalently

$j=\sigma\big(cE+\mathbb P(u\times B)\big),$

with $\mathbb P$ the Leray projection onto divergence-free vector fields. In this formulation, the scalar $\tilde p$ is not an electrostatic potential in the elementary sense but the pressure-like Lagrange multiplier that enforces the solenoidal constraint on $j$ (Guo et al., 16 Jul 2025).

A second meaning arises in nonlocal nanoplasmonics. There, the hydrodynamic constitutive law modifies only the compressional part of the current through a term proportional to $\nabla(\nabla\cdot J)$. If the current is purely solenoidal, $\nabla\cdot J_T=0$, the correction vanishes identically and the transverse sector obeys

$J_T=\sigma E_T.$

In that precise sense, the divergence-free sector retains an exact local Ohmic law inside a broader nonlocal constitutive theory (Mortensen et al., 2012).

A third, older and more classical structure appears in stationary conduction. If

$\nabla\cdot j=0$0

then

$\nabla\cdot j=0$1

Here the solenoidal property is not imposed separately but follows from combining generalized Ohm conduction with steady current conservation (Tachiquin, 2010).

2. Incompressible Navier–Stokes–Maxwell form

The most direct use of the phrase occurs in the singular limit from a two-fluid incompressible Navier–Stokes–Maxwell model to an incompressible Navier–Stokes–Maxwell system with solenoidal Ohm’s law. With

$\nabla\cdot j=0$2

and fixed speed of light $\nabla\cdot j=0$3, the limit $\nabla\cdot j=0$4 removes the inertial, transport, and diffusion terms from the current equation and yields

$\nabla\cdot j=0$5

$\nabla\cdot j=0$6

$\nabla\cdot j=0$7

together with

$\nabla\cdot j=0$8

The paper also writes the Ohm law in projected form as

$\nabla\cdot j=0$9

which makes explicit that the current law is the divergence-free correction of the naive pointwise relation $$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$0 (Guo et al., 16 Jul 2025).

This pressure-corrected structure is constitutively essential rather than decorative. The gradient term is the direct analogue of incompressible hydrodynamic pressure: it removes the non-solenoidal component of $$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$1. The same paper ties the law to an energy identity,

$$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$2

so the Ohm law supplies the $$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$3-control of $$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$4 needed to control the Lorentz force $$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$5. The convergence proof is global in time, strong in Sobolev spaces, and avoids the derivative loss present in earlier work by a frequency-envelope argument (Guo et al., 16 Jul 2025).

3. Kinetic origins and asymptotic derivations

A closely related incompressible law is derived from the diffusive limit of the two-species Vlasov–Maxwell–Landau system with Coulomb potential. The limiting incompressible two-fluid Navier–Stokes–Fourier–Maxwell system contains

$$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$6

equivalently

$$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$7

along with

$$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$8

$$j=\sigma(-\nabla \tilde p+cE+u\times B), \qquad \nabla\cdot j=0,$$9

The paper does not explicitly coin the phrase “solenoidal Ohm’s law,” but it states that in the incompressible Maxwell framework the limiting current is understood as solenoidal, $j=\sigma\big(cE+\mathbb P(u\times B)\big),$0, so the label is appropriate as an interpretation of the limit system (Lei et al., 2023).

The formal kinetic derivation of MHD Ohm’s law from the two-species Vlasov–Maxwell–Boltzmann system gives a complementary route. In the $j=\sigma\big(cE+\mathbb P(u\times B)\big),$1-$j=\sigma\big(cE+\mathbb P(u\times B)\big),$2 formulation, the leading-order equilibrium for the total distribution $j=\sigma\big(cE+\mathbb P(u\times B)\big),$3 is a local Maxwellian $j=\sigma\big(cE+\mathbb P(u\times B)\big),$4, while the charge-disparity variable $j=\sigma\big(cE+\mathbb P(u\times B)\big),$5 solves a linearized kinetic equation driven by $j=\sigma\big(cE+\mathbb P(u\times B)\big),$6. Taking the first velocity moment yields

$j=\sigma\big(cE+\mathbb P(u\times B)\big),$7

equivalently

$j=\sigma\big(cE+\mathbb P(u\times B)\big),$8

The same asymptotic hierarchy gives the charge conservation law

$j=\sigma\big(cE+\mathbb P(u\times B)\big),$9

In quasineutral or reduced-Maxwell limits one further obtains

$\mathbb P$0

and therefore

$\mathbb P$1

This identifies the solenoidal current condition not as an independent axiom, but as the MHD closure of the kinetic charge-balance structure. The paper also states that the Hall term appears as a higher-order correction rather than at leading order (Jang et al., 2012).

4. Transverse-sector Ohmicity in nonlocal electrodynamics

In nanoplasmonics, “solenoidal Ohm’s law” enters through a different mechanism. The local-response approximation uses

$\mathbb P$2

The hydrodynamic nonlocal theory replaces this by

$\mathbb P$3

with

$\mathbb P$4

Because the correction depends only on $\mathbb P$5, it acts directly only on the longitudinal current. Writing

$\mathbb P$6

one gets immediately

$\mathbb P$7

for the transverse or solenoidal component, whereas the longitudinal response remains nonlocal and pressure-stiffened (Mortensen et al., 2012).

The physical importance of this distinction is that longitudinal current is tied to induced charge through

$\mathbb P$8

Local Drude theory therefore concentrates induced charge into a mathematical surface sheet, while the hydrodynamic model spreads charge over a finite near-surface region. The paper associates the nonlocal length scale with $\mathbb P$9 in the lossless scaling sense and estimates the correction strength as

$\tilde p$0

This is why the deviations from local Ohm conduction are strongest for characteristic dimensions in the $\tilde p$1–$\tilde p$2 nm range, in narrow gaps, sharp tips, corrugations, and highly curved surfaces. The same framework predicts size-dependent blue shifts of plasmon resonances, the appearance of confined longitudinal bulk-plasmon resonances above the plasma frequency, and the regularization of field singularities in dimer gaps and hot spots (Mortensen et al., 2012).

5. Reduced-plasma and relativistic generalizations

In drift-reduced plasma fluid models, the nearest analogue of solenoidal Ohm’s law is a parallel Ohm law embedded in a system constrained by

$\tilde p$3

The reduced equations include

$\tilde p$4

together with current continuity and a vorticity equation generated by balancing the divergence of the parallel current with polarization current. The paper’s principal conclusion is that if the electron-pressure gradient is retained in Ohm’s law, then electromagnetic effects through $\tilde p$5 and finite electron inertia $\tilde p$6 must also be kept; otherwise the model supports unphysical parallel propagation or diffusion (Dudson et al., 2021).

Relativistic response theory gives a different caution. Starting from the covariant linear response law

$\tilde p$7

the response tensor must satisfy the gauge-invariance and continuity constraints

$\tilde p$8

This framework is exactly equivalent to the microscopic conductivity law

$\tilde p$9

Here a genuinely solenoidal current is not generic; rather,

$j$0

remains the basic relation, and $j$1 requires additional restrictions such as stationary or charge-free conditions (Starke et al., 2014).

The $j$2 covariant relativistic plasma treatment pushes the same point further. Its generalized Ohm law is a propagation equation for the spatial current $j$3, not initially an algebraic constitutive law. In the cold proton–electron limit it recovers

$j$4

so Hall and Biermann-battery terms emerge as explicit corrections. Under global neutrality, Gauss’s law becomes $j$5, which in an irrotational frame reduces to $j$6; this is the closest explicit solenoidal electric-field statement in that formulation (0711.3573).

A general-relativistic two-component plasma around a rotating compact object adds further terms. For stationary plasma the paper derives

$j$7

with Hall, pressure-gradient, and frame-dragging contributions. It predicts azimuthal current generation from radial current through the gravitomagnetic term, but also concludes that classical Cowling-type antidynamo behavior remains valid for typical neutron-star plasma parameters (Ahmedov, 2010).

6. Static conduction and computational realizations

In stationary inhomogeneous conduction, the solenoidal structure is built directly into the PDE

$j$8

With

$j$9

one has

$\nabla(\nabla\cdot J)$0

The pseudoanalytic-function formulation develops this into a constructive theory of divergence-free current fields in separable media. In two dimensions, the conductivity equation is reformulated as a Vekua equation and expanded in formal powers; in three dimensions it is rewritten quaternionically as

$\nabla(\nabla\cdot J)$1

This makes the paper one of the clearest examples in which a solenoidal current law is realized as an elliptic conservation law rather than an incompressible projection (Tachiquin, 2010).

Numerical formulations do not always encode solenoidality inside Ohm’s law itself. In the space-time finite-element treatment of the vectorial wave equation under Ohm’s law, the constitutive relation is simply

$\nabla(\nabla\cdot J)$2

and, under the Weyl gauge $\nabla(\nabla\cdot J)$3, this yields

$\nabla(\nabla\cdot J)$4

The divergence structure is handled separately through

$\nabla(\nabla\cdot J)$5

together with compatible initial data. The paper is therefore relevant to solenoidal Ohm’s law only indirectly: it keeps divergence compatibility as a Maxwell constraint, not as a distinct constitutive reformulation (Hauser, 2023).

7. Conceptual boundaries and recurrent ambiguities

A persistent ambiguity is that many papers relevant to the subject do not explicitly use the term “solenoidal Ohm’s law.” In axisymmetric dynamo theory, for example, a modified Ohm law

$\nabla(\nabla\cdot J)$6

is introduced inside the usual divergence-free magnetic-flux framework

$\nabla(\nabla\cdot J)$7

but the novelty lies in the added restoring/friction force, not in a divergence-projected constitutive law (Sato et al., 2024).

The same caution applies in tokamak turbulence and collisionless reconnection. The turbulent tokamak current-balance equation is a reduced parallel Ohm law with hyper-resistive, anomalous-resistive, and cross-resistive terms, but it is not presented as a full vectorial solenoidal law (Chavdarovski et al., 2017). The reconnection review likewise does not define a separate solenoidal Ohm’s law; its closest analogue is the spatially uniform out-of-plane reconnection electric field in steady two-dimensional geometry, sustained locally by pressure-tensor and inertia terms in generalized Ohm’s law (Liu et al., 2024).

The resulting usage is therefore genuinely plural. In one class of papers, solenoidal Ohm’s law means an incompressible projected constitutive relation for $\nabla(\nabla\cdot J)$8. In another, it denotes the exact survival of local Ohmic behavior in the divergence-free transverse sector of a nonlocal theory. In a third, it is simply the steady-current conservation form obtained from $\nabla(\nabla\cdot J)$9. This suggests that the phrase is best read as a structural description of how Ohm’s law interacts with divergence constraints, rather than as a single universal equation.