Field-Based Macroscopic QED Overview

• Field-based macroscopic QED is a framework that quantizes electromagnetic fields via macroscopic constitutive responses rather than explicit many-body Hamiltonians.
• It employs Green tensors and operator-valued noise currents to incorporate dispersive, lossy, and inhomogeneous media while preserving causality and fluctuation-dissipation relations.
• The framework extends to nonlocal, time-varying, and moving media, enabling the analysis of phenomena such as quantum friction, vacuum polarization, and gain media dynamics.

Searching arXiv for recent and foundational papers on macroscopic QED and related field-based formulations. arXiv search query: "macroscopic quantum electrodynamics field-based Green tensor noise currents duality time-varying media" Field-based macroscopic quantum electrodynamics denotes a family of formulations in which electromagnetic fields are treated as the primary objects, while macroscopic environments enter through constitutive response functions, Green operators, or effective actions rather than through explicit microscopic many-body matter Hamiltonians. In its standard media-based form, macroscopic quantum electrodynamics (MQED) is a field-based quantization framework for electromagnetic radiation in realistic media: inhomogeneous, dispersive, and lossy macroscopic bodies are incorporated directly through their classical response functions, while the quantum nature of the field enters through operator-valued noise sources (Oue, 25 Feb 2026). In closely related external-field formulations, the “medium” can instead be the polarized quantum vacuum itself, as in Heisenberg–Euler–Schwinger electrodynamics for magnetar-strength backgrounds (Kim et al., 2021). Across these variants, the subject is organized by constitutive response, Green propagation, and operator field structure, and thereby differs from microscopic bound-state QED, cavity-mode truncations, and lattice gauge-theory simulations of vacuum QED (Feist et al., 2020).

1. Definition and conceptual domain

