---
title: Self-Consistent Maxwell–Pauli Framework
url: https://www.emergentmind.com/topics/self-consistent-maxwell-pauli-framework
type: topic
---

# Self-Consistent Maxwell–Pauli Framework

The self-consistent Maxwell–Pauli framework denotes a family of coupled matter–field models in which a two-component Pauli spinor, or Pauli–Kohn–Sham spinors, evolves under minimal electromagnetic coupling while the electromagnetic potentials satisfy Maxwell or Poisson-type equations sourced by charge and current densities constructed from the same spinorial state. Across the literature, the framework appears in one-body travelling-wave problems, semi-relativistic \(1/c\) and \(1/c^2\) expansions of Dirac–Maxwell theory, many-body classical-field models, and ab initio real-space light–matter simulations beyond the dipole approximation [1402.3936] [1403.3854] [2005.06556] [2409.08959]. This suggests that the expression identifies a common self-consistency principle—mutual matter–field back-reaction—rather than a single unique equation set.

## 1. Canonical coupled equations and gauge structure

In the one-body Coulomb-gauge formulation, the matter variable is a Pauli spinor \(\psi=(\psi_1,\psi_2)^T:\mathbb R^3\times\mathbb R\to\mathbb C^2\), with charge \(Q\) and mass \(m\), coupled to electromagnetic potentials \(\phi:\mathbb R^3\times\mathbb R\to\mathbb R\) and \(A:\mathbb R^3\times\mathbb R\to\mathbb R^3\). The fields are
\[
E=-\nabla\phi-\frac1c\partial_tA,\qquad B=\nabla\times A,
\]
and Coulomb gauge imposes
\[
\nabla\cdot A=0,\qquad -\Delta\phi=4\pi Q|\psi|^2.
\]
The Pauli equation is
\[
i\hbar\,\partial_t\psi=H_P(A,\phi)\psi,
\]
with
\[
H_P(A,\phi)=\frac1{2m}\bigl(\sigma\cdot(-i\hbar\nabla+(Q/c)A)\bigr)^2+Q\phi
=\frac1{2m}(-i\hbar\nabla+(Q/c)A)^2-\frac{\hbar Q}{2c}\sigma\cdot B+Q\phi.
\]
The corresponding current density is
\[
J_P[\psi,A](x,t)=-\frac{Q}{m}\,\mathrm{Re}\,\langle \psi,\sigma(\sigma\cdot(-i\hbar\nabla+(Q/c)A))\psi\rangle,
\]
and the transverse vector potential satisfies Maxwell’s equation for the radiation field in Coulomb gauge [1402.3936].

A distinct but closely related second-order weakly relativistic formulation keeps \(\phi\) and \(A\) as variational variables and starts from the Lagrangian density
\[
\mathcal L=\psi^\dagger\bigl(i\hbar\,\partial_t-H_{P2}[\phi,A]\bigr)\psi+\frac1{8\pi}(E^2-B^2),
\]
where \(H_{P2}\) is the Schrödinger–Pauli Hamiltonian through \(O(1/c^2)\). Varying the action with respect to \(\phi\) and \(A\) yields
\[
\nabla\cdot E=4\pi\rho,\qquad
\nabla\times B-\frac1c\partial_tE=\frac{4\pi}{c}J,
\]
together with the homogeneous Maxwell equations. In that formulation, the closed theory is stated in terms of the six functions \(\{\phi,A,\rho,J\}\), with gauge freedom removing one degree of freedom in \(\{\phi,A\}\) and the continuity equation removing one in \(\{\rho,J\}\), leaving exactly six independent field components [2407.09570].

Gauge fixing is not unique across the literature. Coulomb gauge is central in the travelling-wave existence theory and in full minimal-coupling Maxwell–TDDFT [1402.3936] [2409.08959], whereas Lorentz gauge is used in one semi-relativistic Lagrangian derivation and then expanded consistently to \(O(1/c^2)\) [1403.3854]. In the first-order semi-relativistic Pauli–Poisswell model, Coulomb gauge is again standard, although Lorenz gauge may also be adopted before dropping \(O(c^{-2})\) terms [2304.06660].

