Tutorials › Stochastic Processes › The Embedded Jump Chain

Continuous-Time Markov Chains and Generator Matrices · Tutorial 466 of 994

The Embedded Jump Chain

Learn how a generator determines the discrete-time sequence of states visited at jumps, including the special case of absorbing states.

Advanced 10 min read

What You'll Learn

  • Define jump epochs and the embedded sequence of visited states.
  • Derive the probability of jumping from state i to state j as q_ij divided by the total exit rate.
  • Verify that the embedded transition probabilities form stochastic rows when the exit rate is positive.
  • Distinguish genuine jumps from the convention used to extend the chain at absorbing states.
  • Understand why the jump chain does not by itself determine the holding times.

From Continuous Time to a Sequence of Jumps

A continuous-time Markov chain records both which states are visited and when they are visited. Sometimes the order of visited states is the main object of interest. Recording the state at each successive jump removes the clock times and leaves a discrete-time chain. Its transition probabilities are determined by the off-diagonal rates of the generator, normalized by the total rate of leaving the current state.

Let \((X_t)_{t\ge0}\) be a time-homogeneous continuous-time Markov chain on a finite or countable state space \(S\), defined on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), with filtration \((\mathcal{F}_t)_{t\ge0}\). Write its generator as \(Q=(q_{ij})\). For distinct states \(i\) and \(j\), \(q_{ij}\ge0\) is the rate of a jump from \(i\) to \(j\). The generator’s rows sum to zero, so the total exit rate from \(i\) is \(q_i=-q_{ii}=\sum_{j\ne i}q_{ij}\). We assume these rates are finite. The holding time at \(i\), when \(q_i>0\), is exponential with rate \(q_i\).

Set \(\sigma_0=0\). For each successive jump, define \(\sigma_{n+1}=\inf\{t>\sigma_n:X_t\ne X_{\sigma_n}\}\), whenever this time is finite, and let \(Y_n=X_{\sigma_n}\). Thus \(Y_n\) records the state at the \(n\)th jump epoch, with \(Y_0=X_0\). The sequence is considered only for jump epochs that occur before any accumulation of jump times. If the process reaches a state with no possible departure, there is no further genuine jump epoch.

Definition: The sequence of states observed at successive jump epochs is the embedded jump chain. At a state \(i\) with \(q_i>0\), its one-step transition probability to a distinct state \(j\) is \(r_{ij}=q_{ij}/q_i\), and \(r_{ii}=0\). A state with \(q_i=0\) has no next genuine jump, so this formula does not define a row there.

The zero diagonal in the embedded chain is deliberate: a jump means that the state changes. The diagonal entry \(q_{ii}\) of the generator is not a rate of jumping from \(i\) back to \(i\); it is the negative total exit rate. The embedded chain describes jump destinations, not elapsed time.

Why the Rates Normalize to Jump Probabilities

Theorem: Suppose \(i\in S\) has finite positive exit rate \(q_i\), and the chain’s jump from \(i\) is governed by the competing exponential clocks with rates \(q_{ij}\), \(j\ne i\). Then the probability that the next jump is to \(j\) is \(r_{ij}=q_{ij}/q_i\). These probabilities sum to one over \(j\ne i\). At each finite jump epoch, conditional on the states observed so far and the current state \(i\), the next destination has this same distribution; hence the embedded sequence, while it continues, has the Markov transition probabilities \(r_{ij}\).

