Decoherence Mechanisms in Qubits

Every physical qubit is an open quantum system. Its two computational levels are never truly alone: they couple, however weakly, to a vast reservoir of electromagnetic modes, lattice phonons, nuclear and electron spins, charged defects, and amorphous two-level systems in the materials around them. Decoherence is the process by which this coupling quietly converts a pure, controllable quantum state into a statistical mixture, erasing the relative phases that quantum algorithms depend on. It is the single most important practical obstacle to scalable quantum computation, and it is precisely what the DiVincenzo criteria have in mind when they demand coherence times long compared to gate times. Phenomenologically the story of a single qubit reduces to two timescales. The longitudinal relaxation time T1T_1 measures the decay of populations — the qubit dumping a quantum of energy into the bath and falling from ∣1⟩\lvert 1\rangle toward equilibrium. The transverse coherence time T2T_2 measures the decay of the off-diagonal density-matrix elements, the loss of phase memory in a superposition. They are tied together by the universal inequality T2≤2T1T_2 \le 2T_1, with the slack quantified by a pure-dephasing time TϕT_\phi. What makes the subject deep rather than empirical is that both rates follow from a microscopic system–bath Hamiltonian through the Bloch–Redfield or Lindblad master equation, and the answer is startlingly simple: each rate is a value of the environmental noise spectral density S(ω)S(\omega) sampled at one frequency. Relaxation reads SS at the qubit frequency ω01\omega_{01}; pure dephasing reads SS near DC. Everything else — why 1/f1/f noise dephases without relaxing, why sweet spots and spin echoes work, why transmons beat charge qubits — follows from that one sentence.

Key takeaways
  • Decoherence is the loss of quantum coherence as a qubit entangles with uncontrolled environmental degrees of freedom; it is captured by the relaxation time T1T_1 and the dephasing time T2T_2, related by 1/T2=1/(2T1)+1/Tϕ1/T_2 = 1/(2T_1) + 1/T_\phi.
  • T1T_1 (energy relaxation) is driven by noise at the qubit frequency ω01\omega_{01}; TϕT_\phi (pure dephasing) by low-frequency noise that jitters the splitting without exchanging energy — and complete positivity forces T2≤2T1T_2 \le 2T_1, so a reported T2>2T1T_2 > 2T_1 is always an error.
  • The Bloch–Redfield master equation makes this exact: Γ1∝S(ω01)\Gamma_1 \propto S(\omega_{01}) and Γϕ∝S(ω ⁣→ ⁣0)\Gamma_\phi \propto S(\omega\!\to\!0), where S(ω)S(\omega) is the noise power spectral density; ubiquitous 1/f1/f noise gives non-exponential (Gaussian) Ramsey decay.
  • Dominant mechanisms differ by platform: charge, flux, two-level-system, Purcell and quasiparticle loss in superconducting qubits; magnetic-field, motional-heating and laser noise in trapped ions; photon loss in photonics.
  • Coherence is engineered upward at first-order-insensitive sweet spots, by dynamical decoupling, and by materials work — all of which either reshape S(ω)S(\omega) or decouple the qubit from it.
  • Key idea: decoherence rates are values of the environmental noise spectral density read at specific frequencies — relaxation samples S(ω01)S(\omega_{01}) at the qubit frequency, pure dephasing samples S(ω→0)S(\omega\to 0) near DC — so engineering coherence means either reshaping S(ω)S(\omega) or decoupling the qubit from it.
The core insight

Decoherence rates are dictated by the environmental noise spectral density evaluated at specific frequencies: relaxation samples S(ω01)S(\omega_{01}) at the qubit frequency, while pure dephasing samples S(ω→0)S(\omega\to 0) near DC. Engineering coherence therefore means either reshaping S(ω)S(\omega) (better materials, filtering) or decoupling the qubit from it (sweet spots, dynamical decoupling).

Where this lives today

Decoherence is the number every quantum-hardware group is measured by, and the past two decades have been a campaign against it. Superconducting-qubit coherence has climbed from the nanoseconds of the first Cooper-pair box in 1999 to hundreds of microseconds in three-dimensional and tantalum-based transmons, driven by attacking the terms in S(ω)S(\omega): better dielectrics to starve the two-level-system bath, gap engineering to suppress quasiparticles, Purcell filters to close the readout-loss channel, and flux sweet spots to evade 1/f1/f noise. The same filter-function machinery that explains spin echo now runs in reverse as noise spectroscopy — sweeping the spacing of dynamical-decoupling pulses reconstructs Sλ(ω)S_\lambda(\omega) and diagnoses which bath limits a device. And decoherence sets the target for quantum error correction: the surface code only helps once physical error rates, dominated by decoherence during gates and idling, fall below roughly one percent — the threshold the coherence budget of Worked Example 3 makes concrete. The open frontier is the low-frequency noise no one has fully explained: the microscopic origin of 1/f1/f flux and charge noise, and of the individual defects that occasionally swallow a qubit's T1T_1, remains an active and humbling problem.

