Tutorials › Stochastic Processes › The Generator and the Small-Time Expansion

Kolmogorov Equations and Transient Behavior · Tutorial 490 of 994

The Generator and the Small-Time Expansion

Learn how the generator determines short-interval transition probabilities, including the leading probability of staying in the current state.

Advanced 10 min read

What You'll Learn

  • Define the generator through right derivatives of transition probabilities at time zero.
  • Derive the entrywise expansion of the transition matrix for a short time interval.
  • Relate off-diagonal generator entries to short-interval jump probabilities.
  • Interpret diagonal entries through exit rates and the probability of no jump.
  • Distinguish ending in the starting state from making no jump.
  • Check generator row sums and verify expansions in finite-state examples.

Why the Generator Describes Short Intervals

For a continuous-time Markov chain, the transition semigroup records probabilities over every time interval, while the generator records their instantaneous rates of change at time zero. The connection between them is local: over a short interval of length \(h\), the probability of moving from one state to another is approximately a rate multiplied by \(h\). The diagonal entries describe the matching first-order loss from the starting state.

Let \((X_t)_{t\ge0}\) be a time-homogeneous, finite-state continuous-time Markov chain on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), with state space \(S\). Its transition probabilities are \(p_{ij}(t)=\mathbb{P}(X_t=j\mid X_0=i)\), understood as transition-kernel values defined for every starting state \(i\). Write \(P(t)=(p_{ij}(t))\). Each \(P(t)\) is row-stochastic: its entries are nonnegative and every row sums to \(1\).

Definition: If the right-hand limits exist and are finite, the generator \(Q=(q_{ij})\) is defined by \(q_{ij}=\lim_{h\downarrow0}(p_{ij}(h)-\mathbf{1}_{\{i=j\}})/h\). The symbol \(\mathbf{1}_{\{i=j\}}\) is \(1\) when \(i=j\) and \(0\) otherwise. For a finite-state chain, the small-time expansion \(P(h)=I+hQ+o(h)\) means that, for every pair \(i,j\), \(p_{ij}(h)=\mathbf{1}_{\{i=j\}}+h q_{ij}+r_{ij}(h)\), where \(r_{ij}(h)/h\to0\) as \(h\downarrow0\).

Here \(o(h)\) is entrywise; because \(S\) is finite, there are only finitely many remainders to control. The generator definition gives the expansion one entry at a time. The stochastic rows of \(P(h)\) additionally force the familiar rate constraints: off-diagonal entries of \(Q\) are nonnegative, and each row of \(Q\) sums to zero.

Theorem (Generator and small-time expansion): Suppose \(S\) is finite and \(P(h)\) is a row-stochastic transition semigroup whose right derivatives at zero exist and are finite. If \(Q\) is defined by the generator limits above, then \(P(h)=I+hQ+o(h)\) entrywise. Moreover, \(q_{ij}\ge0\) for \(j\ne i\), and \(\sum_{j\in S}q_{ij}=0\) for every \(i\). Thus \(q_{ii}=-q_i\), where \(q_i=\sum_{j\ne i}q_{ij}\) is the exit rate from state \(i\).

Proof. Fix \(i,j\). By the definition of \(q_{ij}\), the quotient \((p_{ij}(h)-\mathbf{1}_{\{i=j\}})/h\) tends to \(q_{ij}\). Subtracting \(q_{ij}\) from this quotient shows that \(r_{ij}(h)=p_{ij}(h)-\mathbf{1}_{\{i=j\}}-h q_{ij}\) satisfies \(r_{ij}(h)/h\to0\). This proves the entrywise expansion. If \(j\ne i\), then \(p_{ij}(h)\ge0\), so \(q_{ij}=\lim_{h\downarrow0}p_{ij}(h)/h\ge0\). Finally, for each \(h\), \(\sum_jp_{ij}(h)=1\). Subtracting \(\sum_j\mathbf{1}_{\{i=j\}}=1\), dividing by \(h\), and taking the limit gives \(\sum_j q_{ij}=0\); the sum has finitely many terms, so its limit is the sum of the limits. Therefore \(q_{ii}=-\sum_{j\ne i}q_{ij}=-q_i\). \(\square\)

Reading the Off-Diagonal and Diagonal Entries

