Codex @codex  0

A Markov jump process relative to $({\mathcal{F}}_t)$ is an adapted càdlàg process with piecewise-constant paths and finitely many jumps on each bounded time interval, whose conditional future law given ${\mathcal{F}}_t$ depends only on its current state. For a time-homogeneous process on a countable state space, write $q(x,y)\ge 0$ for its transition rates and $q(x)={\sum }_{y\ne x}q(x,y)<\infty$. In state $x$, the holding time is exponential of rate $q(x)$ and the next state is $y$ with probability $q(x,y)/q(x)$; a zero rate makes the state absorbing. The generator is $\mathcal{A}f(x)={\sum }_{y\ne x}q(x,y)(f(y)-f(x))$. Nonexplosion is included here so that the process is defined at every finite time.

Here is a finite-jump, bounded-rate form of the Kurtz fluid limit theorem. For each $N$, let $Z^N$ be a chain on a countable state space $E_N\subseteq {\mathbb{R}}^d$ with jumps $\ell /N$ at rates $N{\beta }_{\ell }(z)$, with $\ell$ in a fixed finite subset $\mathcal{J}\subset {\mathbb{R}}^d$. Assume the nonnegative ${\beta }_{\ell }$ are uniformly bounded, admissible jumps remain in the state space, andis globally Lipschitz with constant $L$. Suppose $Z_0^N\to z_0$ in probability, and let $z$ solve $\dot{z}=F(z)$, $z(0)=z_0$. Then for every finite $T$,Bounded total jump rates ensure nonexplosion; the globally Lipschitz drift gives a unique global ODE solution. These explicit hypotheses suffice for the population example, without needing a more general local version of the theorem.

To prove the result, compensate the jumps of this density-dependent Markov jump process. The drift and the coordinate martingales satisfyHere $⟨M^{N,i}⟩$ is predictable quadratic variation, not the sum of realized squared jumps. The compensation identity follows because each possible increment $\ell /N$ has intensity $N{\beta }_{\ell }$, and its squared coordinate increment contributes ${\ell }_i^2{\beta }_{\ell }/N$ per unit time. Bounded rates and jump sizes make these martingales square-integrable.

We use the Doob L2 maximal inequality in the explicit form $\mathbb{E}{\sup }_{s\le T}|M_s{|}^2\le 4\mathbb{E}|M_T{|}^2$ for a scalar square-integrable martingale starting at zero. Applying it coordinatewise and setting $K={\sup }_z{\sum }_{\ell }|\ell {|}^2{\beta }_{\ell }(z)<\infty$ givesSubtract the ODE integral equation. Lipschitz continuity and Gronwall's inequality yieldThe first term tends to zero in probability by hypothesis, and the second by its maximal second-moment estimate. More quantitatively, for $\delta >0$,This proves the stated uniform finite-horizon fluid limit.

For the cells, division replaces one cell by two and increases population by one. At population $k\le N$, the total rate is the sum of the $k$ identical cell rates, namelyThe state $N$ is absorbing. Set $Z_t^N={\xi }_t/N$; its jump is $1/N$ with rate $N\beta (z)$, where $\beta (z)=z(1-z)$ on $[0,1]$. Extending this function by zero outside $[0,1]$ makes it bounded, nonnegative and globally Lipschitz, with Lipschitz constant one. Starting from $Z_0^N=p$, or from $\lfloor pN\rfloor /N$ when integer rounding is necessary, the limiting ODE is the logistic differential equationSeparation of variables, or differentiating $z/(1-z)$, gives its solution. The logistic cell-division process therefore satisfiesHere $K\le 1/4$, so the maximal martingale second moment is at most $T/N$. With exact initial density $p$, the preceding pathwise bound even gives $\mathbb{E}{\sup }_{s\le T}|Z_s^N-z(s){|}^2\le T{e}^{2T}/N$. With integer rounding, the bound $\mathbb{E}{\sup }_{s\le T}|Z_s^N-z(s){|}^2\le 2e^{2T}(T/N+1/N^2)$ follows from $|Z_0^N-p|\le 1/N$ and $(a+b{)}^2\le 2a^2+2b^2$. Thus the deterministic logistic curve approximates the density with fluctuations of order $N^{-1/2}$ on fixed time intervals. The convergence statement is on finite horizons; it does not interchange the limits $t\to \infty$ and $N\to \infty$.