## 2. Source structure: charge, current, polarization, and magnetization

A defining feature of the self-consistent Maxwell–Pauli framework is that the sources are not inserted externally but derived from the same matter fields that evolve under the electromagnetic potentials. At second order in the weakly relativistic expansion, one finds
\[
\rho(\mathbf r,t)=e\,\psi^\dagger\psi-\nabla\cdot\bigl[P_{so}+P_D\bigr],
\]
\[
J(\mathbf r,t)=\frac{e}{2m}\bigl[\psi^\dagger\pi\psi+(\pi\psi)^\dagger\psi\bigr]
+c\,\nabla\times M_{spin}-eF+\partial_t\bigl[P_{so}+P_D\bigr],
\]
with
\[
P_{so}=-\frac{e\hbar}{4m^2c^2}\psi^\dagger(\sigma\times\pi)\psi,\qquad
P_D=-\frac{e\hbar^2}{8m^2c^2}\nabla(\psi^\dagger\psi),
\]
\[
M_{spin}=\frac{e\hbar}{2mc}\psi^\dagger\sigma\psi,\qquad
F=\frac{\hbar^2}{8m^2c^2}\nabla(\psi^\dagger\psi).
\]
Defining the total electric polarization by \(P=P_{so}+P_D\) and the total magnetization \(M\) so that
\[
J=\partial_tP+c\,\nabla\times M,
\]
one obtains the microscopic constitutive relations
\[
\nabla\cdot P=-\rho,\qquad \partial_tP+c\,\nabla\times M=J.
\]
These relations are presented as microscopic analogues of the usual macroscopic constitutive relations [2407.09570].

An equivalent free-plus-bound decomposition appears in the semi-relativistic Lagrangian treatment,
\[
\rho=\rho_{\mathrm{free}}-\nabla\cdot P,\qquad
J=J_{\mathrm{free}}+\nabla\times M+\partial_tP,
\]
where
\[
\rho_{\mathrm{free}}=\Psi^\dagger\Psi,\qquad
J_{\mathrm{free}}=\frac{-i\hbar}{2m}\bigl[\Psi^\dagger\nabla\Psi-(\nabla\Psi^\dagger)\Psi\bigr]-\frac{q}{m}A\,\Psi^\dagger\Psi,
\]
\[
M=M_{\mathrm{spin}}=\frac{q\hbar}{2m}\Psi^\dagger\sigma\Psi,
\]
\[
P=P_{\mathrm{spin}}+P_{\mathrm{Darwin}},
\qquad
P_{\mathrm{spin}}=-\frac{q\hbar}{4m^2c^2}\Psi^\dagger\sigma\Psi\times A,
\qquad
P_{\mathrm{Darwin}}=-\frac{q\hbar^2}{8m^2c^2}\nabla(\Psi^\dagger\Psi).
\]
In this presentation, the Darwin correction is interpreted as smearing due to Zitterbewegung, while the spin-polarization term is interpreted as the Lorentz boost of the rest-frame magnetization [1403.3854].

| Quantity | Expression | Role |
|---|---|---|
| \(P_{so}\) | \(-\frac{e\hbar}{4m^2c^2}\psi^\dagger(\sigma\times\pi)\psi\) | spin–orbit polarization |
| \(P_D\) | \(-\frac{e\hbar^2}{8m^2c^2}\nabla(\psi^\dagger\psi)\) | Darwin polarization |
| \(M_{spin}\) | \(\frac{e\hbar}{2mc}\psi^\dagger\sigma\psi\) | spin magnetization |

