Skip to content
Shadow Theory

Non-Source Projection and Internal Identifiability

Exact Markovian projection: strong lumpability

Jeremy RodgersSection 10 of 14Original PDF ↗Publication DOI ↗

Exact Markovian projection: strong lumpability

The dynamical descent problem of §5 has a sharp stochastic counterpart, in which the question “does the projected process close?” admits a complete and checkable answer.

Definition 9.1 (Setting and conventions).

Let S={1,…,N}S=\{1,\dots,N\} be a finite source state space and let (Xt)t≥0(X_t)_{t\ge0} be a continuous-time, time-homogeneous Markov chain on SS with generator G=(gij)i,j∈SG=(g_{ij})_{i,j\in S} satisfying the standard conventions gij≥0  (i≠j),∑j∈Sgij=0for every i,(31) g_{ij}\ge0\ \ (i\neq j),\qquad \sum_{j\in S}g_{ij}=0\quad\text{for every }i ,\tag{31} so that gii=−∑j≠igij≤0g_{ii}=-\sum_{j\neq i}g_{ij}\le0. Transition probabilities are P(t)=etGP(t)=e^{tG}, and a row-vector law π(t)\pi(t) obeys the (forward) master equation π˙j(t)=∑i∈Sπi(t) gij.(32) \dot\pi_j(t)=\sum_{i\in S}\pi_i(t)\,g_{ij}.\tag{32} For f ⁣:S→Rf\colon S\to\mathbb R the backward generator acts by (Gf)(i)=∑jgijf(j)(Gf)(i)=\sum_{j}g_{ij}f(j).

Let p ⁣:S→Tp\colon S\to T be a surjection onto a finite set TT, partitioning SS into the blocks Ba:=p−1(a),a∈T.B_a:=p^{-1}(a),\qquad a\in T .

Definition 9.2 (Strong lumpability).

The chain is strongly lumpable with respect to pp if, for every initial source distribution π(0)\pi(0), the projected process (p(Xt))t≥0\bigl(p(X_t)\bigr)_{t\ge0} is a time-homogeneous Markov chain on TT whose transition law depends on the initial data only through the initial block distribution πˉ(0)\bar\pi(0).

The quantifier “for every initial distribution” is the whole point, and it is what distinguishes strong lumpability from the weaker, initial-distribution-dependent notion. A projected process can perfectly well be Markov for one carefully chosen π(0)\pi(0) , e.g. one supported on a single state, or a stationary within-block law , while failing to be Markov for another. Such a coincidence gives an observer a Markov description that is valid for their particular preparation and invalid in general. Strong lumpability demands that the projected description be universally valid, and that is the property that corresponds to genuine closure of the readout dynamics in the sense of §5.

Theorem 9.3 (Strong lumpability).

