Prodiabatic Elimination: Higher Order Elimination of Fast Variables with Quantum Noise

1. Problem Statement

In open quantum systems, particularly those involving light-matter coupling (e.g., cavity QED, circuit QED), exact descriptions of dynamics are often intractable. A common approximation is adiabatic elimination, which removes fast degrees of freedom (like cavity photons) to derive an effective model for slow variables (like atomic states).

However, standard adiabatic elimination has significant limitations:

• Leading Order Only: It is typically restricted to the leading order in the timescale separation parameter ($\epsilon$), limiting its accuracy when the separation is not infinite.
• Noise Handling: Extending the method to higher orders is challenging because it requires a systematic treatment of quantum noise (vacuum fluctuations). Existing higher-order methods often rely on master equations, which struggle with multi-time correlation functions (e.g., $g^{(2)}$) and are computationally complex.
• Short-Time Dynamics: Standard elimination often fails to capture short-time behaviors where vacuum noise contributions are critical.

2. Methodology: Prodiabatic Elimination

The authors introduce Prodiabatic Elimination, a systematic approximation technique that extends adiabatic elimination to higher orders while consistently incorporating quantum noise.

Core Framework:

• Quantum Langevin Equations (QLE): Unlike previous higher-order approaches that use master equations, this method is based on QLEs. This naturally facilitates the description of multi-time correlation functions.
• Timescale Separation: The method assumes a "bad-cavity" limit where the cavity decay rate $\kappa$ is much larger than the atomic frequency $\Omega$, coupling $g$, and drive $f$. A small parameter $\epsilon$ is defined such that:
  $$O(\gamma/\kappa) = O(\Omega/\kappa) = \epsilon^2, \quad O(g/\kappa) = O(f/\kappa) = \epsilon$$
• Formal Solution and Expansion:
  1. The cavity field operator $\hat{a}(t)$ is expressed via a formal integral solution involving the atomic operator $\hat{b}(t-\tau)$.
  2. The atomic operator $\hat{b}(t-\tau)$ is expanded in a Taylor series with respect to the delay $\tau$.
  3. Crucially, the expansion retains terms up to order $\epsilon^3$ (and higher for noise), capturing corrections beyond the leading order.
• Noise Inclusion: The derivation explicitly includes input noise operators ($\hat{a}_{in}, \hat{b}_{in}$). The authors derive a specific commutator relation between the integrated noise term $\hat{A}_{in}$ and the atomic operators, showing that noise contributions are essential for short-time dynamics and multi-time correlations.

Key Resulting Formula:
The cavity annihilation operator is approximated as:
$$\hat{a}(t) \approx \hat{a}_{pdb}(t) + \hat{A}_{in}(t) + \text{noise correction terms}$$
where $\hat{a}_{pdb}$ is an operator acting solely on the atom, containing higher-order corrections to the standard adiabatic result. This allows the calculation of multi-time correlation functions $G(\{\tau_m, t_n\})$ using only atomic operators and the Quantum Regression Theorem.

3. Key Contributions

• Systematic Higher-Order Expansion: The method provides a controlled expansion that goes beyond the leading order ($\epsilon^1$), including terms up to $\epsilon^3$ and higher, significantly improving accuracy.
• Consistent Noise Treatment: It is the first approach to systematically include vacuum noise contributions in higher-order elimination for multi-time observables without resorting to complex master equation derivations.
• Computational Efficiency: Despite the higher order, the method retains the simplicity of standard adiabatic elimination, expressing all dynamics in terms of the slow variables (the atom) alone.
• Analytical Expressions: It yields analytical formulas for complex observables, such as the second-order coherence function $g^{(2)}(t)$, which are difficult to obtain with other methods.

4. Results and Demonstrations

The authors validated the method on two distinct physical setups:

A. Driven-Dissipative Jaynes-Cummings Model

• Setup: A two-level atom coupled to a driven cavity.
• Findings:
  • The prodiabatic elimination accurately captures the time-dependent behavior of the atomic population $\langle \hat{\sigma}_z \rangle$ for drive strengths where standard adiabatic elimination fails.
  • $g^{(2)}$ Function: The method provides an analytical expression for the second-order correlation function. It correctly predicts the short-time behavior (the "bunching" or "antibunching" peak) driven by vacuum noise, which standard adiabatic elimination misses.
  • Effective Master Equation: The authors derived an equivalent Lindblad master equation that reproduces the prodiabatic results up to order $\epsilon^3$. This equation includes a modified decay rate and a new jump operator ($\hat{\sigma} + \xi \hat{\sigma}_z$), representing the physical process of photon re-absorption.

B. Stimulated Raman Adiabatic Passage (STIRAP)

• Setup: A three-level $\Lambda$-system in a two-mode cavity used for population transfer between ground states.
• Findings:
  • The method successfully derives the "dark state" of the system with higher-order corrections.
  • In non-ideal STIRAP protocols (where pulses change rapidly), the standard adiabatic elimination fails to capture the delay in the system's response mediated by the cavity.
  • Prodiabatic elimination correctly models the population of the excited state and the delay, showing clear improvement over the standard method and matching numerical simulations closely.

5. Significance and Outlook

• Robustness: The technique is shown to be a robust tool for analyzing open quantum systems where timescale separation exists but is not infinite.
• Experimental Relevance: By retaining the ability to calculate multi-time correlations (crucial for quantum optics experiments like photon counting), the method bridges the gap between theoretical approximations and experimental observables.
• Future Applications: The authors suggest this method will be vital for optimal control protocols in cavity QED, where shaped pulses and cavity filtering interact. It also opens the door to treating thermal noise and developing general procedures for deriving Lindblad master equations for complex open systems.

In summary, Prodiabatic Elimination represents a significant advancement in the theoretical toolkit for quantum optics, offering a balance between the analytical simplicity of adiabatic elimination and the accuracy of full numerical simulations, specifically by solving the long-standing problem of incorporating higher-order noise effects.