Showing 1–28 of 28 results for author: Laing, C

Activity patterns in ring networks of quadratic integrate-and-fire neurons with synaptic and gap junction coupling

Authors: Oleh E. Omel'chenko, Carlo R. Laing

Abstract: We consider a ring network of quadratic integrate-and-fire neurons with nonlocal synaptic and gap junction coupling. The corresponding neural field model supports solutions such as standing and travelling waves, and also lurching waves. We show that many of these solutions satisfy self-consistency equations which can be used to follow them as parameters are varied. We perform numerical bifurcation… ▽ More We consider a ring network of quadratic integrate-and-fire neurons with nonlocal synaptic and gap junction coupling. The corresponding neural field model supports solutions such as standing and travelling waves, and also lurching waves. We show that many of these solutions satisfy self-consistency equations which can be used to follow them as parameters are varied. We perform numerical bifurcation analysis of the neural field model, concentrating on the effects of varying gap junction coupling strength. Our methods are generally applicable to a wide variety of networks of quadratic integrate-and-fire neurons. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Periodic solutions in next generation neural field models

Authors: Carlo R. Laing, Oleh E. Omel'chenko

Abstract: We consider a next generation neural field model which describes the dynamics of a network of theta neurons on a ring. For some parameters the network supports stable time-periodic solutions. Using the fact that the dynamics at each spatial location are described by a complex-valued Riccati equation we derive a self-consistency equation that such periodic solutions must satisfy. We determine the s… ▽ More We consider a next generation neural field model which describes the dynamics of a network of theta neurons on a ring. For some parameters the network supports stable time-periodic solutions. Using the fact that the dynamics at each spatial location are described by a complex-valued Riccati equation we derive a self-consistency equation that such periodic solutions must satisfy. We determine the stability of these solutions, and present numerical results to illustrate the usefulness of this technique. The generality of this approach is demonstrated through its application to several other systems involving delays, two-population architecture and networks of Winfree oscillators. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: 15 pages, 11 figures

Explosive behaviour in networks of Winfree oscillators

Authors: Shawn Means, Carlo Laing

Abstract: We consider directed networks of Winfree oscillators with power law distributed in- and out-degree distributions. Gaussian and power law distributed intrinsic frequencies are considered, and these frequencies are positively correlated with oscillators' in-degrees. The Ott/Antonsen ansatz is used to derive degree-based mean field equations for the expected dynamics of networks, and these are numeri… ▽ More We consider directed networks of Winfree oscillators with power law distributed in- and out-degree distributions. Gaussian and power law distributed intrinsic frequencies are considered, and these frequencies are positively correlated with oscillators' in-degrees. The Ott/Antonsen ansatz is used to derive degree-based mean field equations for the expected dynamics of networks, and these are numerically analysed. In a variety of cases "explosive" transitions between either two different steady states or between a steady state and a periodic solution are found, and these transitions are explained using bifurcation theory. △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: To appear in Chaos, Solitons and Fractals

Theta neuron subject to delayed feedback: a prototypical model for self-sustained pulsing

Authors: Carlo R. Laing, Bernd Krauskopf

Abstract: We consider a single theta neuron with delayed self-feedback in the form of a Dirac delta function in time. Because the dynamics of a theta neuron on its own can be solved explicitly -- it is either excitable or shows self-pulsations -- we are able to derive algebraic expressions for existence and stability of the periodic solutions that arise in the presence of feedback. These periodic solutions… ▽ More We consider a single theta neuron with delayed self-feedback in the form of a Dirac delta function in time. Because the dynamics of a theta neuron on its own can be solved explicitly -- it is either excitable or shows self-pulsations -- we are able to derive algebraic expressions for existence and stability of the periodic solutions that arise in the presence of feedback. These periodic solutions are characterized by one or more equally spaced pulses per delay interval, and there is an increasing amount of multistability with increasing delay time. We present a complete description of where these self-sustained oscillations can be found in parameter space; in particular, we derive explicit expressions for the loci of their saddle-node bifurcations. We conclude that the theta neuron with delayed self-feedback emerges as a prototypical model: it provides an analytical basis for understanding pulsating dynamics observed in other excitable systems subject to delayed self-coupling. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Chimeras with uniformly distributed heterogeneity: two coupled populations

Authors: Carlo R. Laing

