Source author record

Steven H. Strogatz

Steven H. Strogatz 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

26works
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

26 published item(s)

preprint2022arXiv

Asymptotic absorption-time distributions in extinction-prone Markov processes

We characterize absorption-time distributions for birth-death Markov chains with an absorbing boundary. For "extinction-prone" chains (which drift on average toward the absorbing state) the asymptotic distribution is Gaussian, Gumbel, or belongs to a family of skewed distributions. The latter two cases arise when the dynamics slow down dramatically near the boundary. Several models of evolution, epidemics, and chemical reactions fall into these classes; in each case we establish new results for the absorption-time distribution. Applications to African sleeping sickness are discussed.

preprint2022arXiv

Comparing the efficiency of forward and backward contact tracing

Tracing potentially infected contacts of confirmed cases is important when fighting outbreaks of many infectious diseases. The COVID-19 pandemic has motivated researchers to examine how different contact tracing strategies compare in terms of effectiveness (ability to mitigate infections) and cost efficiency (number of prevented infections per isolation). Two important strategies are so-called forward contact tracing (tracing to whom disease spreads) and backward contact tracing (tracing from whom disease spreads). Recently, Kojaku and colleagues reported that backward contact tracing was ``profoundly more effective'' than forward contact tracing, that contact tracing effectiveness ``hinges on reaching the `source' of infection'', and that contact tracing outperformed case isolation in terms of cost efficiency. Here we show that these conclusions are not true in general. They were based in part on simulations that vastly overestimated the effectiveness and efficiency of contact tracing. Our results show that the efficiency of contact tracing strategies is highly contextual; faced with a disease outbreak, the disease dynamics determine whether tracing infection sources or new cases is more impactful. Our results also demonstrate the importance of simulating disease spread and mitigation measures in parallel rather than sequentially.

preprint2022arXiv

Coupled metronomes on a moving platform with Coulomb friction

Using a combination of theory, experiment, and simulation, we revisit the dynamics of two coupled metronomes on a moving platform. Our experiments show that the platform's motion is damped by a dry friction force of Coulomb type, not the viscous linear friction force that has often been assumed in the past. Prompted by this result, we develop a new mathematical model that builds on previously introduced models, but departs from them in its treatment of the friction on the platform. We analyze the model by a two-timescale analysis and derive the slow-flow equations that determine its long-term dynamics. The derivation of the slow flow is challenging, due to the stick-slip motion of the platform in some parameter regimes. Simulations of the slow flow reveal various kinds of long-term behavior including in-phase and antiphase synchronization of identical metronomes, phase locking and phase drift of non-identical metronomes, and metronome suppression and death. In these latter two states, one or both of the metronomes come to swing at such low amplitude that they no longer engage their escapement mechanisms. We find good agreement between our theory, simulations, and experiments, but stress that our exploration is far from exhaustive. Indeed, much still remains to be learned about the dynamics of coupled metronomes, despite their simplicity and familiarity.

preprint2021arXiv

Modeling the interplay between seasonal flu outcomes and individual vaccination decisions

Seasonal influenza presents an ongoing challenge to public health. The rapid evolution of the flu virus necessitates annual vaccination campaigns, but the decision to get vaccinated or not in a given year is largely voluntary, at least in the United States, and many people decide against it. In early attempts to model these yearly flu vaccine decisions, it was often assumed that individuals behave rationally, and do so with perfect information -- assumptions that allowed the techniques of classical economics and game theory to be applied. However, the usual assumptions are contradicted by the emerging empirical evidence about human decision-making behavior in this context. We develop a simple model of coupled disease spread and vaccination dynamics that instead incorporates experimental observations from social psychology to model annual vaccine decision-making more realistically. We investigate population-level effects of these new decision-making assumptions, with the goal of understanding whether the population can self-organize into a state of herd immunity, and if so, under what conditions. Our model agrees with established results while also revealing more subtle population-level behavior, including biennial oscillations about the herd immunity threshold.

preprint2020arXiv

Dense networks that do not synchronize and sparse ones that do

For any network of identical Kuramoto oscillators with identical positive coupling, there is a critical connectivity above which the system is guaranteed to converge to the in-phase synchronous state, for almost all initial conditions. But the precise value of this critical connectivity remains unknown. In 2018, Ling, Xu, and Bandeira proved that if each oscillator is coupled to at least 79.29 percent of all the others, global synchrony is ensured. In 2019, Lu and Steinerberger improved this bound to 78.89 percent. Here, by focusing on circulant networks, we find clues that the critical connectivity may be exactly 75 percent. Our methods yield a slight improvement on the best known lower bound on the critical connectivity, from $68.18\%$ to $68.28\%$. We also consider the opposite end of the connectivity spectrum, where the networks are sparse rather than dense. In this regime, we ask how few edges one needs to add to a ring of $n$ oscillators to turn it into a globally synchronizing network. We prove a partial result: all the twisted states in a ring of size $n=2^m$ can be destabilized by adding just $\mathcal{O}(n \log_2 n)$ edges. To finish the proof, one also needs to rule out all other candidate attractors. We have done this for $n=8$ with computational algebraic geometry, but the problem remains open for larger $n$. Thus, even for systems as simple as Kuramoto oscillators, much remains to be learned about dense networks that do not globally synchronize and sparse ones that do.

preprint2020arXiv

Descendant distributions for the impact of mutant contagion on networks

Contagion, broadly construed, refers to anything that can spread infectiously from peer to peer. Examples include communicable diseases, rumors, misinformation, ideas, innovations, bank failures, and electrical blackouts. Sometimes, as in the 1918 Spanish flu epidemic, a contagion mutates at some point as it spreads through a network. Here, using a simple susceptible-infected (SI) model of contagion, we explore the downstream impact of a single mutation event. Assuming that this mutation occurs at a random node in the contact network, we calculate the distribution of the number of "descendants," $d$, downstream from the initial "Patient Zero" mutant. We find that the tail of the distribution decays as $d^{-2}$ for complete graphs, random graphs, small-world networks, networks with block-like structure, and other infinite-dimensional networks. This prediction agrees with the observed statistics of memes propagating and mutating on Facebook, and is expected to hold for other effectively infinite-dimensional networks, such as the global human contact network. In a wider context, our approach suggests a possible starting point for a mesoscopic theory of contagion. Such a theory would focus on the paths traced by a spreading contagion, thereby furnishing an intermediate level of description between that of individual nodes and the total infected population. We anticipate that contagion pathways will hold valuable lessons, given their role as the conduits through which single mutations, innovations, or failures can sweep through a network as a whole.

preprint2016arXiv

Correlated disorder in the Kuramoto model: Effects on phase coherence, finite-size scaling, and dynamic fluctuations

We consider a mean-field model of coupled phase oscillators with quenched disorder in the natural frequencies and coupling strengths. A fraction $p$ of oscillators are positively coupled, attracting all others, while the remaining fraction $1-p$ are negatively coupled, repelling all others. The frequencies and couplings are deterministically chosen in a manner which correlates them, thereby correlating the two types of disorder in the model. We first explore the effect of this correlation on the system's phase coherence. We find that there is a a critical width $γ_c$ in the frequency distribution below which the system spontaneously synchronizes. Moreover, this $γ_c$ is independent of $p$. Hence, our model and the traditional Kuramoto model (recovered when $p=1$) have the same critical width $γ_c$. We next explore the critical behavior of the system by examining the finite-size scaling and the dynamic fluctuation of the traditional order parameter. We find that the model belongs to the same universality class as the Kuramoto model with deterministically (not randomly) chosen natural frequencies for the case of $p<1$.

preprint2016arXiv

Dynamics of a population of oscillatory and excitable elements

We analyze a variant of a model proposed by Kuramoto, Shinomoto, and Sakaguchi for a large population of coupled oscillatory and excitable elements. Using the Ott-Antonsen ansatz, we reduce the behavior of the population to a two-dimensional dynamical system with three parameters. We present the stability diagram and calculate several of its bifurcation curves analytically, for both excitatory and inhibitory coupling. Our main result is that when the coupling function is broad, the system can display bistability between steady states of constant high and low activity, whereas when the coupling function is narrow and inhibitory, one of the states in the bistable regime can show persistent pulsations in activity.

preprint2016arXiv

Phase coherence induced by correlated disorder

We consider a mean-field model of coupled phase oscillators with quenched disorder in the coupling strengths and natural frequencies. When these two kinds of disorder are uncorrelated (and when the positive and negative couplings are equal in number and strength), it is known that phase coherence cannot occur and synchronization is absent. Here we explore the effects of correlating the disorder. Specifically, we assume that any given oscillator either attracts or repels all the others, and that the sign of the interaction is deterministically correlated with the given oscillator's natural frequency. For symmetrically correlated disorder with zero mean, we find that the system spontaneously synchronizes, once the width of the frequency distribution falls below a critical value. For asymmetrically correlated disorder, the model displays coherent traveling waves: the complex order parameter becomes nonzero and rotates with constant frequency different from the system's mean natural frequency. Thus, in both cases, correlated disorder can trigger phase coherence.

preprint2015arXiv

Evolutionary game dynamics of controlled and automatic decision-making

We integrate dual-process theories of human cognition with evolutionary game theory to study the evolution of automatic and controlled decision-making processes. We introduce a model where agents who make decisions using either automatic or controlled processing compete with each other for survival. Agents using automatic processing act quickly and so are more likely to acquire resources, but agents using controlled processing are better planners and so make more effective use of the resources they have. Using the replicator equation, we characterize the conditions under which automatic or controlled agents dominate, when coexistence is possible, and when bistability occurs. We then extend the replicator equation to consider feedback between the state of the population and the environment. Under conditions where having a greater proportion of controlled agents either enriches the environment or enhances the competitive advantage of automatic agents, we find that limit cycles can occur, leading to persistent oscillations in the population dynamics. Critically, however, these limit cycles only emerge when feedback occurs on a sufficiently long time scale. Our results shed light on the connection between evolution and human cognition, and demonstrate necessary conditions for the rise and fall of rationality.

preprint2015arXiv

Limit Cycles Sparked by Mutation in the Repeated Prisoner's Dilemma

We explore a replicator-mutator model of the repeated Prisoner's Dilemma involving three strategies: always cooperate (ALLC), always defect (ALLD), and tit-for-tat (TFT). The dynamics resulting from single unidirectional mutations are considered, with detailed results presented for the mutations TFT $\rightarrow$ ALLC and ALLD $\rightarrow$ ALLC. For certain combinations of parameters, given by the mutation rate $μ$ and the complexity cost $c$ of playing tit-for-tat, we find that the population settles into limit cycle oscillations, with the relative abundance of ALLC, ALLD, and TFT cycling periodically. Surprisingly, these oscillations can occur for unidirectional mutations between any two strategies. In each case, the limit cycles are created and destroyed by supercritical Hopf and homoclinic bifurcations, organized by a Bogdanov-Takens bifurcation. Our results suggest that stable oscillations are a robust aspect of a world of ALLC, ALLD, and costly TFT; the existence of cycles does not depend on the details of assumptions of how mutation is implemented.

preprint2015arXiv

Nonlinear Dynamics of the Rock-Paper-Scissors Game with Mutations

We analyze the replicator-mutator equations for the Rock-Paper-Scissors game. Various graph-theoretic patterns of mutation are considered, ranging from a single unidirectional mutation pathway between two of the species, to global bidirectional mutation among all the species. Our main result is that the coexistence state, in which all three species exist in equilibrium, can be destabilized by arbitrarily small mutation rates. After it loses stability, the coexistence state gives birth to a stable limit cycle solution created in a supercritical Hopf bifurcation. This attracting periodic solution exists for all the mutation patterns considered, and persists arbitrarily close to the limit of zero mutation rate and a zero-sum game.

preprint2015arXiv

Simple model for the Darwinian transition in early evolution

It has been hypothesized that in the era just before the last universal common ancestor emerged, life on earth was fundamentally collective. Ancient life forms shared their genetic material freely through massive horizontal gene transfer (HGT). At a certain point, however, life made a transition to the modern era of individuality and vertical descent. Here we present a minimal model for this hypothesized "Darwinian transition." The model suggests that HGT-dominated dynamics may have been intermittently interrupted by selection-driven processes during which genotypes became fitter and decreased their inclination toward HGT. Stochastic switching in the population dynamics with three-point (hypernetwork) interactions may have destabilized the HGT-dominated collective state and led to the emergence of vertical descent and the first well-defined species in early evolution. A nonlinear analysis of a stochastic model dynamics covering key features of evolutionary processes (such as selection, mutation, drift and HGT) supports this view. Our findings thus suggest a viable route from early collective evolution to the start of individuality and vertical Darwinian evolution, enabling the emergence of the first species.

preprint2015arXiv

Synchronization as Aggregation: Cluster Kinetics of Pulse-Coupled Oscillators

We consider models of identical pulse-coupled oscillators with global interactions. Previous work showed that under certain conditions such systems always end up in sync, but did not quantify how small clusters of synchronized oscillators progressively coalesce into larger ones. Using tools from the study of aggregation phenomena, we obtain exact results for the time-dependent distribution of cluster sizes as the system evolves from disorder to synchrony.

preprint2014arXiv

Self-organization in soliton modelocked parametric frequency combs

We show that self-organization occurs in the phase dynamics of soliton modelocking in paramet- ric frequency combs. Reduction of the Lugiato-Lefever equation (LLE) to a simpler set of phase equations reveals that this self-organization arises via mechanisms akin to those in the Kuramoto model for synchronization of coupled oscillators. In addition, our simulations show that the phase equations evolve to a broadband phase-locked state, analogous to the soliton formation process in the LLE. Our simplified equations intuitively explain the origin of the pump phase offset in soliton- modelocked parametric frequency combs. They also predict that the phase of the intracavity field undergoes an anti-symmetrization that precedes phase synchronization, and they clarify the role of chaotic states in soliton formation in parametric combs.

preprint2013arXiv

Phase diagram for the Kuramoto model with van Hemmen interactions

We consider a Kuramoto model of coupled oscillators that includes quenched random interactions of the type used by van Hemmen in his model of spin glasses. The phase diagram is obtained analytically for the case of zero noise and a Lorentzian distribution of the oscillators' natural frequencies. Depending on the size of the attractive and random coupling terms, the system displays four states: complete incoherence, partial synchronization, partial antiphase synchronization, and a mix of antiphase and ordinary synchronization.

preprint2012arXiv

Encouraging moderation: Clues from a simple model of ideological conflict

Some of the most pivotal moments in intellectual history occur when a new ideology sweeps through a society, supplanting an established system of beliefs in a rapid revolution of thought. Yet in many cases the new ideology is as extreme as the old. Why is it then that moderate positions so rarely prevail? Here, in the context of a simple model of opinion spreading, we test seven plausible strategies for deradicalizing a society and find that only one of them significantly expands the moderate subpopulation without risking its extinction in the process.

preprint2011arXiv

Conformists and contrarians in a Kuramoto model with identical natural frequencies

We consider a variant of the Kuramoto model, in which all the oscillators are now assumed to have the same natural frequency, but some of them are negatively coupled to the mean field. These "contrarian" oscillators tend to align in antiphase with the mean field, whereas the positively coupled "conformist" oscillators favor an in-phase relationship. The interplay between these effects can lead to rich dynamics. In addition to a splitting of the population into two diametrically opposed factions, the system can also display traveling waves, complete incoherence, and a blurred version of the two-faction state. Exact solutions for these states and their bifurcations are obtained by means of the Watanabe-Strogatz transformation and the Ott-Antonsen ansatz. Curiously, this system of oscillators with identical frequencies turns out to exhibit more complicated dynamics than its counterpart with heterogeneous natural frequencies.

preprint2010arXiv

Analysis of a continuous-time model of structural balance

It is not uncommon for certain social networks to divide into two opposing camps in response to stress. This happens, for example, in networks of political parties during winner-takes-all elections, in networks of companies competing to establish technical standards, and in networks of nations faced with mounting threats of war. A simple model for these two-sided separations is the dynamical system dX/dt = X^2 where X is a matrix of the friendliness or unfriendliness between pairs of nodes in the network. Previous simulations suggested that only two types of behavior were possible for this system: either all relationships become friendly, or two hostile factions emerge. Here we prove that for generic initial conditions, these are indeed the only possible outcomes. Our analysis yields a closed-form expression for faction membership as a function of the initial conditions, and implies that the initial amount of friendliness in large social networks (started from random initial conditions) determines whether they will end up in intractable conflict or global harmony.

preprint2009arXiv

Phase oscillators with global sinusoidal coupling evolve by Mobius group action

Systems of N identical phase oscillators with global sinusoidal coupling are known to display low-dimensional dynamics. Although this phenomenon was first observed about 20 years ago, its underlying cause has remained a puzzle. Here we expose the structure working behind the scenes of these systems, by proving that the governing equations are generated by the action of the Mobius group, a three-parameter subgroup of fractional linear transformations that map the unit disc to itself. When there are no auxiliary state variables, the group action partitions the N-dimensional state space into three-dimensional invariant manifolds (the group orbits). The N-3 constants of motion associated with this foliation are the N-3 functionally independent cross ratios of the oscillator phases. No further reduction is possible, in general; numerical experiments on models of Josephson junction arrays suggest that the invariant manifolds often contain three-dimensional regions of neutrally stable chaos.

preprint2009arXiv

Solvable Model of Spiral Wave Chimeras

Spiral waves are ubiquitous in two-dimensional systems of chemical or biological oscillators coupled locally by diffusion. At the center of such spirals is a phase singularity, a topological defect where the oscillator amplitude drops to zero. But if the coupling is nonlocal, a new kind of spiral can occur, with a circular core consisting of desynchronized oscillators running at full amplitude. Here we provide the first analytical description of such a spiral wave chimera, and use perturbation theory to calculate its rotation speed and the size of its incoherent core.

preprint2009arXiv

The energy landscape of social balance

We model a close-knit community of friends and enemies as a fully connected network with positive and negative signs on its edges. Theories from social psychology suggest that certain sign patterns are more stable than others. This notion of social "balance" allows us to define an energy landscape for such networks. Its structure is complex: numerical experiments reveal a landscape dimpled with local minima of widely varying energy levels. We derive rigorous bounds on the energies of these local minima and prove that they have a modular structure that can be used to classify them.

preprint2009arXiv

The Structure of Phonological Networks Across Multiple Languages

The network characteristics based on the phonological similarities in the lexicons of several languages were examined. These languages differed widely in their history and linguistic structure, but commonalities in the network characteristics were observed. These networks were also found to be different from other networks studied in the literature. The properties of these networks suggest explanations for various aspects of linguistic processing and hint at deeper organization within human language.

preprint2008arXiv

Solvable model for chimera states of coupled oscillators

Networks of identical, symmetrically coupled oscillators can spontaneously split into synchronized and desynchronized sub-populations. Such chimera states were discovered in 2002, but are not well understood theoretically. Here we obtain the first exact results about the stability, dynamics, and bifurcations of chimera states by analyzing a minimal model consisting of two interacting populations of oscillators. Along with a completely synchronous state, the system displays stable chimeras, breathing chimeras, and saddle-node, Hopf and homoclinic bifurcations of chimeras.

preprint2006arXiv

Chimera States in a Ring of Nonlocally Coupled Oscillators

Arrays of identical limit-cycle oscillators have been used to model a wide variety of pattern-forming systems, such as neural networks, convecting fluids, laser arrays, and coupled biochemical oscillators. These systems are known to exhibit rich collective behavior, from synchrony and traveling waves to spatiotemporal chaos and incoherence. Recently, Kuramoto and his colleagues reported a strange new mode of organization--here called the chimera state--in which coherence and incoherence exist side by side in the same system of oscillators. Such states have never been seen in systems with either local or global coupling; they are apparently peculiar to the intermediate case of nonlocal coupling. Here we give an exact solution for the chimera state, for a one-dimensional ring of phase oscillators coupled nonlocally by a cosine kernel. The analysis reveals that the chimera is born in a continuous bifurcation from a spatially modulated drift state, and dies in a saddle-node collision with an unstable version of itself.

preprint2004arXiv

Chimera States for Coupled Oscillators

Arrays of identical oscillators can display a remarkable spatiotemporal pattern in which phase-locked oscillators coexist with drifting ones. Discovered two years ago, such "chimera states" are believed to be impossible for locally or globally coupled systems; they are peculiar to the intermediate case of nonlocal coupling. Here we present an exact solution for this state, for a ring of phase oscillators coupled by a cosine kernel. We show that the stable chimera state bifurcates from a spatially modulated drift state, and dies in a saddle-node bifurcation with an unstable chimera.