Field-based MQED begins from macroscopic Maxwell theory and asks how to quantize the electromagnetic field in arbitrary linear environments without explicitly quantizing every microscopic constituent. In the nanophotonic formulation, the environment is divided into a macroscopic background structure described by $\epsilon(\mathbf r,\omega)$ and $\mu(\mathbf r,\omega)$, together with a small number of microscopic emitters treated explicitly; the medium-assisted field is then represented by bosonic operators $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$, $\lambda\in\{e,m\}$, and all geometry enters through the dyadic Green tensor $\mathbf G(\mathbf r,\mathbf r',\omega)$ (Feist et al., 2020). In the tutorial treatment of passive and active media, the same point is expressed more directly: MQED quantizes the field in arbitrary linear media by using macroscopic constitutive response functions plus operator-valued noise currents, with the field correlation function as the unifying object behind both radiative and mechanical observables (Oue, 25 Feb 2026).

This domain is broader than a single formalism. One branch emphasizes Green-tensor quantization in dispersive and absorptive media (Feist et al., 2020, Horsley, 2022). A second branch generalizes that program to non-local, bianisotropic, and Onsager-violating media by formulating the theory in terms of a conductivity kernel $\mathbf Q(\mathbf r,\mathbf r',\omega)$ rather than only local $\boldsymbol\varepsilon$ and $\boldsymbol\mu$ (Buhmann et al., 2011). A third branch recasts MQED in terms of other field variables, notably the polarization density operator of a finite dielectric object (Forestiere et al., 2022) or the first-order dual electromagnetic field $\mathcal E=[\mathbf E,Z_0\mathbf H]^T$ (Agarwal et al., 29 Mar 2026). A fourth, closely allied branch treats the quantum vacuum itself as a macroscopic nonlinear medium in a prescribed classical background field, deriving effective constitutive response from the one-loop QED effective action rather than from a material susceptibility (Kim et al., 2021).

A recurring misconception is that MQED is merely a mode expansion with loss added phenomenologically. The literature summarized here takes the opposite position. The reservoir or noise sector is not auxiliary bookkeeping but the mechanism by which dissipation, causality, fluctuation relations, and canonical commutators become compatible within a field theory of realistic media (Horsley, 2022, Oue, 25 Feb 2026). This suggests that “field-based” is not a stylistic label but the defining structural feature of the subject.

2. Stationary linear media: Green tensors, reservoirs, and field operators

In passive media, the central operator Maxwell equation can be written as

$$\nabla \times \nabla \times \hat{\mathbf E}^{-}(\mathbf r,\omega) -\frac{\omega^2}{c^2}\,\varepsilon(\mathbf r,\omega)\hat{\mathbf E}^{-}(\mathbf r,\omega) = i\omega\mu_0\,\hat{\mathbf I}_N(\mathbf r,\omega),$$

with $\mu(\mathbf r,\omega)$0 the Langevin noise current and $\mu(\mathbf r,\omega)$1 the Green tensor of the corresponding classical Helmholtz operator. The field then has the exact representation

$\mu(\mathbf r,\omega)$2

and similarly for $\mu(\mathbf r,\omega)$3 (Oue, 25 Feb 2026). In the more general media-assisted formulation, the full field operator is expanded as

$\mu(\mathbf r,\omega)$4

with electric and magnetic reservoir sectors coupled to $\mu(\mathbf r,\omega)$5 and $\mu(\mathbf r,\omega)$6, respectively (Feist et al., 2020).

The reservoir construction is fixed by linear response. In the pedagogical Hamiltonian treatment, a continuum of harmonic oscillator fields $\mu(\mathbf r,\omega)$7 is coupled to the electromagnetic field with strength

$\mu(\mathbf r,\omega)$8

so that eliminating the reservoir reproduces the causal constitutive law and Kramers–Kronig-compatible permittivity (Horsley, 2022). In the Green-tensor formulation for passive dielectrics, the same logic appears as the operator rule

$\mu(\mathbf r,\omega)$9

whose correlations satisfy the fluctuation-dissipation relation and determine the field correlator (Oue, 25 Feb 2026).

The most compact statement of stationary MQED is that the field correlation is governed by the imaginary part of the Green tensor: $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$0 up to the normalization convention employed in the source. This relation controls spontaneous emission, Lamb shifts, Purcell factors, and stress-tensor observables in a single stroke (Oue, 25 Feb 2026). A plausible implication is that the Green tensor, rather than any particular mode basis, is the natural invariant of field-based MQED.

3. General constitutive response, duality, and operator symmetry

The standard local, reciprocal dielectric model is not the endpoint of the theory. For the most general linear absorbing media, the induced current is written in non-local form as

$\hat{\mathbf f}_\lambda(\mathbf r,\omega)$1

with $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$2 an arbitrary conductivity tensor that may be anisotropic, bianisotropic, spatially non-local, and Onsager-violating (Buhmann et al., 2011). In this setting, ordinary transpose symmetries of $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$3 and $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$4 fail, and the appropriate generalized real and imaginary parts are

$\hat{\mathbf f}_\lambda(\mathbf r,\omega)$5

$\hat{\mathbf f}_\lambda(\mathbf r,\omega)$6

With these replacements, the theory reproduces Maxwell’s equations, the fluctuation–dissipation theorem, and the canonical equal-time commutators even in non-reciprocal media (Buhmann et al., 2011).

Duality symmetry becomes nontrivial precisely in this generalized setting. In local bianisotropic media, the constitutive relations are written in terms of $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$7, $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$8, and magnetoelectric tensors $\hat{\mathbf f}_\lambda(\mathbf r,\omega)$9, and continuous electric–magnetic duality survives only when the medium class is broad enough to include non-reciprocal responses (Buhmann et al., 2011). A related but sharper issue appears for microscopic magnetic interactions in magnetodielectrics: the usual coupling $\lambda\in\{e,m\}$0 breaks Heaviside–Larmor duality, whereas local-field treatment restores dual symmetry by coupling the microscopic dipole to the local magnetic field, with magnetic local-field factor $\lambda\in\{e,m\}$1 mirroring the electric Clausius–Mossotti factor $\lambda\in\{e,m\}$2 (Westerberg et al., 2021). In that formulation, macroscopic QED can be made dual symmetric at the operator level rather than only after taking expectation values (Westerberg et al., 2021).

A further generalization replaces the second-order electric-field Helmholtz equation by a first-order Maxwell operator acting on the dual six-component field

$\lambda\in\{e,m\}$3

with $\lambda\in\{e,m\}$4 and $\lambda\in\{e,m\}$5 (Agarwal et al., 29 Mar 2026). The resulting first-order Green operator $\lambda\in\{e,m\}$6 yields propagation, Lorentz reciprocity, and a generalized optical theorem in a form that retains both bulk dissipation and boundary flux. This suggests that operator symmetry in MQED is not exhausted by reciprocity inside the volume; it also has a boundary-theoretic content.

4. Finite objects, open boundaries, and exact reduced descriptions

Finite bodies in vacuum expose a limitation of the standard source-only Langevin picture. For a finite lossy object surrounded by vacuum, the electric field generally decomposes as

$\lambda\in\{e,m\}$7

where $\lambda\in\{e,m\}$8 is the medium-assisted field generated by fluctuating material currents and $\lambda\in\{e,m\}$9 is the homogeneous scattering field associated with radiation incident from infinity (Ciattoni, 2024). The modified Langevin noise formalism proves that both sectors are required for arbitrary finite, linear, inhomogeneous magnetodielectric samples in vacuum, and that the Hamiltonian diagonalizes into two independent bosonic sectors: medium-assisted polaritons and scattering polaritons (Ciattoni, 2024). A common workaround in older formulations was to keep an infinitesimal vacuum loss and take $\mathbf G(\mathbf r,\mathbf r',\omega)$0 only at the end; MLNF removes the need for that limiting prescription (Ciattoni, 2024).

For finite dispersive dielectric objects, an alternative integral formulation takes the polarization density operator $\mathbf G(\mathbf r,\mathbf r',\omega)$1 as the central unknown. The electromagnetic field operators are reconstructed from $\mathbf G(\mathbf r,\mathbf r',\omega)$2 through the vacuum retarded Green function, and the coupled light–matter dynamics reduce to the operator-valued volume integral equation

$\mathbf G(\mathbf r,\mathbf r',\omega)$3

This is a retarded scattering equation over the object volume $\mathbf G(\mathbf r,\mathbf r',\omega)$4, with explicit boundary contributions on $\mathbf G(\mathbf r,\mathbf r',\omega)$5, and it enables direct use of classical computational electromagnetics methods for quantum calculations in open absorbing systems (Forestiere et al., 2022).

In multi-emitter nanophotonics, the same field-based logic can be compressed without approximation. For $\mathbf G(\mathbf r,\mathbf r',\omega)$6 emitters, the full reservoir continuum can be re-expanded into at most $\mathbf G(\mathbf r,\mathbf r',\omega)$7 orthonormal emitter-centered continua $\mathbf G(\mathbf r,\mathbf r',\omega)$8 plus dark modes that do not couple to the emitters (Feist et al., 2020). The nonorthogonal bright modes are defined from the Green-tensor coupling profiles sampled at the emitter positions, with overlap matrix

$\mathbf G(\mathbf r,\mathbf r',\omega)$9

and after orthogonalization the Hamiltonian becomes an exact minimal field basis for the emitter dynamics (Feist et al., 2020). This is one of the clearest demonstrations that field-based MQED need not mean computational intractability.

The first-order Maxwell operator formalism extends the open-system perspective further by keeping the boundary term explicitly: $\mathbf Q(\mathbf r,\mathbf r',\omega)$0 After quantization, the interior field is driven by two independent noise channels: bulk Langevin polarization noise and boundary input-output field operators. The resulting exact commutator closes as

$\mathbf Q(\mathbf r,\mathbf r',\omega)$1

which remains valid even when dielectrics extend to the boundary (Agarwal et al., 29 Mar 2026). This directly addresses waveguide and open-scattering problems that sit awkwardly inside electric-field-only formulations.

5. Nonequilibrium extensions: motion, time variation, and gain

Stationary constitutive response is not the only regime accessible to field-based MQED. In moving media, the same Hamiltonian field theory leads directly to quantum electromagnetic forces. For a uniformly moving dielectric, the constitutive response becomes Doppler shifted; in the one-dimensional treatment,

$\mathbf Q(\mathbf r,\mathbf r',\omega)$2

and the framework reproduces Pendry’s expression for quantum friction between sliding plates (Horsley, 2022). The field-theoretic significance is that quantum stress, field momentum, and center-of-mass motion arise from the same Hamiltonian structure rather than being appended afterwards (Horsley, 2022).

Time-varying dispersive media require a more radical modification. For a time-dependent Drude medium, the naïve substitution $\mathbf Q(\mathbf r,\mathbf r',\omega)$3 inside standard MQED leads to nonphysical polarization currents that become singular for a step change in carrier density (Horsley et al., 2024). The remedy is to keep the field–reservoir coupling fixed and instead modify the reservoir dynamics themselves. In the resulting theory, the homogeneous reservoir solution acquires temporally reflected components, the effective permittivity remains continuous across abrupt modulation, and the non-equilibrium noise-current correlator develops new $\mathbf Q(\mathbf r,\mathbf r',\omega)$4 structures absent from stationary fluctuation theory (Horsley et al., 2024). A plausible implication is that temporal reflection in MQED is encoded not only in wave propagation but in the operator statistics of the material noise sector itself.

Active media modify the same structure in a different direction. In gain-inclusive MQED, the field expansion remains

$\mathbf Q(\mathbf r,\mathbf r',\omega)$5

but gain channels must enter the noise current through creation-like rather than annihilation-like operators: $\mathbf Q(\mathbf r,\mathbf r',\omega)$6 This generates a nonzero reverse-ordered correlator proportional to $\mathbf Q(\mathbf r,\mathbf r',\omega)$7, so the vacuum in active media is not the passive zero-temperature vacuum (Oue, 25 Feb 2026). The theory remains consistent only if all poles of the Green tensor stay in the lower half-plane; otherwise the linear MQED description breaks down (Oue, 25 Feb 2026).

Taken together, motion, temporal modulation, and gain show that MQED is not intrinsically an equilibrium theory. What persists across these extensions is the field-based architecture: constitutive response, Green propagation, and operator noise sources remain the organizing principles (Horsley, 2022, Horsley et al., 2024, Oue, 25 Feb 2026).

6. Vacuum as an effective medium and analog realizations

A broader usage of the field-based perspective treats the vacuum itself as a macroscopic medium induced by a prescribed classical background field. In magnetar electrodynamics, the electromagnetic background is not quantized; instead one computes the quantum response of the Dirac sea through the one-loop effective action of spinor QED in a slowly varying external field (Kim et al., 2021). The relevant regime is

$\mathbf Q(\mathbf r,\mathbf r',\omega)$8

and the effective Lagrangian $\mathbf Q(\mathbf r,\mathbf r',\omega)$9 encodes both vacuum polarization and, when an electric component is present, pair creation (Kim et al., 2021). The vacuum is then treated as an effective nonlinear medium with constitutive response

$\boldsymbol\varepsilon$0

or in tensor form $\boldsymbol\varepsilon$1, $\boldsymbol\varepsilon$2 (Kim et al., 2021). In the weak-field limit this yields the standard birefringent refractive indices

$\boldsymbol\varepsilon$3

while the full treatment extends to $\boldsymbol\varepsilon$4, the regime relevant for magnetars (Kim et al., 2021).

This is field-based macroscopic QED in a literal sense: the chain

$\boldsymbol\varepsilon$5

treats the polarized quantum vacuum as the “medium” (Kim et al., 2021). The same perspective underlies analog constructions. In the linearized QED vacuum under a magnetostatic bias $\boldsymbol\varepsilon$6, the effective constitutive dyadics are uniaxial,

$\boldsymbol\varepsilon$7

with an analogous formula for $\boldsymbol\varepsilon$8, where

$\boldsymbol\varepsilon$9

(Mackay et al., 2011). Because the anisotropy is extremely small, an affine spatial transformation is introduced, after which inverse Bruggeman homogenization specifies a homogenized composite material electromagnetically equivalent to the transformed QED vacuum (Mackay et al., 2011). This does not quantize the field in the analog material, but it demonstrates that QED-vacuum constitutive response can be reformulated as a macroscopic-medium problem.

7. Relation to adjacent QED frameworks

Field-based MQED overlaps with, but is not identical to, several neighboring programs. In dipolar cavity QED, a single cavity mode is kept explicitly and coupled to collective matter polarization through a dipole-gauge Hamiltonian

$\boldsymbol\mu$0

with the $\boldsymbol\mu$1 term essential for gauge consistency and the correct electrostatic limit (Schuler et al., 2020). This is field-based in the sense that the electromagnetic mode remains an explicit quantum degree of freedom, but it is not macroscopic QED proper: there is no dispersive absorbing continuum, no Green-tensor quantization, and no operator-valued noise currents (Schuler et al., 2020).

Nonperturbative waveguide QED likewise remains conceptually adjacent rather than equivalent. There the starting point is a full minimal-coupling Hamiltonian

$\boldsymbol\mu$2

followed by an asymptotic decoupling unitary transformation that produces effective mass renormalization, potential reshaping, and many-body bound states in structured photonic continua (Ashida et al., 2021). The structured reservoir viewpoint is close to MQED, but the formulation is mode-based rather than Green-tensor-based and does not incorporate absorbing media through macroscopic constitutive functions (Ashida et al., 2021).

Hamiltonian lattice QED on quantum hardware is further away. In that setting, compact gauge links $\boldsymbol\mu$3, electric flux operators $\boldsymbol\mu$4, and plaquettes $\boldsymbol\mu$5 define a gauge-invariant lattice theory in temporal gauge, with finite-dimensional $\boldsymbol\mu$6 regularization approaching compact $\boldsymbol\mu$7 as $\boldsymbol\mu$8 (Mou et al., 31 Oct 2025). This is a field-based gauge theory of the vacuum electromagnetic field, but it does not contain the defining ingredients of macroscopic QED in media: no $\boldsymbol\mu$9, no $\mathcal E=[\mathbf E,Z_0\mathbf H]^T$0, no noise polarization or Langevin currents, and no Green-function quantization in structured matter (Mou et al., 31 Oct 2025).

The boundary of the field is therefore conceptually sharp. Macroscopic QED in the strict sense quantizes electromagnetic fields in prescribed macroscopic environments—material or effective—through constitutive response, Green propagation, and fluctuation operators. Neighboring theories may retain field degrees of freedom and illuminate ultrastrong coupling, structured reservoirs, or gauge constraints, but they do not replace MQED’s role as the quantum field theory of realistic macroscopic electrodynamic environments (Schuler et al., 2020, Ashida et al., 2021, Mou et al., 31 Oct 2025).