The physical picture

Picture a qubit superposition α∣0⟩+β∣1⟩\alpha\lvert 0\rangle + \beta\lvert 1\rangle as a tiny clock whose hand points somewhere on the equator of the Bloch sphere. Two physically distinct things can go wrong. Relaxation, the T1T_1 process, is the clock's mainspring running down: the qubit emits its energy quantum into the bath and slides back toward the ground state. Dephasing, the T2T_2 process, is the clock running at a slightly wrong, fluctuating rate: the hand's angle becomes uncertain, and averaging over many noisy repetitions washes the equatorial component to zero even though no energy was lost. The first figure puts both decays side by side with real numbers.

Figure
Two clocks running down. T₁ relaxation slides the Bloch vector back toward the ground state while T₂ shrinks its equatorial length; with T₁ = 80 μs and T₂ = 50 μs the pure-dephasing time follows from 1/T₂ = 1/(2T₁) + 1/Tφ, giving Tφ ≈ 72.7 μs
Two clocks running down: T₁ empties the excited state, T₂ erases the phase |0⟩ ground |1⟩ excited coherence |ρ₀₁| shrinks at 1/T₂ relaxes at 1/T₁ 1 0 T₂ = 50 μs T₁ = 80 μs time t (μs) population e^(−t/T₁) coherence e^(−t/T₂) 1/T₂ = 1/(2T₁) + 1/Tφ ⟹ Tφ ≈ 72.7 μs

The crucial insight is frequency matching. To absorb or emit a quantum ℏω01\hbar\omega_{01}, the bath must have spectral weight at ω01\omega_{01} — so relaxation is set by the noise power at the qubit frequency, S(ω01)S(\omega_{01}). Pure dephasing needs no energy exchange; it is driven by slow wobbles of the transition frequency, so it is set by S(ω)S(\omega) near DC. This is why sluggish 1/f1/f noise — enormous weight at low frequency, almost none at gigahertz — dephases a qubit ferociously while barely relaxing it. The second figure unifies the two timescales: the noise spectral density itself, with both sampling points marked.

Figure
The spectral rule behind both timescales. The environmental noise power S(ω) is large at low frequency (1/f, slope −1) and flat at high frequency; energy relaxation samples S at the qubit frequency ω₀₁/2π ≈ 5 GHz while pure dephasing samples S as ω → 0, which is why slow 1/f noise dephases far more than it relaxes
One spectrum, two rates: Γ₁ reads S at ω₀₁, Γφ reads S at ω → 0 1 Hz 1 kHz 1 MHz 1 GHz frequency ω/2π (log scale) noise PSD S(ω) (log scale) S(ω) = A/ω (1/f, slope −1) white-noise floor ω → 0 samples S → Γφ = 1/Tφ ω₀₁/2π ≈ 5 GHz samples S(ω₀₁) → Γ₁ = 1/T₁

A sweet spot is the noise-engineering masterstroke this picture makes obvious. If the qubit frequency depends on a noisy control parameter λ\lambda — a gate charge, a flux — tune the bias to where ∂ω01/∂λ=0\partial\omega_{01}/\partial\lambda = 0. There the qubit is, to first order, blind to fluctuations of λ\lambda; only the weaker second-order curvature survives. The transmon does this for charge noise at large EJ/ECE_J/E_C, the flux qubit at its symmetry point. Everything this article develops is a way to flatten S(ω)S(\omega) where the qubit reads it, or move the qubit to read a quieter part of the spectrum.

The core equations

Start with the phenomenology. Independent of any microscopic model, the density matrix of a single qubit subject to relaxation and dephasing evolves so that populations relax at rate 1/T11/T_1 and coherences decay at rate 1/T21/T_2 on top of free precession at ω01\omega_{01}. This is the two-timescale skeleton that every experiment fits to.

ρ(t)=(ρ00+ρ11(0) ⁣(1−e−t/T1)ρ01(0) e−t/T2 e−iω01tρ10(0) e−t/T2 e+iω01tρ11(0) e−t/T1)(1)\rho(t)=\begin{pmatrix}\rho_{00}+\rho_{11}^{(0)}\!\left(1-e^{-t/T_1}\right) & \rho_{01}^{(0)}\,e^{-t/T_2}\,e^{-i\omega_{01}t}\\[4pt] \rho_{10}^{(0)}\,e^{-t/T_2}\,e^{+i\omega_{01}t} & \rho_{11}^{(0)}\,e^{-t/T_1}\end{pmatrix} \qquad (1)

where: ρij\rho_{ij} — density-matrix elements, diagonal entries are populations and off-diagonal entries are coherences; T1T_1 — longitudinal (energy) relaxation time; T2T_2 — transverse (coherence) decay time; ω01\omega_{01} — qubit transition frequency. Reading: the diagonal always sums to one, so the map is trace-preserving; the excited population decays with T1T_1, the superposition-encoding coherence with T2T_2. The two timescales are bound together because energy relaxation itself destroys coherence. Splitting the coherence decay into the part inherited from relaxation and the part from genuine phase noise gives the master relation of the subject.

