skip to main content

Title: Non-monotonic transients to synchrony in Kuramoto networks and electrochemical oscillators

We performed numerical simulations with the Kuramoto model and experiments with oscillatory nickel electrodissolution to explore the dynamical features of the transients from random initial conditions to a fully synchronized (one-cluster) state. The numerical simulations revealed that certain networks (e.g., globally coupled or dense Erdős–Rényi random networks) showed relatively simple behavior with monotonic increase of the Kuramoto order parameter from the random initial condition to the fully synchronized state and that the transient times exhibited a unimodal distribution. However, some modular networks with bridge elements were identified which exhibited non-monotonic variation of the order parameter with local maximum and/or minimum. In these networks, the histogram of the transients times became bimodal and the mean transient time scaled well with inverse of the magnitude of the second largest eigenvalue of the network Laplacian matrix. The non-monotonic transients increase the relative standard deviations from about 0.3 to 0.5, i.e., the transient times became more diverse. The non-monotonic transients are related to generation of phase patterns where the modules are synchronized but approximately anti-phase to each other. The predictions of the numerical simulations were demonstrated in a population of coupled oscillatory electrochemical reactions in global, modular, and irregular tree networks. The findings clarify the role of network structure in generation of complex transients that can, for example, play a role in intermittent desynchronization of the circadian clock due to external cues or in deep brain stimulations where long transients are required after a desynchronization stimulus.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Journal of Physics: Complexity
Page Range / eLocation ID:
Article No. 015010
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In networks of nonlinear oscillators, symmetries place hard constraints on the system that can be exploited to predict universal dynamical features and steady states, providing a rare generic organizing principle for far-from-equilibrium systems. However, the robustness of this class of theories to symmetry-disrupting imperfections is untested in free-running (i.e., non-computer-controlled) systems. Here, we develop a model experimental reaction-diffusion network of chemical oscillators to test applications of the theory of dynamical systems with symmeries in the context of self-organizing systems relevant to biology and soft robotics. The network is a ring of four microreactors containing the oscillatory Belousov-Zhabotinsky reaction coupled to nearest neighbors via diffusion. Assuming homogeneity across the oscillators, theory predicts four categories of stable spatiotemporal phase-locked periodic states and four categories of invariant manifolds that guide and structure transitions between phase-locked states. In our experiments, we observed that three of the four phase-locked states were displaced from their idealized positions and, in the ensemble of measurements, appeared as clusters of different shapes and sizes, and that one of the predicted states was absent. We also observed the predicted symmetry-derived synchronous clustered transients that occur when the dynamical trajectories coincide with invariant manifolds. Quantitative agreement between experiment and numerical simulations is found by accounting for the small amount of experimentally determined heterogeneity in intrinsic frequency. We further elucidate how different patterns of heterogeneity impact each attractor differently through a bifurcation analysis. We show that examining bifurcations along invariant manifolds provides a general framework for developing intuition about how chemical-specific dynamics interact with topology in the presence of heterogeneity that can be applied to other oscillators in other topologies. 
    more » « less
  2. Key points

    Cardiac electrophysiology and Ca2+handling change rapidly during the fight‐or‐flight response to meet physiological demands.

    Despite dramatic differences in cardiac electrophysiology, the cardiac fight‐or‐flight response is highly conserved across species.

    In this study, we performed physiological sympathetic nerve stimulation (SNS) while optically mapping cardiac action potentials and intracellular Ca2+transients in innervated mouse and rabbit hearts.

    Despite similar heart rate and Ca2+handling responses between mouse and rabbit hearts, we found notable species differences in spatio‐temporal repolarization dynamics during SNS.

    Species‐specific computational models revealed that these electrophysiological differences allowed for enhanced Ca2+handling (i.e. enhanced inotropy) in each species, suggesting that electrophysiological responses are fine‐tuned across species to produce optimal cardiac fight‐or‐flight responses.


    Sympathetic activation of the heart results in positive chronotropy and inotropy, which together rapidly increase cardiac output. The precise mechanisms that produce the electrophysiological and Ca2+handling changes underlying chronotropic and inotropic responses have been studied in detail in isolated cardiac myocytes. However, few studies have examined the dynamic effects of physiological sympathetic nerve activation on cardiac action potentials (APs) and intracellular Ca2+transients (CaTs) in the intact heart. Here, we performed bilateral sympathetic nerve stimulation (SNS) in fully innervated, Langendorff‐perfused rabbit and mouse hearts. Dual optical mapping with voltage‐ and Ca2+‐sensitive dyes allowed for analysis of spatio‐temporal AP and CaT dynamics. The rabbit heart responded to SNS with a monotonic increase in heart rate (HR), monotonic decreases in AP and CaT duration (APD, CaTD), and a monotonic increase in CaT amplitude. The mouse heart had similar HR and CaT responses; however, a pronounced biphasic APD response occurred, with initial prolongation (50.9 ± 5.1 ms att = 0 svs. 60.6 ± 4.1 ms att = 15 s,P < 0.05) followed by shortening (46.5 ± 9.1 ms att = 60 s,P = NSvs. t = 0). We determined the biphasic APD response in mouse was partly due to dynamic changes in HR during SNS and was exacerbated by β‐adrenergic activation. Simulations with species‐specific cardiac models revealed that transient APD prolongation in mouse allowed for greater and more rapid CaT responses, suggesting more rapid increases in contractility; conversely, the rabbit heart requires APD shortening to produce optimal inotropic responses. Thus, while the cardiac fight‐or‐flight response is highly conserved between species, the underlying mechanisms orchestrating these effects differ significantly.

    more » « less
  3. In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019 ; Lee et al., 2020 ). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution. 
    more » « less
  4. In this article, we devise two related control algorithms to change the degree of synchrony of a population of noise-free, identical, uncoupled neural oscillators using a single control input. The algorithms are based on phase reduction, and use a population-level partial differential equation formulation to change the phase distribution of the neurons as desired. Motivated by the pathological neural synchronization hypothesized to be present in patients suffering from essential and parkinsonian tremor, we take our control objective to be the desynchronization of an initially synchronized neural population. Through numerical simulations, we are able to show that our algorithms work for both Type I and Type II neural populations. To demonstrate the versatility of our control algorithms, we also show that they can be applied to synchronize an initially desynchronized neural population as well. For the systems considered in this paper, the control algorithms can be applied to achieve any desired traveling-wave neural phase distribution, as long as the combination of initial and desired phase distributions is non-degenerate. 
    more » « less
  5. Abstract

    Continuous-time Markov chains are frequently used as stochastic models for chemical reaction networks, especially in the growing field of systems biology. A fundamental problem for these Stochastic Chemical Reaction Networks (SCRNs) is to understand the dependence of the stochastic behavior of these systems on the chemical reaction rate parameters. Towards solving this problem, in this paper we develop theoretical tools called comparison theorems that provide stochastic ordering results for SCRNs. These theorems give sufficient conditions for monotonic dependence on parameters in these network models, which allow us to obtain, under suitable conditions, information about transient and steady-state behavior. These theorems exploit structural properties of SCRNs, beyond those of general continuous-time Markov chains. Furthermore, we derive two theorems to compare stationary distributions and mean first passage times for SCRNs with different parameter values, or with the same parameters and different initial conditions. These tools are developed for SCRNs taking values in a generic (finite or countably infinite) state space and can also be applied for non-mass-action kinetics models. When propensity functions are bounded, our method of proof gives an explicit method for coupling two comparable SCRNs, which can be used to simultaneously simulate their sample paths in a comparable manner. We illustrate our results with applications to models of enzymatic kinetics and epigenetic regulation by chromatin modifications.

    more » « less