Abstract: Chimeras occur in networks of two coupled populations of oscillators when the oscillators in one population synchronise while those in the other are asynchronous. We consider chimeras of this form in networks of planar oscillators for which one parameter associated with the dynamics of an oscillator is randomly chosen from a uniform distribution. A generalisation of the approach in [C.R. Laing, Ph… ▽ More Chimeras occur in networks of two coupled populations of oscillators when the oscillators in one population synchronise while those in the other are asynchronous. We consider chimeras of this form in networks of planar oscillators for which one parameter associated with the dynamics of an oscillator is randomly chosen from a uniform distribution. A generalisation of the approach in [C.R. Laing, Physical Review E, 100, 042211, 2019], which dealt with identical oscillators, is used to investigate the existence and stability of chimeras for these heterogeneous networks in the limit of an infinite number of oscillators. In all cases, making the oscillators more heterogeneous destroys the stable chimera in a saddle-node bifurcation. The results help us understand the robustness of chimeras in networks of general oscillators to heterogeneity. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: To appear in Phys Rev E

Global and Local Reduced Models for Interacting, Heterogeneous Agents

Authors: Thomas N. Thiem, Felix P. Kemeth, Tom Bertalan, Carlo R. Laing, Ioannis G. Kevrekidis

Abstract: Large collections of coupled, heterogeneous agents can manifest complex dynamical behavior presenting difficulties for simulation and analysis. However, if the collective dynamics lie on a low-dimensional manifold then the original agent-based model may be approximated with a simplified surrogate model on and near the low-dimensional space where the dynamics live. This is typically accomplished by… ▽ More Large collections of coupled, heterogeneous agents can manifest complex dynamical behavior presenting difficulties for simulation and analysis. However, if the collective dynamics lie on a low-dimensional manifold then the original agent-based model may be approximated with a simplified surrogate model on and near the low-dimensional space where the dynamics live. This is typically accomplished by deriving coarse variables that summarize the collective dynamics, these may take the form of either a collection of scalars or continuous fields (e.g. densities), which are then used as part of a reduced model. Analytically identifying such simplified models is challenging and has traditionally been accomplished through the use of mean-field reductions or an Ott-Antonsen ansatz, but is often impossible. Here we present a data-driven coarse-graining methodology for discovering such reduced models. We consider two types of reduced models: globally-based models which use global information and predict dynamics using information from the whole ensemble, and locally-based models that use local information, that is, information from just a subset of agents close (close in heterogeneity space, not physical space) to an agent, to predict the dynamics of an agent. For both approaches we are able to learn laws governing the behavior of the reduced system on the low-dimensional manifold directly from time series of states from the agent-based system. A nontrivial conclusion is that the dynamics can be equally well reproduced by an all-to-all coupled as well as by a locally coupled model of the same agents. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 16 pages, 14 figures

The effects of degree distributions in random networks of Type-I neurons

Authors: Carlo R. Laing

Abstract: We consider large networks of theta neurons and use the Ott/Antonsen ansatz to derive degree-based mean field equations governing the expected dynamics of the networks. Assuming random connectivity we investigate the effects of varying the widths of the in- and out-degree distributions on the dynamics of excitatory or inhibitory synaptically coupled networks, and gap junction coupled networks. For… ▽ More We consider large networks of theta neurons and use the Ott/Antonsen ansatz to derive degree-based mean field equations governing the expected dynamics of the networks. Assuming random connectivity we investigate the effects of varying the widths of the in- and out-degree distributions on the dynamics of excitatory or inhibitory synaptically coupled networks, and gap junction coupled networks. For synaptically coupled networks, the dynamics are independent of the out-degree distribution. Broadening the in-degree distribution destroys oscillations in inhibitory networks and decreases the range of bistability in excitatory networks. For gap junction coupled neurons, broadening the degree distribution varies the values of parameters at which there is an onset of collective oscillations. Many of the results are shown to also occur in networks of more realistic neurons. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: To appear in Physical Review E

Learning emergent PDEs in a learned emergent space

Authors: Felix P. Kemeth, Tom Bertalan, Thomas Thiem, Felix Dietrich, Sung Joon Moon, Carlo R. Laing, Ioannis G. Kevrekidis

Abstract: We extract data-driven, intrinsic spatial coordinates from observations of the dynamics of large systems of coupled heterogeneous agents. These coordinates then serve as an emergent space in which to learn predictive models in the form of partial differential equations (PDEs) for the collective description of the coupled-agent system. They play the role of the independent spatial variables in this… ▽ More We extract data-driven, intrinsic spatial coordinates from observations of the dynamics of large systems of coupled heterogeneous agents. These coordinates then serve as an emergent space in which to learn predictive models in the form of partial differential equations (PDEs) for the collective description of the coupled-agent system. They play the role of the independent spatial variables in this PDE (as opposed to the dependent, possibly also data-driven, state variables). This leads to an alternative description of the dynamics, local in these emergent coordinates, thus facilitating an alternative modeling path for complex coupled-agent systems. We illustrate this approach on a system where each agent is a limit cycle oscillator (a so-called Stuart-Landau oscillator); the agents are heterogeneous (they each have a different intrinsic frequency $ωおめが$) and are coupled through the ensemble average of their respective variables. After fast initial transients, we show that the collective dynamics on a slow manifold can be approximated through a learned model based on local "spatial" partial derivatives in the emergent coordinates. The model is then used for prediction in time, as well as to capture collective bifurcations when system parameters vary. The proposed approach thus integrates the automatic, data-driven extraction of emergent space coordinates parametrizing the agent dynamics, with machine-learning assisted identification of an "emergent PDE" description of the dynamics in this parametrization. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Emergent spaces for coupled oscillators

Authors: Thomas N. Thiem, Mahdi Kooshkbaghi, Tom Bertalan, Carlo R. Laing, Ioannis G. Kevrekidis

Abstract: In this paper we present a systematic, data-driven approach to discovering "bespoke" coarse variables based on manifold learning algorithms. We illustrate this methodology with the classic Kuramoto phase oscillator model, and demonstrate how our manifold learning technique can successfully identify a coarse variable that is one-to-one with the established Kuramoto order parameter. We then introduc… ▽ More In this paper we present a systematic, data-driven approach to discovering "bespoke" coarse variables based on manifold learning algorithms. We illustrate this methodology with the classic Kuramoto phase oscillator model, and demonstrate how our manifold learning technique can successfully identify a coarse variable that is one-to-one with the established Kuramoto order parameter. We then introduce an extension of our coarse-graining methodology which enables us to learn evolution equations for the discovered coarse variables via an artificial neural network architecture templated on numerical time integrators (initial value solvers). This approach allows us to learn accurate approximations of time derivatives of state variables from sparse flow data, and hence discover useful approximate differential equation descriptions of their dynamic behavior. We demonstrate this capability by learning ODEs that agree with the known analytical expression for the Kuramoto order parameter dynamics at the continuum limit. We then show how this approach can also be used to learn the dynamics of coarse variables discovered through our manifold learning methodology. In both of these examples, we compare the results of our neural network based method to typical finite differences complemented with geometric harmonics. Finally, we present a series of computational examples illustrating how a variation of our manifold learning methodology can be used to discover sets of "effective" parameters, reduced parameter combinations, for multi-parameter models with complex coupling. We conclude with a discussion of possible extensions of this approach, including the possibility of obtaining data-driven effective partial differential equations for coarse-grained neuronal network behavior. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Moving bumps in theta neuron networks

Authors: Carlo R. Laing, Oleh Omel'chenko

Abstract: We consider large networks of theta neurons on a ring, synaptically coupled with an asymmetric kernel. Such networks support stable "bumps" of activity, which move along the ring if the coupling kernel is asymmetric. We investigate the effects of the kernel asymmetry on the existence, stability and speed of these moving bumps using continuum equations formally describing infinite networks. Dependi… ▽ More We consider large networks of theta neurons on a ring, synaptically coupled with an asymmetric kernel. Such networks support stable "bumps" of activity, which move along the ring if the coupling kernel is asymmetric. We investigate the effects of the kernel asymmetry on the existence, stability and speed of these moving bumps using continuum equations formally describing infinite networks. Depending on the level of heterogeneity within the network we find complex sequences of bifurcations as the amount of asymmetry is varied, in strong contrast to the behaviour of a classical neural field model. △ Less

Submitted 1 April, 2020; originally announced April 2020.

Comments: To appear in Chaos

Journal ref: Chaos 30, 043117 (2020)

The effects of within-neuron degree correlations in networks of spiking neurons

Authors: Carlo R. Laing, Christian Blasche

Abstract: We consider the effects of correlations between the in- and out-degrees of individual neurons on the dynamics of a network of neurons. By using theta neurons, we can derive a set of coupled differential equations for the expected dynamics of neurons with the same in-degree. A Gaussian copula is used to introduce correlations between a neuron's in- and out-degree and numerical bifurcation analysis… ▽ More We consider the effects of correlations between the in- and out-degrees of individual neurons on the dynamics of a network of neurons. By using theta neurons, we can derive a set of coupled differential equations for the expected dynamics of neurons with the same in-degree. A Gaussian copula is used to introduce correlations between a neuron's in- and out-degree and numerical bifurcation analysis is used determine the effects of these correlations on the network's dynamics. For excitatory coupling we find that inducing positive correlations has a similar effect to increasing the coupling strength between neurons, while for inhibitory coupling it has the opposite effect. We also determine the propensity of various two- and three-neuron motifs to occur as correlations are varied and give a plausible explanation for the observed changes in dynamics. △ Less

Submitted 1 April, 2020; originally announced April 2020.

Comments: To appear in Biological Cybernetics

Degree assortativity in networks of spiking neurons

Authors: Christian Blasche, Shawn Means, Carlo R. Laing

Abstract: Degree assortativity refers to the increased or decreased probability of connecting two neurons based on their in- or out-degrees, relative to what would be expected by chance. We investigate the effects of such assortativity in a network of theta neurons. The Ott/Antonsen ansatz is used to derive equations for the expected state of each neuron, and these equations are then coarse-grained in degre… ▽ More Degree assortativity refers to the increased or decreased probability of connecting two neurons based on their in- or out-degrees, relative to what would be expected by chance. We investigate the effects of such assortativity in a network of theta neurons. The Ott/Antonsen ansatz is used to derive equations for the expected state of each neuron, and these equations are then coarse-grained in degree space. We generate families of effective connectivity matrices parametrised by assortativity coefficient and use SVD decompositions of these to efficiently perform numerical bifurcation analysis of the coarse-grained equations. We find that of the four possible types of degree assortativity, two have no effect on the networks' dynamics, while the other two can have a significant effect. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: To appear in Journal of Computational Dynamics

Dynamics and stability of chimera states in two coupled populations of oscillators

Authors: Carlo R. Laing

Abstract: We consider networks formed from two populations of identical oscillators, with uniform strength all-to-all coupling within populations, and also between populations, with a different strength. Such systems are known to support chimera states in which oscillators within one population are perfectly synchronised while in the other the oscillators are incoherent, and have a different mean frequency… ▽ More We consider networks formed from two populations of identical oscillators, with uniform strength all-to-all coupling within populations, and also between populations, with a different strength. Such systems are known to support chimera states in which oscillators within one population are perfectly synchronised while in the other the oscillators are incoherent, and have a different mean frequency from those in the synchronous population. Assuming that the oscillators in the incoherent population always lie on a closed smooth curve $\mathcal{C}$, we derive and analyse the dynamics of the shape of $\mathcal{C}$ and the probability density on $\mathcal{C}$, for four different types of oscillators. We put some previously derived results on a rigorous footing, and analyse two new systems. △ Less

Submitted 26 August, 2019; originally announced August 2019.

Understanding the dynamics of biological and neural oscillator networks through exact mean-field reductions: a review

Authors: Christian Bick, Marc Goodfellow, Carlo R. Laing, Erik A. Martens

Abstract: Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well… ▽ More Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well as the microscopic properties of individual units, shape the emerging collective dynamics is critical to find factors that lead to malfunction. However, many biological systems such as the brain consist of a large number of dynamical units. Hence, their analysis has either relied on simplified heuristic models on a coarse scale, or the analysis comes at a huge computational cost. Here we review recently introduced approaches, known as the Ott-Antonsen and Watanabe-Strogatz reductions, allowing one to simplify the analysis by bridging small and large scales. Thus, reduced model equations are obtained that exactly describe the collective dynamics for each subpopulation in the oscillator network via few collective variables only. The resulting equations are next-generation models: Rather than being heuristic, they exactly link microscopic and macroscopic descriptions and therefore accurately capture microscopic properties of the underlying system. At the same time, they are sufficiently simple to analyze without great computational effort. In the last decade, these reduction methods have become instrumental in understanding how network structure and interactions shape the collective dynamics and the emergence of synchrony. We review this progress based on concrete examples and outline possible limitations. Finally, we discuss how linking the reduced models with experimental data can guide the way towards the development of new treatment approaches, for example, for neurological disease. △ Less

Submitted 28 May, 2020; v1 submitted 14 February, 2019; originally announced February 2019.

Journal ref: Journal of Mathematical Neuroscience, 10(1):9, 2020

Coarse-Grained Descriptions of Dynamics for Networks with both Intrinsic and Structural Heterogeneities

Authors: Tom Bertalan, Yan Wu, Carlo Laing, C. William Gear, Ioannis G. Kevrekidis

Abstract: Finding accurate reduced descriptions for large, complex, dynamically evolving networks is a crucial enabler to their simulation, analysis, and, ultimately, design. Here we propose and illustrate a systematic and powerful approach to obtaining good collective coarse-grained observables-- variables successfully summarizing the detailed state of such networks. Finding such variables can naturally le… ▽ More Finding accurate reduced descriptions for large, complex, dynamically evolving networks is a crucial enabler to their simulation, analysis, and, ultimately, design. Here we propose and illustrate a systematic and powerful approach to obtaining good collective coarse-grained observables-- variables successfully summarizing the detailed state of such networks. Finding such variables can naturally lead to successful reduced dynamic models for the networks. The main premise enabling our approach is the assumption that the behavior of a node in the network depends (after a short initial transient) on the node identity: a set of descriptors that quantify the node properties, whether intrinsic (e.g. parameters in the node evolution equations) or structural (imparted to the node by its connectivity in the particular network structure). The approach creates a natural link with modeling and "computational enabling technology" developed in the context of Uncertainty Quantification. In our case, however, we will not focus on ensembles of different realizations of a problem, each with parameters randomly selected from a distribution. We will instead study many coupled heterogeneous units, each characterized by randomly assigned (heterogeneous) parameter value(s). One could then coin the term Heterogeneity Quantification for this approach, which we illustrate through a model dynamic network consisting of coupled oscillators with one intrinsic heterogeneity (oscillator individual frequency) and one structural heterogeneity (oscillator degree in the undirected network). The computational implementation of the approach, its shortcomings and possible extensions are also discussed. △ Less

Submitted 28 April, 2017; originally announced May 2017.

CHSalign: A Web Server That Builds upon Junction-Explorer and RNAJAG for Pairwise Alignment of RNA Secondary Structures with Coaxial Helical Stacking

Authors: Lei Hua, Yang Song, Namhee Kim, Christian Laing, Jason T. L. Wang, Tamar Schlick

Abstract: RNA junctions are important structural elements of RNA molecules. They are formed when three or more helices come together in three-dimensional space. Recent studies have focused on the annotation and prediction of coaxial helical stacking (CHS) motifs within junctions. Here we exploit such predictions to develop an efficient alignment tool to handle RNA secondary structures with CHS motifs. Speci… ▽ More RNA junctions are important structural elements of RNA molecules. They are formed when three or more helices come together in three-dimensional space. Recent studies have focused on the annotation and prediction of coaxial helical stacking (CHS) motifs within junctions. Here we exploit such predictions to develop an efficient alignment tool to handle RNA secondary structures with CHS motifs. Specifically, we build upon our Junction-Explorer software for predicting coaxial stacking and RNAJAG for modelling junction topologies as tree graphs to incorporate constrained tree matching and dynamic programming algorithms into a new method, called CHSalign, for aligning the secondary structures of RNA molecules containing CHS motifs. Thus, CHSalign is intended to be an efficient alignment tool for RNAs containing similar junctions. Experimental results based on thousands of alignments demonstrate that CHSalign can align two RNA secondary structures containing CHS motifs more accurately than other RNA secondary structure alignment tools. CHSalign yields a high score when aligning two RNA secondary structures with similar CHS motifs or helical arrangement patterns, and a low score otherwise. This new method has been implemented in a web server, and the program is also made freely available, at http://bioinformatics.njit.edu/CHSalign/. △ Less

Submitted 5 September, 2016; originally announced September 2016.

Comments: 45 pages, 5 figures

Dimension reduction in heterogeneous neural networks: generalized Polynomial Chaos (gPC) and ANalysis-Of-VAriance (ANOVA)

Authors: Minseok Choi, Tom Bertalan, Carlo R. Laing, Ioannis G. Kevrekidis

Abstract: We propose, and illustrate via a neural network example, two different approaches to coarse-graining large heterogeneous networks. Both approaches are inspired from, and use tools developed in, methods for uncertainty quantification in systems with multiple uncertain parameters - in our case, the parameters are heterogeneously distributed on the network nodes. The approach shows promise in acceler… ▽ More We propose, and illustrate via a neural network example, two different approaches to coarse-graining large heterogeneous networks. Both approaches are inspired from, and use tools developed in, methods for uncertainty quantification in systems with multiple uncertain parameters - in our case, the parameters are heterogeneously distributed on the network nodes. The approach shows promise in accelerating large scale network simulations as well as coarse-grained fixed point, periodic solution and stability analysis. We also demonstrate that the approach can successfully deal with structural as well as intrinsic heterogeneities. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Modeling Heterogeneity in Networks using Uncertainty Quantification Tools

Authors: Karthikeyan Rajendran, Andreas C. Tsoumanis, Constantinos I. Siettos, Carlo R. Laing, Ioannis G. Kevrekidis

Abstract: Using the dynamics of information propagation on a network as our illustrative example, we present and discuss a systematic approach to quantifying heterogeneity and its propagation that borrows established tools from Uncertainty Quantification. The crucial assumption underlying this mathematical and computational "technology transfer" is that the evolving states of the nodes in a network quickly… ▽ More Using the dynamics of information propagation on a network as our illustrative example, we present and discuss a systematic approach to quantifying heterogeneity and its propagation that borrows established tools from Uncertainty Quantification. The crucial assumption underlying this mathematical and computational "technology transfer" is that the evolving states of the nodes in a network quickly become correlated with the corresponding node "identities": features of the nodes imparted by the network structure (e.g. the node degree, the node clustering coefficient). The node dynamics thus depend on heterogeneous (rather than uncertain) parameters, whose distribution over the network results from the network structure. Knowing these distributions allows us to obtain an efficient coarse-grained representation of the network state in terms of the expansion coefficients in suitable orthogonal polynomials. This representation is closely related to mathematical/computational tools for uncertainty quantification (the Polynomial Chaos approach and its associated numerical techniques). The Polynomial Chaos coefficients provide a set of good collective variables for the observation of dynamics on a network, and subsequently, for the implementation of reduced dynamic models of it. We demonstrate this idea by performing coarse-grained computations of the nonlinear dynamics of information propagation on our illustrative network model using the Equation-Free approach △ Less

Submitted 24 November, 2015; originally announced November 2015.

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

Authors: Mark J. Panaggio, Daniel M. Abrams, Peter Ashwin, Carlo R. Laing

Abstract: 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 or… ▽ More 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. △ Less

Submitted 12 August, 2015; originally announced August 2015.

Comments: 13 pages, 16 figures

MSC Class: 34C15; 34C23

Journal ref: Phys. Rev. E 93, 012218 (2016)

Chimeras in networks with purely local coupling

Authors: Carlo R. Laing

Abstract: Chimera states in spatially extended networks of oscillators have some oscillators synchronised while the remainder are asynchronous. These states have primarily been studied in networks with nonlocal coupling, and more recently in networks with global coupling. Here we present three networks with only local coupling (diffusive, to nearest neighbours) which are numerically found to support chimera… ▽ More Chimera states in spatially extended networks of oscillators have some oscillators synchronised while the remainder are asynchronous. These states have primarily been studied in networks with nonlocal coupling, and more recently in networks with global coupling. Here we present three networks with only local coupling (diffusive, to nearest neighbours) which are numerically found to support chimera states. One of the networks is analysed using a self-consistency argument in the continuum limit, and this is used to find the boundaries of existence of a chimera state in parameter space. △ Less

Submitted 18 June, 2015; originally announced June 2015.

Comments: submitted

Conservation of writhe helicity under anti-parallel reconnection

Authors: Christian E. Laing, Renzo L. Ricca, De Witt L. Sumners

Abstract: Reconnection is a fundamental event in many areas of science, from the interaction of vortices in classical and quantum fluids, and magnetic flux tubes in magnetohydrodynamics and plasma physics, to the recombination in polymer physics and DNA biology. By using fundamental results in topological fluid mechanics, the helicity of a flux tube can be calculated in terms of writhe and twist contributio… ▽ More Reconnection is a fundamental event in many areas of science, from the interaction of vortices in classical and quantum fluids, and magnetic flux tubes in magnetohydrodynamics and plasma physics, to the recombination in polymer physics and DNA biology. By using fundamental results in topological fluid mechanics, the helicity of a flux tube can be calculated in terms of writhe and twist contributions. Here we show that the writhe is conserved under anti-parallel reconnection. Hence, for a pair of interacting flux tubes of equal flux, if the twist of the reconnected tube is the sum of the original twists of the interacting tubes, then helicity is conserved during reconnection. Thus, any deviation from helicity conservation is entirely due to the intrinsic twist inserted or deleted locally at the reconnection site. This result has important implications for helicity and energy considerations in various physical contexts. △ Less

Submitted 14 October, 2014; originally announced October 2014.

Comments: 21 pages, 4 figures

Sonification Aesthetics and Listening for Network Situational Awareness

Authors: Paul Vickers, Christopher Laing, Mohamed Debashi, Tom Fairfax

Abstract: This paper looks at the problem of using sonification to enable network administrators to maintaining situational awareness about their network environment. Network environments generate a lot of data and the need for continuous monitoring means that sonification systems must be designed in such a way as to maximise acceptance while minimising annoyance and listener fatigue. It will be argued that… ▽ More This paper looks at the problem of using sonification to enable network administrators to maintaining situational awareness about their network environment. Network environments generate a lot of data and the need for continuous monitoring means that sonification systems must be designed in such a way as to maximise acceptance while minimising annoyance and listener fatigue. It will be argued that solutions based on the concept of the soundscape offer an ecological advantage over other sonification designs. △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: Workshop paper presented at SoniHED --- Conference on Sonification of Health and Environmental Data, York, UK, 12 September, 2014

Sonification of a Network's Self-Organized Criticality

Authors: Paul Vickers, Chris Laing, Tom Fairfax

Abstract: Communication networks involve the transmission and reception of large volumes of data. Research indicates that network traffic volumes will continue to increase. These traffic volumes will be unprecedented and the behaviour of global information infrastructures when dealing with these data volumes is unknown. It has been shown that complex systems (including computer networks) exhibit self-organi… ▽ More Communication networks involve the transmission and reception of large volumes of data. Research indicates that network traffic volumes will continue to increase. These traffic volumes will be unprecedented and the behaviour of global information infrastructures when dealing with these data volumes is unknown. It has been shown that complex systems (including computer networks) exhibit self-organized criticality under certain conditions. Given the possibility in such systems of a sudden and spontaneous system reset the development of techniques to inform system administrators of this behaviour could be beneficial. This article focuses on the combination of two dissimilar research concepts, namely sonification (a form of auditory display) and self-organized criticality (SOC). A system is described that sonifies in real time an information infrastructure's self-organized criticality to alert the network administrators of both normal and abnormal network traffic and operation. △ Less

Submitted 17 July, 2014; originally announced July 2014.

ACM Class: H.5.2; C.2.3; H.5.5

Chimeras in networks of planar oscillators

Authors: Carlo R. Laing

Abstract: Chimera states occur in networks of coupled oscillators, and are characterized by having some fraction of the oscillators perfectly synchronized, while the remainder are desynchronized. Most chimera states have been observed in networks of phase oscillators with coupling via a sinusoidal function of phase differences, and it is only for such networks that any analysis has been performed. Here we p… ▽ More Chimera states occur in networks of coupled oscillators, and are characterized by having some fraction of the oscillators perfectly synchronized, while the remainder are desynchronized. Most chimera states have been observed in networks of phase oscillators with coupling via a sinusoidal function of phase differences, and it is only for such networks that any analysis has been performed. Here we present the first analysis of chimera states in a network of planar oscillators, each of which is described by both an amplitude and a phase. We find that as the attractivity of the underlying periodic orbit is reduced chimeras are destroyed in saddle-node bifurcations, and supercritical Hopf and homoclinic bifurcations of chimeras also occur. △ Less

Submitted 22 June, 2010; originally announced June 2010.

Comments: To appear, Phys. Rev. E

Solvable Model of Spiral Wave Chimeras

Authors: Erik A. Martens, Carlo R. Laing, Steven H. Strogatz

Abstract: 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.… ▽ More 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. △ Less

Submitted 12 December, 2009; v1 submitted 28 October, 2009; originally announced October 2009.

Comments: 4 pages, 4 figures; added reference, figure, further numerical tests

Journal ref: Phys. Rev. E 79, 026204 (2009)

Chimera states in heterogeneous networks

Authors: Carlo R. Laing

Abstract: Chimera states in networks of coupled oscillators occur when some fraction of the oscillators synchronise with one another, while the remaining oscillators are incoherent. Several groups have studied chimerae in networks of identical oscillators, but here we study these states in a heterogeneous model for which the natural frequencies of the oscillators are chosen from a distribution. We obtain… ▽ More Chimera states in networks of coupled oscillators occur when some fraction of the oscillators synchronise with one another, while the remaining oscillators are incoherent. Several groups have studied chimerae in networks of identical oscillators, but here we study these states in a heterogeneous model for which the natural frequencies of the oscillators are chosen from a distribution. We obtain exact results by reduction to a finite set of differential equations. We find that heterogeneity can destroy chimerae, destroy all states except chimerae, or destabilise chimerae in Hopf bifurcations, depending on the form of the heterogeneity. △ Less

Submitted 17 December, 2008; v1 submitted 23 September, 2008; originally announced September 2008.

Comments: Revised text. To appear, Chaos

Coarse-grained dynamics of an activity bump in a neural field model

Authors: C. R. Laing, T. A. Frewen, I. G. Kevrekidis

Abstract: We study a stochastic nonlocal PDE, arising in the context of modelling spatially distributed neural activity, which is capable of sustaining stationary and moving spatially-localized ``activity bumps''. This system is known to undergo a pitchfork bifurcation in bump speed as a parameter (the strength of adaptation) is changed; yet increasing the noise intensity effectively slowed the motion of… ▽ More We study a stochastic nonlocal PDE, arising in the context of modelling spatially distributed neural activity, which is capable of sustaining stationary and moving spatially-localized ``activity bumps''. This system is known to undergo a pitchfork bifurcation in bump speed as a parameter (the strength of adaptation) is changed; yet increasing the noise intensity effectively slowed the motion of the bump. Here we revisit the system from the point of view of describing the high-dimensional stochastic dynamics in terms of the effective dynamics of a single scalar "coarse" variable. We show that such a reduced description in the form of an effective Langevin equation characterized by a double-well potential is quantitatively successful. The effective potential can be extracted using short, appropriately-initialized bursts of direct simulation. We demonstrate this approach in terms of (a) an experience-based "intelligent" choice of the coarse observable and (b) an observable obtained through data-mining direct simulation results, using a diffusion map approach. △ Less

Submitted 20 December, 2006; originally announced December 2006.

Comments: Corrected aknowledgements

Periodically-forced finite networks of heterogeneous coupled oscillators: a low-dimensional approach

Authors: Carlo R. Laing, Ioannis G. Kevrekidis

Abstract: We study a network of 500 coupled modified van der Pol oscillators. The value of a parameter associated with each oscillator is drawn from a normal distribution, giving a heterogeneous network. For strong enough coupling the oscillators all have the same period, and we consider periodic forcing of the network when it is in this state. By exploiting the correlations that quickly develop between t… ▽ More We study a network of 500 coupled modified van der Pol oscillators. The value of a parameter associated with each oscillator is drawn from a normal distribution, giving a heterogeneous network. For strong enough coupling the oscillators all have the same period, and we consider periodic forcing of the network when it is in this state. By exploiting the correlations that quickly develop between the state of an oscillator and the value of its parameter we obtain an approximate low-dimensional description of the system in terms of the first few coefficients in a polynomial chaos expansion. Standard bifurcation analysis can then be performed on this low-dimensional system, and the results obtained from this predict very well the behaviour of the high-dimensional system for any set of realisations of the random parameter. Situations in which the method begins to fails are also discussed. △ Less

Submitted 14 December, 2006; originally announced December 2006.