Comparison of the standard and dressed-picture master equations for the quantum Rabi model in the ultrastrong coupling regime

This paper numerically compares the standard Gorini-Kossakowski-Sudarshan-Lindblad master equation with the dressed-picture Markovian master equation for the quantum Rabi model in the ultrastrong coupling regime, demonstrating that the latter provides a more consistent description of relaxation and dephasing by explicitly accounting for light-matter hybridization across various initial states and noise spectra.

The Gist

Imagine you are trying to predict how a tiny, dancing pair of partners—a qubit (a two-level atom) and a photon (a particle of light)—move together inside a mirrored box (a cavity).

In the world of quantum physics, this dance is called the Quantum Rabi Model. Usually, when they dance, they are gentle and polite. But in the Ultrastrong Coupling Regime, they are dancing so wildly, so close, and so fast that they essentially become a single, hybrid creature. They are no longer two separate dancers; they are a "super-dancer."

This paper is a guidebook for scientists trying to predict how this "super-dancer" behaves when the room gets noisy, hot, or when the music changes. Specifically, it compares two different rulebooks (mathematical equations) used to predict this behavior.

Here is the breakdown of the paper's story, using simple analogies:

1. The Two Rulebooks (The Master Equations)

To predict how the system loses energy or gets "messy" (dissipation and decoherence), scientists use mathematical formulas called Master Equations. The paper compares two of them:

• The "Standard" Rulebook (GKSL):

  • The Analogy: Imagine a choreographer who watches the qubit and the photon separately. They say, "Okay, the qubit gets tired and slows down. The photon gets absorbed by the walls and disappears."
  • The Problem: In the "Ultrastrong" regime, the qubit and photon are so mixed up that they aren't separate anymore. The Standard Rulebook tries to treat them as individuals even though they are a hybrid. It's like trying to describe a smoothie by listing the calories of the apple and the banana separately, ignoring that they are now blended. This leads to wrong predictions, like saying the system gains energy when it should be losing it.

• The "Dressed" Rulebook (DME):

  • The Analogy: This choreographer looks at the "smoothie" (the hybrid state). They understand that the qubit and photon are fused. When the system interacts with the noisy environment, they calculate how the entire hybrid creature reacts.
  • The Benefit: This is much more accurate for wild, strong dances. However, it is much harder to calculate. It's like trying to solve a puzzle where the pieces keep changing shape as you look at them.

2. The Experiment: Testing the Rules

The authors didn't just talk about theory; they ran massive computer simulations to see which rulebook works better. They tested the "super-dancer" in many different scenarios:

• Different Starting Poses: They started the dance with the light in different states:
  • Coherent State: A steady, rhythmic beat.
  • Schrödinger Cat State: A spooky mix of "being here" and "being there" at the same time (like a cat that is both alive and dead).
  • Squeezed States: A dance where the uncertainty is squeezed into one direction.
  • Thermal State: A chaotic, jittery dance caused by heat.
• Different Dance Intensities: They turned up the coupling strength (how hard they dance together) from a gentle sway to a violent thrash.
• Different Noises: They simulated the environment being "white noise" (static) or "Ohmic noise" (a specific type of friction).

3. What They Found

The results were a mix of "surprisingly similar" and "drastically different":

• When they agree: For some gentle dances or specific measurements, both rulebooks give similar answers. If you just want a rough idea, the "Standard" book is okay.
• When they disagree: For the wild, ultra-strong dances, the "Standard" book often fails.
  • It predicts the dance stops too quickly or too slowly.
  • It fails to predict how much "entanglement" (the quantum bond between the partners) remains.
  • It sometimes predicts the system creates energy out of nothing (which is impossible).
• The "Dressed" book is usually the truth: The "Dressed" approach correctly predicts that the hybrid nature of the system makes it lose its quantum "magic" (coherence) faster than the Standard book thinks.

4. The Special Trick: Making Light from Nothing

In the final section, the authors looked at a cool trick called the Dynamical Casimir Effect.

• The Analogy: Imagine shaking the mirrored box so violently that you rip photons out of the empty vacuum.
• The Twist: They used a "Post-Selection" filter. They said, "We will only count the light if the qubit ends up in the 'sleeping' state."
• The Result: By filtering the results this way, they created a very special, high-quality light beam that is incredibly useful for quantum metrology (measuring things with extreme precision). Even though the two rulebooks disagreed on how the light was made, they agreed that this special, high-precision light could be created.

The Bottom Line

This paper is a user manual for researchers.

It tells us: "If you are studying light and matter dancing together very strongly, do not just use the old, simple rulebook. It will give you the wrong answer. You need to use the more complex 'Dressed' rulebook that accounts for their hybrid nature."

However, it also offers a practical guide: "Here is exactly how to code these complex equations on a computer, and here is a cheat sheet showing you exactly when the simple rulebook is 'good enough' and when it will lead you astray."

In short: When things get crazy strong in the quantum world, you can't treat the parts separately. You have to treat the whole messy, beautiful hybrid as one.

Technical Summary

1. Problem Statement

The paper addresses the challenge of modeling open quantum systems in the Ultrastrong Coupling (USC) regime ($g \gtrsim 0.1\omega$), where the interaction strength between a two-level system (qubit) and a quantized cavity field becomes a significant fraction of the cavity frequency.

