From Pairwise Balance to Stationarity
Stationarity is a statement about the total probability entering each state after one transition. Detailed balance gives a stronger, pairwise statement: for every pair of states, the probability mass flowing from the first to the second equals the mass flowing in the reverse direction. Summing those pairwise equalities over all possible incoming states gives stationarity. The argument is short, but its hypotheses and the direction of the sum matter, especially when the state space is countable.
Let \((X_n)_{n\ge0}\) be a discrete-time Markov chain on a finite or countable state space \(S\), defined on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), with transition matrix \(P=(p_{ij})\). Each \(p_{ij}=\mathbb{P}(X_{n+1}=j\mid X_n=i)\) is nonnegative, and each row sums to \(1\). Write \(\mu_n\) for the row-vector distribution of \(X_n\). The transition rule gives \(\mu_{n+1}=\mu_nP\). A probability distribution \(\pi\) is stationary if \(\pi P=\pi\), so that starting the chain with distribution \(\pi\) leaves its distribution unchanged after one step.
For a countable state space, the key operation is summing over all \(i\in S\). The terms are nonnegative, so the sum is well-defined; the row-stochastic condition then evaluates the resulting sum. No finite-state assumption is needed.
Proof. Fix any state \(j\in S\). The probability mass at \(j\) after one transition from initial distribution \(\pi\) is \((\pi P)(j)=\sum_{i\in S}\pi(i)p_{ij}\). By detailed balance, each summand equals \(\pi(j)p_{ji}\), so \[ (\pi P)(j)=\sum_{i\in S}\pi(j)p_{ji} =\pi(j)\sum_{i\in S}p_{ji} =\pi(j). \] The last sum is the sum of row \(j\) of \(P\), which equals \(1\). Thus \((\pi P)(j)=\pi(j)\) for every \(j\), proving \(\pi P=\pi\). For countable \(S\), the sums contain only nonnegative terms, and the row sum is \(1\), so the same calculation is valid. \(\square\)
This proof uses the incoming transitions to a fixed state \(j\): it sums over the possible states \(i\) the chain could have occupied before arriving at \(j\). After applying detailed balance, that sum becomes the total of row \(j\), not column \(j\). Confusing those two sums can invalidate the argument.
Using Reversible Weights
Often the candidate distribution is not given in normalized form. Instead, one finds nonnegative weights \(w(i)\) satisfying the balance equations and then checks whether their total is finite. The following result states exactly when normalization works.
Proof. Since \(Z\) is positive and finite, \(\pi(i)=w(i)/Z\) is nonnegative and \(\sum_i\pi(i)=\sum_iw(i)/Z=1\). Dividing the weight balance equation by \(Z\) gives \(\pi(i)p_{ij}=\pi(j)p_{ji}\) for every pair \(i,j\). The theorem therefore gives \(\pi P=\pi\). \(\square\)
In a finite chain, any nonzero collection of finite nonnegative weights has a finite positive sum and can be normalized. On an infinite state space, finiteness of the sum is an additional requirement; balance equations alone do not guarantee a probability distribution.
Worked Example: A Two-State Chain
Consider \(S=\{1,2\}\) with \[ P=\begin{bmatrix}\tfrac7{10}&\tfrac3{10}\\[2pt]\tfrac15&\tfrac45\end{bmatrix}. \] The first row sums to \(\tfrac7{10}+\tfrac3{10}=1\), and the second to \(\tfrac15+\tfrac45=1\). Check detailed balance for \(\pi=(\tfrac25,\tfrac35)\), then verify stationarity directly.
For the pair \(1,2\), the two balanced quantities are \[ \pi(1)p_{12}=\tfrac25\cdot\tfrac3{10}=\tfrac3{25}, \qquad \pi(2)p_{21}=\tfrac35\cdot\tfrac15=\tfrac3{25}. \] For a pair with equal indices, the equation \(\pi(i)p_{ii}=\pi(i)p_{ii}\) holds automatically. Thus all detailed-balance equations hold. The proposed distribution is normalized because \(\tfrac25+\tfrac35=1\).
Check by multiplication: The first component of \(\pi P\) is \(\tfrac25\cdot\tfrac7{10}+\tfrac35\cdot\tfrac15=\tfrac7{25}+\tfrac3{25}=\tfrac25\). The second is \(\tfrac25\cdot\tfrac3{10}+\tfrac35\cdot\tfrac45=\tfrac3{25}+\tfrac{12}{25}=\tfrac35\). Hence \(\pi P=(\tfrac25,\tfrac35)=\pi\).
Worked Example: A Three-State Chain
On \(S=\{1,2,3\}\), take \[ P=\begin{bmatrix} \tfrac12&\tfrac12&0\\ \tfrac14&\tfrac12&\tfrac14\\ 0&\tfrac12&\tfrac12 \end{bmatrix}, \qquad \pi=(\tfrac14,\tfrac12,\tfrac14). \] Every row sums to \(1\): the row sums are \(\tfrac12+\tfrac12=1\), \(\tfrac14+\tfrac12+\tfrac14=1\), and \(\tfrac12+\tfrac12=1\). Also \(\tfrac14+\tfrac12+\tfrac14=1\), so \(\pi\) is a probability distribution.
For states \(1\) and \(2\), the two sides of detailed balance are \(\tfrac14\cdot\tfrac12=\tfrac18\) and \(\tfrac12\cdot\tfrac14=\tfrac18\). For states \(2\) and \(3\), they are \(\tfrac12\cdot\tfrac14=\tfrac18\) and \(\tfrac14\cdot\tfrac12=\tfrac18\). For states \(1\) and \(3\), both transition probabilities are \(0\), so both sides equal \(0\). The equations with identical states hold automatically; reversing any pair repeats one of these checks.
For a direct incoming-mass check, the mass entering state \(1\) is \(\tfrac14\cdot\tfrac12+\tfrac12\cdot\tfrac14+\tfrac14\cdot0=\tfrac18+\tfrac18+0=\tfrac14\). The mass entering state \(2\) is \(\tfrac14\cdot\tfrac12+\tfrac12\cdot\tfrac12+\tfrac14\cdot\tfrac12=\tfrac18+\tfrac14+\tfrac18=\tfrac12\). The mass entering state \(3\) is \(\tfrac14\cdot0+\tfrac12\cdot\tfrac14+\tfrac14\cdot\tfrac12=0+\tfrac18+\tfrac18=\tfrac14\). These are precisely the three components of \(\pi\), confirming \(\pi P=\pi\).
Why Countable Sums Cause No Extra Problem
The theorem’s proof applies to countably many states because, for a fixed \(j\), the incoming terms \(\pi(i)p_{ij}\) are nonnegative. Detailed balance replaces each of them by \(\pi(j)p_{ji}\); summing then uses the row of transitions out of \(j\), whose total is \(1\). The sum is therefore finite and equals \(\pi(j)\). The argument does not require summing a conditionally convergent series or rearranging signed terms.
Worked Example: A Countable Birth-and-Death Chain
Let \(S=\{0,1,2,\ldots\}\). From state \(0\), the chain stays at \(0\) with probability \(\tfrac23\) and moves to \(1\) with probability \(\tfrac13\). From each state \(n\ge1\), it moves to \(n-1\) with probability \(\tfrac23\) and to \(n+1\) with probability \(\tfrac13\). At state \(0\), the row sum is \(\tfrac23+\tfrac13=1\); for every \(n\ge1\), the row sum is \(\tfrac23+\tfrac13=1\).
Define \(w(n)=2^{-n}\). For the edge from \(n\) to \(n+1\), where \(n\ge0\), the forward weight is \(w(n)\cdot\tfrac13=2^{-n}/3\), while the reverse weight is \(w(n+1)\cdot\tfrac23=2^{-(n+1)}\cdot\tfrac23=2^{-n}/3\). If \(|i-j|>1\), both transition probabilities are zero; if \(i=j\), the balance equation is an identity. Thus the weights satisfy detailed balance for all pairs.
The geometric sum is \(Z=\sum_{n=0}^{\infty}2^{-n}=2\), so the normalized candidate is \(\pi(n)=2^{-(n+1)}\). Its total is \(\sum_{n=0}^{\infty}2^{-(n+1)}=1\). The proposition proves that \(\pi\) is stationary.
Check the incoming sums directly. At state \(0\), incoming mass is \(\pi(0)\cdot\tfrac23+\pi(1)\cdot\tfrac23=\tfrac12\cdot\tfrac23+\tfrac14\cdot\tfrac23=\tfrac13+\tfrac16=\tfrac12=\pi(0)\). For any \(n\ge1\), incoming mass is \[ \pi(n-1)\tfrac13+\pi(n+1)\tfrac23 =2\pi(n)\tfrac13+\tfrac{\pi(n)}2\tfrac23 =\tfrac23\pi(n)+\tfrac13\pi(n) =\pi(n). \] This verifies each coordinate of \(\pi P=\pi\), including the boundary state.
A Common Mistake: Checking the Wrong Distribution
A proposed distribution must balance the total incoming probability at every state. Checking only that its entries sum to \(1\), or checking only one incoming transition, is not enough. For example, consider \[ P=\begin{bmatrix}\tfrac12&\tfrac12\\[2pt]\tfrac14&\tfrac34\end{bmatrix}. \] Both rows sum to \(1\). A mistaken calculation might claim that \((\tfrac12,\tfrac12)\) is stationary because it is normalized and the chain has only two states. But the mass entering state \(1\) is \(\tfrac12\cdot\tfrac12+\tfrac12\cdot\tfrac14=\tfrac14+\tfrac18=\tfrac38\), not \(\tfrac12\). In fact, detailed balance also fails: \(\tfrac12\cdot\tfrac12=\tfrac14\), whereas \(\tfrac12\cdot\tfrac14=\tfrac18\).
Solving the balance equation gives \(\pi(1)\cdot\tfrac12=\pi(2)\cdot\tfrac14\), hence \(\pi(2)=2\pi(1)\). Together with \(\pi(1)+\pi(2)=1\), this gives \(\pi=(\tfrac13,\tfrac23)\). Substitution confirms the answer: the first incoming mass is \(\tfrac13\cdot\tfrac12+\tfrac23\cdot\tfrac14=\tfrac16+\tfrac16=\tfrac13\), and the second is \(\tfrac13\cdot\tfrac12+\tfrac23\cdot\tfrac34=\tfrac16+\tfrac12=\tfrac23\).
Detailed balance is a sufficient condition for stationarity: its pairwise equalities imply the required total incoming equalities by summation. The proof should always make that summation explicit, use the row-stochastic condition for the correct row, and verify that the proposed \(\pi\) is a probability distribution.
Check Your Understanding
For each question, identify the states being summed over and use exact arithmetic in any numerical check.
- In the detailed-balance proof, why does the sum \(\sum_i p_{ji}\) equal \(1\), rather than a column sum of \(P\)?
- For a two-state chain with \(p_{12}=\tfrac14\) and \(p_{21}=\tfrac12\), find a normalized distribution satisfying detailed balance.
- Suppose nonnegative weights satisfy detailed balance on a countable state space but sum to infinity. Which condition needed for the normalization proposition fails?
- For the three-state example, calculate the incoming mass at state \(2\) directly and compare it with \(\pi(2)\).
- Explain why nonnegativity of the terms makes the countable-state summation in the theorem valid.