Source author record

Peter Ashwin

Peter Ashwin 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

21works
13topics
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

21 published item(s)

preprint2022arXiv

Climate Response and Sensitivity: Timescales and Late Tipping Points

Climate response metrics are used to quantify the Earth's climate response to anthropogenic changes of atmospheric CO2. Equilibrium Climate Sensitivity (ECS) is one such metric that measures the equilibrium response to CO2 doubling. However, both in their estimation and their usage, such metrics make assumptions on the linearity of climate response, although it is known that, especially for larger forcing levels, response can be nonlinear. Such nonlinear responses may become visible immediately in response to a larger perturbation, or may only become apparent after a long transient. In this paper, we illustrate some potential problems and caveats when estimating ECS from transient simulations. We highlight ways that very slow timescales may lead to poor estimation of ECS even if there is seemingly good fit to linear response over moderate timescales. Moreover, such slow timescale might lead to late abrupt responses ("late tipping points") associated with a system's nonlinearities. We illustrate these ideas using simulations on a global energy balance model with dynamic albedo. We also discuss the implications for estimating ECS for global climate models, highlighting that it is likely to remain difficult to make definitive statements about the simulation times needed to reach an equilibrium.

preprint2021arXiv

Quantification and interpretation of the climate variability record

The spectral view of variability is a compelling and adaptable tool for understanding variability of the climate. In Mitchell (1976) seminal paper, it was used to express, on one graph with log scales, a very wide range of climate variations from millions of years to days. The spectral approach is particularly useful for suggesting causal links between forcing variability and climate response variability. However, a substantial degree of variability is intrinsic and the Earth system may respond to external forcing in a complex manner. There has been an enormous amount of work on understanding climate variability over the last decades. Hence in this paper, we address the question: Can we (after 40 years) update the Mitchell (1976) diagram and provide it with a better interpretation? By reviewing both the extended observations available for such a diagram and new methodological developments in the study of the interaction between internal and forced variability over a wide range of timescales, we give a positive answer to this question. In addition, we review alternative approaches to the spectral decomposition and pose some challenges for a more detailed quantification of climate variability.

preprint2020arXiv

The Echo Index and multistability in input-driven recurrent neural networks

A recurrent neural network (RNN) possesses the echo state property (ESP) if, for a given input sequence, it ``forgets'' any internal states of the driven (nonautonomous) system and asymptotically follows a unique, possibly complex trajectory. The lack of ESP is conventionally understood as a lack of reliable behaviour in RNNs. Here, we show that RNNs can reliably perform computations under a more general principle that accounts only for their local behaviour in phase space. To this end, we formulate a generalisation of the ESP and introduce an echo index to characterise the number of simultaneously stable responses of a driven RNN. We show that it is possible for the echo index to change with inputs, highlighting a potential source of computational errors in RNNs due to characteristics of the inputs driving the dynamics.

preprint2019arXiv

Bifurcation of critical sets and relaxation oscillations in singular fast-slow systems

Fast-slow dynamical systems have subsystems that evolve on vastly different timescales, and bifurcations in such systems can arise due to changes in any or all subsystems. We classify bifurcations of the critical set (the equilibria of the fast subsystem) and associated fast dynamics, parametrized by the slow variables. Using a distinguished parameter approach we are able to classify bifurcations for one fast and one slow variable. Some of these bifurcations are associated with the critical set losing manifold structure. We also conjecture a list of generic bifurcations of the critical set for one fast and two slow variables. We further consider how the bifurcations of the critical set can be associated with generic bifurcations of attracting relaxation oscillations under an appropriate singular notion of equivalence.

preprint2019arXiv

Extreme sensitivity and climate tipping points