1T2=12T1+1Tϕ(2)\frac{1}{T_2}=\frac{1}{2T_1}+\frac{1}{T_\phi} \qquad (2)

where: TϕT_\phi — pure-dephasing time, phase loss with no energy exchange. Reading: relaxation feeds the coherence-decay rate with a factor 12\tfrac12, so a perfectly relaxation-limited qubit (Tϕ→∞T_\phi\to\infty) has T2=2T1T_2 = 2T_1; any real dephasing only shortens T2T_2. Now the microscopics. Applying Fermi's golden rule to the coupling, the relaxation rate is the squared transition matrix element times the noise power at the qubit frequency — the only frequency at which the bath can absorb the energy quantum.

Γ1=1T1=1ℏ2 ∣⟨0∣A^∣1⟩∣2 Sλ(ω01)(3)\Gamma_1=\frac{1}{T_1}=\frac{1}{\hbar^2}\,\bigl\lvert\langle 0\rvert \hat{A}\lvert 1\rangle\bigr\rvert^2\,S_\lambda(\omega_{01}) \qquad (3)

where: A^\hat{A} — qubit operator coupling to the noisy variable λ\lambda; Sλ(ω)=∫⟨λ(t)λ(0)⟩ eiωt dtS_\lambda(\omega)=\int \langle\lambda(t)\lambda(0)\rangle\,e^{i\omega t}\,\mathrm dt — power spectral density of the noise; ω01\omega_{01} — qubit frequency. Reading: relaxation is resonant — it cares only about the bath's spectral weight at ω01\omega_{01}, nothing about DC. Pure dephasing is the complementary limit. Slow noise modulates the qubit splitting without flipping the state, so the relevant quantity is the frequency's sensitivity to the noisy parameter times the noise power near zero frequency.

Γϕ=1Tϕ=12ℏ2(∂H^sys∂λ)01 ⁣2 ⁣Sλ(ω→0)(4)\Gamma_\phi=\frac{1}{T_\phi}=\frac{1}{2\hbar^2}\left(\frac{\partial \hat H_{\rm sys}}{\partial\lambda}\right)^{\!2}_{01}\! S_\lambda(\omega\to 0) \qquad (4)

where: (∂H^sys/∂λ)01(\partial \hat H_{\rm sys}/\partial\lambda)_{01} — first-order sensitivity of the qubit splitting to λ\lambda; Sλ(ω→0)S_\lambda(\omega\to 0) — near-DC noise power. Reading: at a sweet spot the prefactor vanishes and first-order dephasing switches off, leaving only weak second-order curvature — the whole justification for sweet-spot design. When the noise is 1/f1/f rather than white, equation (4) fails because Sλ(ω→0)S_\lambda(\omega\to 0) diverges. Carrying the low-frequency weight through the Ramsey phase integral gives a decay that is Gaussian in time, not exponential, with only a weak logarithmic dependence on the averaging time.

⟨σx⟩Ramsey(t)=e−(t/Tϕ1/f)2,(Tϕ1/f)−2=A(∂ω01∂λ) ⁣2∣ln⁡(ωirt)∣(5)\langle\sigma_x\rangle_{\rm Ramsey}(t)=e^{-(t/T_\phi^{1/f})^2},\qquad \bigl(T_\phi^{1/f}\bigr)^{-2}=A\left(\frac{\partial\omega_{01}}{\partial\lambda}\right)^{\!2}\bigl\lvert\ln(\omega_{\rm ir}t)\bigr\rvert \qquad (5)

where: AA — 1/f1/f noise amplitude, Sλ(ω)=A/∣ω∣S_\lambda(\omega)=A/\lvert\omega\rvert; ωir\omega_{\rm ir} — infrared cutoff set by the total measurement time. Reading: the Gaussian envelope is the fingerprint of 1/f1/f dephasing, and its refocusability by echo is what makes dynamical decoupling effective. Two mechanisms dominate real superconducting devices. The first is Purcell decay: in the dispersive regime the qubit inherits a small photonic admixture from its lossy readout resonator and decays at the cavity rate weighted by the square of the coupling-to-detuning ratio.

ΓPurcell=κ(gΔ) ⁣2(6)\Gamma_{\rm Purcell}=\kappa\left(\frac{g}{\Delta}\right)^{\!2} \qquad (6)

where: κ\kappa — photon decay rate of the readout cavity; gg — qubit–cavity coupling; Δ=ω01−ωr\Delta=\omega_{01}-\omega_r — qubit–cavity detuning. Reading: the 1/Δ21/\Delta^2 scaling means detuning sharply lengthens the Purcell-limited T1T_1, at the cost of a weaker dispersive readout — the tension a Purcell filter relieves. The second is quasiparticle tunneling across the Josephson junction, which even at a density as low as xqp∼10−6x_{qp}\sim 10^{-6} can dominate T1T_1 in an otherwise clean device.

