Perturbative Analysis of Dark State Dynamics in Weakly Anharmonic Photon-Emitter Pairs

This paper investigates the origin of dissipation in dark states of weakly anharmonic photon-emitter pairs by applying first and second-order perturbative corrections to the wavefunction and analyzing their impact on system dynamics via the master equation.

The Gist

The Big Picture: The "Silent Room" Problem

Imagine a group of people in a room trying to whisper a secret to each other without anyone outside the room hearing them. In the quantum world, these "people" are tiny energy packets (photons) and "emitters" (like tiny circuits called transmons).

Usually, when these emitters interact, they are noisy. They leak energy into the environment, like a leaky bucket losing water. This noise destroys their ability to hold onto information. However, scientists have found a special trick: if they arrange these emitters just right, they can create a "Dark State."

Think of a Dark State as a perfectly synchronized dance. If two dancers move in perfect opposition (one steps left, the other steps right), their movements cancel each other out from the perspective of the audience. To the outside world, it looks like nothing is happening. No energy is leaked; the secret is safe.

The Problem: The "Rigid" Dance Floor

In the ideal world (which the paper calls the "harmonic" regime), these dancers are perfectly rigid. They follow the rules exactly, and the Dark State is stable.

However, real-world quantum devices (like the transmons mentioned in the paper) aren't perfectly rigid. They have a little bit of "wiggle" or flexibility. The paper calls this anharmonicity.

The paper asks a simple question: What happens to our perfect silent dance when the dancers start to wiggle?

The authors found that even a tiny bit of wiggle breaks the perfect cancellation. The "Dark State" is no longer perfectly dark. It starts to leak energy (dissipate) and eventually collapses. The silence is broken.

The Solution: A "Predictive Map"

The authors wanted to understand exactly how and why this silence breaks, so they could fix it.

Usually, calculating what happens when things get "wiggly" is a nightmare for computers. It's like trying to predict the path of a leaf in a hurricane; the math gets messy, and the computer often crashes or gives wrong answers (a problem the paper calls "numerical instability").

Instead of brute-forcing the math, the authors used a perturbation method.

• The Analogy: Imagine you are trying to predict how a car drives on a bumpy road. Instead of simulating every single rock and pothole, you start with a model of the car on a smooth road. Then, you add a small "correction" for the bumps. You calculate the first bump, then the second, and so on.
• The Paper's Approach: They treated the "wiggle" (anharmonicity) as a tiny disturbance. They calculated the first-order correction (the immediate effect of the wiggle) and the second-order correction (the effect of the wiggle on the wiggle).

What They Discovered

By using this "correction" method, they mapped out the fate of the Dark State:

1. The Leak: The wiggle causes the Dark State to mix with a "Bright State" (a noisy, loud state). It's as if the dancers accidentally step out of sync, and suddenly the audience can hear them.
2. The Burst: When the system starts to relax (lose energy), it doesn't just fade away quietly. Because the Dark State is now slightly connected to the Bright State, the system releases a sudden, intense burst of energy (photons) before settling down. The authors call this a "superradiant burst."
   • Analogy: Imagine a dam that is supposed to hold back water perfectly. A small crack (the wiggle) forms. Instead of a slow drip, the water pressure builds up for a split second and then bursts out in a massive wave before the water level finally drops.
3. Even vs. Odd: They found a quirky rule based on the number of energy packets:
   • If you start with an even number of packets, the system eventually drains all the way to the ground (zero energy).
   • If you start with an odd number, the system gets "stuck" in a middle state and can't drain all the way down.

Why This Matters (According to the Paper)

The authors showed that their "correction" method (the predictive map) works almost exactly as well as the heavy, messy computer simulations, but it is much faster and more stable.

• The Benefit: Because they have this map, they can predict exactly how the "Dark State" will decay.
• The Goal: If you know exactly how the dance breaks down, you might be able to adjust the music (the control parameters) to keep the dancers in sync longer. This helps in keeping quantum information safe for longer periods.

Summary in One Sentence

The paper shows that even a tiny bit of "wiggle" in quantum devices breaks their perfect silence, causing a sudden burst of energy, but the authors have created a simple mathematical tool to predict exactly how this happens so we can potentially fix it.

Technical Summary

Technical Summary: Perturbative Analysis of Dark State Dynamics in Weakly Anharmonic Photon-Emitters Pairs

Problem Statement
Dark states are excited quantum states in open systems that decouple from their environment, preventing the emission or absorption of photons. While these states offer potential for decoherence-free quantum information storage and long-range entanglement, their stability is compromised when the system includes on-site interactions (anharmonicity). In the context of transmon qubits coupled to a waveguide, modeled by the Bose-Hubbard model, the inclusion of anharmonicity ($U$) breaks the perfect decoupling of dark states, inducing dissipation.

