Slow Escape as a Test for Mixing
A chain can have a stationary distribution and still take a long time to resemble it when started from an unfavorable state. One way to detect this delay is to find a set that holds substantial stationary probability but has little probability flow across its boundary. A chain started in that set is unlikely to leave quickly, so its distribution remains concentrated there while stationarity assigns positive mass outside. This discrepancy gives a lower bound on total variation distance.
Let \((X_n)_{n\ge0}\) be a discrete-time Markov chain on a finite state space \(S\), defined on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), with transition matrix \(P=(p_{ij})\). The matrix is row-stochastic: \(p_{ij}\ge0\) and \(\sum_{j\in S}p_{ij}=1\) for every \(i\). Let \(\pi\) be a stationary distribution, so \(\pi P=\pi\). For probability distributions \(\nu,\eta\) on \(S\), use the total variation distance \(\|\nu-\eta\|_{\mathrm{TV}}=\sup_{B\subseteq S}|\nu(B)-\eta(B)|\). Define \(d(n)=\max_{i\in S}\|P^n(i,\cdot)-\pi\|_{\mathrm{TV}}\), and for \(0<\varepsilon<1\) define \(t_{\mathrm{mix}}(\varepsilon)=\inf\{n\ge0:d(n)\le\varepsilon\}\), with the infimum set to \(\infty\) if no such integer exists.
The key comparison starts the chain from the conditional distribution \(\pi_A(i)=\pi_i\mathbf{1}_{\{i\in A\}}/\pi(A)\). If escape by time \(n\) is unlikely under this initial distribution, then the probability of being in \(A\) at time \(n\) remains much larger than \(\pi(A)\). Although this start is a distribution rather than a single state, it still gives a lower bound on the worst-case distance: a mixture of starting states cannot be farther from stationarity than the largest distance among its components.
Proof. Write \(\mu=\pi_A\), and start the chain with distribution \(\mu\). For each \(k\ge1\), consider the event that the first exit from \(A\) occurs on step \(k\). The distribution at time \(k-1\), without imposing any restriction on the path, is \(\mu P^{k-1}\). Since \(\mu(i)\le\pi_i/\pi(A)\) for each \(i\), multiplying by the nonnegative matrix \(P^{k-1}\) gives \((\mu P^{k-1})(i)\le(\pi P^{k-1})(i)/\pi(A)=\pi_i/\pi(A)\). The first-exit event is contained in the event \(X_{k-1}\in A,\ X_k\notin A\). Therefore its probability is at most \[ \sum_{i\in A}\sum_{j\notin A}(\mu P^{k-1})(i)p_{ij} \le \frac{Q(A,A^c)}{\pi(A)} =\Phi(A). \] The events that the first exit occurs on steps \(1,\ldots,n\) are disjoint. Adding their probability bounds shows that the probability of any exit by time \(n\) is at most \(n\Phi(A)\). Consequently \((\mu P^n)(A)\ge1-n\Phi(A)\), including \(n=0\), when the chain starts in \(A\) with probability one. By the definition of total variation, \[ \|\mu P^n-\pi\|_{\mathrm{TV}}\ge(\mu P^n)(A)-\pi(A)\ge1-\pi(A)-n\Phi(A). \] Finally, \(\mu P^n=\sum_{i\in A}\mu(i)P^n(i,\cdot)\). The triangle inequality for total variation gives \(\|\mu P^n-\pi\|_{\mathrm{TV}}\le\sum_{i\in A}\mu(i)\|P^n(i,\cdot)-\pi\|_{\mathrm{TV}}\le d(n)\). This proves the bound. If \(\Phi(A)>0\), the stated strict inequality on \(n\) makes the bound strictly greater than \(\varepsilon\). If \(\Phi(A)=0\), the bound reduces to \(d(n)\ge1-\pi(A)\) for every \(n\). \(\square\)
The theorem requires stationarity of \(\pi\) and positive mass \(\pi(A)\), but it does not require irreducibility, aperiodicity, or reversibility. It is a lower bound, not a convergence theorem: a positive escape ratio does not by itself guarantee that the chain mixes. The bound is most informative when \(A\) has substantial stationary mass and a small escape ratio.
Worked Example: A Four-State Chain with a Narrow Bridge
Consider a chain on \(S=\{1,2,3,4\}\), with the following transition matrix. The two pairs \(\{1,2\}\) and \(\{3,4\}\) are joined by a small transition probability.
Each row sums to \(1\): the second and third rows each sum to \(\tfrac14+\tfrac{37}{50}+\tfrac1{100}=\tfrac{25}{100}+\tfrac{74}{100}+\tfrac1{100}=1\), and the first and fourth rows each sum to \(\tfrac34+\tfrac14=1\). The matrix is symmetric, so its column sums also equal \(1\), and \(\pi=(\tfrac14,\tfrac14,\tfrac14,\tfrac14)\) is stationary: each coordinate of \(\pi P\) is the corresponding column sum multiplied by \(\tfrac14\), hence equals \(\tfrac14\).
Take \(A=\{1,2\}\). Then \(\pi(A)=\tfrac12\). The only transition from \(A\) to \(A^c\) is \(2\to3\), so \(Q(A,A^c)=\pi_2p_{23}=\tfrac14\cdot\tfrac1{100}=\tfrac1{400}\). Thus \(\Phi(A)=(\tfrac1{400})/(\tfrac12)=\tfrac1{200}\), and the theorem gives \(d(n)\ge\tfrac12-\tfrac{n}{200}\). For every integer \(n<50\), this bound is greater than \(\tfrac14\). In particular, at \(n=49\), \(d(49)\ge\tfrac12-\tfrac{49}{200}=\tfrac{100}{200}-\tfrac{49}{200}=\tfrac{51}{200}>\tfrac14\). Therefore \(t_{\mathrm{mix}}(\tfrac14)\ge50\). The estimate uses only the small boundary flow; no exact diagonalization of \(P\) is needed.
Checking the Bound on a Two-State Chain
A small example helps distinguish the stationary mass of a set from its escape ratio. The mass determines how much probability must eventually be found outside the set under stationarity; the ratio controls how quickly a start conditioned to be inside can lose that concentration.
Worked Example: One-Step Escape and the Total Variation Bound
Let \(S=\{1,2\}\) and take \[ P=\begin{bmatrix}\tfrac{19}{20}&\tfrac1{20}\\[2pt]\tfrac1{10}&\tfrac9{10}\end{bmatrix}. \] The first row sums to \(\tfrac{19}{20}+\tfrac1{20}=1\), and the second to \(\tfrac1{10}+\tfrac9{10}=1\). A stationary distribution is \(\pi=(\tfrac23,\tfrac13)\): its first coordinate after multiplication is \(\tfrac23\cdot\tfrac{19}{20}+\tfrac13\cdot\tfrac1{10}=\tfrac{38}{60}+\tfrac2{60}=\tfrac23\); its second is \(\tfrac23\cdot\tfrac1{20}+\tfrac13\cdot\tfrac9{10}=\tfrac1{30}+\tfrac3{10}=\tfrac13\). The entries are nonnegative and sum to \(1\).
For \(A=\{1\}\), \(\pi(A)=\tfrac23\), \(Q(A,A^c)=\tfrac23\cdot\tfrac1{20}=\tfrac1{30}\), and \(\Phi(A)=(\tfrac1{30})/(\tfrac23)=\tfrac1{20}\). At \(n=1\), the theorem gives \(d(1)\ge1-\tfrac23-\tfrac1{20}=\tfrac13-\tfrac1{20}=\tfrac{20}{60}-\tfrac3{60}=\tfrac{17}{60}\).
Check by direct calculation. Starting from the conditional distribution on \(A\) means starting at state \(1\). After one step the probability of \(A\) is \(\tfrac{19}{20}\), so its excess over stationary mass is \(\tfrac{19}{20}-\tfrac23=\tfrac{57}{60}-\tfrac{40}{60}=\tfrac{17}{60}\), exactly the bound. Starting instead from state \(2\), the probability of \(A\) after one step is \(\tfrac1{10}\), whose difference from \(\pi(A)\) has absolute value \(\tfrac23-\tfrac1{10}=\tfrac{20}{30}-\tfrac3{30}=\tfrac{17}{30}\). Thus the worst-case distance is at least \(\tfrac{17}{30}\), consistent with the theorem’s weaker, mixture-based lower bound.
Zero Flow and a Necessary Qualification
When \(\Phi(A)=0\), the theorem is useful without dividing by the escape ratio. The zero-flow case can occur even when \(A\) and \(A^c\) both have positive stationary mass. It says that a chain started according to \(\pi\) conditioned on \(A\) never crosses from \(A\) to \(A^c\) at a step sampled under that distribution. The theorem then prevents the worst-case distance from falling below \(1-\pi(A)\).
Worked Example: An Identity Chain Never Mixes to a Nontrivial Stationary Law
Let \(S=\{1,2\}\) and \(P=\begin{bmatrix}1&0\\0&1\end{bmatrix}\). Both rows sum to \(1\). On a probability space with \(X_0\) distributed as \(\pi=(\tfrac12,\tfrac12)\), define \(X_n=X_0\) for all \(n\); this has transition matrix \(P\), and \(\pi P=\pi\). Set \(A=\{1\}\). Then \(\pi(A)=\tfrac12\) and \(Q(A,A^c)=\pi_1p_{12}=\tfrac12\cdot0=0\), so \(\Phi(A)=0\).
The zero-flow part of the theorem gives \(d(n)\ge1-\tfrac12=\tfrac12\) for every \(n\). In fact, starting from state \(1\), the distribution remains concentrated at \(1\), whose total variation distance from \(\pi\) is \(|1-\tfrac12|=\tfrac12\); starting from state \(2\) gives the same distance. Hence \(d(n)=\tfrac12\) for all \(n\), and \(t_{\mathrm{mix}}(\varepsilon)=\infty\) for every \(\varepsilon<\tfrac12\).
Common mistake. It is invalid to use the positive-flow threshold when \(\Phi(A)=0\). For the identity chain just computed, choosing, for example, \(\varepsilon=\tfrac14\) gives \(1-\pi(A)-\varepsilon=\tfrac14>0\). Writing \(n<(\tfrac14)/0\) is undefined, not a mixing-time conclusion. The correct argument is the separate zero-flow conclusion: \(d(n)\ge1-\pi(A)=\tfrac12>\tfrac14\) for every \(n\), so no finite time achieves that accuracy.
More generally, to use a bottleneck set effectively, check the quantities in order: \(\pi(A)>0\) so the conditional start is defined; \(\pi(A)<1\) so there is stationary mass outside; and then whether \(\Phi(A)\) is positive or zero. In the positive case, compare members of the family \(n\) at the same threshold using the strict inequality in the theorem. In the zero case, use the separate all-times bound. This turns a structural observation about a set into a quantitative obstruction to rapid mixing.
Check Your Understanding
Use the definitions and the proved bound to answer each question. When a question asks for a numerical comparison, retain exact fractions.
- For a stationary distribution \(\pi\) and set \(A\) with \(\pi(A)>0\), write the formula for \(Q(A,A^c)\) and explain the interpretation of \(\Phi(A)\).
- Suppose \(\pi(A)=\tfrac13\) and \(\Phi(A)=\tfrac1{30}\). For which nonnegative integers \(n\) does the theorem guarantee \(d(n)>\tfrac12\)?
- Why does the proof bound the probability of a first exit on step \(k\) by comparing \(\mu P^{k-1}\) with \(\pi/\pi(A)\)? State where stationarity is used.
- If \(\pi(A)=\tfrac25\) and \(\Phi(A)=0\), what lower bound holds for \(d(n)\) at every \(n\), and for which \(\varepsilon\) does it rule out finite mixing time?
- Does the slow-escape lower bound require reversibility or aperiodicity? Identify the assumptions it does require.