Tutorials › Stochastic Processes › A Stationary Chain That Is Not Reversible

Reversibility and Detailed Balance · Tutorial 261 of 994

A Stationary Chain That Is Not Reversible

See how a chain can preserve its stationary distribution while its transitions and path probabilities are not invariant under time reversal.

Advanced 10 min read

What You'll Learn

  • Check stationarity of a finite chain by multiplying a row distribution by its transition matrix.
  • Verify why the uniform distribution is stationary for a doubly stochastic transition matrix.
  • Test detailed balance directly on a directed three-state cycle.
  • Compute the reversed transition matrix from a stationary distribution.
  • Compare the probabilities of a path and its reversal under stationarity.

Stationarity Does Not Require Reversibility

A Markov chain can preserve a distribution even while its transitions have a preferred direction. A three-state cycle makes this distinction especially clear: probability mass is spread uniformly across the states at every time, but each transition moves in the same direction. Reversibility asks for more than preservation of a distribution; it asks whether the stationary chain looks the same when time is reversed.

Work on a probability space \((\Omega,\mathcal{F},\mathbb{P})\). Consider a discrete-time Markov chain \((X_n)_{n\ge0}\) with finite state space \(S\), transition matrix \(P=(p_{ij})\), and filtration \(\mathcal{F}_n=\sigma(X_0,\ldots,X_n)\). The matrix is row-stochastic: \(p_{ij}=\mathbb{P}(X_{n+1}=j\mid X_n=i)\), and each row sums to \(1\). Distributions are row vectors, so the distribution at time \(n+1\) is the distribution at time \(n\) multiplied on the right by \(P\).

Definition: A probability row vector \(\pi\) is stationary for \(P\) if \(\pi P=\pi\). A chain started with \(X_0\) distributed as \(\pi\) is stationary: its distribution remains \(\pi\) at every time. For a stationary distribution with \(\pi_i>0\) for every state, the chain is reversible if it satisfies detailed balance, \(\pi_i p_{ij}=\pi_j p_{ji}\) for all \(i,j\). Equivalently, its stationary transition law agrees with the transition law obtained by reversing time.

Stationarity balances total probability entering and leaving each state. Detailed balance requires a stronger, pair-by-pair balance: the stationary probability flow from \(i\) to \(j\) must equal the flow from \(j\) to \(i\). The difference matters because stationarity can hold through an overall circulation around a cycle, without any pairwise balance.

A Uniform Stationary Distribution on a Directed Cycle

Use the state space \(S=\{0,1,2\}\), with transitions moving deterministically from \(0\) to \(1\), from \(1\) to \(2\), and from \(2\) to \(0\). The transition matrix is

$$ P= \begin{bmatrix} 0&1&0\\ 0&0&1\\ 1&0&0 \end{bmatrix}. $$

Each row sums to \(1\): the three row sums are \(0+1+0=1\), \(0+0+1=1\), and \(1+0+0=1\). Each column also sums to \(1\). This column check is what makes the uniform distribution stationary.

Theorem: If \(P\) is a row-stochastic transition matrix on a finite state space and every column of \(P\) sums to \(1\), then the uniform distribution \(\pi_i=1/|S|\) is stationary.

Proof. Let \(m=|S|\), and let \(\pi\) be the row vector whose \(m\) entries are all \(1/m\). It is a probability distribution because its entries are nonnegative and sum to \(m(1/m)=1\). For each state \(j\), the \(j\)th entry of \(\pi P\) is \((\pi P)_j=\sum_{i\in S}(1/m)p_{ij}=(1/m)\sum_{i\in S}p_{ij}=1/m=\pi_j\), using the column-sum assumption. Thus \(\pi P=\pi\), as required. \(\square\)

Worked Example: Stationary but Not Reversible

For the directed cycle above, test whether the uniform distribution \(\pi=(\tfrac13,\tfrac13,\tfrac13)\) is stationary, and then test detailed balance.

First check that \(\pi\) is a probability row vector: its entries sum to \(\tfrac13+\tfrac13+\tfrac13=1\). Multiplication gives

$$ \pi P = \left( \tfrac13\cdot0+\tfrac13\cdot0+\tfrac13\cdot1,\ \tfrac13\cdot1+\tfrac13\cdot0+\tfrac13\cdot0,\ \tfrac13\cdot0+\tfrac13\cdot1+\tfrac13\cdot0 \right) = \left(\tfrac13,\tfrac13,\tfrac13\right)=\pi. $$

So the distribution is stationary. But detailed balance fails for the pair \(0,1\): the stationary transition probability from \(0\) to \(1\) is \(\pi_0p_{01}=\tfrac13\cdot1=\tfrac13\), while the stationary transition probability from \(1\) to \(0\) is \(\pi_1p_{10}=\tfrac13\cdot0=0\). Since these values differ, this stationary chain is not reversible.

