Unified Effective Field Theory for Nonlinear and Quantum Optics

TL;DR

The paper presents a unified tensor-field theory for light–matter interactions, encoded in a covariant action $\mathcal{S}_{\rm U}$ that couples the electromagnetic field to a polarization multiplet with linear, nonlinear, and topological terms, plus a dissipative sector. Gauge symmetry is preserved via Dirac–BRST quantization within a real-time Keldysh formalism, enabling controlled one-loop renormalization of the Kerr nonlinearity $\chi^{(3)}$ and a ghost-free physical subspace. A tensor-network (MPO/MPS) solver maps the non-equilibrium dynamics of the unified theory to tractable simulations, bridging quantum optics, nonlinear propagation, and topology. The framework is validated across five platforms (GaAs, THz filamentation, silicon photonic lattices, ENZ ITO, and superconducting quarton circuits) through consistent predictions for $g^{(2)}(0)$, $n_2$, topological transitions, ENZ enhancements, and cross-Kerr couplings, highlighting a universal 0-D scaling for quartic nonlinearities. While currently restricted to 1-D geometries and sub-cutoff frequencies, the approach promises extensions to higher dimensions and stronger coupling regimes, enabling room-temperature quantum logic and multi-scale photonics.

Abstract

Predicting phenomena that mix few-photon quantum optics with strong field nonlinear optics is hindered by the use of separate theoretical formalisms for each regime. We close this gap with a unified effective field theory valid for frequencies lower than the material-dependent cutoff set by the band gap, plasma frequency, or similar scale. The action couples the electromagnetic gauge field to vector polarisation modes. An isotropic potential generates the optical susceptibilities, while a higher-dimension axion-like term captures magnetoelectric effects; quantisation on the Schwinger-Keldysh contour with doubled BRST ghosts preserves gauge symmetry in dissipative media. One-loop renormalisation-group equations reproduce the measured dispersion of the third-order susceptibility from terahertz to near-visible frequencies after matching a single datum per material. Real-time dynamics solved with a matrix-product-operator engine yield good agreement with published results for GaAs polariton cavities, epsilon-near-zero indium-tin-oxide films and superconducting quarton circuits. The current formulation is limited to these 1-D geometries and sub-cut-off frequencies; higher-dimensional or above-cut-off phenomena will require additional degrees of freedom or numerical methods.

Figures (5)

• Figure 1: Research workflow: from unified Lagrangian design, through quantization, renormalization, numerical simulation, to experiment.
• Figure 2: Leading one–loop correction to the Kerr vertex in the Keldysh formalism. The diagram represents the polarization self–energy (“bubble”) formed by one retarded and one advanced propagator, $G_{R}$ and $G_{A}$, linked by the Keldysh line $G_{K}$. External insertions $\delta P$ correspond to the operator $P_{q}P_{c}^{3}$. Sewing $q\!\leftrightarrow\!c$ indices across the two Kerr vertices twice yields the product $G_{R}(k)G_{A}(k)$, the only UV–divergent contraction contributing to the renormalization of $\chi^{(3)}$ at one loop; all other index patterns remain finite (see App. \ref{['app:index-census']} and App. \ref{['app:loop']}).
• Figure 3: Tensor–Keldysh numerical architecture used to simulate the unified action $\mathcal{S}_{\rm U}$.
• Figure 4: (a) Time evolution of the mean photon number $\langle n(t)\rangle$, showing relaxation from the initially pumped state toward the steady-state limit. (b) Second–order correlation $g^{(2)}(0;t)$ obtained from the same simulation, exhibiting a transient bunching peak $(g^{(2)}(0)>1)$ before returning to the coherent limit $g^{(2)}(0)\!\to\!1$. Both datasets are computed from the real–time tensor–network (MPO) simulation described in Sec. V.
• Figure 5: (a) Temporal evolution of the mean photon number $\langle n(t)\rangle$ for various damping rates $\eta$ and nonlinear coefficients $\chi_{3}$. (b) Corresponding dynamics of the second–order correlation function $g^{(2)}(0;t)$. Increasing $\eta$ accelerates relaxation and suppresses nonlinear oscillations, while the sign of $\chi_{3}$ determines whether transient photon bunching $(g^{(2)}(0)>1)$ or antibunching $(g^{(2)}(0)<1)$ occurs. All trajectories converge toward the coherent limit $g^{(2)}(0)\!\to\!1$, consistent with the steady–state solution of the unified action $\mathcal{S}_{\mathrm{U}}$.