Skip to content

Noise Spectra

A noise spectrum is the frequency-domain version of a correlation function. It tells which fluctuation frequencies are available to dephase a system, drive transitions, broaden a line, or pass through a control sequence.

In open quantum systems, spectra are not just plotting tools. They enter weak-coupling master equations, detailed-balance relations, dynamical-decoupling estimates, quantum-limited measurement, and experimental noise budgets.

For the conceptual distinction between classical stochastic noise, quantum bath noise, vacuum noise, thermal noise, and technical noise, see Quantum Noise. For the modeling distinction between a bath, reservoir, noise source, and record, see Baths, Reservoirs, and Environments.

The main warning is simple:

classical spectra are usually symmetric;
quantum unsymmetrized spectra generally are not.

That asymmetry is what lets a thermal bath distinguish absorption from emission.

Let B(t)B(t) be a bath operator in the bath Heisenberg picture, and let the bath state be stationary. The two-time correlation depends only on the time difference:

CBB(t)=⟨B(t)B(0)⟩.C_{BB}(t) = \langle B(t)B(0)\rangle.

With the Fourier convention used on this page, the two-sided unsymmetrized spectrum is

SBB(ω)=∫−∞∞dt eiωtCBB(t).S_{BB}(\omega) = \int_{-\infty}^{\infty} dt\, e^{i\omega t} C_{BB}(t).

The inverse convention is

CBB(t)=∫−∞∞dω2π e−iωtSBB(ω).C_{BB}(t) = \int_{-\infty}^{\infty} \frac{d\omega}{2\pi}\, e^{-i\omega t} S_{BB}(\omega).

If BB has units of energy, then SBBS_{BB} has units of energy squared times time. If BB is a frequency fluctuation, then its spectrum has units of angular frequency.

Fourier conventions vary. Some authors place 2π2\pi factors in the forward transform, and some use e−iωte^{-i\omega t} instead of eiωte^{i\omega t}. Before comparing rates, check the definition.

For a quantum operator, the order of B(t)B(t) and B(0)B(0) matters. 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.

This is often what a classical detector reports in a high-temperature or weakly invasive limit. It is also the spectrum that appears in many engineering noise measurements.

The symmetrized spectrum does not by itself determine upward and downward quantum transition rates. Those rates use ordered, unsymmetrized correlations because the bath must either absorb or supply energy.

For a real classical stationary random process ξ(t)\xi(t),

Cξξ(t)=E[ξ(t)ξ(0)],Sξξ(ω)=∫−∞∞dt eiωtCξξ(t).C_{\xi\xi}(t) = \mathbb E[\xi(t)\xi(0)], \qquad S_{\xi\xi}(\omega) = \int_{-\infty}^{\infty} dt\, e^{i\omega t} C_{\xi\xi}(t).

Since Cξξ(t)=Cξξ(−t)C_{\xi\xi}(t)=C_{\xi\xi}(-t) for a real stationary process with ordinary multiplication,

Sξξ(−ω)=Sξξ(ω).S_{\xi\xi}(-\omega)=S_{\xi\xi}(\omega).

This symmetry is not generally true for a quantum unsymmetrized spectrum. It is a common source of wrong temperature factors in relaxation formulas.

Physics derivations usually use two-sided spectra over positive and negative frequencies. Many experimental plots use one-sided spectra over ω≥0\omega\ge0 or f≥0f\ge0.

A common classical convention is

Sξξ(1)(ω)=2Sξξ(ω),ω>0,S_{\xi\xi}^{(1)}(\omega) = 2S_{\xi\xi}(\omega), \qquad \omega>0,

with the same total variance:

⟨ξ2⟩=∫−∞∞dω2πSξξ(ω)=∫0∞dω2πSξξ(1)(ω).\langle \xi^2\rangle = \int_{-\infty}^{\infty} \frac{d\omega}{2\pi} S_{\xi\xi}(\omega) = \int_{0}^{\infty} \frac{d\omega}{2\pi} S_{\xi\xi}^{(1)}(\omega).

The same idea is used with ordinary frequency ff rather than angular frequency ω=2πf\omega=2\pi f, but then the units change. Confusing radians per second with cycles per second creates factors of 2π2\pi.

