Source author record

Alan J. McKane

Alan J. McKane appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2016arXiv

Demographic noise can reverse the direction of deterministic selection

Deterministic evolutionary theory robustly predicts that populations displaying altruistic behaviours will be driven to extinction by mutant cheats that absorb common benefits but do not themselves contribute. Here we show that when demographic stochasticity is accounted for, selection can in fact act in the reverse direction to that predicted deterministically, instead favouring cooperative behaviors that appreciably increase the carrying capacity of the population. Populations that exist in larger numbers experience a selective advantage by being more stochastically robust to invasions than smaller populations, and this advantage can persist even in the presence of reproductive costs. We investigate this general effect in the specific context of public goods production and find conditions for stochastic selection reversal leading to the success of public good producers. This insight, developed here analytically, is missed by both the deterministic analysis as well as standard game theoretic models that enforce a fixed population size. The effect is found to be amplified by space; in this scenario we find that selection reversal occurs within biologically reasonable parameter regimes for microbial populations. Beyond the public good problem, we formulate a general mathematical framework for models that may exhibit stochastic selection reversal. In this context, we describe a stochastic analogue to r-K theory, by which small populations can evolve to higher densities in the absence of disturbance.

preprint2016arXiv

Intrinsic noise in systems with switching environments

We study individual-based dynamics in finite populations, subject to randomly switching environmental conditions. These are inspired by models in which genes transition between on and off states, regulating underlying protein dynamics. Similarly switches between environmental states are relevant in bacterial populations and in models of epidemic spread. Existing piecewise-deterministic Markov process (PDMP) approaches focus on the deterministic limit of the population dynamics while retaining the randomness of the switching. Here we go beyond this approximation and explicitly include effects of intrinsic stochasticity at the level of the linear-noise approximation. Specifically we derive the stationary distributions of a number of model systems, in good agreement with simulations. This improves existing approaches which are limited to the regimes of fast and slow switching.

preprint2016arXiv

Simplification and analysis of a model of social interaction in voting

A recently proposed model of social interaction in voting is investigated by simplifying it down into a version that is more analytically tractable and which allows a mathematical analysis to be performed. This analysis clarifies the interplay of the different elements present in the system --- social influence, heterogeneity and noise --- and leads to a better understanding of its properties. The origin of a regime of bistability is identified. The insight gained in this way gives further intuition into the behaviour of the original model.

preprint2016arXiv

Staged Models for Interdisciplinary Research

Modellers of complex biological or social systems are often faced with an invidious choice: to use simple models with few mechanisms that can be fully analysed, or to construct complicated models that include all the features which are thought relevant. The former ensures rigour, the latter relevance. We discuss a method that combines these two approaches, beginning with a complex model and then modelling the complicated model with simpler models. The resulting "chain" of models ensures some rigour and relevance. We illustrate this process on a complex model of voting intentions, constructing a reduced model which agrees well with the predictions of the full model. Experiments with variations of the simpler model yield additional insights which are hidden by the complexity of the full model. This approach facilitated collaboration between social scientists and physicists -- the complex model was specified based on the social science literature, and the simpler model constrained to agree (in core aspects) with the complicated model.

preprint2016arXiv

Stochastic epidemic dynamics on extremely heterogeneous networks

Networks of contacts capable of spreading infectious diseases are often observed to be highly heterogeneous, with the majority of individuals having fewer contacts than the mean, and a significant minority having relatively very many contacts. We derive a two-dimensional diffusion model for the full temporal behavior of the stochastic susceptible-infectious-recovered (SIR) model on such a network, by making use of a time-scale separation in the deterministic limit of the dynamics. This low-dimensional process is an accurate approximation to the full model in the limit of large populations, even for cases when the time-scale separation is not too pronounced, provided the maximum degree is not of the order of the population size.

preprint2016arXiv

Suppressing escape events in maps of the unit interval with demographic noise

We explore the properties of discrete-time stochastic processes with a bounded state space, whose deterministic limit is given by a map of the unit interval. We find that, in the mesoscopic description of the system, the large jumps between successive iterates of the process can result in probability leaking out of the unit interval, despite the fact that the noise is multiplicative and vanishes at the boundaries. By including higher-order terms in the mesoscopic expansion, we are able to capture the non-Gaussian nature of the noise distribution near the boundaries, but this does not preclude the possibility of a trajectory leaving the interval. We propose a number of prescriptions for treating these escape events, and we compare the results with those obtained for the metastable behavior of the microscopic model, where escape events are not possible. We find that, rather than truncating the noise distribution, censoring this distribution to prevent escape events leads to results which are more consistent with the microscopic model. The addition of higher moments to the noise distribution does not increase the accuracy of the final results, and it can be replaced by the simpler Gaussian noise.

