Skip to content

Fluctuation–Dissipation Relation

The fluctuation–dissipation relation says that, in thermal equilibrium, noise and damping are not independent. The same microscopic degrees of freedom that exert random forces also absorb energy from a driven system.

In open quantum systems this relation is used to check thermal master equations, connect bath spectra to dissipative response, interpret vacuum noise, and avoid treating equilibrium noise as an arbitrary fitting function. The source-to-observable susceptibility is introduced in Green Functions and Response Preview and developed for many-body systems in the Kubo Formula. The Caldeira–Leggett Model is the standard oscillator-bath example where the same spectral density produces both a damping kernel and a thermal force-noise kernel.

The scope here is practical: how the relation appears in bath spectra, damping kernels, transition rates, and weak-coupling open-system modeling. Fluctuation–Dissipation Theorem owns the general many-body KMS and Lehmann proof, energy-normalized convention dictionary, matrix and momentum forms, and static-response caveat. The equilibrium static identities and the distinction between ordinary variance and Kubo–Mori covariance are developed in Fluctuations and Susceptibilities.

At equilibrium:

fluctuations reveal what the bath can do spontaneously;
dissipation reveals how the bath absorbs energy when driven.

The fluctuation–dissipation relation connects these two descriptions. It is not a generic statement about every noisy environment. It relies on equilibrium, stationarity, and linear response.

If a reservoir is driven, inverted, feedback-controlled, measurement-conditioned, or made of multiple baths at different temperatures, the equilibrium relation generally fails.

Let B(t)B(t) be a bath observable in a thermal stationary state. The symmetrized spectrum is

SBBsym(ω)=12∫−∞∞dt eiωt⟨{B(t),B(0)}⟩.S_{BB}^{\mathrm{sym}}(\omega) = \frac12 \int_{-\infty}^{\infty} dt\, e^{i\omega t} \langle \{B(t),B(0)\}\rangle.

Now perturb the bath by a weak classical source coupled to BB:

Hdrive(t)=−f(t)B.H_{\mathrm{drive}}(t) = -f(t)B.

The retarded susceptibility is defined by

χBBR(t)=iℏθ(t)⟨[B(t),B(0)]⟩.\chi_{BB}^R(t) = \frac{i}{\hbar} \theta(t) \langle [B(t),B(0)]\rangle.

It gives the linear response

