Search PubMed⌕ Search
Home/Archive/PubMed

PubMed · 12689112

Stochastic population dynamics: the Poisson approximation.

Hernán G Solari·Mario A Natiello
View original source ↗Export metadata ↓

Abstract

We introduce an approximation to stochastic population dynamics based on almost independent Poisson processes whose parameters obey a set of coupled ordinary differential equations. The approximation applies to systems that evolve in terms of events such as death, birth, contagion, emission, absorption, etc., and we assume that the event-rates satisfy a generalized mass-action law. The dynamics of the populations is then the result of the projection from the space of events into the space of populations that determine the state of the system (phase space). The properties of the Poisson approximation are studied in detail. Especially, error bounds for the moment generating function and the generating function receive particular attention. The deterministic approximation for the population fractions and the Langevin-type approximation for the fluctuations around the mean value are recovered within the framework of the Poisson approximation as particular limit cases. However, the proposed framework allows to treat other limit cases and general situations with small populations that lie outside the scope of the standard approaches. The Poisson approximation can be viewed as a general (numerical) integration scheme for this family of problems in population dynamics.

Explore related subjects

Monte Carlo MethodPoisson DistributionPopulation DynamicsJournal Article

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Hernán G Solari, Mario A Natiello. 2003-03-26. Stochastic population dynamics: the Poisson approximation.. https://doi.org/10.1103/physreve.67.031918

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2025

RadiSeq: a single- and bulk-cell whole-genome DNA sequencing simulator for radiation-damaged cell models.

Objective.To build and validate a simulation framework to perform single-cell and bulk-cell whole genome sequencing simulation of radiation-exposed Monte Carlo (MC) cell models to assist radiation genomics studies.Approach.Sequencing the genomes of radiation-damaged cells can provide useful insight into radiation action for radiobiology research. However, carrying out post-irradiation sequencing experiments can often be challenging, expensive, and time-consuming. Although computational simulations have the potential to provide solutions to these experimental challenges, and aid in designing optimal experiments, the absence of tools currently limits such application. MC toolkits exist to simulate radiation exposures of cell models but there are no tools to simulate single- and bulk-cell sequencing of cell models containing radiation-damaged DNA. Therefore, we aimed to develop a MC simulation framework to address this gap by designing a tool capable of simulating sequencing processes for radiation-damaged cells. Main results.We developed RadiSeq-a multi-threaded whole-genome DNA sequencing simulator written in C++. RadiSeq can be used to simulate Illumina sequencing of radiation-damaged cell models produced by MC simulations. RadiSeq has been validated through comparative analysis, where simulated data were matched against experimentally obtained data, demonstrating reasonable agreement between the two. Additionally, it comes with numerous features designed to closely resemble actual whole-genome sequencing. RadiSeq is also highly customizable with a single input parameter file.Significance.RadiSeq enables the research community to perform complex simulations of radiation-exposed DNA sequencing, supporting the optimization, planning, and validation of costly and time-intensive radiation biology experiments. This framework provides a powerful tool for advancing radiation genomics research.

Monte Carlo Method↗
PubMed2006

Thermal versus mechanical unfolding of ubiquitin.

The authors studied the temperature-induced unfolding of ubiquitin by all-atom Monte Carlo simulations. The unfolding behavior is compared with that seen in previous simulations of the mechanical unfolding of this protein, based on the same model. In mechanical unfolding, secondary-structure elements were found to break in a quite well-defined order. In thermal unfolding, the authors saw somewhat larger event-to-event fluctuations, but the unfolding pathway was still far from random. Two long-lived secondary-structure elements could be identified in the simulations. These two elements have been found experimentally to be the thermally most stable ones. Interestingly, one of these long-lived elements, the first beta-hairpin, was found to break early in the mechanical unfolding simulations. Their combined simulation results thus enable the authors to predict in detail important differences between the thermal and mechanical unfolding behaviors of ubiquitin.

Monte Carlo Method↗