Beyond the adiabatic limit in systems with fast environments: a $\tau$-leaping algorithm

Abstract

We propose a $\tau$-leaping simulation algorithm for stochastic systems subject to fast environmental changes. Similar to conventional $\tau$-leaping the algorithm proceeds in discrete time steps, but as a principal addition it captures environmental noise beyond the adiabatic limit. The key idea is to treat the input rates for the $\tau$-leaping as (clipped) Gaussian random variables with first and second moments constructed from the environmental process. In this way, each step of the algorithm retains environmental stochasticity to sub-leading order in the time scale separation between system and environment. We test the algorithm on several toy examples with discrete and continuous environmental states, and find good performance in the regime of fast environmental dynamics. At the same time, the algorithm requires significantly less computing time than full simulations of the combined system and environment. In this context we also discuss several methods for the simulation of stochastic population dynamics in time-varying environments with continuous states.

22 pages, 6 figures

References in corpus (7)

• Statistical physics of social dynamics
• Fixation in finite populations evolving in fluctuating environments
• Evolution of a Fluctuating Population in a Randomly Switching Environment
• Stochastic amplification in an epidemic model with seasonal forcing
• Noise-induced multistability in chemical systems: Discrete vs Continuum modeling
• Population Dynamics in a Changing Environment: Random versus Periodic Switching
• Efficient analysis of stochastic gene dynamics in the non-adiabatic regime using piecewise deterministic Markov processes

Cited by in corpus (1)

• Evolutionary Dynamics in a Varying Environment: Continuous versus Discrete Noise