In a stochastic reaction network setting we consider the problem of tracking the fate of individual molecules. We show that by using the classical large volume limit results, we may approximate the dynamics of a single tracked molecule in a simple and computationally efficient way. We give examples on how this approach may be used to obtain various characteristics of single-molecule dynamics (for instance, the distribution of the number of infections in a single individual in the course of an epidemic or the activity time of a single enzyme molecule). Moreover, we show how to approximate the overall dynamics of species of interest in the full system with a collection of independent single-molecule trajectories, and give explicit bounds for the approximation error in terms of the reaction rates. This approximation, which is well defined for all times, leads to an efficient and fully parallelizable simulation technique for which we provide some numerical examples.
Individual Molecules Dynamics in Reaction Network Models / Cappelletti, D; Rempala, Ga. - In: SIAM JOURNAL ON APPLIED DYNAMICAL SYSTEMS. - ISSN 1536-0040. - 22:2(2023), pp. 1344-1382. [10.1137/21M1459563]
Individual Molecules Dynamics in Reaction Network Models
Cappelletti, D;
2023
Abstract
In a stochastic reaction network setting we consider the problem of tracking the fate of individual molecules. We show that by using the classical large volume limit results, we may approximate the dynamics of a single tracked molecule in a simple and computationally efficient way. We give examples on how this approach may be used to obtain various characteristics of single-molecule dynamics (for instance, the distribution of the number of infections in a single individual in the course of an epidemic or the activity time of a single enzyme molecule). Moreover, we show how to approximate the overall dynamics of species of interest in the full system with a collection of independent single-molecule trajectories, and give explicit bounds for the approximation error in terms of the reaction rates. This approximation, which is well defined for all times, leads to an efficient and fully parallelizable simulation technique for which we provide some numerical examples.File | Dimensione | Formato | |
---|---|---|---|
revised.pdf
accesso aperto
Tipologia:
2. Post-print / Author's Accepted Manuscript
Licenza:
Pubblico - Tutti i diritti riservati
Dimensione
5.48 MB
Formato
Adobe PDF
|
5.48 MB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
https://hdl.handle.net/11583/2982255