A recurrent point of controversy concerns “spin current.” In the second-order variational formulation, the only sources of Maxwell’s equations are \(\rho\) and \(J\) as defined above. One may define a spin density \(S(\mathbf r,t)=\frac{\hbar}{2}\psi^\dagger\sigma\psi\) and various candidate spin-current tensors \(J_s\), but none arise from \(\delta\mathcal L/\delta A\) or \(\delta\mathcal L/\delta\phi\), and none enter Maxwell’s equations. The local spin-continuity relation has the form
\[
\partial_tS+\nabla\cdot J_s=\tau_{so},
\]
with a spin-torque term, and is therefore not a true conservation law in the same sense as charge continuity [2407.09570].

## 3. Variational travelling waves and solitary structures

A mathematically sharp realization of the self-consistent Maxwell–Pauli framework is given by the travelling-wave problem for the one-body Maxwell–Pauli system. The monochromatic travelling-wave ansatz is
\[
\psi(x,t)=e^{-i\omega t}\,\tilde\psi(x-vt),\qquad
A(x,t)=\tilde A(x-vt),
\]
with constant velocity \(v\) and unknown frequency \(\omega\in\mathbb R\). After setting \(y=x-vt\), the coupled PDE reduces to a stationary system on \(\mathbb R^3\):
\[
-\hbar(\theta+i\,v\cdot\nabla)\tilde\psi
=\frac1{2m}\bigl(\sigma\cdot(-i\hbar\nabla+(Q/c)\tilde A)\bigr)^2\tilde\psi,
\]
\[
\Bigl(\frac{(v\cdot\nabla)^2}{c^2}-\Delta\Bigr)\tilde A=\frac{4\pi}{c}\,P\,J_P[\tilde\psi,\tilde A],
\qquad \nabla\cdot\tilde A=0,
\]
where \(P=1-\nabla\,\mathrm{div}\,\Delta^{-1}\) is the Helmholtz projector onto divergence-free fields, and
\[
\theta=\hbar^{-1}\bigl(\mathrm{EM\mbox{-}energy}(\tilde A,(v\cdot\nabla)\tilde A)-\hbar\omega\bigr).
\]
The electromagnetic energy entering this reduction is
\[
\mathrm{EM\mbox{-}energy}(A,E)=\frac1{8\pi}\int\bigl(|\nabla\times A|^2+|E|^2\bigr)\,dx,
\qquad E=-(v\cdot\nabla)A.
\]

The stationary problem is recast variationally by fixing the \(L^2\)-mass \(\lambda=\|\tilde\psi\|_{L^2}^2>0\) and considering
\[
S_\lambda=\{(\psi,A)\in H^1(\mathbb R^3;\mathbb C^2)\times D^1(\mathbb R^3;\mathbb R^3)\mid \|\psi\|_{L^2}^2=\lambda,\ \nabla\cdot A=0\}.
\]
For the Pauli case, the functional is
\[
\mathcal E^v_P(\psi,A)
=
\frac1{2m}\left\|(\sigma\cdot(-i\hbar\nabla+(Q/c)A))\psi+\frac{mc}{Q}v\psi\right\|_{L^2}^2
-\frac{Q}{c}(\psi,v\cdot A\,\psi)_{L^2}
-\frac{m|v|^2}{2}\lambda
+\frac1{8\pi}\Bigl(\|\nabla\otimes A\|_{L^2}^2-\|(v\cdot\nabla)A/c\|_{L^2}^2\Bigr).
\]
Minimizers under the mass constraint solve the reduced travelling-wave system with a Lagrange multiplier \(\theta\) [1402.3936].

