Spectral gaps and relaxation times
Use the weighted reversible spectrum while separating slow positive modes from nearly periodic negative modes.
On this page
The question for today
What you’ll understand
- Distinguish ordinary and absolute spectral gaps for a reversible finite chain.
- Derive a total-variation upper bound from weighted L² contraction.
- Explain precisely how laziness changes negative and positive eigenvalues.
A reversible transition kernel is self-adjoint in a stationary-weighted inner product, so its nonconstant eigenfunctions describe independently decaying modes. The slowest decay in discrete time depends on eigenvalue magnitude, including eigenvalues near minus one. Looking only at the second largest eigenvalue can therefore predict rapid equilibration for a chain that almost alternates deterministically. This lesson fixes the gap convention, derives a usable distributional bound and shows why adding holding probability can remove an oscillatory obstruction even though holding often slows positive modes.
There are two gaps to keep separate
Assume a finite irreducible reversible chain on at least two states, with stationary law π. List its real eigenvalues in decreasing order, beginning with the simple eigenvalue one. The ordinary gap is one minus the second largest eigenvalue. The absolute gap instead subtracts the largest modulus among all nonconstant eigenvalues. The ordinary gap is positive under irreducibility; the absolute gap is positive exactly when the chain is also aperiodic. A near-minus-one eigenvalue barely affects the ordinary gap but can dominate the actual discrete-time decay.
Here relaxation time means the reciprocal absolute gap, following the convention used by the cited text; for lazy chains it also equals the reciprocal ordinary gap. Naming this convention matters because some literature uses ordinary-gap relaxation time without the absolute qualifier. Weighted variance contracts by the square of the largest nonconstant eigenvalue modulus at each transition. An eigenfunction attaining that modulus achieves equality, so this rate is not merely a convenient estimate: it detects the slowest reversible L² mode.
L² contraction becomes a probability bound
Reversibility makes the density of a t-step law relative to π evolve by the same self-adjoint operator. Subtract its constant component and apply the spectral contraction. The initial point-mass density at x has squared weighted norm one divided by π(x), minus one. Cauchy–Schwarz converts its weighted L² error into an L¹ error, introducing the total-variation factor one half. Taking the smallest stationary mass gives a bound uniform over starting states. This prefactor records the cost of concentrating initially on an improbable state.
The resulting inequality is a general upper bound and can exceed one at early times, in which case the trivial bound one is better. Replacing the geometric factor by an exponential gives a convenient logarithmic mixing estimate involving the inverse absolute gap and the smallest stationary mass. A gap estimate alone is not an exact mixing time, because different modes and prefactors contribute. The finite-state assumptions also avoid unbounded point-mass densities and continuous spectral complications that require separate arguments on countable or general spaces.
Laziness moves every eigenvalue
Replacing P by its average with the identity changes an eigenvalue λ to (1+λ)/2. All new eigenvalues lie between zero and one. A mode near minus one then moves near zero and loses its oscillatory obstruction, while a positive eigenvalue near one moves even closer to one and decays more slowly per transition. Thus the blanket claim that laziness always halves speed, or always accelerates mixing, is false. Its effect depends on which end of the original spectrum controls convergence.
For the two-state chain, the only nonconstant eigenvalue is 1−p−q, so the distinction is completely visible. A large sum of transition probabilities produces alternating convergence; a small sum produces slow monotone convergence. The same stationary law can coexist with very different decay rates if p and q are scaled together. This observation links spectral analysis back to actual transition design: preserving the target distribution is necessary for sampling, but the chosen dynamics determine how long one must wait to obtain an approximately stationary draw.
Work it through
An ordinary gap can hide almost-periodic behavior
Let p=q=19/20. Compute both gaps, the relaxation time and the exact mixing time at tolerance 1/10. Compare the lazy modification.
The stationary law is uniform and the nonconstant eigenvalue is −9/10. Thus the ordinary gap is 19/10, the absolute gap is 1/10 and the absolute-gap relaxation time is ten.
The exact profile is d(t)=(1/2)(9/10)ᵗ. At t=15 it is approximately 0.10295; at t=16 it is approximately 0.09265. The required mixing time is sixteen.
The lazy modification has off-diagonal probabilities 19/40 and eigenvalue 1/20. Its distance after one transition is 1/40, already below 1/10, whereas its time-zero distance is 1/2.
The original large ordinary gap does not certify rapid discrete-time mixing. Laziness removes this particular near-periodic obstruction and reduces the requested mixing time to one.
The living notebook
Explore the model
Interactive lab · markov
Watch a negative mode change sign
Pause and Reset, then set p=0.95,q=0.95. Predict transient P(state 1)=0.95 after one Step and 0.095 after two, oscillating around stationary 0.50. For a fresh comparison Reset and set p=0.47,q=0.47; this representable lazy-style comparison has multiplier 0.06, so one Step gives 0.47 and two give 0.4982. The exact lazy value 0.475 is outside the slider’s 0.01 increments, so this second run is a nearby comparison, not that exact modified kernel.
- P(state 1) at step 0
- 0%
- Stationary P(state 1)
- 60%
- Chain behavior
- Monotone convergence
Exact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forA negative eigenvalue near −1 controls slow alternating convergence; a stationary readout alone cannot reveal its decay rate.
Pause & reason
Check your understanding
Your turn
Try it yourself
For p=0.10,q=0.20, compute πmin, both gaps, exact worst-case d(4), and the sharper L²-derived upper bound shown in the second display.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The nonconstant eigenvalue is positive. Use the larger stationary mass as the exact two-state prefactor and √(πmin⁻¹−1)/2 for the bound.
Compare with the solution
The stationary law is (2/3,1/3), hence πmin=1/3. The eigenvalue is 7/10, so both gaps are 3/10. Exact d(4)=(2/3)(7/10)⁴=2401/15000, approximately 0.160067. The L²-derived bound is (√2/2)(7/10)⁴, approximately 0.169776, which is valid but larger. Positivity of this mode makes the ordinary and absolute gaps agree without requiring a laziness modification.
Keep the idea
Take it with you
- Ordinary gap and absolute gap can differ sharply in discrete time.
- This module’s relaxation time is the reciprocal absolute gap.
- Laziness removes negative eigenvalues while slowing positive modes per transition.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.