Detailed Balance

Detailed balance is the equilibrium consistency condition that relates forward and reverse transition rates. In open quantum systems it explains why a single thermal bath tends to drive an undriven weakly coupled system toward a Gibbs state, and why upward thermal transitions are suppressed at low temperature.

For a two-level system with transition frequency $\omega_0>0$, the most familiar form is

$$\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = e^{-\beta\hbar\omega_0}.$$

The downward rate $\Gamma_\downarrow$ describes relaxation from excited to ground state. The upward rate $\Gamma_\uparrow$ describes thermal excitation. Their ratio is fixed by the bath temperature if the bath is a single equilibrium reservoir and the weak-coupling Markovian assumptions apply.

Detailed balance is stronger than merely having some steady state. It is a statement of equilibrium reversibility: in equilibrium, every elementary transition is balanced by its reverse transition with the correct Gibbs weights.

Classical Rate Balance

For a classical Markov process with states $m,n$ and transition rates $W_{m\leftarrow n}$, detailed balance with respect to an equilibrium distribution $p_n^{\mathrm{eq}}$ means

$$W_{m\leftarrow n}p_n^{\mathrm{eq}} = W_{n\leftarrow m}p_m^{\mathrm{eq}}.$$

If

$$p_n^{\mathrm{eq}} = \frac{e^{-\beta E_n}}{Z},$$

then

$$\frac{W_{m\leftarrow n}}{W_{n\leftarrow m}} = e^{-\beta(E_m-E_n)}.$$

This equation says that transitions raising the system energy are Boltzmann suppressed relative to the reverse transitions.

Quantum detailed balance reduces to this statement for populations after the secular approximation decouples populations from coherences. The corresponding population dynamics is treated in Pauli Rate Equations. The full quantum version also tracks coherences and operator structure.

Qubit Upward and Downward Rates

Let the ground and excited states be $\lvert g\rangle$ and $\lvert e\rangle$, with

$$E_e-E_g=\hbar\omega_0, \qquad \omega_0>0.$$

A finite-temperature two-level master equation has

$$\frac{d\rho}{dt} = \Gamma_\downarrow\mathcal D[\sigma_-]\rho + \Gamma_\uparrow\mathcal D[\sigma_+]\rho,$$

where

$$\sigma_-=\lvert g\rangle\langle e\rvert, \qquad \sigma_+=\lvert e\rangle\langle g\rvert.$$

The excited-state population obeys

$$\dot p_e = -\Gamma_\downarrow p_e + \Gamma_\uparrow p_g, \qquad p_g=1-p_e.$$

At stationarity,

$$\Gamma_\downarrow p_e^{\mathrm{ss}} = \Gamma_\uparrow p_g^{\mathrm{ss}},$$

so

$$\frac{p_e^{\mathrm{ss}}}{p_g^{\mathrm{ss}}} = \frac{\Gamma_\uparrow}{\Gamma_\downarrow}.$$

For a thermal steady state,

$$\frac{p_e^{\beta}}{p_g^{\beta}} = e^{-\beta\hbar\omega_0}.$$

Therefore detailed balance requires

$$\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = e^{-\beta\hbar\omega_0}.$$

At zero temperature, $\beta\to\infty$, this gives $\Gamma_\uparrow/\Gamma_\downarrow\to0$: an ordinary passive bath cannot thermally excite the system.

Bohr-Frequency Convention

In the weak-coupling construction used in this volume, system coupling operators are decomposed as

$$A_\alpha(\omega) = \sum_{\epsilon'-\epsilon=\hbar\omega} \Pi(\epsilon)A_\alpha\Pi(\epsilon').$$

They satisfy

$$[H_S,A_\alpha(\omega)] = -\hbar\omega A_\alpha(\omega).$$

Thus positive $\omega$ labels an operator that lowers the system energy by $\hbar\omega$. It is associated with emission or downward relaxation. The reverse operator is

$$A_\alpha^\dagger(\omega) = A_\alpha(-\omega),$$

which raises the system energy.

With this convention, a thermal bath spectrum satisfies the single-operator relation

$$S(-\omega) = e^{-\beta\hbar\omega}S(\omega), \qquad \omega>0.$$

The upward rate samples the negative-frequency spectrum and is Boltzmann suppressed relative to the downward rate.

Fourier-transform signs differ across books. Before comparing formulas, check whether the source uses $e^{i\omega t}$ or $e^{-i\omega t}$ and whether positive $\omega$ labels absorption by the system or emission into the bath.

KMS Condition

For a bath in thermal equilibrium,

$$\rho_B = \frac{e^{-\beta H_B}}{Z_B},$$

correlation functions obey the Kubo–Martin–Schwinger condition. For bath operators $B_\alpha(t)$ and $B_\beta(0)$, define