Check: the row sums of \(P\) are \(1,1,1\), and the stationary calculation gives all three entries of \(\pi P\) as \(\tfrac13\). The failed detailed-balance equation uses the actual reverse entry \(p_{10}=0\).

The example demonstrates the distinction directly. Each time step moves the state, but if the initial state is uniform, every state is equally likely at every later time. Preservation of the distribution does not imply that the chain retraces its transitions with the same probabilities.

Computing the Reverse Chain

For a stationary chain whose distribution has positive entries, the one-step transition probability of the reversed chain from \(i\) to \(j\) is the conditional probability that the previous state was \(j\), given that the current state is \(i\). The joint probability of being at \(j\) and then moving to \(i\) is \(\pi_jp_{ji}\); the probability of being at \(i\) is \(\pi_i\). Therefore the reversed transition matrix \(P^{\mathrm{rev}}=(p^{\mathrm{rev}}_{ij})\) is given by the following formula.

Reversal formula: If \(\pi P=\pi\) and \(\pi_i>0\) for every state, then \(p^{\mathrm{rev}}_{ij}=\dfrac{\pi_jp_{ji}}{\pi_i}\). Each row of \(P^{\mathrm{rev}}\) sums to \(1\), since \(\sum_j p^{\mathrm{rev}}_{ij}=\sum_j\pi_jp_{ji}/\pi_i=(\pi P)_i/\pi_i=1\).

In particular, detailed balance holds exactly when \(p^{\mathrm{rev}}_{ij}=p_{ij}\) for all \(i,j\): multiplying that equality by \(\pi_i\) gives \(\pi_i p_{ij}=\pi_jp_{ji}\), and dividing the detailed-balance equation by the positive \(\pi_i\) gives the reverse implication.

Worked Example: The Reverse of the Directed Cycle

Compute the reversed matrix for the stationary directed cycle. Every entry of \(\pi\) is \(\tfrac13>0\), so the reversal formula applies. For each pair \(i,j\), it reduces to \(p^{\mathrm{rev}}_{ij}=p_{ji}\). Thus

$$ P^{\mathrm{rev}} = \begin{bmatrix} 0&0&1\\ 1&0&0\\ 0&1&0 \end{bmatrix}. $$

Its row sums are \(0+0+1=1\), \(1+0+0=1\), and \(0+1+0=1\), so it is a valid transition matrix. It moves in the opposite direction from \(P\): for example, \(p_{01}=1\) but \(p^{\mathrm{rev}}_{01}=0\). Hence \(P^{\mathrm{rev}}\ne P\), confirming that the chain is not reversible.

The reversed matrix also preserves the uniform distribution. Direct multiplication gives \((\tfrac13,\tfrac13,\tfrac13)P^{\mathrm{rev}}=(\tfrac13,\tfrac13,\tfrac13)\): each coordinate is the sum of one term equal to \(\tfrac13\) and two terms equal to \(0\). This does not make the original chain reversible; it says the reversed chain has the same stationary distribution.

The reversal formula follows from conditioning, and also identifies what equality of the two transition matrices means. In a reversible chain, the probability of a transition \(i\to j\), when measured under stationarity, matches the probability of the reversed transition \(j\to i\). In the directed cycle, the reversed transition is absent, so equality cannot hold.

Path Probabilities Reveal the Direction

A stationary chain is reversible only if reversing whole paths leaves their probabilities unchanged. To see why detailed balance gives this property, consider a finite path \(x_0,x_1,\ldots,x_n\) under a stationary initial distribution. Its probability is \(\pi_{x_0}\prod_{r=0}^{n-1}p_{x_rx_{r+1}}\). If detailed balance holds, then \(\pi_{x_r}p_{x_rx_{r+1}}=\pi_{x_{r+1}}p_{x_{r+1}x_r}\) at each step. Multiplying these equalities yields \(\pi_{x_0}\prod_{r=0}^{n-1}p_{x_rx_{r+1}}=\pi_{x_n}\prod_{r=0}^{n-1}p_{x_{r+1}x_r}\), because the intermediate factors of \(\pi\) cancel. The right side is the probability of the reversed path \(x_n,\ldots,x_0\) under the same stationary chain. Thus detailed balance implies equality of forward and reversed path probabilities.

Worked Example: A Path and Its Reversal

For the deterministic directed cycle started from \(\pi=(\tfrac13,\tfrac13,\tfrac13)\), compare the path \(0,1,2,0\) with its reversal \(0,2,1,0\). The probability of the first path is

$$ \mathbb{P}(X_0=0,X_1=1,X_2=2,X_3=0) =\pi_0p_{01}p_{12}p_{20} =\tfrac13\cdot1\cdot1\cdot1 =\tfrac13. $$