preprint2016arXiv

The role of demographic stochasticity in a speciation model with sexual reproduction

Recent theoretical studies have shown that demographic stochasticity can greatly increase the tendency of asexually reproducing phenotypically diverse organisms to spontaneously evolve into localised clusters, suggesting a simple mechanism for sympatric speciation. Here we study the role of demographic stochasticity in a model of competing organisms subject to assortative mating. We find that in models with sexual reproduction, noise can also lead to the formation of phenotypic clusters in parameter ranges where deterministic models would lead to a homogeneous distribution. In some cases, noise can have a sizeable effect, rendering the deterministic modelling insufficient to understand the phenotypic distribution.

preprint2015arXiv

Impact of commuting on disease persistence in heterogeneous metapopulations

We use a stochastic metapopulation model to study the combined effects of seasonality and spatial heterogeneity on disease persistence. We find a pronounced effect of enhanced persistence associated with strong heterogeneity, intermediate coupling strength and moderate seasonal forcing. Analytic calculations show that this effect is not related with the phase lag between epidemic bursts in different patches, but rather with the linear stability properties of the attractor that describes the steady state of the system in the large population limit.

preprint2015arXiv

Models of genetic drift as limiting forms of the Lotka-Volterra competition model

The relationship between the Moran model and stochastic Lotka-Volterra competition (SLVC) model is explored via timescale separation arguments. For neutral systems the two are found to be equivalent at long times. For systems with selective pressure, their behavior differs. It is argued that the SLVC is preferable to the Moran model since in the SLVC population size is regulated by competition, rather than arbitrarily fixed as in the Moran model. As a consequence, ambiguities found in the Moran model associated with the introduction of more complex processes, such as selection, are avoided.

preprint2015arXiv

Modes of competition and the fitness of evolved populations

Competition between individuals drives the evolution of whole species. Although the fittest individuals survive the longest and produce the most offspring, in some circumstances the resulting species may not be optimally fit. Here, using theoretical analysis and stochastic simulations of a simple model ecology, we show how the mode of competition can profoundly affect the fitness of evolved species. When individuals compete directly with one another, the adaptive dynamics framework provides accurate predictions for the number and distribution of species, which occupy positions of maximal fitness. By contrast, if competition is mediated by the consumption of a common resource then demographic noise leads to the stabilization of species with near minimal fitness.

preprint2015arXiv

Stationary solutions for metapopulation Moran models with mutation and selection