In the convention

SBB(ω)=∫dt eiωt⟨B(t)B(0)⟩,S_{BB}(\omega) = \int dt\, e^{i\omega t} \langle B(t)B(0)\rangle,

positive and negative frequencies have different physical meanings. For a thermal bath and a system transition of angular frequency ω0>0\omega_0>0:

  • SBB(+ω0)S_{BB}(+\omega_0) controls processes in which the bath absorbs energy from the system;
  • SBB(−ω0)S_{BB}(-\omega_0) controls processes in which the bath supplies energy to the system.

In equilibrium the two are related by detailed balance:

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

for a single Hermitian bath operator with this convention. At low temperature, the negative-frequency thermal excitation spectrum is suppressed.

This is the spectral form of Detailed Balance. The same thermal structure underlies the general Fluctuation–Dissipation Theorem, while the open-system Fluctuation–Dissipation Relation develops bath noise and damping applications.

For a system–bath interaction

HI=∑αAα⊗Bα,H_I = \sum_\alpha A_\alpha\otimes B_\alpha,

weak-coupling Markovian master equations involve bath spectra

Γαβ(ω)=1ℏ2∫−∞∞dt eiωt⟨Bα(t)Bβ(0)⟩.\Gamma_{\alpha\beta}(\omega) = \frac{1}{\hbar^2} \int_{-\infty}^{\infty} dt\, e^{i\omega t} \langle B_\alpha(t)B_\beta(0)\rangle.

The system coupling operators are decomposed into Bohr-frequency components:

[HS,Aα(ω)]=−ℏωAα(ω).[H_S,A_\alpha(\omega)] = -\hbar\omega A_\alpha(\omega).

After a secular approximation, the dissipator has the schematic form

∑ω,α,βΓαβ(ω)(Aβ(ω)ρAα†(ω)−12{Aα†(ω)Aβ(ω),ρ}),\sum_{\omega,\alpha,\beta} \Gamma_{\alpha\beta}(\omega) \left( A_\beta(\omega)\rho A_\alpha^\dagger(\omega) - \frac12 \{A_\alpha^\dagger(\omega)A_\beta(\omega),\rho\} \right),

plus a Lamb-shift Hamiltonian from the imaginary principal-value part of the same bath response.

This formula shows why spectra must be defined with their ordering and sign conventions. Changing the definition of A(ω)A(\omega) or the Fourier transform changes which side of the spectrum is called positive frequency.

For a two-level system with transition frequency ω0\omega_0, transverse bath noise drives energy exchange. Schematically,

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

The longitudinal relaxation time satisfies

1T1=Γ↓+Γ↑.\frac{1}{T_1} = \Gamma_\downarrow+\Gamma_\uparrow.

At zero temperature, Γ↑\Gamma_\uparrow vanishes for an ordinary equilibrium bath, but Γ↓\Gamma_\downarrow need not vanish because the bath can absorb energy. This is the spectral language behind spontaneous emission and amplitude damping.

Longitudinal frequency noise samples the spectrum near zero frequency. For a qubit with

H(t)=ℏ2[ω0+ξ(t)]σz,H(t) = \frac{\hbar}{2} \left[ \omega_0+\xi(t) \right]\sigma_z,

the off-diagonal density-matrix element accumulates the random phase

φ(t)=∫0tds ξ(s).\varphi(t) = \int_0^t ds\,\xi(s).

For Gaussian stationary noise,

⟨e−iφ(t)⟩=exp⁡ ⁣[−12⟨φ(t)2⟩].\langle e^{-i\varphi(t)}\rangle = \exp\!\left[ - \frac12 \langle\varphi(t)^2\rangle \right].

If the noise is effectively white on the timescale of the experiment, then

1Tϕ=12Sξξ(0),\frac{1}{T_\phi} = \frac12 S_{\xi\xi}(0),

with the two-sided convention above. Other Hamiltonian normalizations move factors of 22 between ξ\xi and SξξS_{\xi\xi}.

The key physical point is robust: pure dephasing is sensitive to low-frequency or elastic fluctuations, while relaxation samples noise at transition frequencies.

