T1 Relaxation Time

Prepare a qubit in its excited state and walk away, and you will not find it there when you return. A real two-level system is never truly isolated: it is bolted to a vast, restless environment — the electromagnetic vacuum, thermal photons, lattice phonons, tunnelling defects in amorphous oxides, stray quasiparticles — and any of those channels can carry off the qubit's energy quantum ℏω01\hbar\omega_{01}. The longitudinal relaxation time T1T_1 is the single number that says how fast: the timescale on which the excited-state population drains away toward thermal equilibrium. The name is inherited from nuclear magnetic resonance, where T1T_1 relaxes the magnetization along the static field — the "longitudinal" direction — and the same symbol has followed the physics from spins in a magnet to superconducting circuits on a chip. What makes T1T_1 more than a decay constant is what sets it. Relaxation is a population-transfer process, and Fermi's golden rule ties its rate to one specific thing: how much environmental noise exists at exactly the qubit frequency. 1/T11/T_1 is proportional to the noise spectral density S(ω01)S(\omega_{01}), so lengthening T1T_1 is the concrete task of quieting the bath at one frequency — by filtering, by detuning the qubit from a lossy resonator (the Purcell effect), or by removing the dielectric defects that host the noise. T1T_1 is also prior to the dephasing time T2T_2: because losing energy also erases phase, relaxation puts a hard ceiling T2≤2T1T_2 \le 2T_1 on any coherence. It is the clock every operation races against — the number of good gates scales as T1/tgateT_1/t_{\mathrm{gate}} — which is why long T1T_1 sits at the heart of the DiVincenzo criteria and headlines the datasheet of every qubit ever built.

Key takeaways
  • T1T_1 is the timescale on which a qubit irreversibly loses its energy to the environment, decaying from ∣1⟩|1\rangle toward thermal equilibrium.
  • The excited-state population follows a single exponential, P1(t)=P1(∞)+[P1(0)−P1(∞)] e−t/T1P_1(t)=P_1(\infty)+[P_1(0)-P_1(\infty)]\,e^{-t/T_1}, relaxing toward a thermal floor P1(∞)P_1(\infty), not always to zero.
  • The rate is fixed by Fermi's golden rule, 1/T1∝S(ω01)1/T_1 \propto S(\omega_{01}) — the noise spectral density at the qubit frequency, the one knob T1T_1 engineering turns.
  • Energy loss also destroys phase, so T1T_1 caps coherence through 1/T2=1/(2T1)+1/Tϕ1/T_2 = 1/(2T_1)+1/T_\phi, giving the hard bound T2≤2T1T_2 \le 2T_1.
  • The hardware-independent face of T1T_1 is the amplitude-damping channel, a single Kraus parameter γ=1−e−t/T1\gamma = 1-e^{-t/T_1}.
  • In today's superconducting qubits T1T_1 runs from tens of microseconds to a millisecond, limited by dielectric two-level-system defects, Purcell decay, and quasiparticles.
  • Key idea: T1T_1 is the inverse of the energy-relaxation rate set by the environmental noise spectral density at the qubit frequency, 1/T1=Γ1∝S(ω01)1/T_1 = \Gamma_1 \propto S(\omega_{01}) — it measures how fast a qubit forgets which state it is in, while T2T_2 measures how fast it forgets the phase between states.

Where this lives today

T1T_1 is the number every superconducting-qubit group lives and dies by, and its history is a story of chasing down loss channels one at a time. The first Cooper-pair-box qubits relaxed in nanoseconds; the transmon and 3D cavities pushed T1T_1 into the tens of microseconds by the mid-2010s; the current frontier reaches hundreds of microseconds to a millisecond. The gains came from identifying the dominant bath: amorphous oxide at interfaces hosts a dense forest of two-level-system (TLS) defects whose noise sits right at ω01\omega_{01}, so replacing niobium with tantalum and cleaning the substrate — Place and coworkers at Princeton, 2021 — roughly tripled transmon T1T_1 to the 0.3–0.5 ms range, while low-frequency fluxonium qubits have reached T1T_1 beyond a millisecond. Purcell filters are now standard, decoupling fast readout from qubit decay exactly as Worked Example 2 demands. And T1T_1 is not even constant in time: it fluctuates as TLS defects drift through resonance, and — startlingly — ionizing radiation matters. Cosmic rays and ambient radioactivity dump bursts of quasiparticles that cause correlated, chip-wide dips in T1T_1 (Vepsäläinen and coworkers, 2020; McEwen and coworkers at Google, 2022), a real obstacle for error correction that has sent groups to underground labs and shielded cryostats. The abstract 1/T1∝S(ω01)1/T_1\propto S(\omega_{01}) has become, in the lab, a materials-science program.

The physical picture