Γqp∝ω01π xqp 2Δscℏω01(7)\Gamma_{qp}\propto \frac{\omega_{01}}{\pi}\,x_{qp}\,\sqrt{\frac{2\Delta_{sc}}{\hbar\omega_{01}}} \qquad (7)

where: xqpx_{qp} — density of non-equilibrium quasiparticles normalized to the Cooper-pair density; Δsc\Delta_{sc} — superconducting gap. Reading: the rate is linear in xqpx_{qp} with a kinematic factor from the BCS density of states near the gap edge, so quasiparticle poisoning is fought by gap engineering and shielding, not pulse tricks.

Derivation ABloch-Redfield rates from a system-bath Hamiltonian

The goal is to derive the two master rates Γ1\Gamma_1 and Γϕ\Gamma_\phi from a single microscopic model — a qubit whose flip amplitude and whose energy splitting are each modulated by the same fluctuating bath — and to see why one rate reads the noise spectrum at the qubit frequency while the other reads it at DC. We start from a transverse-plus-longitudinal system–bath coupling, move to the interaction picture, apply the second-order Born–Markov master equation, and trace out the bath so that the noise spectral density SB(ω)S_B(\omega) appears explicitly in each rate.

Step 1. Write the qubit as H^0=ℏω012σz\hat H_0=\tfrac{\hbar\omega_{01}}{2}\sigma_z and couple it to a single fluctuating bath operator B^\hat B through two channels,

H^int=λ⊥ σx B^+λ∥ σz B^,\hat H_{\rm int}=\lambda_\perp\,\sigma_x\,\hat B+\lambda_\parallel\,\sigma_z\,\hat B,

a transverse term and a longitudinal term. Why: σx\sigma_x is the only operator that can flip the qubit between ∣0⟩\lvert0\rangle and ∣1⟩\lvert1\rangle, so it sources relaxation; σz\sigma_z commutes with H^0\hat H_0 and merely wobbles the energy gap, so it sources pure dephasing.

Step 2. Go to the interaction picture, where the transverse coupling acquires the qubit's own phase,

σx(t)=σ+eiω01t+σ−e−iω01t,\sigma_x(t)=\sigma_+e^{i\omega_{01}t}+\sigma_-e^{-i\omega_{01}t},

while σz(t)=σz\sigma_z(t)=\sigma_z is stationary. Why: the transverse term oscillates at ±ω01\pm\omega_{01}, so beaten against the bath correlation function it selects the spectral component at ω01\omega_{01}; the stationary longitudinal term selects the DC component.

Step 3. Insert the coupling into the second-order Born–Markov master equation,

ρ˙=−1ℏ2∫0∞dτ  TrB[H^int(t), [H^int(t−τ), ρ⊗ρB]],\dot\rho=-\frac{1}{\hbar^2}\int_0^\infty \mathrm d\tau\;\mathrm{Tr}_B\bigl[\hat H_{\rm int}(t),\,[\hat H_{\rm int}(t-\tau),\,\rho\otimes\rho_B]\bigr],

and trace over the bath, replacing each pair of bath operators by the correlation function whose Fourier transform is SB(ω)=∫⟨B^(τ)B^(0)⟩eiωτ dτS_B(\omega)=\int\langle\hat B(\tau)\hat B(0)\rangle e^{i\omega\tau}\,\mathrm d\tau. Why: the Born (weak-coupling) and Markov (short-memory) approximations are exactly what let the double commutator collapse into rates rather than a memory kernel.

Step 4. Collect the terms proportional to σ±\sigma_\pm; they transport population between the levels at a rate

Γ1=λ⊥2ℏ2 SB(ω01),\Gamma_1=\frac{\lambda_\perp^2}{\hbar^2}\,S_B(\omega_{01}),

the microscopic content of equation (3). Why: the oscillating transverse coupling filters the bath at ω01\omega_{01}, so relaxation literally measures the bath's ability to exchange one energy quantum — pure resonance.

Step 5. Collect the terms proportional to σz\sigma_z; they leave populations untouched but decay the coherence at

Γϕ=λ∥22ℏ2 SB(0),\Gamma_\phi=\frac{\lambda_\parallel^2}{2\hbar^2}\,S_B(0),

the microscopic content of equation (4), with the characteristic factor 12\tfrac12 relative to the naive expectation absorbed into the definition. Why: the stationary longitudinal coupling filters the bath at zero frequency, so pure dephasing measures only the slow, near-DC fluctuations that jitter the phase without moving energy. *The resulting map reduces to the standard Lindblad form with jump operators σ−\sigma_- (rate Γ1\Gamma_1) and σz\sigma_z (rate Γϕ/2\Gamma_\phi/2), and identifying the two channels reproduces 1/T2=1/(2T1)+1/Tϕ1/T_2=1/(2T_1)+1/T_\phi exactly.*

Go deeper · Derivation B · 1 minWhy T2 cannot exceed twice T1
Derivation BWhy T2 cannot exceed twice T1

