Bath-induced deviations from Gibbs statistics for strongly interacting oscillators

This paper demonstrates that for two strongly interacting quantum oscillators coupled to independent baths, non-secular terms in the Redfield master equation can drive the system into a non-Gibbs steady state with significant deviations from Boltzmann statistics when the oscillators are unequally damped, due to bath-induced coherences between nearly-degenerate levels.

The Gist

Imagine you have two identical pendulums (let's call them Pendulum A and Pendulum B) hanging next to each other. They are connected by a stiff spring, so when one swings, it strongly pulls the other along. This is what physicists call "strongly interacting oscillators."

Now, imagine each pendulum is swinging in a different room. Room A is a bit breezy (it has a little bit of air resistance), while Room B is very windy (it has a lot of air resistance). Both rooms are at the exact same temperature.

The Old Way of Thinking (The "Gibbs" Rule)
For a long time, scientists believed that if you waited long enough, both pendulums would eventually settle into a predictable, calm rhythm based solely on the temperature of the rooms. This is called a "Gibbs state." In this ideal world, the pendulums would act like they are in perfect thermal equilibrium, and their energy levels would follow a standard, boring rulebook.

The New Discovery
This paper says: "Wait a minute. That rulebook isn't always right."

The authors found that because the two pendulums are so tightly connected (strongly interacting) and because they are being slowed down by the air in their rooms at different rates (unequal damping), they don't settle into that standard calm rhythm. Instead, they get stuck in a weird, persistent state that breaks the usual rules.

The "Leaky Bucket" Analogy
Think of the two pendulums as two buckets connected by a pipe.

• Bucket A has a small hole (low damping).
• Bucket B has a huge hole (high damping).
• Both buckets are being filled with water from a faucet at the same rate (the same temperature).

In a normal world, you'd expect the water levels to stabilize based on the faucet's pressure. But because the buckets are connected by a special pipe (the strong interaction) and the holes are different sizes, something strange happens. The water doesn't just sit there. It starts flowing in a continuous loop: water moves from Bucket A to Bucket B, but because Bucket B leaks so fast, the system creates a constant, invisible current.

This "current" is what the paper calls an excitation flux. It's a steady stream of energy flowing from the less-damped oscillator to the more-damped one, driven by a subtle quantum "ghost" connection (called coherence) between the two.

Why Does This Happen?
Usually, scientists ignore the tiny, fast-oscillating details of how these systems interact to keep the math simple. They use a shortcut called the "secular approximation." This shortcut assumes the system will eventually become perfectly calm and follow the standard rules.

However, this paper shows that when you have two strongly connected pendulums with different amounts of friction, those "tiny details" you ignored actually matter. They act like a hidden engine that keeps the system from ever truly settling down into the standard "Gibbs" state.

The Key Takeaways

1. Unequal Friction is the Trigger: If both pendulums had the same amount of air resistance, they would behave normally. The "weird" behavior only happens because one is damped more than the other.
2. Resonance is Key: This effect is strongest when the pendulums are naturally tuned to swing at the same frequency (resonance). If they are tuned to very different frequencies, the effect disappears, and they go back to following the normal rules.
3. A New Steady State: The system reaches a "steady state," but it's not the calm, predictable one we expected. It's a state where the energy levels of the two pendulums are permanently unbalanced, and energy is constantly flowing between them, even though the whole setup is at a constant temperature.

In Summary
The paper demonstrates that when two strongly connected quantum objects are cooled by environments that treat them differently, they don't just "cool down" to a standard temperature. Instead, they enter a unique, non-standard state where energy flows continuously between them, defying the traditional expectations of how heat and equilibrium work.

Technical Summary

Technical Summary: Bath-Induced Deviations from Gibbs Statistics for Strongly Interacting Oscillators

Problem Statement
The study addresses a fundamental limitation in modeling open quantum systems composed of strongly interacting subsystems. While the Redfield quantum master equation is widely used for systems weakly coupled to baths, it does not inherently guarantee positive probabilities. To ensure positivity and thermalization into a Gibbs state, the "secular approximation" is typically applied, which neglects non-secular terms (population-coherence couplings) in the energy basis. However, for strongly interacting subsystems, this approximation may break down or obscure long-time effects. Previous approaches, such as local-Lindblad master equations, can yield thermodynamically inconsistent steady states (e.g., violating the second law) for strongly coupled systems. Conversely, the mean force Gibbs state accounts for system-bath interactions but generally deviates from standard Gibbs statistics. The central question is whether non-secular terms, when preserved in a physically consistent manner, drive strongly interacting oscillators into a steady state that deviates from the canonical Gibbs distribution, even when the baths are at thermal equilibrium.

Methodology
The authors model an open quantum system consisting of two quantum harmonic oscillators ($A$ and $B$) with frequencies $\omega_a$ and $\omega_b$, coupled via a strength $g$. Each oscillator is weakly coupled to an independent thermal bath at the same temperature $T$, but with potentially different damping rates ($\gamma_a \neq \gamma_b$) determined by Ohmic spectral densities.

The dynamics are derived using the Redfield master equation under the Born-Markov approximation. To address the issue of negative probabilities while retaining relevant physics, the authors employ a partial secularization (or coarse-graining) of the Redfield dynamics.

• Transformation: The interacting oscillators are transformed into normal modes ($\hat{c}_1, \hat{c}_2$) with frequencies $\Omega_1$ and $\Omega_2$.
• Approximation: In the regime where oscillators are nearly resonant ($|\Delta| < g$, where $\Delta = \omega_a - \omega_b$), the authors retain the slower oscillating non-secular terms (frequency $\Omega_1 - \Omega_2 = 2g$) while discarding faster oscillating terms. This preserves the coupling between normal mode populations and coherences.
• Analysis: The resulting master equation includes secular terms (Lamb shifts, standard dissipation) and non-secular terms (coherent coupling between modes and cross-dissipation). The authors analytically evaluate the time derivative of the system density operator at the Gibbs state to determine if it is a stationary solution. They further analyze the steady-state occupation numbers and the excitation currents flowing between the baths.

Key Contributions

1. Derivation of a Non-Gibbs Steady State: The paper demonstrates that for two strongly interacting oscillators unequally damped by independent baths at the same temperature, the steady state is not a Gibbs state ($\hat{\rho}_G \propto e^{-\beta \hat{H}_S}$).
2. Identification of Mechanism: The deviation is driven by the interplay between the imaginary part of the bath correlation functions ($S(\omega)$) and the non-secular terms in the master equation. Specifically, bath-induced coherences between nearly-degenerate oscillator levels generate a steady excitation flux.
3. Analytical Conditions: The authors establish that the deviation vanishes in the dispersive regime ($|\Delta| \gg g$) or when damping rates are equal ($\gamma_a = \gamma_b$), recovering the Gibbs state. The deviation is maximized under resonant conditions ($\Delta = 0$) with unequal damping.

Results

• Occupation Deviations: In the resonant regime ($T \gtrsim g$), the steady-state occupation numbers of the oscillators deviate significantly (up to ~10%) from the thermal Boltzmann distribution. One oscillator exhibits an increased occupation while the other decreases, maintaining a balance condition ($\gamma_a \delta \langle \hat{a}^\dagger \hat{a} \rangle + \gamma_b \delta \langle \hat{b}^\dagger \hat{b} \rangle \approx 0$) that conserves the total number of excitations in the system.
• Unbalanced Probability Distributions: The probability distributions of the oscillators over Fock states are unbalanced ($P^{(a)}_\nu \neq P^{(b)}_\nu \neq P^{(th)}_\nu$). The more heavily damped oscillator tends to be more excited than the lightly damped one, a feature absent in the Gibbs state.
• Steady Excitation Flux: A non-zero steady flux of excitations exists between the two baths ($d\langle \hat{N}_a \rangle/dt = -d\langle \hat{N}_b \rangle/dt \neq 0$). This flux is driven by the stationary coherence between the oscillators ($\text{Im}\langle \hat{a}\hat{b}^\dagger \rangle_{ss}$) and flows from the lightly damped oscillator to the highly damped one. This flux vanishes in the dispersive limit where the system relaxes to the Gibbs state.
• Connection to Local-Lindblad Models: The dependence of occupation deviations on steady coherence mirrors findings in local-Lindblad models used to explain cavity-modified vibrational occupations, suggesting a link between these different theoretical approaches in the resonant regime.

Significance and Claims
The paper claims that non-Gibbs steady states are a feasible mechanism to describe apparent modifications of thermal occupations in strongly interacting quantum subsystems that are in thermal equilibrium with their environment. The authors argue that these deviations are not artifacts of incorrect modeling but are physical consequences of preserving non-secular terms in the Redfield dynamics when oscillators are unequally damped.

The work suggests that experimental signatures of this phenomenon could be observed in platforms exhibiting strong light-matter coupling, such as superconductor resonators, plasmonic nanocavities, or Fabry-Perot cavities. The primary signature is the resonant deviation of occupation numbers from Boltzmann statistics and the presence of a steady excitation flux between baths at equal temperature, sustained by bath-induced coherences. The study highlights that standard secular approximations may fail to capture these long-time steady-state features in strongly coupled systems.