The reversed path has probability \(\mathbb{P}(X_0=0,X_1=2,X_2=1,X_3=0)=\pi_0p_{02}p_{21}p_{10}=\tfrac13\cdot0\cdot0\cdot0=0\). The two path probabilities differ, even though the chain is stationary. This gives a path-level confirmation of the failed detailed-balance calculation.

Check: each of the three transitions in the forward path has probability \(1\). The first transition in the reversed path, \(0\to2\), has probability \(p_{02}=0\), so its full path probability is \(0\).

A Less Extreme Directed Cycle

The phenomenon does not depend on transitions being deterministic. A chain can allow movement both ways and still prefer one direction. The following example has a uniform stationary distribution because each column sums to \(1\), while unequal clockwise and counterclockwise probabilities prevent detailed balance.

Worked Example: Biased Movement Around the Cycle

On \(S=\{0,1,2\}\), from each state move one step clockwise with probability \(\tfrac23\), and one step counterclockwise with probability \(\tfrac13\). The transition matrix is

$$ P= \begin{bmatrix} 0&\tfrac23&\tfrac13\\ \tfrac13&0&\tfrac23\\ \tfrac23&\tfrac13&0 \end{bmatrix}. $$

Each row sums to \(0+\tfrac23+\tfrac13=1\), \(\tfrac13+0+\tfrac23=1\), and \(\tfrac23+\tfrac13+0=1\). Each column also sums to \(0+\tfrac13+\tfrac23=1\), \(\tfrac23+0+\tfrac13=1\), and \(\tfrac13+\tfrac23+0=1\). By the proved column-sum result, \(\pi=(\tfrac13,\tfrac13,\tfrac13)\) is stationary; direct multiplication gives \(\pi P=(\tfrac13,\tfrac13,\tfrac13)\).

For states \(0\) and \(1\), detailed balance would require \(\pi_0p_{01}=\pi_1p_{10}\). Instead, the left side is \(\tfrac13\cdot\tfrac23=\tfrac29\), while the right side is \(\tfrac13\cdot\tfrac13=\tfrac19\). The reverse transition matrix is

$$ P^{\mathrm{rev}}= \begin{bmatrix} 0&\tfrac13&\tfrac23\\ \tfrac23&0&\tfrac13\\ \tfrac13&\tfrac23&0 \end{bmatrix}. $$

Its rows sum to \(0+\tfrac13+\tfrac23=1\), \(\tfrac23+0+\tfrac13=1\), and \(\tfrac13+\tfrac23+0=1\). For instance, the first row follows from \(p^{\mathrm{rev}}_{0j}=p_{j0}\), since all stationary probabilities are equal. So the reverse chain has the clockwise and counterclockwise probabilities interchanged, and \(P^{\mathrm{rev}}\ne P\).

Check: the row and column sums above verify both stochasticity and uniform stationarity. The unequal values \(\tfrac29\) and \(\tfrac19\) independently verify failure of detailed balance.

Common Mistake: Confusing Stationarity with Pairwise Balance

A common error is to check \(\pi P=\pi\), observe that the distribution is preserved, and conclude that detailed balance must hold. The first condition adds together all transitions into a state; it does not require each transition to be balanced by its reverse. For the deterministic cycle, the incorrect calculation might write \(\pi_0p_{01}=\tfrac13\) and then claim the reverse side is also \(\tfrac13\) by using \(p_{01}=1\) again. That substitutes the forward transition probability for the reverse probability. The correct reverse term is \(\pi_1p_{10}=\tfrac13\cdot0=0\). Therefore stationarity holds, but detailed balance fails.

When checking a proposed reversible chain, calculate both \(p_{ij}\) and \(p_{ji}\) for the same pair of states. Do not infer pairwise equality from stationarity, and do not confuse a chain and its reverse merely because they share a stationary distribution. The directed examples show that stationary distributions describe how probability is spread across states, while reversibility describes whether the transition dynamics themselves are unchanged by reversing time.

Check Your Understanding

Use the definitions and calculations in the tutorial to answer the following questions.

  1. For the deterministic directed cycle, calculate \(\pi_2p_{20}\) and \(\pi_0p_{02}\). What do these values say about detailed balance for the pair \(2,0\)?
  2. Why does a finite row-stochastic matrix with column sums \(1\) preserve the uniform distribution?
  3. For the biased cycle, what is \(p^{\mathrm{rev}}_{12}\), and how does it compare with \(p_{12}\)?
  4. For a stationary chain with positive stationary probabilities, explain why \(P^{\mathrm{rev}}=P\) is equivalent to detailed balance.
  5. In the deterministic cycle, why can a path have positive probability while its reversal has probability zero?