The goal is to prove that any physically legitimate single-qubit channel combining relaxation and dephasing must obey T2≤2T1T_2\le 2T_1 — that the inequality is not a modelling convenience but a hard consequence of the requirement that a quantum operation be completely positive and trace-preserving. We model the combined amplitude-plus-phase damping map, impose complete positivity through its Choi matrix, and read off the inequality between the decay rates before translating back into times.

Step 1. Write the combined channel on its populations and coherences,

ρ11(t)=ρ11(0)e−t/T1,ρ01(t)=ρ01(0)e−t/T2,\rho_{11}(t)=\rho_{11}^{(0)}e^{-t/T_1},\qquad \rho_{01}(t)=\rho_{01}^{(0)}e^{-t/T_2},

folding amplitude damping (T1T_1) and phase damping (TϕT_\phi) into one map. Why: this is the most general Markovian relaxation-plus-dephasing evolution up to unitary rotation, so a bound proved here binds every real qubit.

Step 2. Demand complete positivity. Requiring the map's Choi matrix to be positive semidefinite constrains how fast the off-diagonal may decay relative to the diagonal: the coherence cannot decay faster than the bound set by the population decay, or the map produces negative probabilities on part of an entangled input. Why: complete positivity, not mere positivity, is the correct requirement because the qubit may be entangled with a spectator system the channel must also respect.

Step 3. Working out that positivity condition yields the rate inequality

Γ2≥12 Γ1,\Gamma_2\ge \tfrac12\,\Gamma_1,

which is automatically satisfied by equation (2) since Γ2=12Γ1+Γϕ\Gamma_2=\tfrac12\Gamma_1+\Gamma_\phi with Γϕ≥0\Gamma_\phi\ge 0. Why: pure dephasing can only add to the coherence-decay rate, never subtract, so the smallest possible Γ2\Gamma_2 is 12Γ1\tfrac12\Gamma_1, reached when there is no pure dephasing at all.

Step 4. Invert to times,

T2=1Γ2≤112Γ1=2T1,T_2=\frac{1}{\Gamma_2}\le\frac{1}{\tfrac12\Gamma_1}=2T_1,

with equality if and only if Γϕ=0\Gamma_\phi=0. Why: this pins the ceiling on coherence — the best a given T1T_1 buys is T2=2T1T_2=2T_1, by pure amplitude damping, and any measured excess signals an experimental error. *Pure amplitude damping saturates the bound at T2=2T1T_2=2T_1, and adding any Γϕ>0\Gamma_\phi>0 strictly lowers T2T_2, confirming the inequality is tight.*

Go deeper · Derivation C · 1 minHow a spin echo refocuses low-frequency noise
Derivation CHow a spin echo refocuses low-frequency noise

The goal is to show quantitatively why a single π\pi-pulse placed halfway through a free-evolution interval — the Hahn echo — cancels the dephasing from slow, quasi-static noise while leaving relaxation untouched, expressed through the sequence's filter function. We write the accumulated phase as the noise weighted by a sequence-dependent filter, compute that filter for the Ramsey and echo sequences, and compare their low-frequency behaviour where 1/f1/f noise lives.

Step 1. For Gaussian frequency noise the coherence decays as e−⟨ϕ2⟩/2e^{-\langle\phi^2\rangle/2}, and the accumulated phase variance is the noise spectrum weighted by a filter function set by the pulse sequence,

⟨ϕ2⟩=(∂ω01∂λ) ⁣2∫0∞dωπ Sλ(ω) ∣F(ωt)∣2ω2.\langle\phi^2\rangle=\left(\frac{\partial\omega_{01}}{\partial\lambda}\right)^{\!2}\int_0^\infty\frac{\mathrm d\omega}{\pi}\,S_\lambda(\omega)\,\frac{\lvert F(\omega t)\rvert^2}{\omega^2}.

Why: because the phase is a linear functional of the noise, every pulse sequence acts as a frequency filter on Sλ(ω)S_\lambda(\omega), and coherence loss is just that filtered noise power.

Step 2. For a free Ramsey experiment the filter is FRamsey(ωt)=sin⁡(ωt/2)F_{\rm Ramsey}(\omega t)=\sin(\omega t/2), so ∣F∣2/ω2→t2/4\lvert F\rvert^2/\omega^2\to t^2/4 as ω→0\omega\to0. Why: the Ramsey filter passes DC with full weight, so slow 1/f1/f noise contributes maximally — exactly the regime where equation (5) gives Gaussian decay.

Step 3. Inserting a π\pi-pulse at the midpoint flips the sign of the accumulated phase for the second half, and the filter becomes Fecho(ωt)∝sin⁡2(ωt/4)F_{\rm echo}(\omega t)\propto \sin^2(\omega t/4), which vanishes as (ωt)2(\omega t)^2 when ω→0\omega\to0. Why: the sign flip makes any frequency offset that is constant over the sequence cancel exactly, so the filter now rejects DC quadratically instead of passing it.