Picture the qubit as a tiny pendulum tied to a vast bath of springs — every mode of the environment it can talk to. If any bath spring resonates near the pendulum's own frequency, energy drains into that mode and does not come back; if the bath is silent there, the pendulum keeps swinging. That is the whole content of 1/T1∝S(ω01)1/T_1 \propto S(\omega_{01}): the rate is how "crowded" the environment is with modes at the qubit frequency. Relaxation is one-directional bookkeeping on populations — probability moves from ∣1⟩|1\rangle to ∣0⟩|0\rangle by emission, and at finite temperature back by absorption — until the two rates reach detailed balance and leave a thermal mixture, not a pure ground state. The first figure pins the one relation the whole article turns on.

Figure
Energy relaxation of a qubit: the excited state loses one quantum ħω₀₁ to a continuum of bath modes at a rate fixed by the noise spectral density there, 1/T₁ = Γ₁ ∝ S(ω₀₁). For a 5 GHz transmon at 20 mK the bath is nearly empty, n_th ≈ 6×10⁻⁶, so emission dominates and the qubit decays toward the ground state
A qubit dumps its energy into the bath at a rate set by the noise there: 1/T₁ ∝ S(ω₀₁) |1⟩ |0⟩ excited ground emission Γ↓ absorption Γ↑ quantum ħω₀₁ into the bath ω S(ω) ω₀₁ S(ω₀₁) 1/T₁ = Γ₁ ∝ S(ω₀₁)

The second thing to hold onto is that T1T_1 is not just a rate but a ceiling. A relaxation event does not only move population — it also collapses any phase the qubit was carrying, because a decayed qubit and a surviving one are distinguishable to the environment. So relaxation alone decoheres a superposition at half the population rate, and no engineering can keep a coherent state alive longer than twice the time the energy takes to leak away: T2≤2T1T_2 \le 2T_1. Pure dephasing TϕT_\phi, which scrambles phase without moving energy, only pulls T2T_2 down further. The second figure is what an experimentalist actually measures — the excited-state decay curve — with the numbers of a good modern transmon and the amplitude-damping and T2T_2-ceiling readings marked.

Figure
Measured excited-state population P₁(t) = e^(−t/T₁) for a superconducting transmon with T₁ = 100 μs and ω₀₁/2π = 5 GHz. After one T₁ the survival is 1/e = 0.368 while the amplitude-damping probability rises to γ = 1 − 1/e = 0.632; the two curves cross at P₁ = γ = 0.5 near t = 69 μs, and the relaxation ceiling on coherence is T₂ ≤ 2T₁ = 200 μs
The punchline: single-exponential decay, T₁ = 100 μs — 1/e survival at t = T₁ 1.0 0.5 0.368 0 0 100 200 300 delay t (μs) population P₁(t) = e^(−t/T₁) γ(t) = 1 − e^(−t/T₁) t = T₁, P₁ = 0.368 γ = 0.632 2T₁ = 200 μs ceiling on T₂

The core equations

The starting point is kinematic: the excited-state population relaxes on the single timescale T1T_1, and because populations obey a linear rate equation the decay is a pure exponential toward a thermal floor. Everything downstream is a computation of, or a consequence of, this one law.

P1(t)=P1(∞)+[P1(0)−P1(∞)] e−t/T1(1)P_1(t) = P_1(\infty) + \big[P_1(0) - P_1(\infty)\big]\,e^{-t/T_1} \qquad (1)

where: P1(t)P_1(t) — excited-state population at time tt; P1(0)P_1(0) — initial population (typically 11 after a π\pi pulse); P1(∞)P_1(\infty) — thermal-equilibrium excited population; T1T_1 — longitudinal (energy) relaxation time. Reading it: the population decays not to zero but to P1(∞)P_1(\infty), on the one timescale T1T_1; at T=0T=0 the floor vanishes and the curve is the familiar e−t/T1e^{-t/T_1}. The rate itself is not free — it is fixed by the environment. Fermi's golden rule connects 1/T11/T_1 to a bath operator A^\hat A that the qubit couples to and to the noise spectral density of that bath evaluated at the qubit frequency, summing the emission (+ω01+\omega_{01}) and absorption (−ω01-\omega_{01}) contributions.

1T1=Γ1=∣⟨0∣A^∣1⟩∣2ℏ2[S(ω01)+S(−ω01)]  ∝  S(ω01)(2)\frac{1}{T_1} = \Gamma_1 = \frac{|\langle 0|\hat A|1\rangle|^2}{\hbar^2}\big[S(\omega_{01}) + S(-\omega_{01})\big] \;\propto\; S(\omega_{01}) \qquad (2)

where: Γ1\Gamma_1 — energy-relaxation rate; A^\hat A — qubit operator coupling to the bath; ∣⟨0∣A^∣1⟩∣|\langle 0|\hat A|1\rangle| — transition matrix element; S(ω)S(\omega) — noise spectral density, the Fourier transform of the bath correlation function. Reading it: only noise at the qubit frequency relaxes the qubit, which is why T1T_1 engineering is spectral — you quiet one frequency, not the whole bath. Because energy loss also erases phase, T1T_1 enters the coherence budget directly. The transverse (total dephasing) rate is the sum of the relaxation contribution — exactly half the population rate — and the pure-dephasing rate.

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