Proof. Fix a current state \(i\). For each \(j\ne i\) with \(q_{ij}>0\), let the clock for a jump to \(j\) be exponential with rate \(q_{ij}\), with the clocks independent. Clocks with rate zero can be omitted. The next jump is to \(j\) exactly when the clock for \(j\) rings before all the other clocks. For \(t\ge0\), the probability density that the clock for \(j\) rings at time \(t\) and every other clock has not rung by then is \(q_{ij}e^{-q_{ij}t}\prod_{k\ne i,j}e^{-q_{ik}t}=q_{ij}e^{-q_it}\). Integrating over \(t\ge0\), the probability that \(j\) is the destination is \(\int_0^\infty q_{ij}e^{-q_it}\,dt=q_{ij}/q_i\), since \(q_i>0\). Independent exponential clocks have probability zero of tying, so there is almost surely a unique first clock. Summing the destination probabilities gives \(\sum_{j\ne i}q_{ij}/q_i=q_i/q_i=1\), using \(q_i=\sum_{j\ne i}q_{ij}\). After a jump, the time-homogeneous Markov jump mechanism applies the same clock rule at the new state; consequently, conditional on the jump history and the current state, the next-state probabilities depend only on that current state. This is the Markov property for the embedded sequence. \(\square\)

The theorem requires \(q_i>0\). If \(q_i=0\), nonnegativity of the off-diagonal rates and their sum \(q_i\) imply \(q_{ij}=0\) for every \(j\ne i\). The state has no possible departure, and dividing by \(q_i\) is not valid. This case needs a separate convention if a stochastic matrix on all of \(S\) is desired.

Computing an Embedded Transition Matrix

Worked Example: Three States with Positive Exit Rates

Consider \(S=\{1,2,3\}\) and the generator \[ Q=\begin{bmatrix} -5&2&3\\ 4&-5&1\\ 2&2&-4 \end{bmatrix}. \] First check the generator rows: \(-5+2+3=0\), \(4-5+1=0\), and \(2+2-4=0\). The exit rates are \(q_1=5\), \(q_2=5\), and \(q_3=4\); each is positive.

For state \(1\), the rates to states \(2\) and \(3\) are \(2\) and \(3\), so \(r_{12}=2/5\) and \(r_{13}=3/5\), while \(r_{11}=0\). For state \(2\), \(r_{21}=4/5\), \(r_{23}=1/5\), and \(r_{22}=0\). For state \(3\), \(r_{31}=2/4=1/2\), \(r_{32}=2/4=1/2\), and \(r_{33}=0\). Therefore \[ R=\begin{bmatrix} 0&\tfrac25&\tfrac35\\ \tfrac45&0&\tfrac15\\ \tfrac12&\tfrac12&0 \end{bmatrix}. \]

Check: The embedded rows sum to \(0+2/5+3/5=1\), \(4/5+0+1/5=1\), and \(1/2+1/2+0=1\). For instance, from state \(1\) the next jump is to state \(3\) with probability \(3/5\), even though the holding time at state \(1\) has rate \(5\). The holding rate and the destination probabilities answer different questions.

More generally, when every state has positive exit rate, the embedded transition matrix on \(S\) is row-stochastic: its entries are nonnegative and each row sums to one. If the state space is countable, the row sum is understood as a sum over all distinct destination states; the finite exit-rate assumption ensures that this sum is \(q_i/q_i=1\).

Worked Example: Changing the Clock Speeds

Use the generator from the preceding example, but multiply each row by a different positive factor: multiply row \(1\) by \(2\), row \(2\) by \(1/2\), and row \(3\) by \(3\). The resulting rate table is \[ \widetilde Q=\begin{bmatrix} -10&4&6\\ 2&-\tfrac52&\tfrac12\\ 6&6&-12 \end{bmatrix}. \] Its row sums are \(-10+4+6=0\), \(2-\tfrac52+\tfrac12=0\), and \(6+6-12=0\), so it satisfies the generator row-sum condition. Its exit rates are \(10\), \(5/2\), and \(12\), respectively.

Normalize each off-diagonal row. In row \(1\), the probabilities are \(4/10=2/5\) and \(6/10=3/5\). In row \(2\), they are \(2/(5/2)=4/5\) and \((1/2)/(5/2)=1/5\). In row \(3\), they are \(6/12=1/2\) and \(6/12=1/2\). The embedded matrix is again \[ \widetilde R=\begin{bmatrix} 0&\tfrac25&\tfrac35\\ \tfrac45&0&\tfrac15\\ \tfrac12&\tfrac12&0 \end{bmatrix}. \]