Step 4. The consequence: static and slowly varying frequency errors are refocused, so T2echo≫T2RamseyT_2^{\rm echo}\gg T_2^{\rm Ramsey} whenever 1/f1/f noise dominates, while noise at ω01\omega_{01} — and hence T1T_1 — is untouched. Iterating to NN pulses (CPMG) pushes the filter passband to ever higher frequency, the basis of dynamical decoupling and noise spectroscopy. Why: echo attacks only the low-frequency dephasing channel a π\pi-pulse can reach; it can never fix relaxation. *The echo filter's quadratic suppression at DC is precisely what immunizes the sequence against quasi-static dephasing while leaving relaxation unchanged.*

Limiting cases

The interplay of T1T_1 and T2T_2 has a handful of clean regimes that between them cover almost every real qubit.

CaseResultMeaning
Pure relaxation, no extra dephasingTϕ→∞  ⇒  T2=2T1T_\phi \to \infty \;\Rightarrow\; T_2 = 2T_1Decoherence is energy relaxation only; coherence decays at half the relaxation rate, the maximum T2T_2 a given T1T_1 allows.
Dephasing-dominatedTϕ≪T1  ⇒  T2≈TϕT_\phi \ll T_1 \;\Rightarrow\; T_2 \approx T_\phiTypical of 1/f1/f-limited superconducting qubits and field-sensitive ion transitions; phase scrambles before energy decays.
White noise (Markovian)S(ω)=const  ⇒  ⟨σx⟩(t)∝e−t/T2S(\omega)=\text{const} \;\Rightarrow\; \langle\sigma_x\rangle(t)\propto e^{-t/T_2}Flat spectrum gives exponential, memoryless decay — the Lindblad regime where Bloch–Redfield is exact.
1/f1/f noiseS(ω)=A/∣ω∣  ⇒  ⟨σx⟩(t)∝e−(t/T2)2S(\omega)=A/\lvert\omega\rvert \;\Rightarrow\; \langle\sigma_x\rangle(t)\propto e^{-(t/T_2)^2}Gaussian Ramsey decay; the soft low-frequency weight is what makes echo and dynamical decoupling effective.
Sweet spot∂ω01/∂λ=0  ⇒  Γϕ∝(∂2ω01/∂λ2)2\partial\omega_{01}/\partial\lambda=0 \;\Rightarrow\; \Gamma_\phi\propto (\partial^2\omega_{01}/\partial\lambda^2)^2First-order dephasing vanishes; only weak curvature survives — the design principle of the transmon and flux qubit.

Worked examples

1 · Extracting the pure-dephasing time (confidence builder)

A superconducting qubit is measured to have T1=80 μsT_1 = 80\ \mu\text{s} and, in a Hahn-echo experiment, T2=50 μsT_2 = 50\ \mu\text{s}. Find the pure-dephasing time TϕT_\phi and check that the numbers are self-consistent.

Step 1. Rearrange the master relation (2) for TϕT_\phi:

1Tϕ=1T2−12T1=150−1160=160−508000=1108000 μs−1.\frac{1}{T_\phi}=\frac{1}{T_2}-\frac{1}{2T_1}=\frac{1}{50}-\frac{1}{160}=\frac{160-50}{8000}=\frac{110}{8000}\ \mu\text{s}^{-1}.

Step 2. Invert: Tϕ=8000/110≈72.7 μsT_\phi = 8000/110 \approx \mathbf{72.7\ \mu s}. The dephasing time is comparable to T1T_1, so this qubit is neither purely relaxation-limited nor catastrophically dephasing-dominated.

Step 3. Consistency check: the ceiling is 2T1=160 μs2T_1 = 160\ \mu\text{s}, and the measured T2=50 μsT_2 = 50\ \mu\text{s} sits comfortably below it, so the result is physical. Tϕ≈72.7 μsT_\phi \approx 72.7\ \mu s, and T2=50<2T1=160 μsT_2 = 50 < 2T_1 = 160\ \mu s, so the measurement is self-consistent.

2 · Purcell-limited relaxation and the value of detuning

A transmon couples to its readout resonator with g/2π=100 MHzg/2\pi = 100\ \text{MHz}, the cavity linewidth is κ/2π=1 MHz\kappa/2\pi = 1\ \text{MHz}, and the qubit sits Δ/2π=1 GHz\Delta/2\pi = 1\ \text{GHz} below the cavity. What Purcell-limited T1T_1 does this imply, and how does it change if the detuning is doubled to 2 GHz2\ \text{GHz}?

Step 1. The coupling-to-detuning ratio is g/Δ=100 MHz/1 GHz=0.1g/\Delta = 100\ \text{MHz}/1\ \text{GHz} = 0.1, so (g/Δ)2=0.01(g/\Delta)^2 = 0.01. The 2π2\pi factors cancel because both are frequencies.