The existence theorem is formulated in terms of the critical speeds
\[
\Theta_{P,\pm}^\lambda=
-\frac{8\pi K_S^3Q^2\lambda}{\hbar}
\pm
\sqrt{\left(\frac{8\pi K_S^3Q^2\lambda}{\hbar}\right)^2+c^2},
\]
where \(K_S\) is the sharp Sobolev constant in \(\|f\|_{L^6}\le K_S\|\nabla f\|_{L^2}\). For every \(v\in\mathbb R^3\) with \(0<|v|<\Theta_{P,+}^\lambda\), there exists \((\omega,\psi,A)\) with \(\|\psi\|^2=\lambda\) such that
\[
\psi(x,t)=e^{-i\omega t}\psi(x-vt),\qquad A(x,t)=A(x-vt)
\]
solves the Maxwell–Pauli system in Coulomb gauge, with \((\psi,A)\in H^1\times D^1\). The proof proceeds by showing that \(\mathcal E^v_j\) is bounded below on \(S_\lambda\) below the critical speed, constructing a minimizing sequence, using concentration–compactness to prevent loss of mass at infinity, translating the sequence so that \(\psi_n\) is tight in \(L^2\), extracting a weakly convergent subsequence in \(H^1\times D^1\), and then using lower semicontinuity to recover a minimizer [1402.3936].

For small velocity, the energy satisfies
\[
E_j(v)=\mathcal E^v_j(\psi,A)+\frac{m|v|^2}{2}\lambda
=\frac12\,m|v|^2\lambda+O(|v|^3),
\]
so the leading-order term is exactly the classical kinetic energy of a particle of mass \(m\). The same analysis applies to the Maxwell–Schrödinger case by dropping the \(\sigma\cdot B\) term. No uniqueness of the minimizer is proved, and stability of the travelling waves is not addressed [1402.3936].

## 4. Semi-relativistic expansions and coherent light-induced couplings

Several self-consistent Maxwell–Pauli models are derived as controlled approximations to the Dirac–Maxwell system. In one \(O(1/c^2)\) derivation, the matter Lagrangian contains the Zeeman, Darwin, and spin–orbit terms, and variation with respect to \(\Psi^\dagger\) yields the extended Pauli equation
\[
i\hbar\partial_t\Psi
=
\Bigl[
mc^2+q\phi+\frac{(-i\hbar\nabla-qA)^2}{2m}
+\frac{q\hbar}{2m}\sigma\cdot B
-\frac{q\hbar^2}{8m^2c^2}\nabla\cdot E
-\frac{q\hbar}{4m^2c^2}\sigma\cdot\bigl(E\times(-i\hbar\nabla-qA)+(-i\hbar\nabla-qA)\times E\bigr)
\Bigr]\Psi.
\]
Maxwell’s equations are expanded consistently to the same order. Writing \(\phi=\phi_0+c^{-2}\phi_2+\cdots\) and \(A=c^{-2}A_2+\cdots\), one finds
\[
\nabla^2\phi_0=-\rho_0/\varepsilon_0,\qquad
\nabla^2\phi_2=-\rho_2/\varepsilon_0+\varepsilon_0^{-1}\partial_t^2\phi_0,\qquad
\nabla^2A_2=-\mu_0J_0,
\]
so that there is no magnetic field at leading order, corresponding to the electric limit [1403.3854].

A related semi-relativistic program starts from the Dirac equation with external and internal fields,
\[
i\hbar\,\partial_t\Psi
=
\Big[c\,\boldsymbol\alpha\cdot(\mathbf p-q\mathbf A_{\rm ext}-q\mathbf A_{\rm int})
+mc^2\beta
+q\Phi_{\rm ext}+q\Phi_{\rm int}\Big]\Psi,
\]
performs a Foldy–Wouthuysen expansion to \(\mathcal O(1/c^2)\), and obtains a two-component Pauli Hamiltonian containing minimal coupling, Zeeman, Darwin, and spin–orbit terms built from \(\mathbf E_{\rm ext}+\mathbf E_{\rm int}\) and \(\mathbf A_{\rm ext}+\mathbf A_{\rm int}\). The charge and current are expanded as
\[
\rho=\rho^{(0)}+\rho^{(2)}+\dots,\qquad
\mathbf j=\mathbf j^{(0)}+\mathbf j^{(2)}+\dots,
\]
with \(\rho^{(2)}\) split into orbital, spin, and field-induced pieces, while \(\mathbf A_{\rm int}^{(2)}\) is generated by \(\mathbf j_{\rm T}^{(0)}\). In this formulation, the internal potentials are
\[
\Phi_{\rm int}=\Phi_{\rm int}^{(0)}+\Phi_{\rm int}^{(2)},\qquad
\mathbf A_{\rm int}=\mathbf A_{\rm int}^{(2)}.
\]
The resulting Hamiltonian is decomposed into external-field Foldy–Wouthuysen terms, internal Breit–Pauli terms, and a coherent light-induced mean field \(H_P^{\rm int-ext}\) [1602.05623].