where: T2T_2 — total dephasing time; TϕT_\phi — pure-dephasing time (frequency noise, no energy exchange); T1T_1 — relaxation time. Reading it: with no pure dephasing (Tϕ→∞T_\phi\to\infty) one hits the relaxation-limited maximum T2=2T1T_2 = 2T_1, so a measured T2>2T1T_2 > 2T_1 is a red flag for a calibration error. The thermal floor is set by how many quanta the bath already holds at the qubit frequency — the Bose occupation — which fixes both the steady-state population and the balance between absorption and emission.

nth=1eℏω01/kBT−1,P1(∞)=nth2nth+1(4)n_{\mathrm{th}} = \frac{1}{e^{\hbar\omega_{01}/k_B T} - 1},\qquad P_1(\infty) = \frac{n_{\mathrm{th}}}{2n_{\mathrm{th}}+1} \qquad (4)

where: nthn_{\mathrm{th}} — mean bath occupation at ω01\omega_{01}; TT — effective bath temperature; kBk_B — Boltzmann constant. Reading it: for a 5 GHz qubit in a 20 mK fridge ℏω01/kBT≈12\hbar\omega_{01}/k_BT\approx 12 and nth≈6×10−6n_{\mathrm{th}}\approx 6\times10^{-6}, so a measured P1(∞)P_1(\infty) near a percent signals a hotter effective environment than the fridge. The hardware-independent face of all this is the amplitude-damping channel: the entire T1T_1 decay is packaged into two Kraus operators and a single damping parameter γ\gamma.

E(ρ)=E0ρE0†+E1ρE1†,E0=(1001−γ), E1=(0γ00),γ=1−e−t/T1(5)\mathcal{E}(\rho) = E_0\rho E_0^\dagger + E_1\rho E_1^\dagger,\quad E_0=\begin{pmatrix}1&0\\0&\sqrt{1-\gamma}\end{pmatrix},\ E_1=\begin{pmatrix}0&\sqrt{\gamma}\\0&0\end{pmatrix},\quad \gamma = 1-e^{-t/T_1} \qquad (5)

where: γ\gamma — damping (decay) probability; E0E_0 — no-jump Kraus operator; E1E_1 — quantum-jump (decay) operator; ρ\rho — qubit density matrix. Reading it: E1E_1 is the emission event ∣1⟩→∣0⟩|1\rangle\to|0\rangle, while E0E_0 is the "no-click" evolution that shrinks the surviving excited amplitude — together they map laboratory relaxation onto the theory of quantum channels. Finally, one dominant and controllable relaxation channel in circuit QED has a closed form: a qubit dispersively coupled to a lossy readout resonator inherits the resonator's loss, scaled by the coupling-to-detuning ratio.

Γ1Purcell=(gΔ)2κ(6)\Gamma_1^{\mathrm{Purcell}} = \left(\frac{g}{\Delta}\right)^2 \kappa \qquad (6)

where: gg — qubit–resonator coupling; Δ=ωq−ωr\Delta=\omega_q-\omega_r — qubit–resonator detuning; κ\kappa — resonator energy decay rate. Reading it: the detuned qubit hybridizes with the lossy mode by the small amount (g/Δ)2(g/\Delta)^2 and leaks at that fraction of κ\kappa — large detuning, or a filter at ω01\omega_{01}, protects T1T_1.

Derivation AFrom the population master equation to the exponential decay law

The goal is to show that the single-exponential relaxation law of equation (1) is not an assumption but the forced solution of a linear rate equation for the excited-state population, and that the same equation delivers the thermal floor P1(∞)P_1(\infty) and identifies 1/T11/T_1 as the sum of the emission and absorption rates. We balance the upward and downward transitions, read off the steady state, and solve the resulting first-order ODE.

Step 1. Balance the population flow. In time dt\mathrm{d}t the excited state loses probability at the emission rate Γ↓\Gamma_\downarrow and gains it from the ground state at the absorption rate Γ↑\Gamma_\uparrow, giving P˙1=−Γ↓P1+Γ↑(1−P1)\dot P_1 = -\Gamma_\downarrow P_1 + \Gamma_\uparrow(1-P_1). Why: relaxation is bookkeeping on diagonal density-matrix elements — nothing but two competing one-way transition rates, with the ground-state population written as 1−P11-P_1.

Step 2. Read the total rate off the ODE. Rewriting, P˙1=−(Γ↓+Γ↑)P1+Γ↑\dot P_1 = -(\Gamma_\downarrow+\Gamma_\uparrow)P_1 + \Gamma_\uparrow, so the coefficient of P1P_1 is the full relaxation rate: 1/T1=Γ↓+Γ↑1/T_1 = \Gamma_\downarrow+\Gamma_\uparrow. Why: it is the sum, not the difference — both emission and absorption push the population toward its steady value, so both add to 1/T11/T_1.