The chain is strongly lumpable with respect to pp if and only if the block-sum condition holds: ∑j∈Bbgij  =  ∑j∈Bbgi′jfor all blocks Ba,Bb and all i,i′∈Ba.(33) \sum_{j\in B_b}g_{ij}\;=\;\sum_{j\in B_b}g_{i'j} \qquad\text{for all blocks }B_a,B_b\ \text{and all }i,i'\in B_a .\tag{33} In that case the common value gˉab:=∑j∈Bbgij(any i∈Ba)(34) \bar g_{ab}:=\sum_{j\in B_b}g_{ij}\qquad(\text{any }i\in B_a)\tag{34} is well defined and Gˉ=(gˉab)\bar G=(\bar g_{ab}) is a valid Markov generator on TT, which generates the projected chain.

Proof. Necessity. Suppose the projected process is Markov for every initial distribution, with a transition law determined by πˉ(0)\bar\pi(0). Fix a block BaB_a and i,i′∈Bai,i'\in B_a, and consider the two initial distributions δi\delta_i and δi′\delta_{i'}. Both project to the same initial block law, namely δa\delta_a on TT. Hence by hypothesis the projected chains started from δi\delta_i and from δi′\delta_{i'} have the same law, and in particular the same instantaneous transition rate out of aa into any block BbB_b with b≠ab\neq a.

Compute that rate from the source. Starting at ii, for b≠ab\neq a, Pr⁡[p(Xh)=b∣X0=i]=∑j∈Bb(ehG)ij=∑j∈Bb(δij+h gij+O(h2))=h∑j∈Bbgij+O(h2),\Pr\bigl[p(X_h)=b \mid X_0=i\bigr]=\sum_{j\in B_b}\bigl(e^{hG}\bigr)_{ij} =\sum_{j\in B_b}\bigl(\delta_{ij}+h\,g_{ij}+O(h^2)\bigr) =h\sum_{j\in B_b}g_{ij}+O(h^2), where δij=0\delta_{ij}=0 for j∈Bbj\in B_b since i∈Bai\in B_a and Ba∩Bb=∅B_a\cap B_b=\varnothing. Dividing by hh and letting h↓0h\downarrow0, the instantaneous a→ba\to b rate from ii is ∑j∈Bbgij\sum_{j\in B_b}g_{ij}, and likewise from i′i' it is ∑j∈Bbgi′j\sum_{j\in B_b}g_{i'j}. Equality of the two projected laws forces ∑j∈Bbgij=∑j∈Bbgi′j(b≠a).\sum_{j\in B_b}g_{ij}=\sum_{j\in B_b}g_{i'j}\qquad(b\neq a). For b=ab=a the condition follows from (31): since row sums vanish, ∑j∈Bagij=−∑b≠a∑j∈Bbgij=−∑b≠a∑j∈Bbgi′j=∑j∈Bagi′j.\sum_{j\in B_a}g_{ij}=-\sum_{b\neq a}\sum_{j\in B_b}g_{ij}=-\sum_{b\neq a}\sum_{j\in B_b}g_{i'j}=\sum_{j\in B_a}g_{i'j}. This is (33).

Sufficiency. Assume (33), so that gˉab\bar g_{ab} in (34) is well defined.

Gˉ\bar G is a generator. For a≠ba\neq b and i∈Bai\in B_a, every j∈Bbj\in B_b satisfies j≠ij\neq i, so gij≥0g_{ij}\ge0 by (31); hence gˉab=∑j∈Bbgij≥0\bar g_{ab}=\sum_{j\in B_b}g_{ij}\ge0. And the row sums vanish: ∑b∈Tgˉab=∑b∈T∑j∈Bbgij=∑j∈Sgij=0,\sum_{b\in T}\bar g_{ab}=\sum_{b\in T}\sum_{j\in B_b}g_{ij}=\sum_{j\in S}g_{ij}=0 , using that the blocks partition SS. So Gˉ\bar G satisfies (31) on TT.

(1) Closure of the block probabilities. Let π(t)\pi(t) solve (32) and set πˉa(t):=∑i∈Baπi(t)\bar\pi_a(t):=\sum_{i\in B_a}\pi_i(t). Then πˉ˙b(t)=∑j∈Bbπ˙j(t)=∑j∈Bb∑i∈Sπi(t) gij=∑a∈T∑i∈Baπi(t)(∑j∈Bbgij)=∑a∈T∑i∈Baπi(t) gˉab=∑a∈Tπˉa(t) gˉab,\begin{aligned} \dot{\bar\pi}_b(t) &=\sum_{j\in B_b}\dot\pi_j(t) =\sum_{j\in B_b}\sum_{i\in S}\pi_i(t)\,g_{ij} =\sum_{a\in T}\sum_{i\in B_a}\pi_i(t)\Bigl(\sum_{j\in B_b}g_{ij}\Bigr)\\ &=\sum_{a\in T}\sum_{i\in B_a}\pi_i(t)\,\bar g_{ab} =\sum_{a\in T}\bar\pi_a(t)\,\bar g_{ab}, \end{aligned} where the fourth equality is exactly (33): the inner sum is the same number gˉab\bar g_{ab} for every i∈Bai\in B_a, so it factors out of the sum over i∈Bai\in B_a regardless of how π(t)\pi(t) distributes mass within the block. Thus πˉ\bar\pi obeys the closed master equation πˉ˙=πˉGˉ\dot{\bar\pi}=\bar\pi\bar G, for every initial π(0)\pi(0).

(2) Process-level Markovianity. Closure of the one-time marginals is not by itself the Markov property, and we now establish the latter. Let F:={f ⁣:S→R  :  f is constant on each block}\mathcal F:=\{f\colon S\to\mathbb R\;:\;f\ \text{is constant on each block}\}, and for f∈Ff\in\mathcal F write fbf_b for its value on BbB_b. For i∈Bai\in B_a, (Gf)(i)=∑j∈Sgijf(j)=∑b∈T∑j∈Bbgijfb=∑b∈Tfb(∑j∈Bbgij)=∑b∈Tgˉabfb,(35) (Gf)(i)=\sum_{j\in S}g_{ij}f(j)=\sum_{b\in T}\sum_{j\in B_b}g_{ij}f_b=\sum_{b\in T}f_b\Bigl(\sum_{j\in B_b}g_{ij}\Bigr)=\sum_{b\in T}\bar g_{ab}f_b ,\tag{35} using (33) in the last step. The right-hand side depends on ii only through its block aa. Hence Gf∈FGf\in\mathcal F: the subspace F\mathcal F of block-constant functions is invariant under GG. Being finite-dimensional and GG-invariant, F\mathcal F is invariant under etGe^{tG} for all t≥0t\ge0 (expand the exponential; each partial sum preserves F\mathcal F, and F\mathcal F is closed).

Apply this with f=1Bb∈Ff=\mathbf 1_{B_b}\in\mathcal F. Invariance gives that i↦(etG1Bb)(i)=Pr⁡[p(Xt)=b∣X0=i]i\mapsto\bigl(e^{tG}\mathbf 1_{B_b}\bigr)(i)=\Pr[p(X_t)=b\mid X_0=i] is constant on blocks: writing Pˉab(t)\bar P_{ab}(t) for its common value on BaB_a, Pr⁡[p(Xt)=b∣X0=i]=Pˉab(t)for every i∈Ba.(36) \Pr\bigl[p(X_t)=b\mid X_0=i\bigr]=\bar P_{ab}(t)\qquad\text{for every }i\in B_a .\tag{36} Moreover (35) identifies the generator of Pˉ(t)\bar P(t): reading (35) in F≅RT\mathcal F\cong\mathbb R^{T} shows G∣FG|_{\mathcal F} acts as Gˉ\bar G, so Pˉ(t)=etGˉ\bar P(t)=e^{t\bar G}.

Now take any times 0≤t1<⋯<tn<tn+10\le t_1<\dots<t_n<t_{n+1}, any blocks b1,…,bn+1b_1,\dots,b_{n+1}, and any initial law π(0)\pi(0). Write H:={p(Xt1)=b1,…,p(Xtn)=bn}H:=\{p(X_{t_1})=b_1,\dots,p(X_{t_n})=b_n\} for the projected past, and assume Pr⁡[H]>0\Pr[H]>0; histories of zero probability impose no condition, since the Markov property is a statement about conditional probabilities, which are defined only on positive-probability events. (Equivalently, one may phrase the argument with regular conditional distributions, which exist automatically here because SS is finite; nothing changes.)

Let ν\nu denote the conditional law of XtnX_{t_n} given HH, i.e. ν(i):=Pr⁡[Xtn=i∣H]\nu(i):=\Pr[X_{t_n}=i\mid H]. It is a probability measure supported in BbnB_{b_n}, and it may depend intricately on the whole past HH and on π(0)\pi(0). The Markov property of XX at time tnt_n, together with (36), gives Pr⁡[p(Xtn+1)=bn+1 ∣ H]=∑i∈Bbnν(i)  Pr⁡[p(Xtn+1−tn)=bn+1 ∣ X0=i]=∑i∈Bbnν(i)  Pˉbnbn+1(tn+1−tn)=Pˉbnbn+1(tn+1−tn),(37) \begin{aligned} \Pr\bigl[p(X_{t_{n+1}})=b_{n+1}\,\big|\,H\bigr] &=\sum_{i\in B_{b_n}}\nu(i)\;\Pr\bigl[p(X_{t_{n+1}-t_n})=b_{n+1}\,\big|\,X_0=i\bigr]\\[2pt] &=\sum_{i\in B_{b_n}}\nu(i)\;\bar P_{b_nb_{n+1}}(t_{n+1}-t_n)\\[2pt] &=\bar P_{b_nb_{n+1}}(t_{n+1}-t_n), \end{aligned}\tag{37} where the second equality is (36) , the summand is independent of ii , and the third follows because ν\nu is a probability measure on BbnB_{b_n}, so its masses sum to 11. The right-hand side of (37) depends on the past only through bnb_n, and on the times only through tn+1−tnt_{n+1}-t_n. This is exactly the Markov property with time-homogeneous transition semigroup Pˉ(t)=etGˉ\bar P(t)=e^{t\bar G}, and it holds for every π(0)\pi(0) , which is the assertion of strong lumpability.

The crux, and the reason marginal closure alone would not suffice, is visible in (37): the unknown within-block conditional law ν\nu , which is precisely the information the readout discards, and which does depend on the whole past , drops out only because Pˉbnbn+1\bar P_{b_n b_{n+1}} does not depend on which i∈Bbni\in B_{b_n} carries the mass. That independence is (36), i.e. (33). ◻

Corollary 9.4 (Failure of lumpability).

Suppose (33) fails: there are a block BaB_a, states i,i′∈Bai,i'\in B_a, and a block BbB_b with ∑j∈Bbgij≠∑j∈Bbgi′j\sum_{j\in B_b}g_{ij}\neq\sum_{j\in B_b}g_{i'j}. Then the two initial distributions δi\delta_i and δi′\delta_{i'} , which have the same projected initial state aa , generate different instantaneous projected transition rates into bb. Consequently the instantaneous readout state is insufficient for a universal Markov description of the projected process.

Proof. The instantaneous rates were computed in the necessity half of Theorem 9.3 to be ∑j∈Bbgij\sum_{j\in B_b}g_{ij} and ∑j∈Bbgi′j\sum_{j\in B_b}g_{i'j}, which differ by hypothesis. A time-homogeneous Markov generator on TT would assign a single rate gˉab\bar g_{ab} to the projected state aa, contradiction. ◻

Remark 56 (What failure of lumpability does and does not establish).

Corollary 9.4 is a negative result: the instantaneous projected state does not support a universal Markov law. It does not by itself establish that the projected process is governed by any particular memory kernel or non-Markovian equation. Producing such a kernel requires an additional reduction construction , a Nakajima–Zwanzig projection with a chosen within-block reference law, or a completion in the sense of Definition 4.3 , and different choices give different kernels.

It is worth being exact about what determines the within-block conditional law ν\nu, since this is the discarded information. The generator GG together with an initial distribution π(0)\pi(0) determines it: the pair (G,π(0))(G,\pi(0)) fixes the full law of the source process, hence the conditional law of XtX_t given any positive-probability projected history. The generator alone does not: without an initial law , equivalently, without a specification of how the system was prepared , there is no such thing as “the” within-block distribution, and different preparations consistent with the same projected state yield different ν\nu, which is exactly why Corollary 9.4 bites. Compare §5.2(2): in the deterministic setting likewise, (S,p,Φt)(S,p,\Phi_t) does not determine the within-fiber distribution, which must be supplied. What is theorem here is the impossibility; what is modelling is the replacement.

Example 9.6 (A lumpable and a non-lumpable partition).

Let S={1,2,3}S=\{1,2,3\} with blocks Ba={1,2}B_a=\{1,2\} and Bb={3}B_b=\{3\}, so T={a,b}T=\{a,b\}.

(i) Lumpable. Take G=(−3211−2140−4),G=\begin{pmatrix} -3 & 2 & 1\\ 1 & -2 & 1\\ 4 & 0 & -4\end{pmatrix}, which satisfies (31) (nonnegative off-diagonals, zero row sums). Check (33) for the pair i=1,i′=2∈Bai=1,i'=2\in B_a: ∑j∈Bbg1j=g13=1,∑j∈Bbg2j=g23=1.\sum_{j\in B_b}g_{1j}=g_{13}=1,\qquad \sum_{j\in B_b}g_{2j}=g_{23}=1 . These agree, and the BaB_a column sums then agree automatically by the row-sum identity (g11+g12=−3+2=−1=1−2=g21+g22g_{11}+g_{12}=-3+2=-1=1-2=g_{21}+g_{22}). So the chain is strongly lumpable, with projected generator Gˉ=(gˉaagˉabgˉbagˉbb)=(−114−4),\bar G=\begin{pmatrix}\bar g_{aa} & \bar g_{ab}\\ \bar g_{ba} & \bar g_{bb}\end{pmatrix} =\begin{pmatrix}-1 & 1\\ 4 & -4\end{pmatrix}, a valid generator. Note the internal states 11 and 22 exchange probability at rates 22 and 11 and are in no sense equivalent , the readout is genuinely noninjective, its fiber Ba={1,2}B_a=\{1,2\} is a genuine two-element fiber (it is a set of two source states, not a vector space, and no dimension is being asserted of it), and the internal dynamics within the fiber is nontrivial , yet the projected process is exactly Markov for every initial condition. This is the stochastic analogue of Example 5.5 and of Corollary 5.3: closure of the projected law is no evidence of source completeness.

(ii) Non-lumpable. Now perturb a single entry: G′=(−3215−8340−4),G'=\begin{pmatrix} -3 & 2 & 1\\ 5 & -8 & 3\\ 4 & 0 & -4\end{pmatrix}, again a valid generator. Test (33): ∑j∈Bbg1j′=g13′=1,∑j∈Bbg2j′=g23′=3,\sum_{j\in B_b}g'_{1j}=g'_{13}=1,\qquad \sum_{j\in B_b}g'_{2j}=g'_{23}=3, and 1≠31\neq3. By Corollary 9.4, the projected process is not Markov for all initial distributions: started at δ1\delta_1 the instantaneous a→ba\to b rate is 11; started at δ2\delta_2 it is 33. Both have the same projected initial state aa. A readout observer preparing what they take to be “the state aa” will see different escape statistics depending on a source variable they cannot access, and no single Markov generator on {a,b}\{a,b\} reproduces both.

Observe finally the sense in which this is only a negative result, per Remark 9.5: it tells us no Markov law on TT works universally. It does not tell us what does govern the projected process , that depends on the within-block law, which is precisely the discarded information.

Back to top ↑