For distinct states \(i\) and \(j\), the expansion says \(p_{ij}(h)=h q_{ij}+o(h)\). Thus \(q_{ij}\) is the first-order rate of transition from \(i\) to \(j\). The total rate of leaving \(i\) is \(q_i=\sum_{j\ne i}q_{ij}\), and the diagonal expansion is \(p_{ii}(h)=1-hq_i+o(h)\).

There is an important distinction in that last expression. \(p_{ii}(h)\) is the probability of being in \(i\) at the end of the interval; it is not, in general, the probability of making no jump during the interval. A chain can leave \(i\) and return before time \(h\). For a finite-state CTMC with exponential holding times, the probability of no jump from \(i\) during the interval is \(e^{-q_i h}\), which has the same first-order expansion \(1-q_i h+o(h)\). Return paths affect endpoint probabilities only at smaller order than \(h\).

Short-interval interpretation: In the standard finite-state CTMC construction, if \(q_i>0\), the holding time in \(i\) is exponential with rate \(q_i\), and the next state is \(j\ne i\) with probability \(q_{ij}/q_i\). Hence the probability that the first jump by time \(h\) is to \(j\) is \(h q_{ij}+O(h^2)\). The probability of two or more jumps by time \(h\) is \(O(h^2)\). If \(q_i=0\), the chain has no jump from \(i\).

Proof. Let \(\Lambda=\max_{i\in S}q_i\), which is finite because \(S\) is finite. When \(q_i>0\), the probability that the first jump from \(i\) occurs by time \(h\) and goes to \(j\) is \(\int_0^h e^{-q_i s}q_{ij}\,ds=q_{ij}(1-e^{-q_i h})/q_i\). For \(x\ge0\), \(0\le x-(1-e^{-x})=\int_0^x(1-e^{-u})\,du\le\int_0^x u\,du=x^2/2\), where \(1-e^{-u}\le u\). Taking \(x=q_i h\) shows that this first-jump probability differs from \(q_{ij}h\) by at most \(q_{ij}q_i h^2/2\). The probability of a first jump by time \(h\) is at most \(\Lambda h\), since \(1-e^{-q_i h}\le q_i h\le\Lambda h\). Given a first jump, the conditional probability of a second jump before time \(h\) is at most \(\Lambda h\), by the same bound at the state then occupied. Thus the probability of at least two jumps is at most \(\Lambda^2h^2\). For \(j\ne i\), the event that the endpoint is \(j\) and the event that the first jump is to \(j\) can differ only on paths with at least two jumps. This proves the claimed first-order interpretation. When \(q_i=0\), all \(q_{ij}=0\) for \(j\ne i\), and the holding time is infinite, so the assertion also holds in that case. \(\square\)

Worked Short-Time Expansions

Worked Example: A Two-State Chain

On \(S=\{0,1\}\), let the rate from \(0\) to \(1\) be \(2\), and the rate from \(1\) to \(0\) be \(3\). The generator is \(Q=\begin{bmatrix}-2&2\\3&-3\end{bmatrix}\). Before using it, check the generator rows: \(-2+2=0\) and \(3-3=0\).

Applying the expansion entry by entry gives \(p_{00}(h)=1-2h+o(h)\), \(p_{01}(h)=2h+o(h)\), \(p_{10}(h)=3h+o(h)\), and \(p_{11}(h)=1-3h+o(h)\). Starting from state \(0\), the first-order probability of leaving during the interval is \(2h+o(h)\), while the first-order probability of a jump to \(1\) is also \(2h+o(h)\), since \(1\) is the only possible destination.

Check: The first-order coefficients in row \(0\) sum to \(-2+2=0\); those in row \(1\) sum to \(3-3=0\). This agrees with the zero row sums of \(Q\), as required for the first-order terms to preserve total probability.

Worked Example: Three States with Different Exit Rates

Consider \(S=\{0,1,2\}\) and \(Q=\begin{bmatrix}-2&2&0\\1&-4&3\\0&5&-5\end{bmatrix}\). The row sums are \(-2+2+0=0\), \(1-4+3=0\), and \(0+5-5=0\); all off-diagonal entries are nonnegative. From state \(1\), the exit rate is \(1+3=4\).