$$C_{\alpha\beta}(t) = \operatorname{Tr}_B \left[ B_\alpha(t)B_\beta(0)\rho_B \right].$$

With the spectrum convention

$$\Gamma_{\alpha\beta}(\omega) = \int_{-\infty}^{\infty} dt\, e^{i\omega t} C_{\alpha\beta}(t),$$

the KMS relation implies, schematically,

$$\Gamma_{\alpha\beta}(-\omega) = e^{-\beta\hbar\omega} \Gamma_{\beta\alpha}(\omega),$$

up to the ordering and Fourier-convention choices of the source.

For one Hermitian bath operator this becomes the common detailed-balance spectrum relation:

$$S(-\omega) = e^{-\beta\hbar\omega}S(\omega).$$

This is the microscopic origin of the qubit rate ratio.

Gibbs Stationarity

For a system Hamiltonian

$$H_S = \sum_n E_n\Pi_n,$$

the Gibbs state is

$$\rho_\beta = \frac{e^{-\beta H_S}}{\operatorname{Tr}(e^{-\beta H_S})}.$$

A thermal weak-coupling generator should satisfy

$$\mathcal L(\rho_\beta)=0$$

or the corresponding statement with a properly renormalized Hamiltonian when Lamb shifts or weak-coupling corrections are included.

Detailed balance is stronger than stationarity. A generator may have $\rho_\beta$ as a steady state while still supporting irreversible circulating currents in a degenerate or driven sector. Quantum detailed balance rules out such equilibrium currents by imposing a symmetry of the generator with respect to a thermal inner product.

At the level of this chapter, the practical message is:

single equilibrium bath + weak coupling + secularization
  -> KMS spectra
  -> detailed-balance rates
  -> Gibbs steady state

Each arrow uses assumptions.

Harmonic Oscillator Bath Example

For a harmonic oscillator coupled to a thermal reservoir, a common Lindblad equation is

$$\frac{d\rho}{dt} = - \frac{i}{\hbar}[\hbar\omega a^\dagger a,\rho] + \kappa(\bar n+1)\mathcal D[a]\rho + \kappa\bar n\,\mathcal D[a^\dagger]\rho.$$

The downward rate is proportional to $\bar n+1$: stimulated plus spontaneous emission. The upward rate is proportional to $\bar n$: thermal absorption.

For a Bose occupation

$$\bar n = \frac{1}{e^{\beta\hbar\omega}-1},$$

the ratio is

$$\frac{\kappa\bar n}{\kappa(\bar n+1)} = e^{-\beta\hbar\omega}.$$

This is detailed balance for adjacent oscillator levels.

Fluctuation–Dissipation Relation

Detailed balance is closely related to fluctuation–dissipation relations. For a stationary bath operator $B(t)$, define the symmetrized noise spectrum

$$S_{BB}^{\mathrm{sym}}(\omega) = \frac12 \int_{-\infty}^{\infty} dt\, e^{i\omega t} \langle \{B(t),B(0)\}\rangle$$

and the retarded susceptibility

$$\chi_{BB}^R(t) = \frac{i}{\hbar} \theta(t) \langle [B(t),B(0)]\rangle.$$

In a common convention, equilibrium gives

$$S_{BB}^{\mathrm{sym}}(\omega) = \hbar \coth\!\left( \frac{\beta\hbar\omega}{2} \right) \operatorname{Im}\chi_{BB}^R(\omega).$$

This relation says that equilibrium noise and linear response are not independent. The same KMS condition that fixes upward and downward transition rates also constrains the balance between fluctuations and dissipation.

For a broader discussion of ordered versus symmetrized bath spectra and their use in rates, see Noise Spectra. Fluctuation–Dissipation Theorem owns the general KMS and Lehmann derivation of the equilibrium response relation. Fluctuation–Dissipation Relation owns its bath-noise and damping applications.

When Detailed Balance Fails

Detailed balance is not expected in every Markovian master equation.

It can fail for:

• multiple baths at different temperatures;
• driven systems in rotating frames;
• feedback-controlled dynamics;
• measurement-conditioned dynamics;
• nonthermal reservoirs;
• inverted media;
• approximate generators before secularization;
• strong coupling where the steady state is not the bare Gibbs state.

Failure of detailed balance does not automatically mean the equation is mathematically invalid. It means the model is not an equilibrium single-bath thermal relaxation model.

Practical Checks

When reviewing a proposed thermal master equation, check:

• What Hamiltonian defines the energy gaps?
• Which sign convention defines positive Bohr frequency?
• Which operators lower and raise the system energy?
• Are upward and downward rates related by $e^{-\beta\hbar\omega}$?
• Does the steady state match the expected Gibbs state?
• Are Lamb shifts or renormalized energies included consistently?
• Are degeneracies treated by a proper secular block?
• Is there more than one bath or a drive that should produce nonequilibrium currents?
• Are symmetrized noise spectra being confused with transition-rate spectra?