Step 3. Set P˙1=0\dot P_1=0 for the steady state, P1(∞)=Γ↑/(Γ↓+Γ↑)P_1(\infty)=\Gamma_\uparrow/(\Gamma_\downarrow+\Gamma_\uparrow). Detailed balance fixes the ratio through the Bose factor, Γ↑/Γ↓=nth/(nth+1)\Gamma_\uparrow/\Gamma_\downarrow = n_{\mathrm{th}}/(n_{\mathrm{th}}+1), giving P1(∞)=nth/(2nth+1)P_1(\infty)=n_{\mathrm{th}}/(2n_{\mathrm{th}}+1). Why: the ratio of up- to down-rates is pinned by temperature, which is what makes T1T_1 measurements double as thermometry.

Step 4. Solve the linear ODE. With constant coefficients, P1(t)=P1(∞)+[P1(0)−P1(∞)]e−t/T1P_1(t)=P_1(\infty)+[P_1(0)-P_1(\infty)]e^{-t/T_1}; at T=0T=0, P1(∞)=0P_1(\infty)=0, recovering pure spontaneous decay e−t/T1e^{-t/T_1}. Why: a first-order linear equation can only give exponential approach to its fixed point, so the single-exponential shape is a theorem — a non-exponential decay is evidence of physics beyond this simple bath. *The solution recovers P1(∞)=nth/(2nth+1)P_1(\infty)=n_{\mathrm{th}}/(2n_{\mathrm{th}}+1) and, at T=0T=0, the pure amplitude-damping law e−t/T1e^{-t/T_1}, consistent with the Lindblad master equation for a two-level system coupled to a thermal bath.*

Go deeper · Derivation B · 1 minThe relaxation rate from Fermi's golden rule
Derivation BThe relaxation rate from Fermi's golden rule

The goal is to derive 1/T1∝S(ω01)1/T_1 \propto S(\omega_{01}) from first principles — to show that the qubit's decay rate is literally the environmental noise spectral density at the transition frequency, times a matrix element — by treating the environment as a continuum of weakly coupled modes. We apply the golden rule, package the bath's mode sum as a spectral density, and use detailed balance to fix emission against absorption.

Step 1. Write the interaction as H^int=A^⊗B^\hat H_{\mathrm{int}}=\hat A\otimes\hat B, with A^\hat A a qubit operator and B^\hat B a bath operator. Fermi's golden rule gives the decay rate into all bath final states, Γ1→0=2πℏ2∣⟨0∣A^∣1⟩∣2∑f∣⟨f∣B^∣i⟩∣2 δ(ω01−ωfi)\Gamma_{1\to0}=\frac{2\pi}{\hbar^2}|\langle 0|\hat A|1\rangle|^2\sum_f|\langle f|\hat B|i\rangle|^2\,\delta(\omega_{01}-\omega_{fi}). Why: weak coupling to a dense continuum is the golden rule's home turf, and the delta function enforces energy conservation — the qubit decays only by handing its quantum to a bath mode of matching frequency.

Step 2. Recognize the bath sum as a spectral density. The weighted sum over final states is the Fourier transform of the bath correlation function, S(ω)=∫−∞∞⟨B^(t)B^(0)⟩ eiωt dtS(\omega)=\int_{-\infty}^{\infty}\langle\hat B(t)\hat B(0)\rangle\,e^{i\omega t}\,\mathrm{d}t, evaluated at ω=ω01\omega=\omega_{01}. Why: the microscopic clutter of bath modes collapses into one measurable function of frequency, and the qubit only ever "sees" its value at ω01\omega_{01}.

Step 3. Separate emission and absorption. Positive frequency releases energy, Γ↓∝S(+ω01)\Gamma_\downarrow\propto S(+\omega_{01}); negative absorbs it, Γ↑∝S(−ω01)\Gamma_\uparrow\propto S(-\omega_{01}). Summing, 1/T1=∣⟨0∣A^∣1⟩∣2ℏ2[S(ω01)+S(−ω01)]1/T_1=\frac{|\langle 0|\hat A|1\rangle|^2}{\hbar^2}[S(\omega_{01})+S(-\omega_{01})]. Why: a quantum bath has an asymmetric spectrum — it gives and takes energy at different rates — and the two signed frequencies bookkeep the directions cleanly.

Step 4. Impose detailed balance. The quantum fluctuation–dissipation theorem fixes S(−ω01)/S(+ω01)=e−ℏω01/kBTS(-\omega_{01})/S(+\omega_{01})=e^{-\hbar\omega_{01}/k_BT}, reproducing the thermal steady state and, at T=0T=0, killing absorption so 1/T1=∣⟨0∣A^∣1⟩∣2S(ω01)/ℏ21/T_1=|\langle 0|\hat A|1\rangle|^2 S(\omega_{01})/\hbar^2. Why: dissipation and fluctuation are two faces of one spectral density — the same asymmetry that lets the bath absorb is what forces a thermal, not pure, ground state. *The derivation reproduces 1/T1∝S(ω01)1/T_1\propto S(\omega_{01}) and the detailed-balance ratio, in agreement with Bloch–Redfield theory and the standard circuit-QED noise-spectroscopy relation Γ1=∣⟨0∣A^∣1⟩∣2S(ω01)/ℏ2\Gamma_1=|\langle 0|\hat A|1\rangle|^2 S(\omega_{01})/\hbar^2.*