Control sequences reshape how a system samples noise. For a toggling function y(t)=±1y(t)=\pm1 that changes sign under ideal refocusing pulses, the accumulated phase from classical dephasing noise is

φ(T)=∫0Tdt y(t)ξ(t).\varphi(T) = \int_0^T dt\,y(t)\xi(t).

Define the filter amplitude

Y(ω,T)=∫0Tdt y(t)eiωt.Y(\omega,T) = \int_0^T dt\, y(t)e^{i\omega t}.

For Gaussian dephasing noise,

χ(T)=12∫−∞∞dω2π Sξξ(ω)∣Y(ω,T)∣2,\chi(T) = \frac12 \int_{-\infty}^{\infty} \frac{d\omega}{2\pi}\, S_{\xi\xi}(\omega) |Y(\omega,T)|^2,

and the coherence factor is e−χ(T)e^{-\chi(T)}.

This formula explains why spin echo suppresses very slow noise: the sign-changing control makes Y(ω,T)Y(\omega,T) small near ω=0\omega=0. The canonical pulse-sequence discussion is Dynamical Decoupling.

Several idealized spectra recur across open-system modeling.

SpectrumTime-domain sourceTypical use
whiteshort correlation timeMarkovian dephasing or diffusion
Lorentzianexponential correlationfinite-bandwidth noise, random telegraph approximations
Ohmicoscillator bath with J(ω)∝ωJ(\omega)\propto\omega at low frequencydamping, Brownian motion, circuit environments
sub-Ohmicenhanced low-frequency bath weightslow environments, strong memory effects
super-Ohmicsuppressed low-frequency bath weightphonon-like environments in some dimensions
1/f1/fbroad distribution of slow fluctuatorssolid-state dephasing and drift

For an exponential classical correlation

C(t)=σ2e−∣t∣/τc,C(t)=\sigma^2e^{-|t|/\tau_c},

the two-sided spectrum is Lorentzian:

S(ω)=2σ2τc1+ω2τc2.S(\omega) = \frac{2\sigma^2\tau_c} {1+\omega^2\tau_c^2}.

The correlation time τc\tau_c sets the bandwidth. In the limit τc\tau_c becomes very short while the area remains fixed, the noise approaches a white-noise approximation.

The phrase spectral density is used in two related but distinct ways. The mode-density convention is treated directly in Spectral Densities.

For an oscillator bath, one often defines a mode spectral density such as

J(ω)=∑k∣gk∣2δ(ω−ωk),ω>0.J(\omega) = \sum_k |g_k|^2 \delta(\omega-\omega_k), \qquad \omega>0.

This object describes the density of bath modes weighted by coupling strengths. A finite-temperature ordered noise spectrum also includes occupation factors:

SBB(+ω)∝J(ω) [n(ω)+1],SBB(−ω)∝J(ω) n(ω),S_{BB}(+\omega) \propto J(\omega)\,[n(\omega)+1], \qquad S_{BB}(-\omega) \propto J(\omega)\,n(\omega),

for ω>0\omega>0, in a common bosonic convention.

Thus J(ω)J(\omega) and SBB(ω)S_{BB}(\omega) are not always interchangeable. The spectrum relevant for rates includes both coupling density and bath state.

For the canonical coordinate-coupled oscillator-bath model where Ohmic and Drude spectral densities become damping kernels, see Caldeira–Leggett Model.

For the canonical two-level-system version, where Ohmic, sub-Ohmic, and super-Ohmic spectra control dissipative tunneling and dephasing, see Spin-Boson Model.

A Markov approximation needs the bath correlation to decay on a timescale short compared with system evolution. Spectrally, a short correlation time usually means a broad, smooth spectrum over the system’s relevant transition frequencies.

Useful checks:

  • Is S(ω)S(\omega) smooth near each Bohr frequency?
  • Does the spectrum have sharp resonances, gaps, or band edges?
  • Is there strong low-frequency weight that produces slow drift?
  • Are there finite-size recurrences in the bath?
  • Is the system being driven so that sidebands sample additional frequencies?