Within \(H_P^{\rm int-ext}\), four coherent spin–light–induced mechanisms are identified:

1. **Zeeman-like mechanism (A1)**: \(-\frac{q\hbar}{2m}\,\sigma\cdot[\nabla\wedge\mathbf A_{\rm field}^{(2)}]\), corresponding to \(\mathbf A_{\rm ext}\to \mathbf j_{\rm field}\to \mathbf A_{\rm field}^{(2)}\to\) an effective light-induced magnetic field.

2. **Spin–orbit-like mechanism (A2)**: \(-\frac{q^2\hbar}{4m^2c^2}\,\sigma\cdot(\nabla\Phi^{(0)}\wedge\mathbf A_{\rm ext})\), in which the external vector potential modifies orbital motion in the internal Coulomb field.

3. **“\(A\cdot A\)” spin–other–orbit mechanism (B1)**: \(\frac{q^2}{m}\,\mathbf A_{\rm ext}\cdot\mathbf A_{\rm spin}^{(2)}\), where the internal spin current produces a vector potential that couples parametrically to \(\mathbf A_{\rm ext}\).

4. **Darwin-type spin–orbit mechanism (B2)**: \(q\,\Phi_{\rm field}^{(2)}\), with
\[
\rho_{\rm field}^{(2)}
=
-\frac{q\hbar}{4m^2c^2}\nabla\cdot\bigl[(\phi^\dagger\boldsymbol\sigma\phi)\wedge \mathbf A_{\rm ext}\bigr],
\]
so that a light pulse and local spin density generate an induced second-order Hartree potential [1602.05623].

Taken together, these semi-relativistic derivations show that self-consistency extends beyond the bare Pauli Zeeman term to include Darwin corrections, spin–orbit structures, internal Breit–Pauli mean fields, and light-induced internal sources [1403.3854] [1602.05623].

## 5. Many-body dynamics, energetic stability, and asymptotic limits

The many-body Maxwell–Pauli equations model \(N\ge 1\) non-relativistic electrons with spin interacting with their self-generated classical electromagnetic field and \(K\ge 0\) static nuclei. In the atomic units of Kieffer’s thesis, and in Coulomb gauge \(\nabla\cdot A=0\), the \(2^N\)-component antisymmetric spinor \(\psi=\psi(t,\mathbf x_1,\dots,\mathbf x_N)\) satisfies
\[
i\,\partial_t\psi=H_{\rm P}(A)\psi,\qquad
H_{\rm P}(A)=\sum_{j=1}^N[\sigma_j\cdot(\mathbf p_j+A(\mathbf x_j))]^2+V(\mathbf R,Z),
\]
with the usual electron–electron, electron–nucleus, and nucleus–nucleus Coulomb terms in \(V(\mathbf R,Z)\). Maxwell’s equations are
\[
\nabla\cdot E=4\pi\rho,\qquad
\nabla\times E=-\alpha\,\partial_tB,
\]
\[
\nabla\cdot B=0,\qquad
\nabla\times B=4\pi\alpha\,J+\alpha\,\partial_tE,
\]
with
\[
E=-\nabla\phi-\alpha\,\partial_tA,\qquad B=\nabla\times A,
\]
and Coulomb gauge gives
\[
-\Delta\phi=-4\pi\alpha|\psi|^2,\qquad
\square A=4\pi\alpha\,\mathbb P[J_{\rm P}[\psi,A]],
\qquad \square=\alpha^2\partial_t^2-\Delta.
\]
The conserved total energy is
\[
E_{\rm MP}[\psi,A,\partial_tA]=T_{\rm P}[\psi,A]+V[\psi]+F[A,\partial_tA].
\]
No further renormalization is needed [2005.06556].