Go deeper · Derivation C · 1 minAmplitude damping and the T₂ ≤ 2T₁ bound
Derivation CAmplitude damping and the T₂ ≤ 2T₁ bound

The goal is to prove that energy relaxation alone decoheres a superposition at exactly half the population-decay rate, which fixes the hard ceiling T2≤2T1T_2\le 2T_1, and to show where the pure-dephasing term of equation (3) adds on top. We push a coherence through the amplitude-damping channel, see how the off-diagonal element shrinks, translate that to a rate, and append independent phase noise.

Step 1. Track the coherence under damping. Applying the channel of equation (5) to a general ρ\rho, the off-diagonal element transforms as ρ01→1−γ ρ01\rho_{01}\to\sqrt{1-\gamma}\,\rho_{01}, while the excited population transforms as ρ11→(1−γ)ρ11\rho_{11}\to(1-\gamma)\rho_{11}. Why: E0E_0 carries 1−γ\sqrt{1-\gamma} on the excited amplitude, so a coherence (linear in that amplitude) decays with the square root, whereas a population (quadratic) decays with the full factor.

Step 2. Convert to rates. Since 1−γ=e−t/T11-\gamma=e^{-t/T_1}, the coherence factor is 1−γ=e−t/2T1\sqrt{1-\gamma}=e^{-t/2T_1}, so relaxation decays ∣ρ01∣|\rho_{01}| at rate 1/(2T1)1/(2T_1) — exactly half the population rate. Why: this factor of two is the whole origin of the 2T12T_1 bound — a kinematic consequence of amplitude being the square root of probability, not an extra assumption.

Step 3. Add pure dephasing. Independent, energy-conserving frequency noise multiplies the coherence by e−t/Tϕe^{-t/T_\phi}, so ∣ρ01(t)∣=∣ρ01(0)∣ e−t/2T1 e−t/Tϕ|\rho_{01}(t)|=|\rho_{01}(0)|\,e^{-t/2T_1}\,e^{-t/T_\phi}. Why: the two mechanisms are statistically independent, so their envelopes multiply and their rates add.

Step 4. Read off T2T_2. The combined exponent defines 1/T2=1/(2T1)+1/Tϕ1/T_2=1/(2T_1)+1/T_\phi; with Tϕ→∞T_\phi\to\infty this saturates at T2=2T1T_2=2T_1, and since 1/Tϕ≥01/T_\phi\ge 0 one always has T2≤2T1T_2\le 2T_1. Why: the pure-dephasing rate is non-negative, so it can only shorten T2T_2 — there is no beating the relaxation ceiling, which is why a reported T2>2T1T_2>2T_1 must be an artifact. *The result gives T2=2T1T_2=2T_1 at Tϕ→∞T_\phi\to\infty and T2≤2T1T_2\le 2T_1 in general; the factor-of-two coherence-versus-population relation matches the standard Bloch–Redfield and NMR definitions of T1T_1 and T2T_2.*

Limiting cases

How T1T_1 and its consequences behave in the regimes that matter for real devices.

CaseResultMeaning
Zero temperature (nth→0n_{\mathrm{th}}\to 0)P1(t)=P1(0) e−t/T1P_1(t)=P_1(0)\,e^{-t/T_1}Pure spontaneous emission; the steady state is the ground state ∣0⟩\vert 0\rangle and the channel is pure amplitude damping with rate Γ1=1/T1\Gamma_1=1/T_1.
Finite temperatureP1(∞)=nth2nth+1P_1(\infty)=\dfrac{n_{\mathrm{th}}}{2n_{\mathrm{th}}+1}, 1T1∝(2nth+1)\dfrac{1}{T_1}\propto(2n_{\mathrm{th}}+1)Decay relaxes toward a nonzero thermal floor and speeds up by (2nth+1)(2n_{\mathrm{th}}+1); a resolved floor measures the effective bath temperature.
No pure dephasing (Tϕ→∞T_\phi\to\infty)T2=2T1T_2 = 2T_1The relaxation-limited maximum coherence time; any pure dephasing only lowers T2T_2, so T2≤2T1T_2\le 2T_1 always.
Purcell-dominatedΓ1=(gΔ)2κ\Gamma_1=\left(\dfrac{g}{\Delta}\right)^2\kappaA qubit dispersively coupled to a resonator of linewidth κ\kappa inherits decay set by coupling and detuning; large Δ\Delta protects T1T_1.
Short time (t≪T1t\ll T_1)P1(t)≈P1(0) (1−t/T1)P_1(t)\approx P_1(0)\,(1-t/T_1), γ≈t/T1\gamma\approx t/T_1Linear decay; the per-shot error of an idle qubit grows as t/T1t/T_1, the form used in error budgets.