The row of short-time probabilities starting at \(1\) is \((p_{10}(h),p_{11}(h),p_{12}(h))=(h+o(h),\,1-4h+o(h),\,3h+o(h))\). In particular, the leading probability of a jump to \(0\) is \(h\), and the leading probability of a jump to \(2\) is \(3h\). Adding the two off-diagonal probabilities gives \(4h+o(h)\), the leading probability of leaving state \(1\).

Check: The first-order coefficients in this row add to \(1-4+3=0\). The total exit rate computed from the two possible destinations, \(1+3=4\), equals \(-q_{11}=4\). The row therefore gives consistent off-diagonal jump rates and diagonal loss.

Worked Example: Returning Before the Interval Ends

Take two states \(a,b\), with rate \(1\) from \(a\) to \(b\) and rate \(2\) from \(b\) to \(a\). The generator is \(Q=\begin{bmatrix}-1&1\\2&-2\end{bmatrix}\), whose row sums are \(-1+1=0\) and \(2-2=0\). Starting at \(a\), the probability of making no jump by time \(h\) is \(e^{-h}\).

There is also a positive probability of making exactly two jumps and ending at \(a\): the first jump occurs at time \(s\), the second at time \(u\), and no third jump occurs before \(h\), where \(0<s<u<h\). The probability of these paths is \(\int_0^h\int_s^h e^{-s}\,2e^{-2(u-s)}e^{-(h-u)}\,du\,ds\). Combining the exponential factors gives \(2e^{-h}\int_0^h\int_s^h e^{s-u}\,du\,ds\). The inner integral is \(1-e^{s-h}\), and therefore the expression equals \(2e^{-h}\int_0^h(1-e^{s-h})\,ds=2e^{-h}(h-1+e^{-h})\).

Using \(e^{-h}=1-h+h^2/2+O(h^3)\), the exactly-two-jump probability is \(h^2+O(h^3)\). It is positive for \(h>0\), so the endpoint probability \(p_{aa}(h)\) is not simply the no-jump probability \(e^{-h}\): it also includes paths that return to \(a\). Nonetheless, both have first-order expansion \(1-h+o(h)\), in agreement with \(q_{aa}=-1\).

Check: The two-jump contribution is of order \(h^2\), so it does not change the coefficient of \(h\). The expansion \(p_{aa}(h)=1-h+o(h)\) is consistent with the exit rate \(q_a=1\), while the exact no-jump probability is \(e^{-h}\).

A Common Mistake and Why the Expansion Matters

Common mistake. From the exit rate \(q_i\), one might write \(p_{ii}(h)=e^{-q_i h}\) exactly. That formula gives the probability of no jump from \(i\), not necessarily the probability of ending in \(i\). In the two-state example above, the wrong calculation would set \(p_{aa}(h)=e^{-h}\). The correction is to include paths that leave and return: \(p_{aa}(h)=1-h+o(h)\), and the exactly-two-jump paths alone contribute \(h^2+O(h^3)\). The exact no-jump formula and the endpoint probability agree to first order, but need not agree at a fixed positive \(h\).

The small-time expansion is the bridge from rate specifications to transition probabilities over short intervals. It also explains why each generator row must sum to zero: to first order, probability lost from the diagonal is exactly distributed among the off-diagonal destinations. The expansion is local; it does not by itself give a formula for \(P(t)\) at arbitrary times. That calculation requires additional arguments about the semigroup.

Check Your Understanding

Use the generator definition and the row-sum condition to answer the questions.

  1. If \(q_{ij}=5\) for \(j\ne i\) and all other off-diagonal rates in row \(i\) are zero, what are \(q_i\) and \(q_{ii}\)?
  2. For a row with \(q_{i1}=2\), \(q_{i2}=4\), and no other off-diagonal rates, write the first-order expansion for \(p_{ii}(h)\) and the two listed transition probabilities.
  3. Why does the generator definition give an \(o(h)\) remainder for each entry of \(P(h)\)? Why does finiteness of \(S\) matter when summing the row expansions?
  4. In the two-state return example, why is \(p_{aa}(h)\) not exactly the probability of no jump?
  5. What does a zero exit rate \(q_i=0\) imply for the row of \(Q\) and, in the standard finite-state CTMC construction, for motion from state \(i\)?