If the spectrum has narrow structure on the scale of the system linewidth, a memoryless Lindblad model may miss nonexponential decay, coherent exchange with a mode, or non-Markovian information return. Lorentzian structured spectra are a common entry point for Pseudomode Methods.

For Gaussian baths, another common numerical strategy is to expand the bath correlation function into exponentials and evolve the corresponding Hierarchical Equations of Motion.

  • Using a symmetrized spectrum to compute quantum upward and downward rates without checking ordering.
  • Forgetting whether the source uses angular frequency ω\omega or ordinary frequency ff.
  • Mixing one-sided and two-sided spectra.
  • Assuming S(+ω)=S(−ω)S(+\omega)=S(-\omega) for a quantum bath at low temperature.
  • Treating a mode density J(ω)J(\omega) as the full thermal noise spectrum.
  • Ignoring cross-spectra when several bath operators couple to the system.
  • Reading a white-noise approximation as valid at arbitrarily high frequency.
  • Inferring Markovianity from a smooth-looking plot without comparing timescales.

Show that, with the convention on this page,

⟨ξ2⟩=∫−∞∞dω2πSξξ(ω).\langle \xi^2\rangle = \int_{-\infty}^{\infty} \frac{d\omega}{2\pi} S_{\xi\xi}(\omega).
Solution

Set t=0t=0 in the inverse transform:

Cξξ(0)=∫−∞∞dω2πSξξ(ω).C_{\xi\xi}(0) = \int_{-\infty}^{\infty} \frac{d\omega}{2\pi} S_{\xi\xi}(\omega).

But

Cξξ(0)=E[ξ(0)2]=⟨ξ2⟩.C_{\xi\xi}(0)=\mathbb E[\xi(0)^2]=\langle\xi^2\rangle.

Compute the spectrum of

C(t)=σ2e−∣t∣/τc.C(t)=\sigma^2e^{-|t|/\tau_c}.
Solution

Use the evenness of C(t)C(t):

S(ω)=2σ2∫0∞dt e−t/τccos⁡(ωt).S(\omega) = 2\sigma^2 \int_0^\infty dt\, e^{-t/\tau_c} \cos(\omega t).

The integral is

∫0∞dt e−t/τccos⁡(ωt)=τc1+ω2τc2.\int_0^\infty dt\, e^{-t/\tau_c} \cos(\omega t) = \frac{\tau_c}{1+\omega^2\tau_c^2}.

Therefore

S(ω)=2σ2τc1+ω2τc2.S(\omega) = \frac{2\sigma^2\tau_c} {1+\omega^2\tau_c^2}.

For a thermal bath with

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

what is the ratio Γ↑/Γ↓\Gamma_\uparrow/\Gamma_\downarrow for a two-level system of transition frequency ω0\omega_0?

Solution

With the convention used here,

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

Thus

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

Upward transitions are Boltzmann suppressed in thermal equilibrium.

For

H(t)=ℏ2[ω0+ξ(t)]σz,H(t) = \frac{\hbar}{2} \left[ \omega_0+\xi(t) \right]\sigma_z,

assume Gaussian white frequency noise with two-sided spectrum Sξξ(ω)=S0S_{\xi\xi}(\omega)=S_0. Show that the coherence decays with rate S0/2S_0/2.

Solution

The random phase is

φ(t)=∫0tds ξ(s).\varphi(t) = \int_0^t ds\,\xi(s).

For white noise,

⟨φ(t)2⟩=∫0tds∫0tds′ ⟨ξ(s)ξ(s′)⟩=S0t.\langle \varphi(t)^2\rangle = \int_0^t ds \int_0^t ds'\, \langle \xi(s)\xi(s')\rangle = S_0t.

Gaussian averaging gives

⟨e−iφ(t)⟩=exp⁡ ⁣[−12S0t].\langle e^{-i\varphi(t)}\rangle = \exp\!\left[ - \frac12S_0t \right].

Therefore

1Tϕ=S02.\frac{1}{T_\phi} = \frac{S_0}{2}.
  • 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).
  • H. J. Carmichael, Statistical Methods in Quantum Optics 1, Springer (1999).
  • 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).
  • R. Kubo, “The fluctuation-dissipation theorem,” Reports on Progress in Physics 29, 255–284 (1966).