The existence theory is tied to energetic stability. The absolute ground-state energy \(E_{\rm MP}^{\rm G}(N,K,Z,\alpha)\) is bounded below for general \(N,K\) provided
\[
\alpha\le 0.06,\qquad \alpha^2\max_k Z_k\le 0.041,
\]
yielding energetic stability of second kind,
\[
E_{\rm MP}^{\rm G}(N,K,Z,\alpha)\ge -C(\alpha,Z)(N+K).
\]
For \(N=K=1\), there is a critical nuclear charge \(Z_c\) (\(\sim 10^4\)) such that
\[
E_{\rm MP}^{\rm G}(1,1,Z,\alpha)>-\infty\quad\Leftrightarrow\quad Z<Z_c.
\]
Assuming \(\alpha\) and \(\alpha^2\max Z_k\) are small enough so that \(E_{\rm MP}^{\rm G}>-\infty\), any initial data
\[
(\psi_0,a_0,\dot a_0)\in H^1(\mathbb R^{3N};\mathbb C^{2^N})\times H^1(\mathbb R^3;\mathbb R^3)\times L^2(\mathbb R^3;\mathbb R^3),
\]
with \(\|\psi_0\|_2=1\) and divergence-free initial vector potential and time derivative, generate a global-in-time, finite-energy, weak solution
\[
(\psi,A,\partial_tA)\in C_{\rm weak}\bigl(\mathbb R_+;H^1\oplus H^1\oplus L^2\bigr),
\]
which conserves both \(\|\psi\|_2=1\) and the total energy [2005.06556].

A lower-order semi-relativistic approximation is the Pauli–Poisswell system, described as the first order in \(1/c\) semi-relativistic approximation of the Dirac–Maxwell equation for 4-spinors coupled to self-consistent electromagnetic fields. It consists of
\[
i\hbar\,\partial_t\psi
=
\frac1{2m}(-i\hbar\nabla-A)^2\psi
+\phi\,\psi
-\frac{\hbar}{2mc}(\sigma\cdot B)\psi
+O(c^{-2}),
\]
coupled to the magnetostatic Poisson equations
\[
-\Delta\phi=4\pi\rho,\qquad
-\Delta A=4\pi J,\qquad
\nabla\cdot A=0,
\]
with
\[
\rho=|\psi|^2,\qquad
J=\mathrm{Re}\{\psi^*(-i\hbar\nabla-A)\psi\}+\nabla\times(\psi^*\sigma\psi).
\]
The Hamiltonian
\[
H[\psi,A,\phi]
=
\int_{\mathbb R^3}
\Bigl[
\frac1{2m}|(-i\hbar\nabla-A)\psi|^2
+\phi|\psi|^2
-\frac{\hbar}{2mc}\psi^*(\sigma\cdot B)\psi
\Bigr]dx
+\frac1{8\pi}\int_{\mathbb R^3}\bigl(|\nabla\phi|^2+|B|^2\bigr)\,dx
\]
is conserved. WKB analysis yields a mathematically rigorous semiclassical limit \(\hbar\to 0\) from the Pauli–Poisswell equation to the magnetic Euler–Poisswell equation, with local wellposedness for Euler–Poisswell that is global unless a finite time blow-up occurs [2304.06660].

## 6. Ab initio full minimal coupling and beyond-dipole simulation

