Pith. sign in

REVIEW 1 major objections 5 minor 91 references

Dynamics of fluctuating populations in multi-state switching environments

T0 review · 1 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read In feast-famine cycles with intermediate resource levels, a slow strain's fixation probability can rise above or drop below its value under matched binary switching, and a piecewise-deterministic Markov process approximation predicts the…

desk verdict Solid, careful extension of the binary-switching program to multi-state feast-famine environments; the core multi-state results are well supported by simulations, while the headline comparison to binary environments is tied to one matching convention. read the letter →

arxiv 2608.12208 v1 pith:QTQSSZKJ submitted 2026-08-12 q-bio.PE cond-mat.stat-mechnlin.AOphysics.bio-ph

Mauro Mobilia This is my paper · ORCID
classification q-bio.PEcond-mat.stat-mechnlin.AOphysics.bio-ph MSC 92D1592D2560J27
keywords multi-stateswitchingenvironmentsfeast-faminecyclesfixationprobabilitystochasticpopulationdynamicscarryingcapacityfluctuationspiecewise-deterministicMarkovprocessMoranmodelenvironmentalnoise
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks whether a feast-famine environment with several intermediate resource levels leads to different evolutionary outcomes than the usual two-state (feast/famine) model. It claims that adding intermediate states, modeled as a (2n+1)-state Markov chain of carrying capacities, changes the slow strain's fixation probability in a parameter-dependent way: with a low median carrying capacity K_0, the slow strain fixates more often in the multi-state environment than in a properly matched binary one; with a high K_0, the opposite holds, and the switching bias shifts the crossover. Why it matters: if true, the standard binary coarse-graining of nutrient fluctuations can systematically misestimate fixation and mean fixation time in real feast-famine cycles, and the paper supplies a fast approximation that predicts these quantities without full simulation.

What carries the argument

The central object is the (2n+1)-state switching model, in which the carrying capacity moves up and down a ladder of evenly spaced values between famine and feast, driven by a continuous-time Markov chain. The key approximation is the N-PDMP: a piecewise-deterministic Markov process in which the population size evolves deterministically by a logistic equation between random environmental switches, and whose stationary density is used to average the Moran fixation probability. The load-bearing identity is the mean-transition-time matching condition of Eq. (8), which fixes effective binary switching rates so that the multi-state and binary environments have the same average famine-to-feast and feast-to-famine travel times; this calibration makes the comparison between multi-state and binary fixation probabilities meaningful. The approximation works by replacing the intractable individual-based process with the deterministic-plus-random N-PDMP, then averaging the Moran formula over its stationary density with a rescaled switching rate.

What would settle it

Recompute the multi-state and binary fixation probabilities with the binary environment calibrated to match the autocorrelation time or the stationary variance of the carrying capacity instead of the mean transition time of Eq. (8), and check whether the regions where multi-state fixation exceeds binary fixation, such as low median carrying capacity, persist; if they vanish or flip sign, the reported enhancement of slow-strain fixation in gradual environments is an artifact of the mean-time calibration convention.

Watch

Extended reading notes

Core claim

In multi-state switching environments, the fixation probability of the slower strain is not simply a smoothed version of the binary result. The paper shows that for a (2n+1)-state environment, the slow-strain fixation probability can be larger or smaller than its binary counterpart depending on the median carrying capacity, the environmental bias, and the switching rate, and that the piecewise-deterministic-Markov-process (N-PDMP) approximation quantitatively captures it across slow, intermediate, and fast switching in three- and five-state models. In the slow-switching limit, the fixation probability is a weighted average of Moran fixation probabilities over the 2n+1 carrying capacities; in the fast-switching limit it is governed by an effective carrying capacity that is the harmonic mean of the K_i over the stationary distribution. Intermediate-state effects are strongest when the median carrying capacity is low, because the population then spends enough time at small sizes that demographic noise boosts the slow strain's chances, and they are weaker or reversed when the median carrying capacity is high.

Load-bearing premise

The comparison between multi-state and binary environments assumes that matching only the average time for the carrying capacity to move from famine to feast makes the two environments equivalent, even though the multi-state dwell time is a sum of several exponential waiting times and is not exponentially distributed.

Editorial extensions

