Benchmarking Atomic Ionization Driven by Strong Quantum Light

The recently available high-intensity quantum light pulses provide novel tools for controlling light-matter interactions. However, the rigor of the theoretical frameworks currently used to describe the interaction of strong quantum light with atoms and molecules remains unverified. Here, we establish a rigorous benchmark by solving the fully quantized time-dependent Schrödinger equation for an atom exposed to bright squeezed vacuum light. Our \textit{ab initio} simulations reveal a critical limitation of the widely used $Q$-representation: although it accurately reproduces the total photoelectron spectrum after tracing over photon states, it completely fails to capture the electron-photon joint energy spectrum. To overcome this limitation, we develop a general theoretical framework based on the Feynman path integral that properly incorporates the electron-photon quantum entanglement. Our results provide both quantitative benchmarks and fundamental theoretical insights for the emerging field of strong-field quantum optics.

💡 Research Summary

The authors present a rigorous benchmark for strong‑field ionization driven by bright squeezed‑vacuum (BSV) light, a non‑classical quantum‑light source now experimentally accessible at intensities comparable to atomic Coulomb fields (10¹⁴ W cm⁻²). They solve the fully quantized time‑dependent Schrödinger equation (TDSE) for a one‑dimensional hydrogen atom coupled to a quantized photon mode, treating both the electron and the field quantum mechanically. The wavefunction is expanded in a direct‑product basis of electron position |x⟩ and photon Fock states |n⟩, and propagated with a Crank‑Nicolson scheme. After the pulse, electron and photon observables are extracted using window‑operator techniques.

The key observable is the joint electron‑photon energy spectrum, i.e., the photoelectron energy distribution conditioned on a specific final photon number. The full quantum simulation reveals rich structure: ATI peaks appear as stripes spaced by the photon energy ω, tilted due to the photon‑number‑dependent ponderomotive shift Uₚ(n). Two clear cut‑offs at E = 2Uₚ(n) and E = 10Uₚ(n) correspond to direct and rescattering ionization, respectively, in agreement with semiclassical strong‑field theory. Moreover, fine horizontal fringes show photon‑number‑dependent modulations, and a striking even‑odd alternation: spectra correlated with even photon numbers differ markedly from those with odd numbers.

To assess the commonly used Q‑representation (which replaces the quantized field by a statistical mixture of coherent states weighted by the Q‑function), the authors compute the same joint spectra. While the total photoelectron spectrum (integrated over photon number) from the Q‑method matches the full quantum result, the joint spectra do not: the Q‑approach yields a smooth, continuous distribution lacking the even‑odd modulation and the photon‑number‑resolved fringes. Consequently, the Q‑representation reproduces averaged observables but fails to capture electron‑photon entanglement and the dynamical redistribution of photon statistics.

Photon‑number distributions further illustrate this failure. In the full TDSE, the initial BSV state (strictly even‑photon parity) evolves into a mixed state with both even and odd components, and the local dominance of even versus odd numbers varies across the distribution, directly reflecting entanglement with the electron. By contrast, the Q‑representation predicts a time‑independent photon‑number distribution identical to the initial one, showing no parity mixing.

To resolve the discrepancy, the authors develop a path‑integral formulation using coherent‑state bases. By integrating over photon coherent trajectories, they obtain an approximate state |Ψ(t)⟩ ≈ ∫ d²α ⟨α|ϕ_γ⟩ |ψₑ(t;α)⟩⊗|α⟩, where each coherent amplitude α drives an independent electron wavefunction ψₑ(t;α). The resulting density matrix retains off‑diagonal terms R(α*,β)=⟨α|ϕ_γ⟩⟨ϕ_γ|β⟩ e^{-(|α|²+|β|²)/2}, which encode electron‑photon entanglement. This “R‑representation” reproduces both the joint spectra and the time‑dependent photon‑number distributions observed in the full TDSE. Neglecting the off‑diagonal terms reduces the R‑representation to the Q‑representation, explaining why the latter loses the dynamical photon statistics.

Physically, the authors argue that the dominant mechanism for photon‑field evolution is not the small back‑action of the electron on the coherent amplitude (which scales with the weak coupling ϵ_V) but rather the entanglement‑mediated redistribution of photon number. Hence, even though the electron’s influence on α is minor, the entangled structure of the state drives the observed parity mixing and spectral modulations.

In summary, the paper provides (i) a high‑fidelity benchmark dataset for strong‑field ionization with quantum light, (ii) a clear demonstration that the Q‑representation is adequate only for observables that ignore photon‑electron correlations, and (iii) a practical, analytically tractable R‑representation that captures the essential quantum correlations while remaining computationally efficient. These results establish a solid foundation for future theoretical and experimental work in strong‑field quantum optics, including quantum control of high‑harmonic generation, non‑sequential double ionization, and photon‑resolved spectroscopy of ultrafast processes.