A climate state close to a tipping point will have a degenerate linear response to perturbations, which can be associated with extreme values of the equilibrium climate sensitivity (ECS). In this paper we contrast linearized (`instantaneous') with fully nonlinear geometric (`two-point') notions of ECS, in both presence and absence of tipping points. For a stochastic energy balance model of the global mean surface temperature with two stable regimes, we confirm that tipping events cause the appearance of extremes in both notions of ECS. Moreover, multiple regimes with different mean sensitivities are visible in the two-point ECS. We confirm some of our findings in a physics-based multi-box model of the climate system.

preprint2016arXiv

Chaos in generically coupled phase oscillator networks with nonpairwise interactions

The Kuramoto-Sakaguchi system of coupled phase oscillators, where interaction between oscillators is determined by a single harmonic of phase differences of pairs of oscillators, has very simple emergent dynamics in the case of identical oscillators that are globally coupled: there is a variational structure that means the only attractors are full synchrony (in-phase) or splay phase (rotating wave/full asynchrony) oscillations and the bifurcation between these states is highly degenerate. Here we show that nonpairwise coupling - including three and four-way interactions of the oscillator phases - that appears generically at the next order in normal-form based calculations, can give rise to complex emergent dynamics in symmetric phase oscillator networks. In particular, we show that chaos can appear in the smallest possible dimension of four coupled phase oscillators for a range of parameter values.

preprint2016arXiv

Chaotic Weak Chimeras and their Persistence in Coupled Populations of Phase Oscillators

Nontrivial collective behavior may emerge from the interactive dynamics of many oscillatory units. Chimera states are chaotic patterns of spatially localized coherent and incoherent oscillations. The recently-introduced notion of a weak chimera gives a rigorously testable characterization of chimera states for finite-dimensional phase oscillator networks. In this paper we give some persistence results for dynamically invariant sets under perturbations and apply them to coupled populations of phase oscillators with generalized coupling. In contrast to the weak chimeras with nonpositive maximal Lyapunov exponents constructed so far, we show that weak chimeras that are chaotic can exist in the limit of vanishing coupling between coupled populations of phase oscillators. We present numerical evidence that positive Lyapunov exponents can persist for a positive measure set of this inter-population coupling strength.

preprint2016arXiv

Hopf normal form with $S_N$ symmetry and reduction to systems of nonlinearly coupled phase oscillators

Coupled oscillator models where $N$ oscillators are identical and symmetrically coupled to all others with full permutation symmetry $S_N$ are found in a variety of applications. Much, but not all, work on phase descriptions of such systems consider the special case of pairwise coupling between oscillators. In this paper, we show this is restrictive - and we characterise generic multi-way interactions between oscillators that are typically present, except at the very lowest order near a Hopf bifurcation where the oscillations emerge. We examine a network of identical weakly coupled dynamical systems that are close to a supercritical Hopf bifurcation by considering two parameters, $ε$ (the strength of coupling) and $λ$ (an unfolding parameter for the Hopf bifurcation). For small enough $λ>0$ there is an attractor that is the product of $N$ stable limit cycles; this persists as a normally hyperbolic invariant torus for sufficiently small $ε>0$. Using equivariant normal form theory, we derive a generic normal form for a system of coupled phase oscillators with $S_N$ symmetry. For fixed $N$ and taking the limit $0<ε\llλ\ll 1$, we show that the attracting dynamics of the system on the torus can be well approximated by a coupled phase oscillator system that, to lowest order, is the well-known Kuramoto-Sakaguchi system of coupled oscillators. The next order of approximation genericlly includes terms with up to four interacting phases, regardless of $N$. Using a normalization that maintains nontrivial interactions in the limit $N\rightarrow \infty$, we show that the additional terms can lead to new phenomena in terms of coexistence of two-cluster states with the same phase difference but different cluster size.

preprint2016arXiv

Identical phase oscillator networks: bifurcations, symmetry and reversibility for generalized coupling

For a system of coupled identical phase oscillators with full permutation symmetry, any broken symmetries in dynamical behaviour must come from spontaneous symmetry breaking, i.e. from the nonlinear dynamics of the system. The dynamics of phase differences for such a system depends only on the coupling (phase interaction) function $g(φ)$ and the number of oscillators $N$. This paper briefly reviews some results for such systems in the case of general coupling $g$ before exploring two cases in detail: (a) general two harmonic form: $g(φ)=q\sin(φ-α)+r\sin(2φ-β)$ and $N$ small (b) the coupling $g$ is odd or even. We extend previously published bifurcation analyses to the general two harmonic case, and show for even $g$ that the dynamics of phase differences has a number of time-reversal symmetries. For the case of even $g$ with one harmonic it is known the system has $N-2$ constants of the motion. This is true for $N=4$ and any $g$, while for $N=4$ and more than two harmonics in $g$, we show the system must have fewer independent constants of the motion.

preprint2015arXiv

Chimera states in networks of phase oscillators: the case of two small populations

Chimera states are dynamical patterns in networks of coupled oscillators in which regions of synchronous and asynchronous oscillation coexist. Although these states are typically observed in large ensembles of oscillators and analyzed in the continuum limit, chimeras may also occur in systems with finite (and small) numbers of oscillators. Focusing on networks of $2N$ phase oscillators that are organized in two groups, we find that chimera states, corresponding to attracting periodic orbits, appear with as few as two oscillators per group and demonstrate that for $N>2$ the bifurcations that create them are analogous to those observed in the continuum limit. These findings suggest that chimeras, which bear striking similarities to dynamical patterns in nature, are observable and robust in small networks that are relevant to a variety of real-world systems.

preprint2015arXiv

Designing heteroclinic and excitable networks in phase space using two populations of coupled cells

We give a constructive method for realizing an arbitrary directed graph (with no one-cycles) as a heteroclinic or an excitable dynamic network in the phase space of a system of coupled cells of two types. In each case, the system is expressed as a system of first order differential equations. One of the cell types (the $p$-cells) interacts by mutual inhibition and classifies which vertex (state) we are currently close to, while the other cell type (the $y$-cells) excites the $p$-cells selectively and becomes active only when there is a transition between vertices. We exhibit open sets of parameter values such that these dynamical networks exist and demonstrate via numerical simulation that they can be attractors for suitably chosen parameters.

preprint2015arXiv

Mathematical frameworks for oscillatory network dynamics in neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances when this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear - for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

preprint2014arXiv

Multi-cluster dynamics in coupled phase oscillator networks

In this paper we examine robust clustering behaviour with multiple nontrivial clusters for identically and globally coupled phase oscillators. These systems are such that the dynamics is completely determined by the number of oscillators N and a single scalar function $g(φ)$ (the coupling function). Previous work has shown that (a) any clustering can stably appear via choice of a suitable coupling function and (b) open sets of coupling functions can generate heteroclinic network attractors between cluster states of saddle type, though there seem to be no examples where saddles with more than two nontrivial clusters are involved. In this work we clarify the relationship between the coupling function and the dynamics. We focus on cases where the clusters are inequivalent in the sense of not being related by a temporal symmetry, and demonstrate that there are coupling functions that give robust heteroclinic networks between periodic states involving three or more nontrivial clusters. We consider an example for N=6 oscillators where the clustering is into three inequivalent clusters. We also discuss some aspects of the bifurcation structure for periodic multi-cluster states and show that the transverse stability of inequivalent clusters can, to a large extent, be varied independently of the tangential stability.

preprint2014arXiv

On designing heteroclinic networks from graphs

Robust heteroclinic networks are invariant sets that can appear as attractors in symmetrically coupled or otherwise constrained dynamical systems. These networks may have a very complicated structure that is poorly understood and determined to a large extent by the constraints and dimension of the system. As these networks are of great interest as dynamical models of biological and cognitive processes, it is useful to understand how particular graphs can be realised as robust heteroclinic networks that are attracting. This paper presents two methods of realizing arbitrarily complex directed graphs as robust heteroclinic networks for flows generated by ODEs---we say the ODEs {\em realise} the graphs as heteroclinic networks between equilibria that represent the vertices. Suppose we have a directed graph on $n_v$ vertices with $n_e$ edges. The "simplex realisation" embeds the graph as an invariant set of a flow on an $(n_v-1)$-simplex. This method realises the graph as long as it is one- and two-cycle free. The "cylinder realisation" embeds a graph as an invariant set of a flow on a $(n_e+1)$-dimensional space. This method realises the graph as long as it is one-cycle free. In both cases we find the graph as an invariant set within an attractor, and discuss some illustrative examples, including the influence of noise and parameters on the dynamics. In particular we show that the resulting heteroclinic network may or may not display "memory" of the vertices visited.

preprint2014arXiv

Weak chimeras in minimal networks of coupled phase oscillators

We suggest a definition for a type of chimera state that appears in networks of indistinguishable phase oscillators. Defining a "weak chimera" as a type of invariant set showing partial frequency synchronization, we show that this means they cannot appear in phase oscillator networks that are either globally coupled or too small. We exhibit various networks of four, six and ten indistinguishable oscillators where weak chimeras exist with various dynamics and stabilities. We examine the role of Kuramoto-Sakaguchi coupling in giving degenerate (neutrally stable) families of weak chimera states in these example networks.

preprint2013arXiv

Tipping points in open systems: bifurcation, noise-induced and rate-dependent examples in the climate system

Tipping points associated with bifurcations (B-tipping) or induced by noise (N-tipping) are recognized mechanisms that may potentially lead to sudden climate change. We focus here a novel class of tipping points, where a sufficiently rapid change to an input or parameter of a system may cause the system to "tip" or move away from a branch of attractors. Such rate-dependent tipping, or R-tipping, need not be associated with either bifurcations or noise. We present an example of all three types of tipping in a simple global energy balance model of the climate system, illustrating the possibility of dangerous rates of change even in the absence of noise and of bifurcations in the underlying quasi-static system.

preprint2012arXiv

A model for motor-mediated bidirectional transport along an antipolar microtubule bundle

Long-distance bidirectional transport of organelles depends on the motor proteins kinesin and dynein. Using quantitative data obtained from a fungal model system, we previously developed ASEP-models of bidirectional motion of motors along unipolar microtubules (MTs) near the cell ends of the elongated hyphal cells (herein referred as "unipolar section"). However, recent quantitative live cell imaging in this system has demonstrated that long-range motility of motors and their endosomal cargo mainly occurs along extended antipolar microtubule bundles within the central part of the cell (herein referred to as "bipolar section"). Dynein and kinesin-3 motors coordinate their activity to move early endosomes (EEs) in a bidirectional fashion, with dynein mediating retrograde motility along the unipolar section near the cell poles, whereas kinesin-3 is responsible for bidirectional motions along the antipolar section. Here we extend our modelling approach to simulate bidirectional motility along an antipolar microtubule bundle. In our model, cargos (particles) change direction on each MT with a turning rate $Ω$ and the MTs are linked to each other at the minus ends where particles can hop between MTs with a rate $q_1$ (obstacle-induced switching rate) or $q_2$ (end-induced switching rate). By numerical simulations and mean-field approximations, we investigate the distribution of particles along the MTs for different overall densities $Θ$. We find that even if $Θ$ is low, the system can exhibit shocks in the density profiles near plus and minus ends caused by queueing of particles. We also discuss how the switching rates $q_{1,2}$ influence the type of motor that dominates the active transport in the bundle.

preprint2011arXiv

Bidirectional transport and pulsing states in a multi-lane ASEP model

In this paper, we introduce an ASEP-like transport model for bidirectional motion of particles on a multi-lane lattice. The model is motivated by {\em in vivo} experiments on organelle motility along a microtubule (MT), consisting of thirteen protofilaments, where particles are propelled by molecular motors (dynein and kinesin). In the model, organelles (particles) can switch directions of motion due to "tug-of-war" events between counteracting motors. Collisions of particles on the same lane can be cleared by switching to adjacent protofilaments (lane changes). We analyze transport properties of the model with no-flux boundary conditions at one end of a MT ("plus-end" or tip). We show that the ability of lane changes can affect the transport efficiency and the particle-direction change rate obtained from experiments is close to optimal in order to achieve efficient motor and organelle transport in a living cell. In particular, we find a nonlinear scaling of the mean {\em tip size} (the number of particles accumulated at the tip) with injection rate and an associated phase transition leading to {\em pulsing states} characterized by periodic filling and emptying of the system.

preprint2011arXiv

Chaos in Symmetric Phase Oscillator Networks

Phase-coupled oscillators serve as paradigmatic models of networks of weakly interacting oscillatory units in physics and biology. The order parameter which quantifies synchronization was so far found to be chaotic only in systems with inhomogeneities. Here we show that even symmetric systems of identical oscillators may not only exhibit chaotic dynamics, but also chaotically fluctuating order parameters. Our findings imply that neither inhomogeneities nor amplitude variations are necessary to obtain chaos, i.e., nonlinear interactions of phases give rise to the necessary instabilities.

preprint2011arXiv

On local attraction properties and a stability index for heteroclinic connections

Some invariant sets may attract a nearby set of initial conditions but nonetheless repel a complementary nearby set of initial conditions. For a given invariant set $X\subset\R^n$ with a basin of attraction $N$, we define a stability index $σ(x)$ of a point $x\in X$ that characterizes the local extent of the basin. Let $B_ε$ denote a ball of radius $ε$ about $x$. If $σ(x)>0$, then the measure of $B_ε\setminus N$ relative the measure of the ball is $O(ε^{|σ(x)|})$, while if $σ(x)<0$, then the measure of $B_ε\cap N$ relative the measure of the ball is of the same order. We show that this index is constant along trajectories, and relate this orbit invariant to other notions of stability such as Milnor attraction, essential asymptotic stability and asymptotic stability relative to a positive measure set. We adapt the definition to local basins of attraction (i.e. where $N$ is defined as the set of initial conditions that are in the basin and whose trajectories remain local to $X$). This stability index is particularly useful for discussing the stability of robust heteroclinic cycles, where several authors have studied the appearance of cusps of instability near cycles that are Milnor attractors. We study simple (robust heteroclinic) cycles in $\R^4$ and show that the local stability indices (and hence local stability properties) can be calculated in terms of the eigenvalues of the linearization of the vector field at steady states on the cycle. In doing this, we extend previous results of Krupa and Melbourne (1995,2004) and give criteria for simple heteroclinic cycles in $\R^4$ to be Milnor attractors.

preprint2009arXiv

Within-burst synchrony changes for coupled elliptic bursters

We study the appearance of a novel phenomenon for linearly coupled identical bursters: synchronized bursts where there are changes of spike synchrony within each burst. The examples we study are for normal form elliptic bursters where there is a periodic slow passage through a Bautin (codimension two degenerate Andronov-Hopf) bifurcation. This burster has a subcritical Andronov-Hopf bifurcation at the onset of repetitive spiking while end of burst occurs via a fold limit cycle bifurcation. We study synchronization behavior of two and three Bautin-type elliptic bursters for a linear direct coupling scheme. Burst synchronization is known to be prevalent behavior among such coupled bursters, while spike synchronization is more dependent on the details of the coupling. We note that higher order terms in the normal form that do not affect the behavior of a single burster can be responsible for changes in synchrony pattern; more precisely, we find within-burst synchrony changes associated with a turning point in the spiking frequency.