If this is right

  • If the N-PDMP approximation of Eq. (19) retains its accuracy for more than two intermediate states, fixation probabilities and mean fixation times in environments with many states can be computed without full individual-based simulation, at a fraction of the cost.
  • The sign of the difference between multi-state and binary fixation probabilities is not fixed: a low median carrying capacity makes the slow strain more likely to fixate in multi-state environments, while a high median carrying capacity makes it less likely, and the environmental bias shifts the crossover point.
  • Increasing the number of intermediate states generally lowers the slow strain's fixation probability at fixed median carrying capacity, but parameter sets exist where the multi-state environment beats the binary one in some switching regimes and loses in others, so binary-model conclusions can miss regime-dependent reversals.
  • The unconditional mean fixation time scales with the inverse selection strength, with a prefactor set by the number of states, the environmental bias, and the median carrying capacity; the N-PDMP approximation captures its dependence on the switching rate, giving direct predictions for time-to-fixation in gradual environments.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If a different 'same environment' calibration is used, such as matching the autocorrelation time or the stationary variance of the carrying capacity instead of the mean famine-to-feast transition time of Eq. (8), the boundaries of the regions where multi-state fixation exceeds binary fixation, including the low-median-capacity enhancement regions, may shift or disappear.
  • The mechanism behind the low-median-capacity enhancement, namely repeated population bottlenecks at intermediate carrying capacities amplifying demographic noise, predicts that a chemostat with a gradual nutrient ramp should fixate the slower strain more often than an abrupt two-state chemostat with the same extreme concentrations and the same mean cycle period, which is a directly testable experi
  • Taking the number of intermediate states to infinity with suitably scaled rates would turn the discrete ladder of carrying capacities into a continuum, providing a bridge to established continuous-environmental-noise models and a test of whether the discrete-state effects reported here survive in the continuum limit.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 5 minor

Summary. The paper studies a well-mixed population of two competing strains, a slow-growing S and a faster F, under a carrying capacity K(t) that switches among 2n+1 environmental states with linearly spaced values K_i. The environment is a Markov chain with nearest-neighbour switching, characterized by rates ν_{2n+1}, bias δ, and asymmetry ε. Using 10^5-realization Gillespie simulations and an N-PDMP approximation that averages the Moran fixation probability over the PDMP stationary density with a rescaled switching rate (Eq. 19), the authors compute the population size distribution, average population size, fixation probability of S, and unconditional mean fixation time. They compare these quantities with an effective binary switching model calibrated by matching mean transition times (Eq. 8). The central findings are that multi-state environments produce a broader PSD, that the slow strain's fixation probability can be higher or lower than in the matched binary environment depending on K_0, δ, n, and ν, and that the N-PDMP approximation accurately captures the fixation probability across slow, intermediate, and fast switching (Figs. 5 and 6).

Significance. If the comparative claims are robust, this is a valuable and experimentally motivated extension of binary switching models to multi-state feast-famine cycles. The paper's strengths are its systematic validation against large-scale Gillespie simulations, the reusable N-PDMP approximation of Eq. (19) for fixation probability and mean fixation time, the analytical slow/fast switching limits in Eqs. (16) and (17), and the open data and code. The main reservation concerns the calibration convention used for the binary comparator, which affects the quantitative comparisons but not the internal validation of the multi-state results.

major comments (1)
  1. [II A, Eq. (8); III B 2 and Fig. 6(e)] The effective binary rates in Eq. (8) are defined by matching only the mean transition times between K− and K+, but the multi-state dwell time from famine to feast is a sum of n exponential waiting times (Erlang) whereas the binary model has exponential dwell times. Since the binary curves in Figs. 5 and 6 are plotted against ν_{2n+1} through eν(ν_{2n+1}), an alternative calibration—for example matching the autocorrelation time or the variance of K—shifts the binary curves horizontally. For the symmetric 5-state case the autocorrelation time changes eν by about 24% relative to Eq. (8). This convention-dependence directly affects quantitative statements such as the crossing between φ_5 and φ_2 in Fig. 6(e) and the location of intermediate-regime thresholds. The manuscript states the matching convention but does not test its robustness. I ask the authors to add a short analysis (or an explicit caveat) demonstrating that the qualitative conclusions of Sec. III B are insensitive to physically reasonable alternative matching rules, or to reframe the comparative claims as convention-specific.
minor comments (5)
  1. [Fig. 8 caption] The caption states 'Here, n=1' for the five-state model; this should read n=2.
  2. [Figs. 7 and 8 captions] The captions refer to the N-PDMP approximation 'see Eq. (19)' for the mean fixation time; the correct reference is Eq. (B2).
  3. [IV, extension paragraph] The text mentions '2nenvironmental states' as a possible extension; this is presumably a typo for '(2n+1)-state' or '2n+1 environmental states'.
  4. [III A, near Eq. (15)] There is a duplicated phrase 'the number of the number of intermediate states' in the sentence explaining the slower convergence for n=2.
  5. [Appendix D] The asserted insensitivity of PSD, φ_{2n+1}, and τ_{2n+1} to ε is supported by a small set of parameter combinations and no error analysis; a sentence explaining why the ε-dependence cancels (or a systematic scan) would make the claim more convincing.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the multi-state fixation predictions are checked against independent Gillespie simulations, and the binary comparison is fixed by an explicit mean-transition-time convention rather than by a parameter fitted to the predicted quantity.

full rationale