Step 2. Apply equation (6): ΓPurcell=κ (g/Δ)2=(2π×106 s−1)(0.01)=6.28×104 s−1\Gamma_{\rm Purcell}=\kappa\,(g/\Delta)^2 = (2\pi\times 10^6\ \text{s}^{-1})(0.01) = 6.28\times 10^4\ \text{s}^{-1}. Invert to get the Purcell-limited relaxation time T1Purcell=1/ΓPurcell≈15.9 μsT_1^{\rm Purcell}=1/\Gamma_{\rm Purcell} \approx \mathbf{15.9\ \mu s}.

Step 3. Because ΓPurcell∝1/Δ2\Gamma_{\rm Purcell}\propto 1/\Delta^2, doubling the detuning cuts the rate by four and quadruples the time: T1Purcell→63.7 μsT_1^{\rm Purcell}\to \mathbf{63.7\ \mu s}. Detuning from 1 to 2 GHz lifts the Purcell-limited T1T_1 from 15.9 to 63.7 μs, at the price of a four-times weaker dispersive shift and slower readout — the trade-off a Purcell filter removes.

3 · A coherence budget for gates

A qubit has T2=100 μsT_2 = 100\ \mu\text{s} and single-qubit gates take tg=25 nst_g = 25\ \text{ns}. Roughly how many coherent operations fit inside a coherence time, and what error per gate does that imply?

Step 1. The number of gates before coherence is lost is N≈T2/tg=(100×10−6)/(25×10−9)=4000N \approx T_2/t_g = (100\times 10^{-6})/(25\times 10^{-9}) = \mathbf{4000}.

Step 2. The corresponding decoherence error accrued during one gate is ε≈tg/T2=25 ns/100 μs=2.5×10−4\varepsilon \approx t_g/T_2 = 25\ \text{ns}/100\ \mu\text{s} = \mathbf{2.5\times 10^{-4}}.

Step 3. This is the quantity the DiVincenzo criteria demand be small, and it sits below the surface-code fault-tolerance threshold (a few ×10−3\times 10^{-3}). About 4000 coherent operations per T2T_2, a decoherence floor of 2.5×10−42.5\times 10^{-4} per gate — which is why longer T2T_2 translates directly into deeper circuits.

Physics of the portrait

Decoherence is one of those ideas that was implicit in quantum mechanics for decades before anyone named it. The relaxation times came first, from magnetic resonance; the master equations came next, from the theory of open systems; and the recognition that entanglement with an environment is what turns quantum into classical — and what limits every qubit — was the work of a small group of theorists who insisted on taking the environment seriously rather than wishing it away.

Origins · 1 minWojciech H. Zurek and how this was discovered
WH
Wojciech H. Zurek
b. 1951 · decoherence, einselection, and the quantum-to-classical transition

Wojciech Hubert Zurek, a Polish-American theoretical physicist at Los Alamos, did more than anyone to turn decoherence from a footnote into a research program. Building on Hans Dieter Zeh's 1970 insight that a macroscopic system is never isolated and is therefore continuously monitored by its surroundings, Zurek argued through the early 1980s that the environment does not merely disturb a quantum system — it actively selects a preferred set of stable "pointer states" and destroys superpositions between them, a process he named einselection (environment-induced superselection). This is why we never see a chair in two places at once, and, in the language of this article, why a qubit's off-diagonal density-matrix elements decay: the environment is measuring the qubit whether we want it to or not, at a rate set by the noise spectral density. Zurek later extended the picture to quantum Darwinism, and his 2003 Reviews of Modern Physics synthesis remains the standard reference. For the quantum engineer his message is bracing: decoherence is not a defect to be patched but the generic fate of any system coupled to a bath, and fighting it means engineering the coupling itself.

History
1946
— Bloch introduces the relaxation times T1 and T2 to describe nuclear magnetic resonance
1957
— Redfield derives the master equation for a spin weakly coupled to a bath
1970
— Zeh argues that environmental entanglement, not wavefunction collapse, produces classical behaviour
1981
— Zurek formalizes einselection and the pointer basis, naming the decoherence program
1999
— Nakamura, Pashkin and Tsai observe coherent oscillations in a superconducting charge qubit

Going further

Texts. Nielsen & Chuang, Quantum Computation and Quantum Information, Ch. 8 (quantum operations and noise) · Preskill, Lecture Notes on Quantum Computation, Ch. 3 (open systems and decoherence) · Watrous, The Theory of Quantum Information, Ch. 2 (channels and CPTP maps) · Wilde, Quantum Information Theory, Ch. 4 (quantum channels) · Breuer & Petruccione, The Theory of Open Quantum Systems, Ch. 3 (Markovian master equations) · Krantz et al., "A Quantum Engineer's Guide to Superconducting Qubits", Appl. Phys. Rev. 6, 021318 (2019) · Zurek, "Decoherence, einselection, and the quantum origins of the classical", Rev. Mod. Phys. 75, 715 (2003). Key results: Worked Example 1 gives Tϕ≈72.7 μT_\phi\approx 72.7\ \mus with T2<2T1T_2<2T_1; Example 2 gives T1Purcell=15.9→63.7 μT_1^{\rm Purcell}=15.9\to 63.7\ \mus under Δ→2Δ\Delta\to 2\Delta; Example 3 gives N≈4000N\approx 4000 gates and ε≈2.5×10−4\varepsilon\approx 2.5\times10^{-4}; the Ramsey and echo filter scalings go as ∝t2/4\propto t^2/4 and ∝(ωt)2\propto(\omega t)^2 respectively. No platform-specific record coherence numbers are quoted, as state-of-the-art values evolve.