Check: The three rows sum to \(2/5+3/5=1\), \(4/5+1/5=1\), and \(1/2+1/2=1\). For each state, multiplying all its outgoing rates by the same positive factor changes the total exit rate by that factor, so each ratio \(q_{ij}/q_i\) stays the same. The jump chain alone therefore cannot identify these changed holding rates.

States with No Possible Departure

When \(q_i=0\), the process remains at \(i\) forever once it reaches \(i\). The sequence of genuine jump destinations ends there: there is no next jump to record. In some calculations, it is convenient to extend the sequence by declaring that it stays at \(i\) thereafter. This produces an ordinary discrete-time Markov chain with an absorbing row at \(i\), but the added repetitions are a bookkeeping convention, not actual jumps.

Worked Example: A Jump Chain That Terminates

Let \(S=\{a,b,c\}\) and \[ Q=\begin{bmatrix} -3&1&2\\ 0&0&0\\ 2&0&-2 \end{bmatrix}. \] The generator rows sum to zero: \(-3+1+2=0\), \(0+0+0=0\), and \(2+0-2=0\). The exit rates are \(q_a=3\), \(q_b=0\), and \(q_c=2\). Thus \(b\) is absorbing.

At \(a\), the next genuine jump goes to \(b\) with probability \(1/3\) and to \(c\) with probability \(2/3\). At \(c\), the next genuine jump goes to \(a\) with probability \(2/2=1\). There is no genuine next jump from \(b\), so the genuine embedded sequence ends if it reaches \(b\).

If we adopt the absorbing-row convention, the extended transition matrix is \[ R=\begin{bmatrix} 0&\tfrac13&\tfrac23\\ 0&1&0\\ 1&0&0 \end{bmatrix}. \] Its rows sum to \(0+1/3+2/3=1\), \(0+1+0=1\), and \(1+0+0=1\). The second row says that the extended discrete-time sequence remains at \(b\); it does not assert that the continuous-time process makes a jump from \(b\) to itself.

Common Mistake: Treating the Generator Diagonal as a Jump

A frequent error is to divide every entry of a generator row by the exit rate and treat the result as a transition row. In the first example, this would assign the row \((-5/5,2/5,3/5)=(-1,2/5,3/5)\) to state \(1\). Its entries sum to \(-1+2/5+3/5=0\), and the negative entry cannot be a probability. The correction is to use only off-diagonal rates as destination rates, set the embedded diagonal to zero, and divide those rates by \(q_1\). This gives \((0,2/5,3/5)\), a valid probability row.

The embedded chain is useful when a question concerns the order of visited states—for example, which state is visited next, or which sequence of states occurs before absorption. It is not enough for questions about the time spent in each state. The generator contains both kinds of information: its off-diagonal rates determine the jump destinations after normalization, while the total exit rate determines the holding-time distribution. Keeping those roles separate prevents a jump probability from being mistaken for a rate.

Check Your Understanding

For each question, distinguish the generator’s rates from the embedded chain’s probabilities, and check that any probability row you compute sums to one.

  1. At a state \(i\), the only outgoing rates are \(q_{ij}=3\) and \(q_{ik}=2\). What are the embedded probabilities of the two destinations?
  2. A generator row has off-diagonal rates \(q_{i1}=1\), \(q_{i2}=4\), and all other off-diagonal rates zero. What must its diagonal entry be, and what is its embedded row?
  3. Why is the formula \(r_{ii}=q_{ii}/q_i\) not a valid way to define a self-transition probability?
  4. If all outgoing rates from a state are multiplied by \(5\), which features of the jump mechanism change, and which embedded probabilities remain unchanged?
  5. For a state with \(q_i=0\), explain the difference between the genuine jump sequence and the convention that appends an absorbing row to a transition matrix.