skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Chaotic heteroclinic networks as models of switching behavior in biological systems
Key features of biological activity can often be captured by transitions between a finite number of semi-stable states that correspond to behaviors or decisions. We present here a broad class of dynamical systems that are ideal for modeling such activity. The models we propose are chaotic heteroclinic networks with nontrivial intersections of stable and unstable manifolds. Due to the sensitive dependence on initial conditions, transitions between states are seemingly random. Dwell times, exit distributions, and other transition statistics can be built into the model through geometric design and can be controlled by tunable parameters. To test our model’s ability to simulate realistic biological phenomena, we turned to one of the most studied organisms, C. elegans, well known for its limited behavioral states. We reconstructed experimental data from two laboratories, demonstrating the model’s ability to quantitatively reproduce dwell times and transition statistics under a variety of conditions. Stochastic switching between dominant states in complex dynamical systems has been extensively studied and is often modeled as Markov chains. As an alternative, we propose here a new paradigm, namely, chaotic heteroclinic networks generated by deterministic rules (without the necessity for noise). Chaotic heteroclinic networks can be used to model systems with arbitrary architecture and size without a commensurate increase in phase dimension. They are highly flexible and able to capture a wide range of transition characteristics that can be adjusted through control parameters.  more » « less
Award ID(s):
1901009
PAR ID:
10635835
Author(s) / Creator(s):
;
Publisher / Repository:
AIP Publishing
Date Published:
Journal Name:
Chaos: An Interdisciplinary Journal of Nonlinear Science
Volume:
32
Issue:
12
ISSN:
1054-1500
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Reversible to irreversible (R-IR) transitions arise in numerous periodically driven collectively interacting systems that, after a certain number of driving cycles, organize into a reversible state where the particle trajectories repeat during every or every few cycles. On the irreversible side of the transition, the motion is chaotic. R-IR transitions were first systematically studied for periodically sheared dilute colloids, and have now been found in a wide variety of both soft and hard matter periodically driven systems, including amorphous solids, crystals, vortices in type-II superconductors, and magnetic textures. It has been shown that in several of these systems, the transition to a reversible state is an absorbing phase transition with a critical divergence in the organization timescale at the transition. The same systems are capable of storing multiple memories and may exhibit return point memory. We give an overview of R-IR transitions including recent advances in the field and discuss how the general framework of R-IR transitions could be applied to a much broader class of nonequilibrium systems in which periodic driving occurs, including not only soft and hard condensed matter systems, but also astrophysics, biological systems, and social systems. In particular, some likely candidate systems are commensurate-incommensurate states, systems exhibiting hysteresis or avalanches, nonequilibrium pattern forming states, and other systems with absorbing phase transitions. Periodic driving could be applied to hard condensed matter systems to see if organization into reversible states occurs for metal-insulator transitions, semiconductors, electron glasses, electron nematics, cold atom systems, or Bose-Einstein condensates. R-IR transitions could also be examined in dynamical systems where synchronization or phase locking occurs. We also discuss the possibility of using complex periodic driving, such as changing drive directions or using multiple frequencies, to determine whether these systems can still organize to reversible states or retain complex multiple memories. Finally, we describe features of classical and quantum time crystals that could suggest the occurrence of R-IR transitions in these systems. 
    more » « less
  2. This work presents a framework for studying temporal networks using zigzag persistence, a tool from the field of Topological Data Analysis (TDA). The resulting approach is general and applicable to a wide variety of time-varying graphs. For example, these graphs may correspond to a system modeled as a network with edges whose weights are functions of time, or they may represent a time series of a complex dynamical system. We use simplicial complexes to represent snapshots of the temporal networks that can then be analyzed using zigzag persistence. We show two applications of our method to dynamic networks: an analysis of commuting trends on multiple temporal scales, e.g., daily and weekly, in the Great Britain transportation network, and the detection of periodic/chaotic transitions due to intermittency in dynamical systems represented by temporal ordinal partition networks. Our findings show that the resulting zero- and one-dimensional zigzag persistence diagrams can detect changes in the networks’ shapes that are missed by traditional connectivity and centrality graph statistics. 
    more » « less
  3. Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials. 
    more » « less
  4. Abstract Similar activity patterns may arise from model neural networks with distinct coupling properties and individual unit dynamics. These similar patterns may, however, respond differently to parameter variations and specifically to tuning of inputs that represent control signals. In this work, we analyze the responses resulting from modulation of a localized input in each of three classes of model neural networks that have been recognized in the literature for their capacity to produce robust three-phase rhythms: coupled fast-slow oscillators, near-heteroclinic oscillators, and threshold-linear networks. Triphasic rhythms, in which each phase consists of a prolonged activation of a corresponding subgroup of neurons followed by a fast transition to another phase, represent a fundamental activity pattern observed across a range of central pattern generators underlying behaviors critical to survival, including respiration, locomotion, and feeding. To perform our analysis, we extend the recently developed local timing response curve (lTRC), which allows us to characterize the timing effects due to perturbations, and we complement our lTRC approach with model-specific dynamical systems analysis. Interestingly, we observe disparate effects of similar perturbations across distinct model classes. Thus, this work provides an analytical framework for studying control of oscillations in nonlinear dynamical systems and may help guide model selection in future efforts to study systems exhibiting triphasic rhythmic activity. 
    more » « less
  5. Dynamic systems which underlie controlled systems are expected to increase in complexity as robots, devices, and connected networks become more intelligent. While classical stable systems converge to a stable point (a sink), another type of stability is to consider a stable path rather than a single point. Such stable paths can be made of saddle points that draw in trajectories from certain regions, and then push the trajectory toward the next saddle point. These chains of saddles are called stable heteroclinic channels (SHCs) and can be used in robotic control to represent time sequences. While we have previously shown that each saddle is visualizable as a trajectory waypoint in phase space, how to increase the fidelity of the trajectory was unclear. In this paper, we hypothesized that the waypoints can be individually modified to locally vary fidelity. Specifically, we expected that increasing the saddle value (ratio of saddle eigenvalues) causes the trajectory to slow to more closely approach a particular saddle. Combined with other parameters that control speed and magnitude, a system expressed with an SHC can be modified locally, point by point, without disrupting the rest of the path, supporting their use in motion primitives. While some combinations can enable a trajectory to better reach into corners, other combinations can rotate, distort, and round the trajectory surrounding the modified saddle. Of the system parameters, the saddle value provides the most predictable tunability across 3 orders of magnitude. 
    more » « less