δ⟨B(t)⟩=∫−∞∞dt′ χBBR(t−t′)f(t′).\delta\langle B(t)\rangle = \int_{-\infty}^{\infty} dt'\, \chi_{BB}^R(t-t')f(t').

The imaginary part of χBBR(ω)\chi_{BB}^R(\omega) is the dissipative part of the response, up to sign conventions. With the convention above, the common equilibrium form is

SBBsym(ω)=ℏcoth⁡ ⁣(βℏω2)Im⁡χBBR(ω).S_{BB}^{\mathrm{sym}}(\omega) = \hbar \coth\!\left( \frac{\beta\hbar\omega}{2} \right) \operatorname{Im}\chi_{BB}^R(\omega).

Some books define the retarded susceptibility with an extra minus sign. Then the displayed formula changes sign accordingly. The physical spectrum must remain nonnegative.

The ordered spectrum is

SBB(ω)=∫−∞∞dt eiωt⟨B(t)B(0)⟩.S_{BB}(\omega) = \int_{-\infty}^{\infty} dt\, e^{i\omega t} \langle B(t)B(0)\rangle.

For a Hermitian operator,

SBBsym(ω)=12[SBB(+ω)+SBB(−ω)].S_{BB}^{\mathrm{sym}}(\omega) = \frac12 \left[ S_{BB}(+\omega) + S_{BB}(-\omega) \right].

Thermal equilibrium also gives the Kubo–Martin–Schwinger relation

SBB(−ω)=e−βℏωSBB(+ω),ω>0.S_{BB}(-\omega) = e^{-\beta\hbar\omega} S_{BB}(+\omega), \qquad \omega>0.

The dissipative response is tied to the antisymmetric part:

SBB(+ω)−SBB(−ω)=2ℏ Im⁡χBBR(ω)S_{BB}(+\omega)-S_{BB}(-\omega) = 2\hbar\, \operatorname{Im}\chi_{BB}^R(\omega)

with the susceptibility convention used on this page. Combining this identity with the KMS relation gives

SBB(+ω)+SBB(−ω)SBB(+ω)−SBB(−ω)=coth⁡ ⁣(βℏω2),\frac{ S_{BB}(+\omega)+S_{BB}(-\omega) }{ S_{BB}(+\omega)-S_{BB}(-\omega) } = \coth\!\left( \frac{\beta\hbar\omega}{2} \right),

which is the fluctuation–dissipation relation in spectral form.

When

βℏ∣ω∣≪1,\beta\hbar|\omega|\ll1,

the thermal factor becomes

ℏcoth⁡ ⁣(βℏω2)≈2kBTω.\hbar \coth\!\left( \frac{\beta\hbar\omega}{2} \right) \approx \frac{2k_BT}{\omega}.

The relation becomes the classical fluctuation–dissipation form

SBBsym(ω)≈2kBTωIm⁡χBBR(ω).S_{BB}^{\mathrm{sym}}(\omega) \approx \frac{2k_BT}{\omega} \operatorname{Im}\chi_{BB}^R(\omega).

This is why thermal noise is proportional to temperature in classical regimes. The factor multiplying TT depends on what variable is being measured and on the response function convention.

When

βℏ∣ω∣≫1,\beta\hbar|\omega|\gg1,

the factor approaches

coth⁡ ⁣(βℏω2)⟶sgn⁡(ω).\coth\!\left( \frac{\beta\hbar\omega}{2} \right) \longrightarrow \operatorname{sgn}(\omega).

The symmetrized spectrum does not vanish at zero temperature:

SBBsym(ω)⟶ℏ∣Im⁡χBBR(ω)∣.S_{BB}^{\mathrm{sym}}(\omega) \longrightarrow \hbar |\operatorname{Im}\chi_{BB}^R(\omega)|.

This is the spectral language of zero-point fluctuations. However, zero-temperature symmetrized noise should not be confused with the ability of the bath to excite a system. The ordered negative-frequency spectrum is suppressed for an ordinary ground-state bath.

For a two-level system of transition frequency ω0>0\omega_0>0, weakly coupled to a thermal bath,

Γ↓∝SBB(+ω0),Γ↑∝SBB(−ω0).\Gamma_\downarrow \propto S_{BB}(+\omega_0), \qquad \Gamma_\uparrow \propto S_{BB}(-\omega_0).

The KMS relation gives

Γ↑Γ↓=e−βℏω0.\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = e^{-\beta\hbar\omega_0}.

This is the rate-level form of Detailed Balance. It ensures that a properly derived single-bath thermal weak-coupling generator relaxes toward a Gibbs state, subject to the usual Born, Markov, and secular assumptions.

The fluctuation–dissipation relation adds a response interpretation: the same bath spectral asymmetry that fixes upward and downward rates also fixes the relation between equilibrium noise and dissipative susceptibility.

For a resistor with resistance RR, the two-sided angular-frequency symmetrized voltage noise is often written

SVVsym(ω)=ℏωRcoth⁡ ⁣(βℏω2),S_{VV}^{\mathrm{sym}}(\omega) = \hbar\omega R \coth\!\left( \frac{\beta\hbar\omega}{2} \right),

with this convention. In the classical high-temperature limit,

SVVsym(ω)≈2RkBT.S_{VV}^{\mathrm{sym}}(\omega) \approx 2Rk_BT.

In one-sided ordinary-frequency notation this same result is commonly written with the familiar factor 4RkBT4Rk_BT. The physics is the same; the spectral convention is different.

At low temperature, the symmetrized voltage noise retains a quantum contribution proportional to ℏR∣ω∣\hbar R|\omega|. Whether that noise can excite a detector depends on the ordered spectrum and the detector’s coupling.

For a coordinate xx driven by a weak force ff, the retarded susceptibility satisfies

δ⟨x(ω)⟩=χxxR(ω)f(ω).\delta\langle x(\omega)\rangle = \chi_{xx}^R(\omega)f(\omega).

The fluctuation–dissipation relation gives

Sxxsym(ω)=ℏcoth⁡ ⁣(βℏω2)Im⁡χxxR(ω),S_{xx}^{\mathrm{sym}}(\omega) = \hbar \coth\!\left( \frac{\beta\hbar\omega}{2} \right) \operatorname{Im}\chi_{xx}^R(\omega),

up to the same sign convention for χR\chi^R. Peaks in equilibrium position noise occur where the oscillator also has dissipative response to a force.

This is the oscillator version of the same principle used in bath modeling: equilibrium fluctuations are weighted by the bath’s ability to absorb energy.

The fluctuation–dissipation relation does not mean every noise model must include a matching damping term. It means that an equilibrium thermal bath with linear response has constrained noise and dissipation.

It does not apply without additional work to:

  • nonthermal driven reservoirs;
  • active or inverted media;
  • feedback-controlled baths;
  • measurement-conditioned trajectories;
  • classical technical noise from electronics or drift;
  • multiple reservoirs at different temperatures;
  • aging, glassy, or nonstationary environments;
  • strong-coupling regimes where the bare system and bath cannot be cleanly separated.

In those cases one can still define spectra and response functions, but they are not tied by the equilibrium formula.

  • Using symmetrized zero-temperature noise as if it automatically excites a quantum system.
  • Applying the equilibrium relation to a driven or feedback-controlled environment.
  • Forgetting that the ordered spectrum, not the symmetrized spectrum alone, sets upward and downward transition rates.
  • Comparing formulas without checking the sign convention for χR\chi^R.
  • Mixing one-sided ordinary-frequency spectra with two-sided angular-frequency spectra.
  • Treating 1/f1/f technical noise as if it must satisfy an equilibrium fluctuation–dissipation relation.
  • Assuming dissipation can be added phenomenologically while leaving equilibrium noise arbitrary.

Use

S(−ω)=e−βℏωS(+ω)S(-\omega) = e^{-\beta\hbar\omega}S(+\omega)

to derive the upward-to-downward rate ratio for a two-level system.

Solution

With the convention on this page,

Γ↓∝S(+ω0),Γ↑∝S(−ω0).\Gamma_\downarrow\propto S(+\omega_0), \qquad \Gamma_\uparrow\propto S(-\omega_0).

Therefore

Γ↑Γ↓=S(−ω0)S(+ω0)=e−βℏω0.\frac{\Gamma_\uparrow}{\Gamma_\downarrow} = \frac{S(-\omega_0)}{S(+\omega_0)} = e^{-\beta\hbar\omega_0}.

Let x=βℏωx=\beta\hbar\omega and assume S(−ω)=e−xS(+ω)S(-\omega)=e^{-x}S(+\omega). Show that

S(+ω)+S(−ω)S(+ω)−S(−ω)=coth⁡ ⁣(x2).\frac{S(+\omega)+S(-\omega)} {S(+\omega)-S(-\omega)} = \coth\!\left(\frac{x}{2}\right).
Solution

Substitute S(−ω)=e−xS(+ω)S(-\omega)=e^{-x}S(+\omega):

S(+ω)+S(−ω)S(+ω)−S(−ω)=1+e−x1−e−x.\frac{S(+\omega)+S(-\omega)} {S(+\omega)-S(-\omega)} = \frac{1+e^{-x}}{1-e^{-x}}.

Multiplying numerator and denominator by ex/2e^{x/2} gives

ex/2+e−x/2ex/2−e−x/2=coth⁡ ⁣(x2).\frac{e^{x/2}+e^{-x/2}} {e^{x/2}-e^{-x/2}} = \coth\!\left(\frac{x}{2}\right).

Starting from

SVVsym(ω)=ℏωRcoth⁡ ⁣(βℏω2),S_{VV}^{\mathrm{sym}}(\omega) = \hbar\omega R \coth\!\left( \frac{\beta\hbar\omega}{2} \right),

show the high-temperature two-sided angular-frequency limit.

Solution

For βℏ∣ω∣≪1\beta\hbar|\omega|\ll1,

coth⁡ ⁣(βℏω2)≈2βℏω.\coth\!\left( \frac{\beta\hbar\omega}{2} \right) \approx \frac{2}{\beta\hbar\omega}.

Thus

SVVsym(ω)≈ℏωR2βℏω=2RkBT.S_{VV}^{\mathrm{sym}}(\omega) \approx \hbar\omega R \frac{2}{\beta\hbar\omega} = 2Rk_BT.

At zero temperature an equilibrium bath has no negative-frequency ordered noise for ω>0\omega>0 in this convention, but the symmetrized spectrum is nonzero. Explain why these statements are compatible.

Solution

The ordered negative-frequency spectrum describes the bath supplying energy to a system. At zero temperature an ordinary equilibrium bath cannot supply energy at positive transition frequency, so that part is suppressed.

The symmetrized spectrum is the average of ordered positive- and negative-frequency spectra:

Ssym(ω)=12[S(+ω)+S(−ω)].S^{\mathrm{sym}}(\omega) = \frac12 \left[ S(+\omega)+S(-\omega) \right].

At zero temperature S(−ω)S(-\omega) may vanish for ω>0\omega>0, but S(+ω)S(+\omega) need not. The bath can still absorb energy, and the symmetrized spectrum records zero-point fluctuations associated with that response.

  • H. B. Callen and T. A. Welton, “Irreversibility and generalized noise,” Physical Review 83, 34–40 (1951).
  • R. Kubo, “Statistical-mechanical theory of irreversible processes. I,” Journal of the Physical Society of Japan 12, 570–586 (1957).
  • R. Kubo, “The fluctuation-dissipation theorem,” Reports on Progress in Physics 29, 255–284 (1966).
  • 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 (2004).
  • A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf, “Introduction to quantum noise, measurement, and amplification,” Reviews of Modern Physics 82, 1155–1208 (2010).
  • U. Weiss, Quantum Dissipative Systems, World Scientific (2012).