[Paper Review] Multiscale Modelling of Birth-Death Processes

Tom Kimpson, Domenic P. J. Germano+2 more|arXiv (Cornell University)|Jan 19, 2026

Gene Regulatory Network Analysis0 citations

TL;DR

The paper develops a principled framework for selecting regime-switching thresholds in the Jump–Switch–Flow (JSF) hybrid method to accurately estimate extinction probabilities in multiscale birth-death processes, with a practical heuristic for threshold choice validated on stochastic Lotka–Volterra models.

ABSTRACT

Many biological systems exhibit multiscale dynamics, where some species occur in high copy numbers while others remain rare. This heterogeneity necessitates hybrid modelling approaches: deterministic models are computationally efficient but inaccurate for low-count species, while fully stochastic simulations are accurate but prohibitively expensive. Hybrid methods like the Jump-Switch-Flow (JSF) algorithm address this by simulating low-count species stochastically and high-count species deterministically. However, selecting regime-switching thresholds to control errors for specific observables remains an open challenge. We develop a principled framework for threshold selection targeting extinction probability. We formalise JSF as a piecewise-deterministic Markov process and derive backward equations for extinction under exact and hybrid dynamics. Near extinction boundaries, complex nonlinear dynamics reduce to tractable time-inhomogeneous linear birth-death processes. This structure yields a rigorous error decomposition based on early and late excursions. Isolating the dominant error term motivates a fast, actionable heuristic. We demonstrate via Monte Carlo studies on a stochastic Lotka-Volterra model that our heuristic reliably upper-bounds empirical errors in extinction probability. This enables users to select the smallest threshold that satisfies a target error tolerance. This work paves the way for principled, efficient multiscale modelling and simulation in stochastic biological systems.

Motivation & Objective

• Motivate multiscale modelling in biology where some species are high-count and others are rare.
• Develop a systematic method to choose the JSF threshold to balance accuracy and efficiency, focused on extinction probabilities.
• Formulate extinction analysis using a birth-death framework near extinction boundaries.
• Derive error decompositions and practical heuristics for hybrid approximation errors in extinction probabilities.

Proposed method

• Formalise JSF as a piecewise-deterministic Markov process (PDMP).
• Derive backward equations for extinction under exact and hybrid dynamics.
• Show near-extinction dynamics reduce to time-inhomogeneous linear birth–death processes.
• Decompose error into early and late excursion components and derive an actionable heuristic.
• Provide Monte Carlo validation on stochastic Lotka–Volterra models to bound JSF errors.
• Propose a practical algorithm to compute a threshold that ensures target error tolerance.

Figure 1: Illustration of the critical time $t_{c}$ and point of no return $s_{*}$ . (a) The net growth rate $r(t)=\lambda(t)-\mu(t)$ transitions from negative (death-dominated) to positive (birth-dominated) at $t_{c}$ . The point of no return $s_{*}$ marks the latest time an up-crossing can occur a

Experimental results

Research questions

• RQ1 How can we quantify and bound the error in extinction probability introduced by the JSF hybrid method?
• RQ2 What is a principled way to choose the regime-switching threshold to achieve a desired accuracy-efficiency trade-off?
• RQ3 How do near-extinction dynamics simplify to birth–death processes, enabling analytical extinction computations?
• RQ4 Can a simple heuristic reliably bound the JSF error for practical threshold selection?

Key findings

• Extinction probabilities can be analyzed exactly via branching-process theory once near-extinction dynamics are reduced to birth–death processes.
• An error decomposition separates early and late excursions, enabling targeted threshold optimization.
• A rigorous, though impractical, bound on JSF error is derived and analyzed critically for its limitations.
• A practical heuristic bounds the error mainly by the extinction probability evaluated at the point of no return, enabling straightforward threshold selection.
• Monte Carlo studies on a stochastic Lotka–Volterra model show the heuristic reliably upper-bounds empirical extinction errors.
• The framework enables selecting the smallest threshold that satisfies a prescribed error tolerance, improving computational efficiency while maintaining accuracy.

Figure 2: Comparison of full stochastic and hybrid simulations for the predator–prey model with decoupled predator decay. Parameters: $\alpha=1.10$ , $\beta=0.05$ , $\gamma=0.4$ , initial conditions $(x_{1,0},x_{2,0})=(10,50)$ . Top: 1000 exact SSA trajectories (solid black: mean). Bottom: 1000 JSF