Worked examples

1 · A 5 GHz transmon with T₁ = 100 μs (fully worked)

A superconducting transmon has transition frequency ω01/2π=5\omega_{01}/2\pi = 5 GHz and a measured T1=100 μT_1 = 100\ \mus. The dilution fridge sits at 20 mK. Find the equilibrium excited population, the amplitude-damping probability after one T1T_1, the ceiling on T2T_2, and roughly how many 20 ns gates fit inside the coherence window. Then explain why the measured residual population is nearer 1%.

Step 1. Thermal occupation at the fridge temperature. Form ℏω01/kBT=hf/kBT=(6.626×10−34)(5×109)/[(1.381×10−23)(0.020)]=12.0\hbar\omega_{01}/k_BT = hf/k_BT = (6.626\times10^{-34})(5\times10^{9})/[(1.381\times10^{-23})(0.020)] = \mathbf{12.0}. Then nth=1/(e12−1)=6.2×10−6n_{\mathrm{th}} = 1/(e^{12}-1) = \mathbf{6.2\times10^{-6}}, so the ideal equilibrium population is P1(∞)≈nth=6×10−6P_1(\infty)\approx n_{\mathrm{th}} = \mathbf{6\times10^{-6}} — utterly negligible.

Step 2. Amplitude damping after one T1T_1. At t=T1t=T_1, γ=1−e−1=0.632\gamma = 1-e^{-1} = \mathbf{0.632}, and the survival is P1(T1)=e−1=0.368P_1(T_1)=e^{-1}=\mathbf{0.368}. So after 100 μ100\ \mus an initially excited qubit has a 63% chance of having emitted its quantum — this is the 1/e1/e point marked in the second figure.

Step 3. Coherence ceiling and gate budget. The relaxation-limited maximum is T2≤2T1=200 μsT_2 \le 2T_1 = \mathbf{200\ \mu s}; reaching it requires pure dephasing Tϕ≫T1T_\phi \gg T_1. With a gate time tgate=20t_{\mathrm{gate}}=20 ns, the number of gates before relaxation matters is T1/tgate=100 μs/20 ns=5000T_1/t_{\mathrm{gate}} = 100\ \mu\mathrm{s}/20\ \mathrm{ns} = \mathbf{5000} — the rough depth ceiling this T1T_1 imposes on a circuit.

Step 4. Why the measured floor is nearer a percent. The 6×10⁻⁶ prediction assumes the qubit sees 20 mK, but non-equilibrium quasiparticles and stray thermal photons give an effective temperature closer to 50 mK. Repeating Step 1 at 50 mK: ℏω01/kBT=4.80\hbar\omega_{01}/k_BT = 4.80, nth=1/(e4.80−1)=8.3×10−3n_{\mathrm{th}} = 1/(e^{4.80}-1) = 8.3\times10^{-3}, so P1(∞)=nth/(2nth+1)=0.82%P_1(\infty) = n_{\mathrm{th}}/(2n_{\mathrm{th}}+1) = \mathbf{0.82\%}. The floor is not the fridge temperature — it is a thermometer for the qubit's true environment, here about 50 mK.

2 · Purcell-limited T₁ and how far to detune

A transmon couples to its readout resonator with g/2π=100g/2\pi = 100 MHz at a detuning Δ/2π=1\Delta/2\pi = 1 GHz; the resonator loses photons at κ/2π=1\kappa/2\pi = 1 MHz. Find the Purcell-limited T1T_1. Then find the detuning needed to push the Purcell ceiling above 1 ms.

Step 1. Purcell rate at the given detuning. The hybridization fraction is (g/Δ)2=(100/1000)2=0.01(g/\Delta)^2 = (100/1000)^2 = 0.01, so Γ1Purcell=0.01×(2π×106 s−1)=6.3×104 s−1\Gamma_1^{\mathrm{Purcell}} = 0.01\times(2\pi\times10^6\ \mathrm{s}^{-1}) = 6.3\times10^{4}\ \mathrm{s}^{-1}.

Step 2. Convert to a time. T1Purcell=1/Γ1Purcell=16 μsT_1^{\mathrm{Purcell}} = 1/\Gamma_1^{\mathrm{Purcell}} = \mathbf{16\ \mu s}. If nothing else limited the qubit, this single channel would cap T1T_1 near 16 μs — modest, and a real problem for readout design.