A contemporary computational realization of the self-consistent Maxwell–Pauli idea is the full minimal-coupling Maxwell–TDDFT framework implemented in Octopus. It is described as the first ab initio, non-relativistic QED method that couples light and matter self-consistently beyond the electric dipole approximation and without multipolar truncations. The matter subsystem uses a Pauli–Kohn–Sham Hamiltonian
\[
\hat H_{\mathrm{f.m.c.}}
=
\sum_{j=1}^{N_{\mathrm{occ}}}
\Biggl[
\frac1{2m}\Bigl(-i\hbar\nabla_j+\frac{|e|}{c}A(r_j,t)\Bigr)^2
+\frac{|e|\hbar}{2mc}B(r_j,t)\cdot s_j
+V_H[n](r_j,t)+V_{xc}[n](r_j,t)+V_{nuc}(r_j)
\Biggr],
\]
with transverse \(A(r,t)\) in Coulomb gauge. The Maxwell subsystem evolves the transverse field according to
\[
\nabla\times\nabla\times A(r,t)-\frac1{c^2}\frac{\partial^2A(r,t)}{\partial t^2}=\mu_0J(r,t),
\]
or, equivalently, in the Riemann–Silberstein representation. The quantum current is split into
\[
J(r,t)=J_p(r,t)+J_d(r,t)+J_m(r,t),
\]
with paramagnetic, diamagnetic, and magnetization-current parts, and \(J_m=-\nabla\times M\) [2409.08959].

The multi-system cycle couples matter and Maxwell solvers on different real-space grids. Per global step \(\Delta t_e\), the code regrids \(A\) from the Maxwell grid to the matter grid, propagates the matter subsystem with a high-order exponential midpoint or Taylor propagator (4th order), evaluates \(n(r,t)\) and the total current \(J(r,t)\), regrids \(J\) onto the Maxwell grid, propagates the Maxwell subsystem for \(N_m\) steps of \(\Delta t_m=\Delta t_e/N_m\) with \(N_m\approx 20\text{–}30\), solves the Poisson equation for \(\nabla\times A=B\) and gauge-corrects to maintain \(\nabla\cdot A=0\), then regrids the updated \(A\) back to the matter grid. Matter grids typically use \(O(0.2\text{–}0.3\,\mathrm{bohr})\) spacing, while the Maxwell grid is coarser, \(0.4\text{–}0.6\,\mathrm{bohr}\). Perfectly matched layers absorb outgoing radiation, and Dirichlet conditions can introduce plane waves [2409.08959].

| System | Setup | Reported result |
|---|---|---|
| Cherenkov radiation of an electronic wavepacket | \(v=p/m>c/n\) in a refractive medium with \(n\approx 100\) | back-reaction changes spatially resolved density by \(\sim 3\%\) versus \(\sim 0.1\%\) in dipole-only back-coupling |
| Non-chiral benzene under XUV | \(\omega\approx 270\,\mathrm{eV}\), polarization and \(k\) both in-plane | magneto-optical spectrum has a distinct resonance around \(290\,\mathrm{eV}\), absent in pure dipole spectra |
| \(Na_{55}\) dimer | gap \(0.7\,\mathrm{nm}\); full Maxwell–matter back-coupling | plasmon resonance red-shifts by \(\sim 70\,\mathrm{meV}\); hot-spot phase shifts by \(\sim 0.06\pi\,\mathrm{rad}\) |

The framework is presented as origin-independent, with all multipole orders included non-perturbatively and self-consistently together with retardation and radiation reaction. It also retains exact treatment of paramagnetic and diamagnetic currents and spin–orbit (Pauli) coupling. The stated limitations are higher computational cost due to 3D Maxwell propagation, gauge issues with non-local pseudopotentials, and discretization constraints associated with Courant-type stability [2409.08959].

## 7. Alternative formulations, interpretive variants, and open directions