The paper's central results are the fixation probabilities φ_{2n+1} and mean fixation times τ_{2n+1} obtained from Gillespie simulations of the individual-based model (Appendix C), with analytic approximations (Eqs. (16), (17), (19), (B2)) constructed from the Moran formula (12), the stationary distribution (9), and the N-PDMP dynamics (14). None of these formulas uses the simulated fixation probabilities as an input; the N-PDMP density in Eq. (19) is generated from the switching process and logistic growth, then compared with full stochastic simulations in Figs. 5 and 6. The effective binary rates (Eq. (8)) are derived by equating the mean transition times K∓→K±, an explicit modeling convention stated in Sec. II A, not a parameter fitted to the fixation data; the comparison is therefore defined by this convention rather than being a hidden fit. The author's prior binary-state papers are cited for background methods and for the binary N-PDMP approach, but the multi-state extension is derived from the environmental generator and independently validated, so the self-citations are not load-bearing. No step reduces by construction to its inputs, and no prediction is equivalent to a fitted parameter.

Assumptions & free parameters 0 free parameters · 8 assumptions · 0 invented entities

The central claims rest on standard stochastic population genetics (Moran model), Markov chain theory, and a set of explicit modeling choices. The main choices that could shape the conclusions are the linear spacing of carrying capacities (Eq. 6) and the mean-transition-time matching to the binary model (Eq. 8); these are flagged in the text.

assumptions (8)
  • standard math Moran fixation probability and mean fixation time formulas (Eqs. 12 and B1) are valid for the constant-size population approximation.
    Used as building blocks for the quenched/annealed and PDMP approximations.
  • standard math The finite-state Markov chain for environmental switching has a unique stationary distribution given by Eq. (9).
    Standard result for irreducible birth-death chains; used to initialize simulations and compute averages.
  • domain assumption Population growth is logistic with per capita death rate N/K (Eq. 2).
    This is the model's ecological core; it links carrying capacity to demographic noise, following Refs. [11,13,18,19].
  • domain assumption Environmental switching is a nearest-neighbour Markov chain with rates nu and omega (Eqs. 4-5).
    Defines the multi-state colored noise; a coarse-grained representation of gradual resource changes.
  • ad hoc to paper Carrying capacities are linearly spaced between K- and K+ according to Eq. (6).
    Chosen for simplicity; the paper notes other distributions are possible, so conclusions about amplitude effects may depend on this choice.
  • ad hoc to paper Effective binary switching rates are defined by matching mean transition times between K- and K+ (Eq. 8).
    This is the calibration convention for all binary-vs-multi-state comparisons; it ignores higher-order dwell-time statistics.
  • domain assumption The parameters nu, delta, epsilon satisfy omega=(1+epsilon)nu and |delta|<1 (Eq. 7).
    Restricts the switching rate matrix to a two-parameter family; epsilon is shown numerically to have little effect.
  • domain assumption Environmental noise and carrying capacity are initialized at stationarity, with K drawn from pi (Appendix C).
    Avoids transient biases; standard in this line of work.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Dynamics of fluctuating populations in multi-state switching environments." pith.science (2026). https://pith.science/paper/QTQSSZKJ

@misc{pith2026260812208,
  author       = {Pith},
  title        = {Pith review of: Dynamics of fluctuating populations in multi-state switching environments},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QTQSSZKJ}},
  note         = {Machine review of arXiv:2608.12208}
}
read the original abstract

Microbial populations generally evolve in fluctuating environments under time-varying conditions. These are often described by binary switching models, sometimes seen as coarse-grained feast-famine cycles, in which resource availability switches abruptly between abundant and scarce conditions. However, experimental studies suggest that feast-famine environments actually exhibit more complex temporal dynamics. Here, we study how two strains, one growing slightly slower than the other, compete for the same resources in fluctuating environments comprising a finite number of intermediate states, each having its own carrying capacity. Environmental switching between these states and their carrying capacities represents gradual changes in nutrient availability. This class of multi-state stochastic switching models can be interpreted as a coarse-grained description of feast--famine cycles and allows us to investigate strain competition under the gradual recovery and depletion of resources. By computational and analytical means, we characterise the population dynamics in these multi-state fluctuating environments. In particular, we study how the switching rates and distribution of carrying capacities affect the population-size statistics, fixation probability, and mean fixation time. By comparing these results with their counterparts in binary environments, we clarify how the frequency and amplitude of environmental fluctuations influence population dynamics in coarse-grained feast-famine cycles.

Figures

