Skip to content

Decoherence Timescales

A decoherence timescale is a time extracted from the decay of an interference signal, an off-diagonal density-matrix element, a pointer-state overlap, or a visibility envelope. It is not a universal property of a device or object by itself. It is a property of a model, a basis or pair of alternatives, a protocol, an environment, and a threshold convention.

The practical question is:

how long until the relevant interference is negligible?\text{how long until the relevant interference is negligible?}

The answer may be a Markovian exponential time T2T_2, a Ramsey inhomogeneous time T2∗T_2^*, a pure-dephasing time TϕT_\phi, a spatial decoherence time tdec(Δx)t_{\mathrm{dec}}(\Delta x), or a protocol-specific envelope time. A trustworthy estimate always states which one is meant.

Choose the coherence whose loss matters. For a two-level system in a chosen basis,

ρ=(ρ00ρ01ρ10ρ11),\rho = \begin{pmatrix} \rho_{00}&\rho_{01}\\ \rho_{10}&\rho_{11} \end{pmatrix},

the relevant coherence may be ρ01\rho_{01}. A simple exponential model writes

ρ01(t)=e−iω0te−t/T2ρ01(0).\rho_{01}(t) = e^{-i\omega_0t} e^{-t/T_2} \rho_{01}(0).

Then T2T_2 is the time at which the magnitude of this off-diagonal term has fallen by e−1e^{-1}:

∣ρ01(T2)∣=e−1∣ρ01(0)∣.|\rho_{01}(T_2)| = e^{-1} |\rho_{01}(0)|.

For a position superposition, the relevant object might instead be an off-diagonal position-space element:

ρ(x,x′;t).\rho(x,x';t).

For a matter-wave interferometer, it may be fringe visibility. For a detector pointer, it may be the overlap of environmental record states. The first step in estimating a timescale is therefore not plugging numbers into a formula; it is identifying the coherence or visibility whose decay is physically relevant.

The same experiment can involve several distinct times.

SymbolMeaningTypical diagnostic
T1T_1population relaxation timeenergy populations approach equilibrium
T2T_2transverse coherence decay timeoff-diagonal energy-basis coherence decays
TϕT_\phipure-dephasing timeextra phase decay after separating T1T_1
T2∗T_2^*inhomogeneous Ramsey decay timeensemble or slow-noise dephasing without echo
T2,echoT_{2,\mathrm{echo}}echo coherence timedecay after a spin-echo pulse
tdec(Δx)t_{\mathrm{dec}}(\Delta x)spatial decoherence timeposition coherence at separation Δx\Delta x decays

These are fitted or inferred quantities, not microscopic mechanisms by themselves. Saying “the decoherence time is 10 μs10\,\mu\mathrm{s}” is incomplete unless the basis, protocol, and envelope are specified.

For a weakly coupled qubit with energy relaxation and independent pure dephasing, the standard bookkeeping relation is

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

This relation is one of the most useful quick estimates in spectroscopy and quantum-device work. It says:

  • population relaxation reduces coherence at half the population-relaxation rate;
  • pure dephasing adds directly to transverse coherence decay;
  • T2T_2 cannot exceed 2T12T_1 in this simple model unless the assumptions or fitted definitions differ.

Equivalently, with relaxation rate Γ1=1/T1\Gamma_1=1/T_1 and pure-dephasing rate Γϕ=1/Tϕ\Gamma_\phi=1/T_\phi,

Γ2=Γ12+Γϕ,T2=Γ2−1.\Gamma_2 = \frac{\Gamma_1}{2} + \Gamma_\phi, \qquad T_2=\Gamma_2^{-1}.

This is a model relation, not a law of nature. It assumes a clear two-level subspace, weak coupling, approximately Markovian noise, a well-defined energy basis, and separable relaxation and dephasing channels. For the generator-level derivations, see Pure Dephasing Master Equation and Amplitude Damping Master Equation.

At zero temperature, amplitude damping gives

ρ˙11=−Γ1ρ11,ρ˙01=−Γ12ρ01\dot\rho_{11} = -\Gamma_1\rho_{11}, \qquad \dot\rho_{01} = -\frac{\Gamma_1}{2}\rho_{01}

in the interaction picture. The excited population decays as

ρ11(t)=e−Γ1tρ11(0),\rho_{11}(t) = e^{-\Gamma_1t}\rho_{11}(0),

while coherence involving the excited state decays as

ρ01(t)=e−Γ1t/2ρ01(0).\rho_{01}(t) = e^{-\Gamma_1t/2}\rho_{01}(0).

Pure dephasing with rate Γϕ\Gamma_\phi contributes

ρ01(t)⟼e−Γϕtρ01(t)\rho_{01}(t) \longmapsto e^{-\Gamma_\phi t}\rho_{01}(t)

without changing ρ00\rho_{00} or ρ11\rho_{11}. Multiplying the two independent coherence factors gives

ρ01(t)=e−(Γ1/2+Γϕ)tρ01(0),\rho_{01}(t) = e^{-(\Gamma_1/2+\Gamma_\phi)t} \rho_{01}(0),

so

1T2=Γ12+Γϕ.\frac{1}{T_2} = \frac{\Gamma_1}{2} + \Gamma_\phi.

The half-rate is easy to forget because populations and amplitudes decay differently. A finite-time amplitude-damping channel makes the same point: the population survival factor is e−t/T1e^{-t/T_1}, while the coherence factor is its square root, e−t/(2T1)e^{-t/(2T_1)}.

At finite temperature, upward and downward transition rates both contribute to T1T_1:

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

The transverse decay rate in the same simple Markovian model is

1T2=Γ↑+Γ↓2+Γϕ.\frac{1}{T_2} = \frac{\Gamma_\uparrow+\Gamma_\downarrow}{2} + \Gamma_\phi.

Thermal excitation, spontaneous emission, engineered reservoirs, and leakage can all change what should be included in the measured relaxation rate. Before using the formula, identify the modeled subspace and decide whether leakage outside it is part of T1T_1, a separate loss channel, or a failure of the two-level description.

Many experiments report a Ramsey decay time T2∗T_2^* rather than the Markovian free-induction T2T_2. The star often signals inhomogeneous dephasing: slow detuning noise, static disorder across an ensemble, calibration drift, or shot-to-shot variation.

For a qubit with stochastic detuning ξ(t)\xi(t),

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

the Ramsey coherence envelope is

W(T)=⟨e−i∫0Tξ(t) dt⟩.W(T) = \left\langle e^{-i\int_0^T\xi(t)\,dt} \right\rangle.

If ξ\xi is quasi-static during each run and Gaussian distributed between runs with variance σ2\sigma^2, then

W(T)=e−σ2T2/2.W(T) = e^{-\sigma^2T^2/2}.

This is a Gaussian envelope, not an exponential Markovian one. The threshold time defined by ∣W(T∗)∣=e−1|W(T_*)|=e^{-1} is

T∗=2σ.T_*=\frac{\sqrt2}{\sigma}.

A spin echo can refocus static or slow components by reversing the sign of phase accumulation halfway through the sequence. This distinction is especially important for One-Over-F Noise, where the spectrum is concentrated near low frequencies. Echo therefore often gives a longer time than Ramsey:

T2,echo>T2∗T_{2,\mathrm{echo}} \gt T_2^*

when low-frequency noise dominates. Dynamical decoupling extends this logic by using pulse sequences to filter different noise bands. For the filter-function formalism and a reproducible notebook contract, see Dynamical Decoupling and Decoherence Timescale Estimation.

For classical Gaussian frequency noise, a common filter-function convention writes

W(T)=e−χ(T).W(T)=e^{-\chi(T)}.

With modulation function y(t)y(t) and

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

one two-sided spectral convention gives

χ(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.

The coherence time may then be defined by the threshold

χ(T∗)=1.\chi(T_*)=1.

This framework makes the protocol dependence explicit. Ramsey, echo, CPMG, and other sequences have different Y(ω,T)Y(\omega,T) and therefore different sensitivity to the same noise spectrum. It also makes convention checks essential:

  • Is SξξS_{\xi\xi} one-sided or two-sided?
  • Are frequencies angular frequencies ω\omega or ordinary frequencies ff?
  • Does ξ(t)\xi(t) have units of angular frequency?
  • What ultraviolet and infrared cutoffs are used?
  • Is the reported time a threshold time or a fit parameter?

For the bath side of this language, see Noise Spectra.

For spatial superpositions, a common reduced description is

ρ(x,x′;t)≈e−F(x−x′)tρ(x,x′;0),\rho(x,x';t) \approx e^{-F(x-x')t} \rho(x,x';0),

where F(Δx)F(\Delta x) is a separation-dependent decoherence rate. For small separations compared with the relevant environmental wavelength, one often obtains

F(Δx)≃Λ∣Δx∣2.F(\Delta x) \simeq \Lambda|\Delta x|^2.

The e−1e^{-1} time for coherence between two positions separated by Δx\Delta x is then

tdec(Δx)=1Λ∣Δx∣2.t_{\mathrm{dec}}(\Delta x) = \frac{1}{\Lambda|\Delta x|^2}.

The quadratic scaling means that doubling the separation reduces the decoherence time by a factor of four in this regime. For larger separations, a single scattering event may resolve the alternatives, and F(Δx)F(\Delta x) can saturate near the scattering rate. The correct interpolation depends on wavelengths, cross sections, scatterer density, velocity distribution, geometry, and what information the environment retains.

This is the central lesson of environment-induced spatial decoherence: macroscopic separation can make environmental records highly distinguishable, so position coherence can disappear much faster than energy relaxation or mechanical damping.

If environmental records arrive as independent fragments, coherences multiply. Suppose each fragment reduces the coherence between alternatives aa and bb by an average overlap ηab\eta_{ab}, and fragments arrive at rate RR. A simple Poisson estimate gives

Dab(t)=exp⁡[Rt(ηab−1)].D_{ab}(t) = \exp[ Rt(\eta_{ab}-1) ].

The approximate decoherence rate is

Γab∼R(1−Re⁡ηab),\Gamma_{ab} \sim R \left( 1-\operatorname{Re}\eta_{ab} \right),

so the e−1e^{-1} time is roughly

tdec∼1Γab.t_{\mathrm{dec}} \sim \frac{1}{\Gamma_{ab}}.

This estimate is useful for intuition, but it hides angular averages, phase shifts, spectra, and correlations between fragments. It should be replaced by a scattering calculation or a microscopic master equation when numerical accuracy matters.

A useful hierarchy is:

tgate,tmeas,tdec,T1,tthermal,trec.t_{\mathrm{gate}}, \quad t_{\mathrm{meas}}, \quad t_{\mathrm{dec}}, \quad T_1, \quad t_{\mathrm{thermal}}, \quad t_{\mathrm{rec}}.

Here tgatet_{\mathrm{gate}} may be a control or interaction time, tmeast_{\mathrm{meas}} a measurement duration, tdect_{\mathrm{dec}} a coherence-loss time for selected alternatives, T1T_1 an energy-relaxation time, tthermalt_{\mathrm{thermal}} a thermalization time, and trect_{\mathrm{rec}} a recurrence or recoherence time. Their ordering depends on the system.

Examples:

  • In a good qubit, gates must be much shorter than the relevant coherence envelope.
  • In a measurement apparatus, pointer decoherence should be much faster than record readout.
  • In macroscopic spatial decoherence, tdect_{\mathrm{dec}} can be far shorter than mechanical damping.
  • In quantum sensing, the optimal interrogation time can be close to a protocol-specific coherence time but also depends on readout, duty cycle, and signal waveform.

The phrase “decoheres before it relaxes” usually means T2≪T1T_2\ll T_1 or tdec≪T1t_{\mathrm{dec}}\ll T_1. It does not mean that energy exchange is impossible, only that the interference signal vanishes on a shorter timescale than population relaxation.

For a practical calculation:

  1. Specify the system subspace and the alternatives whose interference matters.
  2. Choose the representation of coherence: matrix element, visibility, environmental overlap, or filter-function envelope.
  3. Identify whether the process is relaxation, pure dephasing, spatial decoherence, measurement-induced dephasing, or a mixture.
  4. Choose a model: Markovian master equation, finite-time channel, stochastic phase model, scattering calculation, or experimental fit.
  5. State the threshold or fit convention: e−1e^{-1} time, half-visibility time, exponential fit, Gaussian fit, stretched exponential, or protocol-specific time.
  6. Report the assumptions: basis, temperature, noise spectrum, pulse sequence, cutoffs, averaging, and whether the environment record is ignored or conditioned on.

This workflow is more important than a memorized formula. The wrong formula can easily give a precise-looking but physically meaningless time.

  • Reporting a single “decoherence time” without basis, protocol, or envelope.
  • Treating T2T_2, T2∗T_2^*, echo T2T_2, and dynamically decoupled coherence times as interchangeable.
  • Inferring TϕT_\phi from T1T_1 and T2T_2 outside the Markovian qubit bookkeeping model.
  • Forgetting that relaxation contributes 1/(2T1)1/(2T_1), not 1/T11/T_1, to transverse coherence decay.
  • Comparing a Gaussian Ramsey time directly to an exponential T2T_2.
  • Ignoring one-sided versus two-sided spectrum conventions.
  • Using a spatial-decoherence formula without specifying the separation Δx\Delta x.
  • Treating a long T1T_1 as proof of long coherence when low-frequency dephasing dominates.
  • Calling a fitted time a microscopic mechanism.
  • Forgetting that conditioned trajectories and unconditional ensemble averages can have different coherence behavior.
  • A. Abragam, The Principles of Nuclear Magnetism, Oxford University Press (1961).
  • C. P. Slichter, Principles of Magnetic Resonance, Springer, 3rd ed. (1990).
  • E. Joos and H. D. Zeh, “The emergence of classical properties through interaction with the environment,” Zeitschrift für Physik B 59, 223-243 (1985).
  • W. H. Zurek, “Decoherence, einselection, and the quantum origins of the classical,” Reviews of Modern Physics 75, 715-775 (2003).
  • M. Schlosshauer, Decoherence and the Quantum-to-Classical Transition, Springer (2007).
  • K. Hornberger, “Introduction to decoherence theory,” in Entanglement and Decoherence, Lecture Notes in Physics 768, Springer (2009).
  • H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).
  • L. Cywiński, R. M. Lutchyn, C. P. Nave, and S. Das Sarma, “How to enhance dephasing time in superconducting qubits,” Physical Review B 77, 174509 (2008).
  1. Pure-dephasing bookkeeping. A qubit has T1=30 μsT_1=30\,\mu\mathrm{s} and an exponential transverse decay time T2=20 μsT_2=20\,\mu\mathrm{s}. In the standard Markovian bookkeeping model, what is TϕT_\phi?
Solution

Use

1Tϕ=1T2−12T1.\frac{1}{T_\phi} = \frac{1}{T_2} - \frac{1}{2T_1}.

Then

1Tϕ=120 μs−160 μs=260 μs.\frac{1}{T_\phi} = \frac{1}{20\,\mu\mathrm{s}} - \frac{1}{60\,\mu\mathrm{s}} = \frac{2}{60\,\mu\mathrm{s}}.

Thus

Tϕ=30 μs.T_\phi=30\,\mu\mathrm{s}.
  1. Impossible inferred rate. A fit gives T1=10 μsT_1=10\,\mu\mathrm{s} and T2=25 μsT_2=25\,\mu\mathrm{s}. What happens if you apply the standard relation?
Solution

The inferred pure-dephasing rate would be

1Tϕ=125 μs−120 μs<0.\frac{1}{T_\phi} = \frac{1}{25\,\mu\mathrm{s}} - \frac{1}{20\,\mu\mathrm{s}} \lt0.

A negative pure-dephasing rate is unphysical in this model. The likely conclusion is not “negative dephasing,” but that at least one assumption or label is wrong: the fitted T2T_2 may not be a free exponential transverse decay time, the T1T_1 measurement may refer to a different subspace, echo or decoupling may be involved, or the weak-coupling Markovian bookkeeping model is not appropriate.

  1. Spatial scaling. If F(Δx)=Λ∣Δx∣2F(\Delta x)=\Lambda|\Delta x|^2, how does the decoherence time change when the separation is tripled?
Solution

The e−1e^{-1} time is

tdec(Δx)=1Λ∣Δx∣2.t_{\mathrm{dec}}(\Delta x) = \frac{1}{\Lambda|\Delta x|^2}.

For 3Δx3\Delta x,

tdec(3Δx)=19Λ∣Δx∣2=19tdec(Δx).t_{\mathrm{dec}}(3\Delta x) = \frac{1}{9\Lambda|\Delta x|^2} = \frac{1}{9} t_{\mathrm{dec}}(\Delta x).

Tripling the separation reduces the decoherence time by a factor of nine in the quadratic regime.

  1. Gaussian versus exponential envelope. A Ramsey envelope is W(T)=e−(T/TG)2W(T)=e^{-(T/T_G)^2}. At what time does it reach e−1e^{-1}, and why should TGT_G not be casually called an exponential T2T_2?
Solution

The threshold equation is

e−(T/TG)2=e−1,e^{-(T/T_G)^2}=e^{-1},

so T=TGT=T_G. The threshold happens to equal TGT_G, but the envelope shape is Gaussian, not exponential. An exponential T2T_2 means e−T/T2e^{-T/T_2}. The two envelopes have different short-time behavior, different fitted parameters over finite windows, and usually different physical origins.

  1. Fragment rate estimate. Environmental fragments arrive at rate RR and each fragment has real overlap η\eta between two alternatives. Estimate the e−1e^{-1} decoherence time using the Poisson model.
Solution

The coherence factor is

D(t)=exp⁡[Rt(η−1)]=exp⁡[−R(1−η)t].D(t) = \exp[ Rt(\eta-1) ] = \exp[ -R(1-\eta)t ].

The e−1e^{-1} time is therefore

tdec=1R(1−η).t_{\mathrm{dec}} = \frac{1}{R(1-\eta)}.

If η\eta is close to one, many fragments are needed. If η\eta is close to zero, one fragment almost fully distinguishes the alternatives and the timescale is near R−1R^{-1}.