Step 3. Detune to reach 1 ms. Writing T1Purcell=(Δ/g)2/κT_1^{\mathrm{Purcell}} = (\Delta/g)^2/\kappa, solve at T1=1T_1=1 ms: (Δ/g)2=κT1=(2π×106)(10−3)=6.3×103(\Delta/g)^2 = \kappa T_1 = (2\pi\times10^6)(10^{-3}) = 6.3\times10^{3}, so Δ/g=79\Delta/g = 79 and Δ/2π≈7.9 GHz\Delta/2\pi \approx \mathbf{7.9\ GHz}. A nearly eightfold-larger detuning buys a Purcell ceiling above 1 ms — but it weakens dispersive readout, which is why real devices add a Purcell filter to block κ\kappa at ω01\omega_{01} instead of detuning alone.

Physics of the portrait

The symbols T1T_1 and T2T_2 did not begin in a quantum computer. They were written down in 1946 to describe nuclear spins relaxing in a magnet, by the physicist who turned magnetic resonance into a quantitative science — and have survived intact all the way to the superconducting chip.

Origins · 1 minFelix Bloch and how this was discovered
FB
Felix Bloch
1905–1983 · nuclear magnetic resonance and the theory of relaxation

Felix Bloch was a Swiss-born physicist who, after foundational work on electrons in crystals — Bloch waves and Bloch's theorem are his — turned to the magnetism of atomic nuclei. In 1946 at Stanford he introduced the phenomenological equations that still bear his name, describing how a nuclear magnetization, tipped away from equilibrium, returns to it. He named two distinct timescales. The longitudinal relaxation time T1T_1 ("spin-lattice") governs recovery of the magnetization along the static field — energy exchange with the lattice — while the transverse time T2T_2 ("spin-spin") governs decay of the magnetization across it — loss of phase coherence among precessing spins. That very distinction, longitudinal energy relaxation versus transverse dephasing, is the one this article rests on, transplanted from spins to qubits with the physics unchanged. Bloch shared the 1952 Nobel Prize with Edward Purcell — whose group detected nuclear resonance in bulk matter the same year, and whose name attaches to the cavity-enhanced emission that limits T1T_1 in circuit QED. The microscopic theory of what sets these times followed from Bloembergen, Purcell and Pound and then Wangsness, Bloch and Redfield, whose master equation is the direct ancestor of the golden-rule rate 1/T1∝S(ω01)1/T_1\propto S(\omega_{01}).

History
1946
— Bloch introduces the phenomenological Bloch equations, defining the longitudinal T₁ and transverse T₂ relaxation times
1946
— Purcell, Torrey and Pound detect nuclear magnetic resonance in bulk matter, and Purcell predicts cavity-enhanced spontaneous emission
1948
— Bloembergen, Purcell and Pound relate nuclear relaxation rates to the spectral density of fluctuating fields
1957
— Wangsness, Bloch and Redfield cast relaxation as a microscopic master equation, the Bloch–Redfield theory
2005
— Ithier and coworkers apply the noise-spectral-density picture to decoherence in a superconducting qubit

Going further

Texts. Nielsen and Chuang, Quantum Computation and Quantum Information, §8.3 (quantum operations and the amplitude-damping channel) · Preskill, Lecture Notes on Quantum Computation, Ch. 3 (decoherence, master equations, and golden-rule relaxation) · Wilde, Quantum Information Theory, Ch. on quantum channels (amplitude damping) · Watrous, The Theory of Quantum Information (channels and CPTP maps) · Kitaev, Shen and Vyalyi, Classical and Quantum Computation (noise models) · Krantz et al., A Quantum Engineer's Guide to Superconducting Qubits, Appl. Phys. Rev. 6, 021318 (2019) (T₁, T₂, and noise spectroscopy in practice). Key results to keep in mind: the thermal steady state of the population rate equation is P1(∞)=nth/(2nth+1)P_1(\infty)=n_{\mathrm{th}}/(2n_{\mathrm{th}}+1); the detailed-balance ratio is S(−ω)/S(+ω)=e−ℏω/kBTS(-\omega)/S(+\omega)=e^{-\hbar\omega/k_BT}; the coherence relation 1/T2=1/(2T1)+1/Tϕ1/T_2=1/(2T_1)+1/T_\phi follows from 1−γ=e−t/2T1\sqrt{1-\gamma}=e^{-t/2T_1}; the Worked Example 1 numbers are nth≈6×10−6n_{\mathrm{th}}\approx 6\times10^{-6} at 20 mK, P1(∞)≈0.82%P_1(\infty)\approx 0.82\% at 50 mK, and γ(T1)=0.632\gamma(T_1)=0.632; the Worked Example 2 Purcell values are T1Purcell≈16 μT_1^{\mathrm{Purcell}}\approx 16\ \mus and Δ/2π≈7.9\Delta/2\pi\approx 7.9 GHz for a 1 ms ceiling. Device T1T_1 values are representative modern figures, not universal constants, since they are strongly material- and design-dependent.

Prerequisites. DiVincenzo Criteria · Density Matrix and Open Systems · Fermi's Golden Rule. Leads to. The Transmon Qubit · Circuit Quantum Electrodynamics · Dispersive Readout · Sideband Cooling.