A distinct static-field usage of Maxwell–Pauli language appears in the Schrödinger–Pauli theory reformulated by the Quantal Newtonian first law,
\[
F_{\rm ext}(\mathbf r)+F_{\rm int}(\mathbf r)=0.
\]
Here the external field is the sum of electrostatic and Lorentz fields, while the internal field is the sum of electron–interaction, differential-density, kinetic, and internal magnetic fields, each defined from expectation values of Hermitian operators. Because the total field is a known functional of the wave function, the binding potential is written as
\[
v[\Psi](\mathbf r)=\int_\infty^{\mathbf r}\mathbf F(\mathbf r')\cdot d\mathbf r',
\]
leading to the generalized self-consistent Schrödinger–Pauli equation
\[
[\hat T+\hat W+\hat V[\Psi]]\Psi=E[\Psi]\Psi.
\]
This formulation further supports a local effective-potential mapping within quantal density-functional theory [1909.09692].

Another interpretive variant introduces an additional vector field \(\mathbf G(\mathbf r,t)\) into a unified action for \(\Psi\), the electromagnetic field, and a Klein–Gordon-type \(\mathbf G\)-field. Variation produces a Maxwell–Pauli–\(\mathbf G\) system in which the Pauli equation contains the nonlinear term \(-e\,\sigma\cdot\mathbf G\), and \(\mathbf G\) satisfies
\[
\frac1{c^2}\partial_t^2\mathbf G-\nabla^2\mathbf G+\mu^2\mathbf G
=
4\pi e\,\Psi^\dagger\sigma\Psi.
\]
Within that framework, spontaneous emission is said to emerge from the atom’s own dipole radiation field, and spin density \(\mathbf s=\Psi^\dagger\sigma\Psi\) together with magnetic moment density \(\mathbf m=-(e\hbar/2m_ec)\Psi^\dagger\sigma\Psi\) are given a classical field interpretation, with \(\boldsymbol\mu=-(e/m_e)\mathbf S\) and \(g_{\rm spin}=2\) [2203.09466].

A further construction treats \(N\) electrons and \(L\) photons through a single configuration-space wave function \(\Psi(t;q_e,q_{ph})\). In that approach, the expectation values of photon-field operators satisfy the inhomogeneous Maxwell equations with source terms given by the electron charge and current densities, and the algebraic structure of bosonic creation and annihilation operators is recovered in first-quantized guise without second quantization of the classical Maxwell field [2002.11106].

The framework also carries explicit unresolved issues. In the travelling-wave existence theory, no uniqueness of minimizers is proved and stability is not addressed [1402.3936]. In the many-body Maxwell–Pauli equations, open problems include local well-posedness without smallness assumptions, blow-up beyond the stability threshold, inclusion of full Coulomb self-interaction, semiclassical and non-relativistic limits, and extension to quantized fields [2005.06556]. Possible extensions identified in the travelling-wave and ab initio literatures include many-body mean-field limits, Maxwell–Dirac, pseudo-relativistic solitary waves, twisted light and orbital angular momentum beams, strong-field physics, inelastic light scattering, and cavity QED or polaritonic chemistry [1402.3936] [2409.08959]. At the level of source variables, the second-order variational theory explicitly rejects an additional “spin current” as an independent electromagnetic variable, retaining \(\{\phi,A;\rho,J\}\) as the minimal set [2407.09570].

Taken together, these developments define the self-consistent Maxwell–Pauli framework as a technically diverse but conceptually unified program: Pauli spinors or spinor Kohn–Sham orbitals generate electromagnetic sources, those sources determine fields through Maxwell or Poisson equations, and the resulting fields feed back into the matter Hamiltonian through minimal coupling, Zeeman, Darwin, spin–orbit, and, in semi-relativistic settings, higher-order internal-field terms.

Source: https://www.emergentmind.com/topics/self-consistent-maxwell-pauli-framework