Figures reproduced from arXiv: 2608.12208 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Illustration of the the 3-state ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Histograms of the quasi-stationary PSD, [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Average population size [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of the [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: (c). As in ternary environments, ϕ5 at given ν5 de￾creases with δ; see [PITH_FULL_IMAGE:figures/full_fig_p012_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Unconditional mean fixation time in the three-state switching model as a function of [PITH_FULL_IMAGE:figures/full_fig_p018_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Unconditional mean fixation time in the five-state switching model as a function of [PITH_FULL_IMAGE:figures/full_fig_p019_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Supplementary figure showing in dark blue the PSD of the five-state switching model ( [PITH_FULL_IMAGE:figures/full_fig_p020_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Supplementary figure showing in blue [PITH_FULL_IMAGE:figures/full_fig_p020_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

91 extracted references · 80 canonical work pages

  1. [1]

    quenched approximation

    Fixation under slow and fast environmental switching Whenν 2n+1 ≪s(slow switching), environmental switching is much slower than selection dynamics. It is therefore is likely that no switches occur on the evolution- ary timescalet∼1/s(see Eq. (10) and Appendix B). In this regime, the population is thus subject to the carrying capacityK i with a probability...

  2. [2]

    It is generally difficult to obtain analytical results for the fixation in the interme- diate switching regime

    Fixation under intermediate environmental switching Under intermediate environmental switching, where ν2n+1 ∼s(|δ|<1), evolutionary and environmental dynamics take place on similar time scales, andNhas a nontrivial PSD; see 2(b),(c). It is generally difficult to obtain analytical results for the fixation in the interme- diate switching regime. Here, in ad...

  3. [3]

    Thus, in the three-state switching model K(t)∈{K −,K 0,K +}; see Fig

    Fixation in the three-state switching model Whenn= 1, the ternary environmental cycle con- sists of the feast and famine states of carrying capac- ityK ±, and the median state whose carrying capac- ity isK 0. Thus, in the three-state switching model K(t)∈{K −,K 0,K +}; see Fig. 1(a,left). At the start of each simulation, the value ofKis ran- domly allocat...

  4. [4]

    The conditions forϕ 0 3 > ϕ0 2 in the examples of Fig

    when ϕM(K0)> ϕ 0 2, andϕ 0 3≤ϕ 0 2 otherwise. The conditions forϕ 0 3 > ϕ0 2 in the examples of Fig. 4(a)-(c) are there- foreK 0 ≲141.36 (δ= 0),K 0 ≲119.61 (δ= 0.2), and K0 ≲166.50 (δ=−0.2). The horizontal dashed lines in Fig. 5 show thatϕ 3(ν3 ≪s)≈ϕ 0 3, and henceϕ 0 3 is a good approximation in the slow switching regime, withϕ 3(ν3≪s)≈ϕ 0 3 < ϕ2(ν3≪s) w...

  5. [5]

    5.Sfixation probability in the three-state model as a function of the switching rateν 3

    Fixation in the five-state switching model Whenn= 2, the five-fold environmental cycle con- sists of the feast, famine and median states of re- spective carrying capacitiesK ±2 =K ± andK 0, and the intermediate mild/harsh states of carrying capacity 11 10 3 10 2 10 1 100 101 3 0.11 0.19 0.27 0.35S Fixation Probability K0 = 250, = = 0, s = 0.02, x0 = 0.5 P...

  6. [6]

    feast” and “famine

    Background: Competition in binary-state switching models As a background for this study, we review the main fea- tures of the competition dynamics in binary-state switch- ing models; see Refs. [18–21] for further details. In two-state switching models, with cyclically alter- nating mild and harsh conditions, the time-varying car- rying capacityK(t)∈ {K −,...

  7. [7]

    Multi-state environmental variability In this work, environmental variability is encoded in the continuous-time Markov processξ(t) =i∈ {−n,...,0,...,n}that is a (2n+ 1)-state coloured noise (non-zero correlation time). Here, the environmental noiseξ(t) is initialised in its stationary distribution (see below) and is therefore a (2n+ 1)-state generalisatio...

  8. [8]

    Master equation The competition dynamics in the (2n+1)-state switch- ing models is a continuous-time multivariate Markov pro- cess – defined by the transition ratesT ± S/F given by Eq. (2) that satisfies the master equation for the joint probabilityp 2n+1(NS,NF,ξ,t)≡p 2n+1(⃗N,ξ,t) of find- ing the population in configuration ⃗N= (N S,NF ) and environmenta...

Show all 91 references
  1. [9]

    uMFT in static environment: the Moran approximation The uMFT can be computed exactly for a population of constant size ¯Kevolving according to the classical Moran process [8, 53]; see Sec. II C. When the initial fraction ofSindividuals isx 0 =k/ ¯K, its expression, here denote...

  2. [10]

    III B for theSfixation probability, the Moran uMFT, given by Eq

    uMFT under environmental switching Similarly to what is done in Sec. III B for theSfixation probability, the Moran uMFT, given by Eq. (B1), can be used to obtain suitable approximations of the uMFT under (2n+ 1)-state environmental switching, denoted byτ 2n+1 and defined as th...

  3. [11]

    Roughgarden.Theory of Population Genetics and Evo- lutionary Ecology: An Introduction

    J. Roughgarden.Theory of Population Genetics and Evo- lutionary Ecology: An Introduction. Macmillan, New York, USA, 1979

    1979

  4. [12]

    M. E. Hibbing, C. Fuqua, M. R. Parsek, and S. B. Pe- terson. Bacterial competition: Surviving and thriving in the microbial jungle.Nat. Rev. Microbiol., 8(1):15–25, 2010. 21

    2010

  5. [13]

    Widder andet al.Challenges in microbial ecology: building predictive understanding of community function and dynamics.ISME J., 10:2557–2568, 2016

    S. Widder andet al.Challenges in microbial ecology: building predictive understanding of community function and dynamics.ISME J., 10:2557–2568, 2016

    2016

  6. [14]

    Nguyen, J

    J. Nguyen, J. Lara-Guti´ errez, and R. Stocker. Environ- mental fluctuations and their effects on microbial com- munities, populations and individuals.FEMS Microbiol. Rev., 45:fuaa068, 2021

    2021

  7. [15]

    S. A. Smits andet al.Seasonal cycling in the gut micro- biome of the hadza hunter-gatherers of tanzania.Science, 357:802––805, 2017

    2017

  8. [16]

    Cignarella andet al.Intermittent fasting confers pro- tection in cns autoimmunity by altering the gut micro- biota.Cell Metab., 27:1222–1235, 2018

    F. Cignarella andet al.Intermittent fasting confers pro- tection in cns autoimmunity by altering the gut micro- biota.Cell Metab., 27:1222–1235, 2018

    2018

  9. [17]

    Rodr´ ıguez-Verdugo, C

    A. Rodr´ ıguez-Verdugo, C. Vulin, and M. Ackermann. The rate of environmental fluctuations shapes ecological dynamics in a two-species microbial system.Ecol. Lett., 22(5):838–846, 2019

    2019

  10. [18]

    Abdul-Rahman, D

    F. Abdul-Rahman, D. Tranchina, and D. Gresham. Fluc- tuating environments maintain genetic diversity through neutral fitness effects and balancing selection.Mol. Biol. Evol., 38(10):4362–4375, 2021

    2021

  11. [19]

    Ewens.Mathematical Population Genetics

    W.J. Ewens.Mathematical Population Genetics. Springer, New York, 2004

    2004

  12. [20]

    J. F. Crow and M. Kimura.An Introduction to Popu- lation Genetics Theory. Blackburn Press, Caldwell, NJ, USA, 2009

    2009

  13. [21]

    Dobramysl, M

    U. Dobramysl, M. Mobilia, M. Pleimling, and U. C. T¨ auber. Stochastic population dynamics in spatially ex- tended predator–prey systems.J. Phys. A: Math. Theor., 51:063001, 2018

    2018

  14. [22]

    Shibasaki, M

    S. Shibasaki, M. Mobilia, and S. Mitri. Exclusion of the fittest predicts microbial community diversity in fluctuat- ing environments.J. R. Soc. Interface, 18(183):20210613, 2021

    2021

  15. [23]

    Melbinger, J

    A. Melbinger, J. Cremer, and E. Frey. Evolutionary game theory in growing populations.Phys. Rev. Lett., 105(17):178101, 2010

    2010

  16. [24]

    Cremer, A

    J. Cremer, A. Melbinger, and E. Frey. Frey, evolutionary and population dynamics: A coupled approach.Phys. Rev. E, 84:051921, 2011

    2011

  17. [25]

    Cremer, A

    J. Cremer, A. Melbinger, and E. Frey. Growth dynamics and the evolution of cooperation in microbial popula- tions.Scientific Reports, 2:281, 2012

    2012

  18. [26]

    Melbinger, J

    A. Melbinger, J. Cremer, and E Frey. The emergence of cooperation from a single mutant during microbial life cycles.J. R. Soc. Interface, 12:20150171, 2015

    2015

  19. [27]

    J. S. Chuang, O. Rivoire, and S. Leibler. Simp- son’s paradox in a synthetic microbial system.Science, 323:20150171, 2009

    2009

  20. [28]

    Verdon, O

    N. Verdon, O. Popescu, S. Titmuss, and R. J. Allen. Habitat fragmentation enhances microbial collective de- fence.J. R. Soc. Interface, 22:20240611, 2025

    2025

  21. [29]

    Wienand, E

    K. Wienand, E. Frey, and M. Mobilia. Evolution of a Fluctuating Population in a Randomly Switching Envi- ronment.Phys. Rev. Lett., 119(15):158301, 2017

    2017

  22. [30]

    Wienand, E

    K. Wienand, E. Frey, and M. Mobilia. Eco-evolutionary dynamics of a population with randomly switching car- rying capacity.J. R. Soc. Interface, 15(145):20180343, 2018

    2018

  23. [31]

    West and M

    R. West and M. Mobilia. Fixation properties of rock- paper-scissors games in fluctuating populations.J. Theor. Biol., 491:110135, 2020

    2020

  24. [32]

    Taitelbaum, R

    A. Taitelbaum, R. West, M. Assaf, and M. Mobilia. Pop- ulation Dynamics in a Changing Environment: Ran- dom versus Periodic Switching.Phys. Rev. Letters, 125(4):048105, 2020

    2020

  25. [33]

    Mahrt, A

    N. Mahrt, A. Tietze, S. K¨ unzel, S. Franzenburg, C. Bar- bosa, G. Jansen, and H. Schulenburg. Bottleneck size and selection level reproducibly impact evolution of antibiotic resistance.Nat. Ecol. Evol., 5(9):1233–1242, 2021

    2021

  26. [34]

    Taitelbaum, R

    A. Taitelbaum, R. West, M. Mobilia, and M. Assaf. Evo- lutionary dynamics in a varying environment: Continu- ous versus discrete noise.Phys. Rev. Res., 5(2):L022004, 2023

    2023

  27. [35]

    Hern´ andez-Navarro, M

    L. Hern´ andez-Navarro, M. Asker, A. M. Rucklidge, and M. Mobilia. Coupled environmental and demographic fluctuations shape the evolution of cooperative antimi- crobial resistance.J. R. Soc. Interface, 20(208):20230393, 2023

    2023

  28. [36]

    Asker, L

    M. Asker, L. Hern´ andez-Navarro, A. M. Rucklidge, and M. Mobilia. Coexistence of Competing Microbial Strains under Twofold Environmental Variability and Demo- graphic Fluctuations.New J. Phys., 25(12):123010, 2023

    2023

  29. [37]

    Hern´ andez-Navarro, M

    L. Hern´ andez-Navarro, M. Asker, and M. Mobilia. Eco- evolutionary dynamics of cooperative antimicrobial resis- tance in a population of fluctuating volume and size.J. Phys. A: Math. Theor., 57(26):265003, 2024

    2024

  30. [38]

    Asker, M

    M. Asker, M. Swailem, U. C. T¨ auber, and M. Mo- bilia. Fixation and extinction in time-fluctuating spa- tially structured metapopulations.Phys. Rev. Res., 7:043205, 2025

    2025

  31. [39]

    Hern´ andez-Navarro, K

    L. Hern´ andez-Navarro, K. Distefano, U. C. T¨ auber, and Mobilia M. Slow spatial migration can help eradicate cooperative antimicrobial resistance in time-varying en- vironments.PLoS Comput. Biol., 22(3):e1013997, 2026

    2026

  32. [40]

    L. M. Wahl, P. J. Gerrish, and I. Saika-Voivod. Evaluat- ing the impact of population bottlenecks in experimental evolution.Genetics, 162:961, 2002

    2002

  33. [41]

    Patwa and L

    Z. Patwa and L. M. Wahl. Adaptation Rates of Lytic Viruses Depend Critically on Whether Host Cells Survive the Bottleneck.Evolution, 64(4):1166–1172, 2010

    2010

  34. [42]

    M. A. Brockhurst. Population bottlenecks promote co- operation in bacterial biofilms.PLoS One, 2:e634, 2007

    2007

  35. [43]

    Coates, B

    J. Coates, B. R. Park, D. Le, E. S ¸im¸ sek, W. Chaudhry, and M. Kim. Antibiotic-induced population fluctuations and stochastic clearance of bacteria.eLife, 7:e32976, 2018

    2018

  36. [44]

    Ashcroft, P

    P. Ashcroft, P. M. Altrock, and T. Galla. Fixa- tion in finite populations evolving in fluctuating en- vironments.Journal of The Royal Society Interface, 11(100):20140663, 2014

    2014

  37. [45]

    Raatz and A

    M. Raatz and A. Traulsen. Promoting extinction or min- imizing growth? The impact of treatment on trait tra- jectories in evolving populations.Evolution, 77(6):1408– 1421, 2023

    2023

  38. [46]

    Fruet, E

    C. Fruet, E. L. M¨ uller, C. Loverdo, and A.-F. Bitbol. Spatial structure facilitates evolutionary rescue by drug resistance.PLoS Comput Biol, 21(4):e1012861, 2025

    2025

  39. [47]

    Q. He, M. Mobilia, and U. C. Tauber. Coexistence in the two-dimensional may-leonard model with random rates. Eur. Phys. J. B, 82:97, 2011

    2011

  40. [48]

    Szolnoki and M

    A. Szolnoki and M. Perc. Vortices determine the dynam- ics of biodiversity in cyclical interactions with protection spillovers.New J. Phys., 17:113033, 2015

    2015

  41. [49]

    Esmaeili, B

    S. Esmaeili, B. L. Brown, and M. Pleimling. Perturbing cyclic predator-prey systems: How a six-species coarsen- ing system with nontrivial in-domain dynamics responds to sudden changes.Phys. Rev. E, 98:062105, 2018. 22

    2018

  42. [50]

    Thattai and A

    M. Thattai and A. van Oudenaarden. Stochastic gene expression in fluctuating environments.Genetics, 167(1):523, 2004

    2004

  43. [51]

    Kussell, R

    E. Kussell, R. Kishony, N. Q. Balaban, and S. Leibler. Bacterial persistence: A model of survival in changing environments.Genetics, 169(4):1807, 2005

    2005

  44. [52]

    M. Acar, J. T Mettetal, and A. van Oudenaarden. Stochastic switching as a survival strategy in fluctuat- ing environments.Nat. Genet., 40(4):471–475, 2008

    2008

  45. [53]

    Lambert and E

    G. Lambert and E. Kussell. Memory and fintess op- timization of bacteria under fluctuating environments. PLoS Genet., 10(9):e1004556, 2014

    2014

  46. [54]

    Hufton, Y

    P. Hufton, Y. T. Lin, T. Galla, and McKane A. J. Intrin- sic noise in systems with switching environments.Phys. Rev. E, 93:052119, 2016

    2016

  47. [55]

    Traulsen and C

    A. Traulsen and C. Hauert. Stochastic evolutionary game dynamics.Reviews of Nonlinear Dynamics and Complex- ity, 2:25–61, 2009

    2009

  48. [56]

    Kalyuzhny, R

    M. Kalyuzhny, R. Kadmon, and N. M. Shnerb. A neutral theory with environmental stochasticity explains static and dynamic properties of ecological communities.Ecol- ogy Letters, 18(6):572–580, 2015

    2015

  49. [57]

    Meyer and N

    I. Meyer and N. M. Shnerb. Evolutionary dynamics in fluctuating environment.Phys. Rev. Res., 2:023308, 2020

    2020

  50. [58]

    Srinivasan and S

    S. Srinivasan and S. Kjelleberg. Cycles of famine and feast: The starvation and outgrowth strategies of a ma- rine vibrio.J. Biosci., 23:501, 1998

    1998

  51. [59]

    Merritt and S

    J. Merritt and S. Kuehn. Frequency- and amplitude- dependent microbial population dynamics during cycles of feast and famine.Phys. Rev. Lett., 121:098101, 2018

    2018

  52. [60]

    Himeoka and N

    Y. Himeoka and N. Mitarai. Dynamics of bacterial pop- ulations under the feast-famine cycles.Phys. Rev. Res., 2:013372, 2020

    2020

  53. [61]

    Niimi, C

    R. Niimi, C. Furusawa, and Himeoka Y. Population dy- namics of generalist and specialist strategies under feast– famine cycles.PLoS Comput Biol, 22(4):e1014265, 2026

    2026

  54. [62]

    Sanchez and J

    A. Sanchez and J. Gore. Feedback between population and evolutionary dynamics determines the fate of social microbial populations.PLoS Biol, 11:e101547, 2013

    2013

  55. [63]

    L. J. S. Allen.An Introduction to Stochastic Processes with Applications to Biology. Pearson Education, 1 edi- tion, April 2003

    2003

  56. [64]

    Moran.The Statistical Processes of Evolutionary Theory

    P.A.P. Moran.The Statistical Processes of Evolutionary Theory. Oxford, UK: Clarendon, 1962

    1962

  57. [65]

    R. A. Blythe and A. J. McKane. Stochastic models of evolution in genetics, ecology and linguistics.J. Stat. Mech., P07018, 2007

    2007

  58. [66]

    D. T. Gillespie. A general method for numerically simu- lating the stochastic time evolution of coupled chemical reactions.J. Comput. Phys., 202:403, 1976

    1976

  59. [67]

    Antal and I

    T. Antal and I. Scheuring. Fixation of strategies for an evolutionary game in finite populations.Bull. Math. Biol., 68(8):1923–1944, 2006

    1923

  60. [68]

    Marrec, I

    L. Marrec, I. Lamberti, and A.-F. Bitbol. Toward a Uni- versal Model for Spatially Structured Populations.Phys. Rev. Lett., 127(21):218102, 2021

    2021

  61. [69]

    Abbara, L

    A. Abbara, L. Pagani, C. Garc´ ıa-Pareja, and A.-F. Bit- bol. Mutant fate in spatially structured populations on graphs: Connecting models to experiments.PLOS Com- put Biol, 20(9):e1012424, 2024

    2024

  62. [70]

    Moawad, A

    A. Moawad, A. Abbara, and A.-F. Bitbol. Evolution of cooperation in deme-structured populations on graphs. Phys. Rev. E, 109(2):024307, 2024

    2024

  63. [71]

    Yagoobi, N

    S. Yagoobi, N. Sharma, and A. Traulsen. Categoriz- ing update mechanisms for graph-structured metapop- ulations.J. R. Soc. Interface, 20(200):20220769, 2023

    2023

  64. [72]

    H. L. Smith and P. E. Waltman.The theory of the chemo- stat: dynamics of microbial competition. Cambridge Uni- versity Press, 1995

    1995

  65. [73]

    S. Li, N. Abdulkadir, F. Schattenberg, U. Nunes da Rocha, V. Grimm, S. M¨ uller, and Z. Liu. Stabilizing microbial communities by looped mass transfer.Proc. Natl. Acad. Sci. U.S.A., 119:e2117814119, 2022

    2022

  66. [74]

    I. Bena. Dichotomous Markov Noise: Exact results for out-of-equilibrium systems.International Journal of Modern Physics B, 20:2825, 2006

    2006

  67. [75]

    Ridolfi, P

    L. Ridolfi, P. D’Odorico, and F. Laio.Noise-Induced Phe- nomena in the Envionmental Sciences. Cambridge Uni- versity Press, Cambridge, U.K., 2011

    2011

  68. [76]

    Horsthemke and R

    W. Horsthemke and R. Lefever.Noise-Induced Transi- tions. Springer, Berlin, 2006

    2006

  69. [77]

    Melbinger and M

    A. Melbinger and M. Vergassola. The impact of envi- ronmental fluctuations on evolutionary fitness functions. Sci. Rep., 5:15211, 2015

    2015

  70. [78]

    A Blythe and A

    R. A Blythe and A. J McKane. Stochastic models of evolution in genetics, ecology and linguistics.J. Stat. Mech.: Theory Exp., 2007(07):P07018, 2007

    2007

  71. [79]

    Wkb theory of large devia- tions in stochastic populations.J

    M Assaf and B Meerson. Wkb theory of large devia- tions in stochastic populations.J. Phys. A: Math. Theor., 50:263001, 2017

    2017

  72. [80]

    Assaf and M

    M. Assaf and M. Mobilia. Large fluctuations and fixation in evolutionary games.J. Stat. Mech., page P09009, 2010

    2010

  73. [81]

    Assaf and M

    M. Assaf and M. Mobilia. Fixation of a deleterious allele under mutation pressure and finite selection intensity.J. Theor. Biol., 275:93, 2011

    2011

  74. [82]

    G. F. Gause. Experimental studies on the struggle for existence.Science, 131:1292, 1960

    1960

  75. [83]

    G. Hardin. The competitive exclusion principle.J. Exp. Biol., 9(4):389, 1992

    1992

  76. [84]

    M. H. A. Davis. Piecewise-Deterministic Markov Pro- cesses: A General Class of Non-Diffusion Stochastic Mod- els.Journal of the Royal Statistical Society. Series B (Methodological), 46(3):353–388, 1984

    1984

  77. [85]

    M. Mobilia. Polarization and consensus in a voter model under time-fluctuating influences.Physics, 5:517, 2023

    2023

  78. [86]

    Meyer, A

    I. Meyer, A. Taitelbaum, M. Assaf, and N. M. Shnerb. Population dynamics in a time-varying environment with fat-tailed correlations.Phys. Rev. E, 110:L012401, 2024

    2024

  79. [87]

    Assaf, M

    M. Assaf, M. Mobilia, and E. Roberts. Cooperation dilemma in finite populations under fluctuating environ- ments.Physical Review Letters, 111(23):238101, 2013

    2013

  80. [89]

    E. A. Yurtsev, H. X. Chao, M. S. Datta, T. Artemova, and J. Gore. Bacterial cheating drives the population dynamics of cooperative antibiotic resistance plasmids. Mol. Syst. Biol., 9:683, 2013

    2013

  81. [90]

    Abbara and A.-F

    A. Abbara and A.-F. Bitbol. Frequent asymmetric mi- grations suppress natural selection in spatially structured populations.PNAS Nexus, 2(11):pgad392, 2023

    2023

  82. [91]

    M. Mobilia. Supporting data and codes. Research Data Leeds Repository, 2026. URL:https://doi.org/10. 5518/1899

    2026

  83. [100]

    3 (dotted lines) The number of independent samplesR ind is large whenR ≫1 and the switching rate is not too low

    In the limitν 2n+1→0, the average population size is given by⟨N⟩ 2n+1≈⟨K⟩ 2n+1; see Fig. 3 (dotted lines) The number of independent samplesR ind is large whenR ≫1 and the switching rate is not too low. Therefore, within the approach outlined above, theN-PDMP-based approximatio...

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.