Ask the tutor: Ask me to walk the golden-rule derivation of 1/T1∝S(ω01)1/T_1\propto S(\omega_{01}) line by line, or to work out the effective temperature behind a measured 1% excited-state floor, and I will take it step by step.

What is the difference between T₁ and T₂, in one breath?

T1T_1 is about energy; T2T_2 is about phase. T1T_1 is how long an excited qubit takes to lose its energy quantum and fall back toward the ground state — it acts on populations, the diagonal of the density matrix. T2T_2 is how long a superposition keeps a well-defined phase between ∣0⟩|0\rangle and ∣1⟩|1\rangle — it acts on the off-diagonal coherences. They are linked but not identical: because any energy-loss event also destroys phase, relaxation contributes to T2T_2 through 1/T2=1/(2T1)+1/Tϕ1/T_2=1/(2T_1)+1/T_\phi, giving the ceiling T2≤2T1T_2\le 2T_1. Pure dephasing TϕT_\phi — random modulation of the qubit frequency with no energy exchange — is the extra ingredient that pulls T2T_2 below that ceiling.

Why does only noise at the qubit frequency matter for T₁?

Because relaxation is a resonant energy exchange. To decay, the qubit must hand its quantum ℏω01\hbar\omega_{01} to some environmental mode, and energy conservation — the delta function in Fermi's golden rule — means the receiving mode must sit at the qubit frequency. Noise at other frequencies cannot accept the quantum, so it does not contribute to T1T_1. That is why the rate is 1/T1∝S(ω01)1/T_1\propto S(\omega_{01}) evaluated at that one frequency, and why T1T_1 engineering is spectral surgery: you filter or detune to open a quiet spot in the bath exactly at ω01\omega_{01}, leaving every other frequency untouched. Low-frequency noise, by contrast, mostly causes pure dephasing.

Does the excited state really not decay all the way to zero?

Only at absolute zero does it. At finite temperature the environment holds thermal quanta at the qubit frequency, so absorption competes with emission and the population settles at a nonzero floor P1(∞)=nth/(2nth+1)P_1(\infty)=n_{\mathrm{th}}/(2n_{\mathrm{th}}+1) set by detailed balance. For a 5 GHz qubit in a 20 mK fridge this floor is theoretically about 6×10−66\times10^{-6} — essentially zero. But measured floors in real qubits are often near a percent, because stray infrared photons and non-equilibrium quasiparticles push the effective temperature closer to 50 mK. Resolving that floor is a genuine thermometer for the qubit's environment, which is why careful experiments fit the offset rather than forcing it to zero.

What is the amplitude-damping channel, and why do information theorists like it?

It is the hardware-independent description of T1T_1 decay as a quantum channel. Rather than tracking a specific bath, it summarizes the entire effect of relaxation with two Kraus operators and a single number γ=1−e−t/T1\gamma=1-e^{-t/T_1}, the probability the qubit has emitted its quantum. One operator, E1E_1, is the actual quantum jump ∣1⟩→∣0⟩|1\rangle\to|0\rangle; the other, E0E_0, is the "no-jump" evolution that quietly shrinks the surviving excited amplitude by 1−γ\sqrt{1-\gamma}. Information theorists like it because it is exactly solvable, it is the canonical lossy channel, and its parameters feed directly into channel capacities and error-correction thresholds — you reason about relaxation abstractly, then plug in γ(t)\gamma(t) from a measured T1T_1.

What is the Purcell effect and why does it hurt T₁?

The Purcell effect is the modification of a system's spontaneous-emission rate by its electromagnetic environment. In circuit QED a qubit is deliberately coupled to a readout resonator so you can measure it, but that resonator is lossy — it leaks photons at rate κ\kappa — and so opens a decay path. Even detuned by Δ\Delta, the qubit hybridizes with the resonator by the small fraction (g/Δ)2(g/\Delta)^2 and inherits that fraction of the loss, giving a rate (g/Δ)2κ(g/\Delta)^2\kappa. Large detuning suppresses it but weakens readout, so the standard fix is a Purcell filter: a frequency-selective element that blocks the resonator's loss at the qubit frequency while passing it at the readout frequency.

How is T₁ actually measured in the lab?

By inversion recovery, the direct implementation of equation (1). You prepare ∣1⟩|1\rangle with a calibrated π\pi pulse, wait a variable delay τ\tau, then read out the excited-state population, repeating over many shots at each delay. Sweeping τ\tau traces the decay curve, and fitting P1(τ)=A e−τ/T1+P1(∞)P_1(\tau)=A\,e^{-\tau/T_1}+P_1(\infty) returns T1T_1 as the time constant, with amplitude AA absorbing state-preparation error and offset P1(∞)P_1(\infty) capturing the thermal floor. Because individual defects drift, T1T_1 is not perfectly stable, so careful characterization repeats the measurement over hours to report a distribution rather than a single number.

Before this

Where this leads