Theoretical analysis of these systems is complicated by their non-Hermitian nature. The effective Hamiltonian governing the dynamics possesses complex eigenvalues and requires the use of biorthogonal left and right eigenvectors. As anharmonicity increases, the system approaches exceptional points where eigenvectors coalesce, leading to numerical instabilities and difficulties in tracking state populations, particularly in higher excitation manifolds beyond the qubit model. Existing methods often struggle to provide a stable, analytical description of these dynamics without losing information through renormalization corrections.

Methodology
The authors propose a perturbative approach to analyze the origin of dissipation in dark states by treating weak anharmonicity as a perturbation to a system of two dissipatively coupled harmonic oscillators.

1. Model Formulation: The system is described using an effective non-Hermitian Bose-Hubbard Hamiltonian. The Hamiltonian includes harmonic terms, a collective decay term (dissipative coupling via a waveguide), and an anharmonic interaction term ($U$).
2. Perturbative Expansion: The density matrix $\hat{\rho}$ and state vectors $|\psi\rangle$ are expanded in a power series with respect to the perturbation strength $\epsilon$ (representing anharmonicity).
   • The zeroth-order solution corresponds to the unperturbed dark state in the harmonic limit ($U=0$).
   • First-order corrections are calculated for both the energy and the wavefunction using non-degenerate perturbation theory adapted for non-Hermitian systems.
   • Second-order corrections are derived, incorporating the effects of jump operators and state mixing.
3. Dynamics Tracking: The corrected density matrices (up to second order) are applied to the master equation to track the time-dependent evolution of the system.
4. Validation: The perturbative results are compared against exact numerical diagonalization of the non-Hermitian Hamiltonian and standard numerical evolution of the master equation. To handle the biorthogonal nature of the eigenvectors near exceptional points, the authors employ a simultaneous Gram-Schmidt orthogonalization method for left and right eigenvectors during initialization.

Key Contributions and Results

• Mechanism of Dissipation: The analysis reveals that anharmonicity induces dissipation in dark states primarily through an induced coupling between the dark state and a neighboring dissipative (bright) state that shares similar collective parity.
• Analytical Corrections:
  • First Order: The first-order correction to the wavefunction introduces a complex-valued correlation between the dark state and the bright state. This creates a transition pathway allowing population to leak from the dark state into dissipative channels.
  • Second Order: The second-order correction to the energy reveals a shift in the imaginary component, directly quantifying the induced dissipation rate. The correction scales with the excitation number $N$ as $\lambda_c(N) \propto -i N(N-1)U^2/16\gamma$.
• Dynamical Behavior:
  • The perturbative method successfully reproduces the numerical relaxation dynamics of the dark state into the ground state with high accuracy for weak anharmonicity ($U/\gamma < 0.5$).
  • The dynamics exhibit an initial burst of photon emission (superradiant-type signature) before settling into a near-exponential decay. This burst arises because the population initially mixes into the bright state (which is maximally dissipative) faster than the state decays.
  • Parity Dependence: The final relaxation state depends on the parity of the initial excitation number. Even excitation numbers relax to the ground state, while odd excitation numbers relax and remain trapped in the $N=1$ dark state, as the first excitation manifold is unaffected by anharmonicity in this model.
• Stability: The perturbative approach provides a numerically stable and trace-preserving method to track population dynamics, circumventing the instabilities associated with biorthogonal wavefunction coalescence near exceptional points.

Significance and Claims
The paper claims that this perturbative framework offers a robust alternative to direct numerical diagonalization for studying non-Hermitian open quantum systems with weak anharmonicity. By providing analytical expressions for wavefunction and energy corrections, the method allows for the construction of a "map" of transition pathways and correlations between sub- and superradiant states.

The authors state that this approach is particularly useful for:

1. Understanding the breakdown of decoherence-free spaces in higher excitation manifolds (beyond the qubit limit).
2. Predicting population dynamics and coherence times for entangled states in transmon arrays.
3. Informing quantum control schemes in superconducting circuits by identifying how anharmonicity affects state fidelity.

The work is modest in its scope, acknowledging that the analytical corrections deviate from numerical results as anharmonicity increases toward and beyond the exceptional point, consistent with the limitations of perturbation theory. The authors suggest future work could involve scaling the system to larger arrays and investigating other dephasing factors, but the current paper focuses strictly on validating the perturbative method against established numerical techniques for a two-transmon system.