Prerequisites. DiVincenzo Criteria · Density matrices and open systems · Bloch sphere and qubit dynamics. Leads to. The Transmon Qubit · Dispersive Readout · Sideband Cooling · Quantum error correction.

Ask the tutor: Give me your qubit's T1T_1, T2T_2, and pulse sequence, and I will tell you which bath is limiting you and whether echo, a sweet spot, or a Purcell filter is your best next move.

What is the difference between T1 and T2, in one breath?

T1T_1 is about energy, T2T_2 is about phase. T1T_1 measures how long the qubit stays excited before dumping its energy quantum into the environment and relaxing toward the ground state — a change in the populations, the diagonal of the density matrix. T2T_2 measures how long a superposition keeps a definite phase before the relative phase between ∣0⟩\lvert0\rangle and ∣1⟩\lvert1\rangle is scrambled — a decay of the off-diagonal coherences. Energy loss inevitably costs phase too, which is why T1T_1 feeds into T2T_2 through 1/T2=1/(2T1)+1/Tϕ1/T_2=1/(2T_1)+1/T_\phi; the term TϕT_\phi is the phase you lose without losing any energy.

Why exactly is there a factor of 2 in 1/T2 = 1/(2T1) + 1/Tφ?

Because coherence is an amplitude, and amplitudes decay at half the rate of the probabilities they come from. A coherence ρ01\rho_{01} is built from the amplitudes of being in ∣0⟩\lvert0\rangle and in ∣1⟩\lvert1\rangle; when relaxation depletes the excited population at rate Γ1\Gamma_1, the coherence tracks the geometric mean of the two survival probabilities, which decays at half the exponent — hence Γ1/2\Gamma_1/2. Pure dephasing then adds its rate Γϕ\Gamma_\phi directly, since it acts on the phase and not the populations, giving Γ2=12Γ1+Γϕ\Gamma_2=\tfrac12\Gamma_1+\Gamma_\phi.

If a bath has noise everywhere, why does relaxation only care about ω₀₁?

Because relaxation is an energy-conserving transition, and energy conservation is a resonance condition. For the qubit to fall from ∣1⟩\lvert1\rangle to ∣0⟩\lvert0\rangle it must hand a quantum ℏω01\hbar\omega_{01} to the bath, which can only accept it with spectral weight at exactly that frequency. Fermi's golden rule makes this quantitative: the rate is the matrix element squared times S(ω01)S(\omega_{01}), and nothing else. Noise at other frequencies is not idle — near DC it drives dephasing — but it cannot cause the energy-exchanging transition relaxation requires.

Can dynamical decoupling fix everything, or just dephasing?

Just dephasing. A spin echo or CPMG sequence flips the qubit so slow frequency errors accumulated in the first half are cancelled in the second; its filter function rejects the low-frequency noise where dephasing lives. But relaxation is caused by noise at ω01\omega_{01}, far above the frequencies a pulse sequence can reach, and no amount of π\pi-pulsing undoes an energy quantum that has already left. So dynamical decoupling can push T2T_2 toward its 2T12T_1 ceiling but never past it — for that you must improve T1T_1 itself.

Is decoherence the same thing as wavefunction collapse?

No, though they are often conflated. Decoherence is a unitary, entirely standard quantum process: the qubit becomes entangled with its environment, and once we trace out the environment the qubit's reduced density matrix loses its off-diagonal terms and looks like a classical mixture. Nothing has "collapsed" — the global state is still pure, the coherence has simply leaked into correlations we cannot track. This is Zurek's point: decoherence explains why we see apparently classical, definite outcomes without invoking any new collapse postulate. Whether that fully resolves the measurement problem is interpretation, but the physics of the qubit's lost coherence is unambiguous and computable.

Do photonic qubits decohere the same way superconducting ones do?

Not really — the dominant channel is different. A flying photon does not sit near a bath of two-level defects jittering its energy, so it barely "dephases" in the T2T_2 sense. Its enemy is loss: absorption, scattering, imperfect coupling into fibres and chips, and finite detector efficiency all remove the photon entirely rather than randomize its phase. In linear-optical quantum computing this shows up as failed heralding and reduced success probability, not phase errors, which is why the field invests in low-loss waveguides, high-efficiency nanowire detectors, and loss-tolerant codes rather than dynamical decoupling. The unifying view is still spectral: loss is a relaxation channel where the "bath" is the set of modes the photon can escape into.

Before this

Where this leads