These are the thermal parts of the broader Approximation Checklist.

Common Mistakes

Reversing the Boltzmann factor

For $\omega_0>0$, upward transitions are suppressed:

$$\Gamma_\uparrow/\Gamma_\downarrow = e^{-\beta\hbar\omega_0}.$$

The inverse ratio would describe an inverted bath or a sign-convention mismatch.

Using zero-temperature amplitude damping at finite temperature

Zero-temperature amplitude damping has no upward transition. If $\Gamma_\uparrow$ is appreciable, use a finite-temperature model. See Amplitude-Damping Channel for the finite-temperature warning.

Confusing stationarity with detailed balance

$\mathcal L(\rho_\beta)=0$ says the Gibbs state is stationary. Detailed balance additionally constrains reverse processes and rules out equilibrium currents.

Ignoring degeneracies

Degenerate or nearly degenerate transitions require careful secularization. Treating each matrix element as an independent classical rate can destroy coherences or violate positivity.

Expecting detailed balance with multiple baths

Two baths at different temperatures usually drive heat currents and nonequilibrium steady states. A single Gibbs detailed-balance condition no longer applies.

Exercises

Two-level steady state

For

$$\dot p_e = -\Gamma_\downarrow p_e + \Gamma_\uparrow(1-p_e),$$

find $p_e^{\mathrm{ss}}$ and the condition for a Gibbs steady state.

Solution

Set $\dot p_e=0$:

$$\Gamma_\downarrow p_e^{\mathrm{ss}} = \Gamma_\uparrow(1-p_e^{\mathrm{ss}}).$$

Thus

$$p_e^{\mathrm{ss}} = \frac{\Gamma_\uparrow} {\Gamma_\uparrow+\Gamma_\downarrow}.$$

The ratio of excited to ground populations is

$$\frac{p_e^{\mathrm{ss}}}{p_g^{\mathrm{ss}}} = \frac{\Gamma_\uparrow}{\Gamma_\downarrow}.$$

For a Gibbs state with spacing $\hbar\omega_0$,

$$\frac{p_e^\beta}{p_g^\beta} = e^{-\beta\hbar\omega_0},$$

so detailed balance requires

$$\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = e^{-\beta\hbar\omega_0}.$$

Oscillator rate ratio

Show that for

$$\bar n=\frac{1}{e^{\beta\hbar\omega}-1},$$

one has

$$\frac{\bar n}{\bar n+1} = e^{-\beta\hbar\omega}.$$

Solution

Compute

$$\bar n+1 = \frac{1}{e^{\beta\hbar\omega}-1}+1 = \frac{e^{\beta\hbar\omega}} {e^{\beta\hbar\omega}-1}.$$

Therefore

$$\frac{\bar n}{\bar n+1} = \frac{1} {e^{\beta\hbar\omega}} = e^{-\beta\hbar\omega}.$$

KMS to upward suppression

Assume a single bath operator has spectrum satisfying

$$S(-\omega)=e^{-\beta\hbar\omega}S(\omega), \qquad \omega>0.$$

If $S(\omega_0)$ gives the downward rate and $S(-\omega_0)$ gives the upward rate, find $\Gamma_\uparrow/\Gamma_\downarrow$.

Solution

With the stated convention,

$$\Gamma_\downarrow\propto S(\omega_0), \qquad \Gamma_\uparrow\propto S(-\omega_0).$$

Therefore

$$\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = \frac{S(-\omega_0)}{S(\omega_0)} = e^{-\beta\hbar\omega_0}.$$

Zero-temperature limit

What happens to $\Gamma_\uparrow/\Gamma_\downarrow$ as $T\to0$ for an ordinary thermal bath?

Solution

As $T\to0$, $\beta=1/(k_BT)\to\infty$. For $\omega_0>0$,

$$e^{-\beta\hbar\omega_0}\to0.$$

Thus

$$\frac{\Gamma_\uparrow}{\Gamma_\downarrow}\to0.$$

An ordinary zero-temperature bath can absorb energy but cannot thermally excite the system.

References

• R. Kubo, “Statistical-Mechanical Theory of Irreversible Processes. I,” Journal of the Physical Society of Japan 12, 570-586 (1957).
• P. C. Martin and J. Schwinger, “Theory of Many-Particle Systems. I,” Physical Review 115, 1342-1373 (1959).
• E. B. Davies, Quantum Theory of Open Systems, Academic Press (1976).
• R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, Springer (1987).
• H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).
• C. W. Gardiner and P. Zoller, Quantum Noise, Springer, 3rd ed. (2004).
• Á. Rivas and S. F. Huelga, Open Quantum Systems: An Introduction, Springer (2012).