We construct an individual-based metapopulation model of population genetics featuring migration, mutation, selection and genetic drift. In the case of a single `island', the model reduces to the Moran model. Using the diffusion approximation and timescale separation arguments, an effective one-variable description of the model is developed. The effective description bears similarities to the well-mixed Moran model with effective parameters which depend on the network structure and island sizes, and is amenable to analysis. Predictions from the reduced theory match the results from stochastic simulations across a range of parameters. The nature of the fast-variable elimination technique we adopt is further studied by applying it to a linear system, where it provides a precise description of the slow-dynamics in the limit of large timescale separation.

preprint2015arXiv

The statistics of fixation times for systems with recruitment

We investigate the statistics of the time taken for a system driven by recruitment to reach fixation. Our model describes a series of experiments where a population is confronted with two identical options, resulting in the system fixating on one of the options. For a specific population size, we show that the time distribution behaves like an inverse Gaussian with an exponential decay. Varying the population size reveals that the timescale of the decay depends on the population size and allows the critical population number, below which fixation occurs, to be estimated from experimental data.

preprint2015arXiv

WKB versus generalized van Kampen system-size expansion: the stochastic logistic equation

Stochastic fluctuations are central to the understanding of extinction dynamics. In the context of population models they allow for the description of the transition from the vicinity of a non-trivial fixed point of the deterministic dynamics to a trivial fixed point, where the population has become extinct. To characterize analytically the fluctuations of a given stochastic population model, one can operate within the so-called linear-noise approximation. Here the fluctuations are taken to be Gaussian, and the phenomenon of extinction is so rare as to be negligible, for all practical purposes. When the size of the population becomes small, non-Gaussian fluctuations are instead found. Two analytical schemes are in principle available to determine the nature of the distribution of associated fluctuations: the generalized van Kampen system-size expansion beyond the conventional order of approximation and a WKB-like method. Here we investigate the accuracy of these two different approximation schemes, with reference to a simple stochastic process that converges to the logistic equation in the deterministic limit.

preprint2014arXiv

Fast-mode elimination in stochastic metapopulation models

We investigate the stochastic dynamics of entities which are confined to a set of islands, between which they migrate. They are assumed to be one of two types, and in addition to migration, they also reproduce and die. Systems which fall into this class are common in biology and social science, occurring in ecology, population genetics, epidemiology, biochemistry, linguistics, opinion dynamics, and other areas. In all these cases the governing equations are intractable, consisting as they do of multidimensional Fokker-Planck equations or, equivalently, coupled nonlinear stochastic differential equations with multiplicative noise. We develop a methodology which exploits a separation in time scales between fast and slow variables to reduce these equations so that they resemble those for a single island, which are amenable to analysis. The technique is generally applicable, but we choose to discuss it in the context of population genetics, in part because of the extra features that appear due to selection. The idea behind the method is simple, its application systematic, and the results in very good agreement with simulations of the full model for a range of parameter values.

preprint2014arXiv

Intrinsic noise and two-dimensional maps: Quasicycles, quasiperiodicity, and chaos

We develop a formalism to describe the discrete-time dynamics of systems containing an arbitrary number of interacting species. The individual-based model, which forms our starting point, is described by a Markov chain, which in the limit of large system sizes is shown to be very well-approximated by a Fokker-Planck-like equation, or equivalently by a set of stochastic difference equations. This formalism is applied to the specific case of two species: one predator species and its prey species. Quasi-cycles --- stochastic cycles sustained and amplified by the demographic noise --- previously found in continuous-time predator-prey models are shown to exist, and their behavior predicted from a linear noise analysis is shown to be in very good agreement with simulations. The effects of the noise on other attractors in the corresponding deterministic map, such as periodic cycles, quasiperiodicity and chaos, are also investigated.

preprint2014arXiv

Noise-Induced Bistable States and Their Mean Switching Time in Foraging Colonies

We investigate a type of bistability where noise not only causes transitions between stable states, but also constructs the states themselves. We focus on the experimentally well-studied system of ants choosing between two food sources to illustrate the essential points, but the ideas are more general. The mean time for switching between the two bistable states of the system is calculated. This suggests a procedure for estimating, in a real system, the critical population size above which bistability ceases to occur.

preprint2014arXiv

The theory of individual based discrete-time processes

A general theory is developed to study individual based models which are discrete in time. We begin by constructing a Markov chain model that converges to a one-dimensional map in the infinite population limit. Stochastic fluctuations are hence intrinsic to the system and can induce qualitative changes to the dynamics predicted from the deterministic map. From the Chapman-Kolmogorov equation for the discrete-time Markov process, we derive the analogues of the Fokker-Planck equation and the Langevin equation, which are routinely employed for continuous time processes. In particular, a stochastic difference equation is derived which accurately reproduces the results found from the Markov chain model. Stochastic corrections to the deterministic map can be quantified by linearizing the fluctuations around the attractor of the map. The proposed scheme is tested on stochastic models which have the logistic and Ricker maps as their deterministic limits.

preprint2013arXiv

Demographic noise and resilience in a semi-arid ecosystem model

The scarcity of water characterising drylands forces vegetation to adopt appropriate survival strategies. Some of these generate water-vegetation feedback mechanisms that can lead to spatial self-organisation of vegetation, as it has been shown with models representing plants by a density of biomass, varying continuously in time and space. However, although plants are usually quite plastic they also display discrete qualities and stochastic behaviour. These features may give rise to demographic noise, which in certain cases can influence the qualitative dynamics of ecosystem models. In the present work we explore the effects of demographic noise on the resilience of a model semi-arid ecosystem. We introduce a spatial stochastic eco-hydrological hybrid model in which plants are modelled as discrete entities subject to stochastic dynamical rules, while the dynamics of surface and soil water are described by continuous variables. The model has a deterministic approximation very similar to previous continuous models of arid and semi-arid ecosystems. By means of numerical simulations we show that demographic noise can have important effects on the extinction and recovery dynamics of the system. In particular we find that the stochastic model escapes extinction under a wide range of conditions for which the corresponding deterministic approximation predicts absorption into desert states.

preprint2013arXiv

Intrinsic noise and discrete-time processes

A general formalism is developed to construct a Markov chain model that converges to a one-dimensional map in the infinite population limit. Stochastic fluctuations are therefore internal to the system and not externally specified. For finite populations an approximate Gaussian scheme is devised to describe the stochastic fluctuations in the non-chaotic regime. More generally, the stochastic dynamics can be captured using a stochastic difference equation, derived through an approximation to the Markov chain. The scheme is demonstrated using the logistic map as a case study.

preprint2013arXiv

Multi-Compartment Linear Noise Approximation

The ability to quantify stochastic fluctuations present in biochemical and other systems is becoming increasing important. Analytical descriptions of these fluctuations are attractive, as stochastic simulations are computationally expensive. Building on previous work, the linear noise approximation is developed for biochemical models with many compartments e.g. cells. The procedure is then implemented in the software package COPASI. This technique is illustrated with two simple examples and is then applied to a more realistic biochemical model. Expressions for the noise, given in the form of covariances matrices, are presented.

preprint2013arXiv

Synchronisation of stochastic oscillators in biochemical systems

A formalism is developed which describes the extent to which stochastic oscillations in biochemical models are synchronised. It is based on the calculation of the complex coherence function within the linear noise approximation. The method is illustrated on a simple example and then applied to study the synchronisation of chemical concentrations in social amoeba. The degree to which variation of rate constants in different cells and the volume of the cells affects synchronisation of the oscillations is explored, and the phase lag calculated. In all cases the analytical results are shown to be in good agreement with those obtained through numerical simulations.

preprint2013arXiv

The linear noise approximation for reaction-diffusion systems on networks

Stochastic reaction-diffusion models can be analytically studied on complex networks using the linear noise approximation. This is illustrated through the use of a specific stochastic model, which displays traveling waves in its deterministic limit. The role of stochastic fluctuations is investigated and shown to drive the emergence of stochastic waves beyond the region of the instability predicted from the deterministic theory. Simulations are performed to test the theoretical results and are analyzed via a generalized Fourier transform algorithm. This transform is defined using the eigenvectors of the discrete Laplacian defined on the network. A peak in the numerical power spectrum of the fluctuations is observed in quantitative agreement with the theoretical predictions.

preprint2012arXiv

Demographic noise and piecewise deterministic Markov processes

We explore a class of hybrid (piecewise deterministic) systems characterized by a large number of individuals inhabiting an environment whose state is described by a set of continuous variables. We use analytical and numerical methods from nonequilibrium statistical mechanics to study the influence that intrinsic noise has on the qualitative behavior of the system. We discuss the application of these concepts to the case of semiarid ecosystems. Using a system-size expansion we calculate the power spectrum of the fluctuations in the system. This predicts the existence of noise-induced oscillations.

preprint2012arXiv

Demographic noise can lead to the spontaneous formation of species

When a collection of phenotypically diverse organisms compete with each other for limited resources, with competition being strongest amongst the most similar, the population can evolve into tightly localised clusters. This process can be thought of as a simple model of the emergence of species. Past studies have neglected the effects of demographic noise and studied the population on a macroscopic scale, where species formation is found to depend upon the shape of the curve describing the decline of competition strength with phenotypic distance. In the following, we will show how including the effects of demographic noise leads to a radically different conclusion. Two situations are identified: a weak-noise regime in which the population exhibits patterns of fluctuation around the macroscopic description, and a strong-noise regime where species appear spontaneously even in the case that all organisms have equal fitness.

preprint2012arXiv

Jamming and pattern formation in models of segregation

We investigate the Schelling model of social segregation, formulated as an intrinsically non-equilibrium system, in which the agents occupy districts (or patches) rather than sites on a grid. We show that this allows the equations governing the dynamical behaviour of the model to be derived. Analysis of these equations reveals a jamming transition in the regime of low-vacancy density, and inclusion of a spatial dimension in the model leads to a pattern forming instability. Both of these phenomena exhibit unusual characteristics which may be studied through our approach.

preprint2012arXiv

Noise-induced metastability in biochemical networks

Intra-cellular biochemical reactions exhibit a rich dynamical phenomenology which cannot be explained within the framework of mean-field rate equations and additive noise. Here, we show that the presence of metastable states and radically different timescales are general features of a broad class of autocatalyic reaction networks, and moreover, that this fact may be exploited to gain analytical results. The latter point is demonstrated by a treatment of the paradigmatic Togashi-Kaneko reaction, which has resisted theoretical analysis for the last decade.

preprint2012arXiv

Spontaneous genetic clustering in populations of competing organisms

We introduce and analyse an individual-based evolutionary model, in which a population of genetically diverse organisms compete with each other for limited resources. Through theoretical analysis and stochastic simulations, we show that the model exhibits a pattern-forming instability which is highly amplified by the effects of demographic noise, leading to the spontaneous formation of genotypic clusters. This mechanism supports the thesis that stochasticity has a central role in the formation and coherence of species.

preprint2012arXiv

Stochastic pattern formation and spontaneous polarisation: the linear noise approximation and beyond

We review the mathematical formalism underlying the modelling of stochasticity in biological systems. Beginning with a description of the system in terms of its basic constituents, we derive the mesoscopic equations governing the dynamics which generalise the more familiar macroscopic equations. We apply this formalism to the analysis of two specific noise-induced phenomena observed in biologically-inspired models. In the first example, we show how the stochastic amplification of a Turing instability gives rise to spatial and temporal patterns which may be understood within the linear noise approximation. The second example concerns the spontaneous emergence of cell polarity, where we make analytic progress by exploiting a separation of time-scales.

preprint2012arXiv

The mode of host-parasite interaction shapes coevolutionary dynamics and the fate of host cooperation

Antagonistic coevolution between hosts and parasites can have a major impact on host population structures, and hence on the evolution of social traits. Using stochastic modelling techniques in the context of bacteria-virus interactions, we investigate the impact of coevolution across a continuum of host-parasite genetic specificity (specifically, where genotypes have the same infectivity/resistance ranges (matching alleles, MA) to highly variable ranges (gene-for-gene, GFG)) on population genetic structure, and on the social behaviour of the host. We find that host cooperation is more likely to be maintained towards the MA end of the continuum, as the more frequent bottlenecks associated with an MA-like interaction can prevent defector invasion, and can even allow migrant cooperators to invade populations of defectors.

preprint2011arXiv

A unified framework for Schelling's model of segregation

Schelling's model of segregation is one of the first and most influential models in the field of social simulation. There are many variations of the model which have been proposed and simulated over the last forty years, though the present state of the literature on the subject is somewhat fragmented and lacking comprehensive analytical treatments. In this article a unified mathematical framework for Schelling's model and its many variants is developed. This methodology is useful in two regards: firstly, it provides a tool with which to understand the differences observed between models; secondly, phenomena which appear in several model variations may be understood in more depth through analytic studies of simpler versions.

preprint2011arXiv

Stability of growing vesicles

We investigate the stability of growing vesicles using the formalism of nonequilibrium thermodynamics. The vesicles are growing due to the accretion of lipids to the bilayer which forms the vesicle membrane. The thermodynamic description is based on the hydrodynamics of a water{/}lipid mixture together with a model of the vesicle as a discontinuous system in the sense of linear nonequilibrium thermodynamics. This formulation allows the forces and fluxes relevant to the dynamic stability of the vesicle to be identified. The method is used to analyze the stability of a spherical vesicle against arbitrary axisymmetric perturbations. It is found that there are generically two critical radii at which changes of stability occur. In the case where the perturbation takes the form of a single zonal harmonic, only one of these radii is physical and is given by the ratio $2 L_p / L_γ$, where $L_p$ is the hydraulic conductivity and $L_γ$ is the Onsager coefficient related to changes in membrane area due to lipid accretion. The stability of such perturbations is related to the value of $l$ corresponding to the particular zonal harmonic: those with lower $l$ are more unstable than those with higher $l$. Possible extensions of the current work and the need for experimental input are discussed.

preprint2011arXiv

Stochastic waves in a Brusselator model with nonlocal interaction

We show that intrinsic noise can induce spatio-temporal phenomena such as Turing patterns and travelling waves in a Brusselator model with nonlocal interaction terms. In order to predict and to characterize these quasi-waves we analyze the nonlocal model using a system-size expansion. The resulting theory is used to calculate the power spectra of the quasi-waves analytically, and the outcome is tested successfully against simulations. We discuss the possibility that nonlocal models in other areas, such as epidemic spread or social dynamics, may contain similar stochastically-induced patterns.

preprint2011arXiv

WKB calculation of an epidemic outbreak distribution

We calculate both the exponential and pre-factor contributions in a WKB approximation of the master equation for a stochastic SIR model with highly oscillatory dynamics. Fixing the basic parameters of the model we investigate how the outbreak distribution changes with the population size. We show that this distribution rapidly becomes highly non-Gaussian, acquiring large tails indicating the presence of rare, but large outbreaks, as the population is made smaller. The analytic results are found to be in excellent agreement with simulations until the systems become so small that the dynamics are dominated by fade-out of the disease.

preprint2010arXiv

A spatial model of autocatalytic reactions

Biological cells with all of their surface structure and complex interior stripped away are essentially vesicles - membranes composed of lipid bilayers which form closed sacs. Vesicles are thought to be relevant as models of primitive protocells, and they could have provided the ideal environment for pre-biotic reactions to occur. In this paper, we investigate the stochastic dynamics of a set of autocatalytic reactions, within a spatially bounded domain, so as to mimic a primordial cell. The discreteness of the constituents of the autocatalytic reactions gives rise to large sustained oscillations, even when the number of constituents is quite large. These oscillations are spatio-temporal in nature, unlike those found in previous studies, which consisted only of temporal oscillations. We speculate that these oscillations may have a role in seeding membrane instabilities which lead to vesicle division. In this way synchronization could be achieved between protocell growth and the reproduction rate of the constituents (the protogenetic material) in simple protocells.

preprint2010arXiv

Diffusion in a crowded environment

We analyze a pair of diffusion equations which are derived in the infinite system--size limit from a microscopic, individual--based, stochastic model. Deviations from the conventional Fickian picture are found which ultimately relate to the depletion of resources on which the particles rely. The macroscopic equations are studied both analytically and numerically, and are shown to yield anomalous diffusion which does not follow a power law with time, as is frequently assumed when fitting data for such phenomena. These anomalies are here understood within a consistent dynamical picture which applies to a wide range of physical and biological systems, underlining the need for clearly defined mechanisms which are systematically analyzed to give definite predictions.

preprint2010arXiv

Dynamical description of vesicle growth and shape change

We systematize and extend the description of vesicle growth and shape change using linear nonequilibrium thermodynamics. By restricting the study to shape changes from spheres to axisymmetric ellipsoids, we are able to give a consistent formulation which includes the lateral tension of the vesicle membrane. This allows us to generalize and correct a previous calculation. Our present calculations suggest that, for small growing vesicles, a prolate ellipsoidal shape should be favored over oblate ellipsoids, whereas for large growing vesicles oblates should be favored over prolates. The validity of this prediction is examined in the light of the various assumptions made in its derivation.

preprint2010arXiv

Evolutionary dynamics, intrinsic noise and cycles of co-operation

We use analytical techniques based on an expansion in the inverse system size to study the stochastic evolutionary dynamics of finite populations of players interacting in a repeated prisoner's dilemma game. We show that a mechanism of amplification of demographic noise can give rise to coherent oscillations in parameter regimes where deterministic descriptions converge to fixed points with complex eigenvalues. These quasi-cycles between co-operation and defection have previously been observed in computer simulations; here we provide a systematic and comprehensive analytical characterization of their properties. We are able to predict their power spectra as a function of the mutation rate and other model parameters, and to compare the relative magnitude of the cycles induced by different types of underlying microscopic dynamics. We also extend our analysis to the iterated prisoner's dilemma game with a win-stay lose-shift strategy, appropriate in situations where players are subject to errors of the trembling-hand type.

preprint2009arXiv

Limit cycles, complex Floquet multipliers and intrinsic noise

We study the effects of intrinsic noise on chemical reaction systems, which in the deterministic limit approach a limit cycle in an oscillatory manner. Previous studies of systems with an oscillatory approach to a fixed point have shown that the noise can transform the oscillatory decay into sustained coherent oscillations with a large amplitude. We show that a similar effect occurs when the stable attractors are limit cycles. We compute the correlation functions and spectral properties of the fluctuations in suitably co-moving Frenet frames for several model systems including driven and coupled Brusselators, and the Willamowski-Roessler system. Analytical results are confirmed convincingly in numerical simulations. The effect is quite general, and occurs whenever the Floquet multipliers governing the stability of the limit cycle are complex, with the amplitude of the oscillations increasing as the instability boundary is approached.