• The Limitation of Standard Approaches: In this regime, the standard Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) master equation, derived in the "bare" (uncoupled) basis, often yields unphysical predictions. This is because strong light-matter interaction hybridizes the bare atom and field states, making the assumption that dissipation acts on uncoupled eigenstates invalid.
• The Need for Dressed Pictures: A more rigorous approach requires the Dressed Master Equation (DME), derived in the eigenbasis of the interacting Hamiltonian (dressed states). However, DMEs are analytically intractable and numerically complex, often requiring diagonalization of large Hamiltonians and evaluation of transition rates in the dressed basis.
• The Gap: There is a lack of accessible, step-by-step numerical implementations comparing these two approaches across various initial states and coupling strengths to determine exactly where the standard GKSL approximation fails.

2. Methodology

The authors provide a comprehensive tutorial and numerical framework for solving both the standard GKSL and the DME (specifically the formulation by Beaudoin, Gambetta, and Blais) for the Quantum Rabi Model.

• Numerical Implementation:

  • GKSL: The density operator is expanded in the bare basis $\{|g, f_n\rangle, |e, f_n\rangle\}$. The authors derive a system of coupled ordinary differential equations (ODEs) for the expansion coefficients, utilizing a Runge-Kutta-Verner method for integration.
  • DME: The density operator is expanded in the dressed eigenbasis $\{|n\rangle\}$ of the Rabi Hamiltonian. This requires diagonalizing the Hamiltonian to find eigenvalues ($\Lambda_n$) and eigenvectors. The dissipative terms are constructed using transition rates dependent on the spectral density of the reservoir and matrix elements of operators (e.g., $\sigma_z, \sigma_x, X$) in the dressed basis.
  • Basis Transformation: To calculate observables in the bare basis (like photon number), the authors provide explicit transformation formulas to convert the dressed-state density matrix coefficients back to the bare-state coefficients.

• Parameters and Conditions:

  • Coupling Strength ($g$): Varied from $0.05\omega$ to $0.8\omega$.
  • Initial States: Coherent, odd Schrödinger cat, squeezed coherent, squeezed vacuum, and thermal states.
  • Reservoir Spectral Densities: Two types were tested: White noise (constant rates) and Ohmic noise (frequency-dependent rates with cutoffs).
  • Dynamic Scenarios: Multiphoton Rabi oscillations (3, 5, 7-photon processes) and non-stationary dynamics involving time-dependent modulation of the qubit frequency (analogous to the dynamical Casimir effect) with post-selection protocols.

3. Key Contributions

1. Step-by-Step Tutorial: The paper offers explicit formulas and algorithmic steps for implementing both GKSL and DME solvers, lowering the barrier for researchers to perform these complex numerical simulations.
2. Systematic Comparison: It provides a direct, quantitative comparison of the two master equations across a wide range of physical observables, including:
   • Qubit excited-state population ($P_e$).
   • Mean photon number ($\langle n \rangle$).
   • Mandel $Q$-factor (photon statistics).
   • Subsystem purities ($\Pi_q, \Pi_c$) and Negativity (entanglement measure).
   • Photon-number probability distributions.
3. Post-Selection and Metrology: The study investigates photon generation from vacuum via parameter modulation. It evaluates the Quantum Fisher Information (QFI) of the post-selected cavity states to determine their metrological advantage over classical light.

4. Key Results

• Discrepancies in Dynamics:
  • Coherence Decay: The DME generally predicts faster decay of coherences and entanglement (negativity) compared to the GKSL equation, particularly in the dispersive regime. The GKSL equation often underestimates the damping of high-frequency oscillations.
  • Dependence on Observables: For some quantities (e.g., mean photon number in specific resonance regimes), the differences between GKSL and DME are small. However, for entanglement measures and purity, the discrepancies can be substantial.
  • Spectral Density Sensitivity: The results highlight that accurate modeling requires precise knowledge of the reservoir's spectral density (White vs. Ohmic), as different densities yield different quantitative predictions, especially at long times.
• Multiphoton Resonances: In 3-, 5-, and 7-photon resonance regimes, the qualitative behavior of the three models (GKSL, DME-White, DME-Ohmic) is similar, but quantitative deviations increase with the coupling strength $g$.
• Vacuum Generation & Metrology:
  • For short times ($\omega t \lesssim 2 \times 10^4$), the GKSL and DME predictions are nearly indistinguishable, validating previous studies that used only GKSL.
  • At longer times, DME and GKSL diverge.
  • Post-Selection: Conditioning on the qubit remaining in the ground state ($|g\rangle$) successfully engineers non-trivial cavity states. These states exhibit metrological advantage (QFI > mean photon number) and possess a highly non-monotonic photon-number distribution with a pronounced tail, making them useful for quantum metrology.

5. Significance

• Methodological Guidance: The work serves as a practical guide for researchers entering the USC field, clarifying when the simpler GKSL equation is sufficient and when the computationally expensive DME is strictly necessary.
• Physical Insight: It confirms that while the standard GKSL equation is robust for certain observables (like mean energy) in specific regimes, it fails to capture the correct decoherence rates and entanglement dynamics in the ultrastrong coupling regime.
• Experimental Relevance: By demonstrating that post-selected states generated via modulation possess metrological advantages, the paper suggests new pathways for generating non-classical light in circuit QED and cavity QED experiments, provided the correct master equation is used to model the system's evolution.

In conclusion, the paper establishes that while the standard GKSL master equation offers a reasonable approximation for some macroscopic observables in the USC regime, the dressed-picture master equation is essential for accurately describing decoherence, entanglement, and the long-time dynamics of open quantum systems in this regime.