Showing papers in "arXiv: Dynamical Systems in 2017"

Deep learning for universal linear embeddings of nonlinear dynamics.

Abstract: Identifying coordinate transformations that make strongly nonlinear dynamics approximately linear is a central challenge in modern dynamical systems. These transformations have the potential to enable prediction, estimation, and control of nonlinear systems using standard linear theory. The Koopman operator has emerged as a leading data-driven embedding, as eigenfunctions of this operator provide intrinsic coordinates that globally linearize the dynamics. However, identifying and representing these eigenfunctions has proven to be mathematically and computationally challenging. This work leverages the power of deep learning to discover representations of Koopman eigenfunctions from trajectory data of dynamical systems. Our network is parsimonious and interpretable by construction, embedding the dynamics on a low-dimensional manifold that is of the intrinsic rank of the dynamics and parameterized by the Koopman eigenfunctions. In particular, we identify nonlinear coordinates on which the dynamics are globally linear using a modified auto-encoder. We also generalize Koopman representations to include a ubiquitous class of systems that exhibit continuous spectra, ranging from the simple pendulum to nonlinear optics and broadband turbulence. Our framework parametrizes the continuous frequency using an auxiliary network, enabling a compact and efficient embedding at the intrinsic rank, while connecting our models to half a century of asymptotics. In this way, we benefit from the power and generality of deep learning, while retaining the physical interpretability of Koopman embeddings.

Extended dynamic mode decomposition with dictionary learning: a data-driven adaptive spectral decomposition of the Koopman operator

Abstract: Numerical approximation methods for the Koopman operator have advanced considerably in the last few years. In particular, data-driven approaches such as dynamic mode decomposition (DMD) and its generalization, the extended-DMD (EDMD), are becoming increasingly popular in practical applications. The EDMD improves upon the classical DMD by the inclusion of a flexible choice of dictionary of observables that spans a finite dimensional subspace on which the Koopman operator can be approximated. This enhances the accuracy of the solution reconstruction and broadens the applicability of the Koopman formalism. Although the convergence of the EDMD has been established, applying the method in practice requires a careful choice of the observables to improve convergence with just a finite number of terms. This is especially difficult for high dimensional and highly nonlinear systems. In this paper, we employ ideas from machine learning to improve upon the EDMD method. We develop an iterative approximation algorithm which couples the EDMD with a trainable dictionary represented by an artificial neural network. Using the Duffing oscillator and the Kuramoto Sivashinsky PDE as examples, we show that our algorithm can effectively and efficiently adapt the trainable dictionary to the problem at hand to achieve good reconstruction accuracy without the need to choose a fixed dictionary a priori. Furthermore, to obtain a given accuracy we require fewer dictionary terms than EDMD with fixed dictionaries. This alleviates an important shortcoming of the EDMD algorithm and enhances the applicability of the Koopman framework to practical problems.

Linearly-Recurrent Autoencoder Networks for Learning Dynamics

Abstract: This paper describes a method for learning low-dimensional approximations of nonlinear dynamical systems, based on neural-network approximations of the underlying Koopman operator. Extended Dynamic Mode Decomposition (EDMD) provides a useful data-driven approximation of the Koopman operator for analyzing dynamical systems. This paper addresses a fundamental problem associated with EDMD: a trade-off between representational capacity of the dictionary and over-fitting due to insufficient data. A new neural network architecture combining an autoencoder with linear recurrent dynamics in the encoded state is used to learn a low-dimensional and highly informative Koopman-invariant subspace of observables. A method is also presented for balanced model reduction of over-specified EDMD systems in feature space. Nonlinear reconstruction using partially linear multi-kernel regression aims to improve reconstruction accuracy from the low-dimensional state when the data has complex but intrinsically low-dimensional structure. The techniques demonstrate the ability to identify Koopman eigenfunctions of the unforced Duffing equation, create accurate low-dimensional models of an unstable cylinder wake flow, and make short-time predictions of the chaotic Kuramoto-Sivashinsky equation.

Data-driven model reduction and transfer operator approximation

Abstract: In this review paper, we will present different data-driven dimension reduction techniques for dynamical systems that are based on transfer operator theory as well as methods to approximate transfer operators and their eigenvalues, eigenfunctions, and eigenmodes. The goal is to point out similarities and differences between methods developed independently by the dynamical systems, fluid dynamics, and molecular dynamics communities such as time-lagged independent component analysis (TICA), dynamic mode decomposition (DMD), and their respective generalizations. As a result, extensions and best practices developed for one particular method can be carried over to other related methods.

Dimension, comparison, and almost finiteness

Abstract: We develop a dynamical version of some of the theory surrounding the Toms-Winter conjecture for simple separable nuclear C*-algebras and study its connections to the C*-algebra side via the crossed product. We introduce an analogue of hyperfiniteness for free actions of amenable groups on compact spaces and show that it plays the role of Z-stability in the Toms-Winter conjecture in its relation to dynamical comparison, and also that it implies Z-stability of the crossed product. This property, which we call almost finiteness, generalizes Matui's notion of the same name from the zero-dimensional setting. We also introduce a notion of tower dimension as partial analogue of nuclear dimension and study its relation to dynamical comparison and almost finiteness, as well as to the dynamical asymptotic dimension and amenability dimension of Guentner, Willett, and Yu.

Eigendecompositions of Transfer Operators in Reproducing Kernel Hilbert Spaces

Abstract: Transfer operators such as the Perron--Frobenius or Koopman operator play an important role in the global analysis of complex dynamical systems. The eigenfunctions of these operators can be used to detect metastable sets, to project the dynamics onto the dominant slow processes, or to separate superimposed signals. We extend transfer operator theory to reproducing kernel Hilbert spaces and show that these operators are related to Hilbert space representations of conditional distributions, known as conditional mean embeddings in the machine learning community. Moreover, numerical methods to compute empirical estimates of these embeddings are akin to data-driven methods for the approximation of transfer operators such as extended dynamic mode decomposition and its variants. One main benefit of the presented kernel-based approaches is that these methods can be applied to any domain where a similarity measure given by a kernel is available. We illustrate the results with the aid of guiding examples and highlight potential applications in molecular dynamics as well as video and text data analysis.

Asymptotics of spectral gaps of quasi-periodic Schrödinger operators

Abstract: For non-critical almost Mathieu operators with Diophantine frequency, we establish exponential asymptotics on the size of spectral gaps, and show that the spectrum is homogeneous. We also prove the homogeneity of the spectrum for Sch\"odinger operators with (measure-theoretically) typical quasi-periodic analytic potentials and fixed strong Diophantine frequency. As applications, we show the discrete version of Deift's conjecture \cite{Deift, Deift17} for subcritical analytic quasi-periodic initial data and solve a series of open problems of Damanik-Goldstein et al \cite{BDGL, DGL1, dgsv, Go} and Kotani \cite{Kot97}.

Rate distortion theory, metric mean dimension and measure theoretic entropy

Abstract: We prove a variational principle for the metric mean dimension analog to the one in [LT]. Instead of using the rate distortion function we use the function $h_\mu(\epsilon,T,\delta)$ that is closely related to the entropy $h_\mu(T)$ of $\mu$. Our formulation has the advantage of being, in the authors opinion, more natural when doing computations. As a corollary we obtain a proof of the standard variational principle. We also obtain some relations between the rate distortion function with our function $\tilde{h}_\mu(\epsilon,T,\delta)$, a modification of $h_\mu(\epsilon,T,\delta)$ when replacing the dynamical metrics with the average dynamical metrics. Using our methods we also reprove the main result in [LT]. We will explain how to construct homeomorphisms on closed manifolds with maximal metric mean dimension. We end this paper with some questions that naturally arise from this work.

Interaction network, state space, and control in social dynamics

Abstract: In the present chapter, we study the emergence of global patterns in large groups in first- and second-order multiagent systems, focusing on two ingredients that influence the dynamics: the interaction network and the state space. The state space determines the types of equilibrium that can be reached by the system. Meanwhile, convergence to specific equilibria depends on the connectivity of the interaction network and on the interaction potential. When the system does not satisfy the necessary conditions for convergence to the desired equilibrium, control can be exerted, both on finite-dimensional systems and on their mean-field limit.

A Geometric Approach to Dynamical Model-Order Reduction

Abstract: Any model order reduced dynamical system that evolves a modal decomposition to approximate the discretized solution of a stochastic PDE can be related to a vector field tangent to the manifold of fixed rank matrices The Dynamically Orthogonal (DO) approximation is the canonical reduced order model for which the corresponding vector field is the orthogonal projection of the original system dynamics onto the tangent spaces of this manifold The embedded geometry of the fixed rank matrix manifold is thoroughly analyzed The curvature of the manifold is characterized and related to the smallest singular value through the study of the Weingarten map Differentiability results for the orthogonal projection onto embedded manifolds are reviewed and used to derive an explicit dynamical system for tracking the truncated Singular Value Decomposition (SVD) of a time-dependent matrix It is demonstrated that the error made by the DO approximation remains controlled under the minimal condition that the original solution stays close to the low rank manifold, which translates into an explicit dependence of this error on the gap between singular values The DO approximation is also justified as the dynamical system that applies instantaneously the SVD truncation to optimally constrain the rank of the reduced solution Riemannian matrix optimization is investigated in this extrinsic framework to provide algorithms that adaptively update the best low rank approximation of a smoothly varying matrix The related gradient flow provides a dynamical system that converges to the truncated SVD of an input matrix for almost every initial data

A landing theorem for entire functions with bounded post-singular sets.

Abstract: The Douady-Hubbard landing theorem for periodic external rays is one of the cornerstones of the study of polynomial dynamics. It states that, for a complex polynomial with bounded postcritical set, every periodic external ray lands at a repelling or parabolic periodic point, and conversely every repelling or parabolic point is the landing point of at least one periodic external ray. We prove an analogue of this theorem for an entire function with bounded postsingular set. If the function has finite order of growth, then it is known that the escaping set contains certain curves called "periodic hairs"; we show that every periodic hair lands at a repelling or parabolic periodic point, and conversely every repelling or parabolic periodic point is the landing point of at least one periodic hair. For a postsingularly bounded entire function of infinite order, such hairs may not exist. Therefore we introduce certain dynamically natural connected sets, called "dreadlocks". We show that every periodic dreadlock lands at a repelling or parabolic periodic point, and conversely every repelling or parabolic periodic point is the landing point of at least one periodic dreadlock. More generally, we prove that every point of a hyperbolic set is the landing point of a dreadlock.

Controllability of Conjunctive Boolean Networks with Application to Gene Regulation

Abstract: A Boolean network is a finite state discrete time dynamical system. At each step, each variable takes a value from a binary set. The value update rule for each variable is a local function which depends only on a selected subset of variables. Boolean networks have been used in modeling gene regulatory networks. We focus in this paper on a special class of Boolean networks, namely the conjunctive Boolean networks (CBNs), whose value update rule is comprised of only logic AND operations. It is known that any trajectory of a Boolean network will enter a periodic orbit. Periodic orbits of a CBN have been completely understood. In this paper, we investigate the orbit-controllability and state-controllability of a CBN: We ask the question of how one can steer a CBN to enter any periodic orbit or to reach any final state, from any initial state. We establish necessary and sufficient conditions for a CBN to be orbit-controllable and state-controllable. Furthermore, explicit control laws are presented along the analysis.

Arnold diffusion in arbitrary degrees of freedom and normally hyperbolic invariant cylinders

Abstract: We prove a form of Arnold diffusion in the a priori stable case. Let H0(p) + $\epsilon$H1($\theta$, p, t), $\theta$ $\in$ T n , p $\in$ B n , t $\in$ T = R/T be a nearly integrable system of arbitrary degrees of freedom n 2 with a strictly convex H0. We show that for a "generic" $\epsilon$H1, there exists an orbit ($\theta$, p)(t) satisfying p(t) -- p(0) {\textgreater} l(H1) {\textgreater} 0, where l(H1) is independent of $\epsilon$. The diffusion orbit travels along a co-dimension one resonance , and the only obstruction to our construction is a finite set of additional resonances. For the proof we use a combination geometric and variational methods, and manage to adapt tools which have recently been developed in the a priori unstable case.

Model error in data assimilation

Abstract: This chapter provides various perspective on an important challenge in data assimilation: model error. While the overall goal is to understand the implication of model error of any type in data assimilation, we emphasize on the effect of model error from unresolved scales. In particular, connection to related subjects under different names in applied mathematics, such as the Mori-Zwanzig formalism and the averaging method, were discussed with the hope that the existing methods can be more accessible and eventually be used appropriately. We will classify existing methods into two groups: the statistical methods for those who directly estimate the low-order model error statistics; and the stochastic parameterizations for those who implicitly estimate all statistics by imposing stochastic models beyond the traditional unbiased white noise Gaussian processes. We will provide theory to justify why stochastic parameterization, as one of the main theme in this book, is an adequate tool for mitigating model error in data assimilation. Finally, we will also discuss challenges in lifting this approach in general applications and provide an alternative nonparametric approach.

Koopman Operator Spectrum for Random Dynamical Systems

Abstract: In this paper we consider the Koopman operator associated with the discrete and the continuous time random dynamical system (RDS). We provide results that characterize the spectrum and the eigenfunctions of the stochastic Koopman operator associated with different types of linear RDS. Then we consider the RDS for which the associated Koopman operator family is a semigroup, especially those for which the generator can be determined. We define a stochastic Hankel DMD (sHankel-DMD) algorithm for numerical approximations of the spectral objects (eigenvalues, eigenfunctions) of the stochastic Koopman operator and prove its convergence. We apply the methodology to a variety of examples, revealing objects in spectral expansions of the stochastic Koopman operator and enabling model reduction.

$C^1$ density of stable ergodicity

Abstract: We prove a $C^1$ version of a conjecture by Pugh and Shub: among partially hyperbolic volume-preserving $C^r$ diffeomorphisms, $r>1$, the stably ergodic ones are $C^1$-dense. To establish these results, we develop new perturbation tools for the $C^1$ topology: linearization of horseshoes while preserving entropy, and creation of "superblenders" from hyperbolic sets with large entropy.

Dimension and measures on sub-self-affine sets

Abstract: We show that in a typical sub-self-affine set, the Hausdorff and the Minkowski dimensions coincide and equal the zero of an appropriate topological pressure. This gives a partial positive answer to the question of Falconer. We also study the properties of the topological pressure and the existence and the uniqueness of natural measures supported on a sub-self-affine set.

Unique equilibrium states for geodesic flows in nonpositive curvature

Abstract: We study geodesic flows over compact rank 1 manifolds and prove that sufficiently regular potential functions have unique equilibrium states if the singular set does not carry full pressure. In dimension 2, this proves uniqueness for scalar multiples of the geometric potential on the interval $(-\infty,1)$, which is optimal. In higher dimensions, we obtain the same result on a neighborhood of 0, and give examples where uniqueness holds on all of $\mathbb{R}$. For general potential functions $\varphi$, we prove that the pressure gap holds whenever $\varphi$ is locally constant on a neighborhood of the singular set, which allows us to give examples for which uniqueness holds on a $C^0$-open and dense set of Holder potentials.

LaSalle Invariance Principle for Discrete-time Dynamical Systems: A Concise and Self-contained Tutorial

Abstract: LaSalle invariance principle was originally proposed in the 1950's and has become a fundamental mathematical tool in the area of dynamical systems and control. In both theoretical research and engineering practice, discrete-time dynamical systems have been at least as extensively studied as continuous-time systems. For example, model predictive control is typically studied in discrete-time via Lyapunov methods. However, there is a peculiar absence in the standard literature of standard treatments of Lyapunov functions and LaSalle invariance principle for discrete-time nonlinear systems. Most of the textbooks on nonlinear dynamical systems focus only on continuous-time systems. In Chapter 1 of the book by LaSalle [11], the author establishes the LaSalle invariance principle for difference equation systems. However, all the useful lemmas in [11] are given in the form of exercises with no proof provided. In this document, we provide the proofs of all the lemmas proposed in [11] that are needed to derive the main theorem on the LaSalle invariance principle for discrete-time dynamical systems. We organize all the materials in a self-contained manner. We first introduce some basic concepts and definitions in Section 1, such as dynamical systems, invariant sets, and limit sets. In Section 2 we present and prove some useful lemmas on the properties of invariant sets and limit sets. Finally, we establish the original LaSalle invariance principle for discrete-time dynamical systems and a simple extension in Section~3. In Section 4, we provide some references on extensions of LaSalle invariance principles for further reading. This document is intended for educational and tutorial purposes and contains lemmas that might be useful as a reference for researchers.

Hopf bifurcation in a delayed reaction-diffusion-advection population model

Abstract: In this paper, we investigate a reaction-diffusion-advection model with time delay effect. The stability/instability of the spatially nonhomogeneous positive steady state and the associated Hopf bifurcation are investigated when the given parameter of the model is near the principle eigenvalue of an elliptic operator. Our result implies that time delay can make the spatially nonhomogeneous positive steady state unstable for a reaction-diffusion-advection model, and the model can exhibit oscillatory pattern through Hopf bifurcation.

On three types of dynamics, and the notion of attractor

Abstract: We propose a theoretical framework for an explanation of the numerically discovered phenomenon of the attractor-repeller merger. We identify regimes which are observed in dynamical systems with attractors as defined in a work by Ruelle and show that these attractors can be of three different types. The first two types correspond to th ewell-known types of chaotic behavior - conservative and dissipative, while the attractors of the third type, the reversible cores, provide a new type of chaos, the so-called mixed dynamics, characterized by the inseparability of dissipative and conservative regimes. We prove that every elliptic orbit of a generic non-conservative time-reversible system is a reversible core. We also prove that a generic reversible system with an elliptic orbit is universal, i.e., it displays dynamics of maximum possible richness and complexity.

Explicit Backbone Curves from Spectral Submanifolds of Forced-Damped Nonlinear Mechanical Systems

Abstract: Spectral submanifolds (SSMs) have recently been shown to provide exact and unique reduced-order models for nonlinear unforced mechanical vibrations. Here we extend these results to periodically or quasiperiodically forced mechanical systems, obtaining analytic expressions for forced responses and backbone curves on modal (i.e. two-dimensional) time dependent SSMs. A judicious choice of the parameterization of these SSMs allows us to simplify the reduced dynamics considerably. We demonstrate our analytical formulae on three numerical examples and compare them to results obtained from available normal form methods.

From rate distortion theory to metric mean dimension: variational principle

Abstract: The purpose of this paper is to point out a new connection between information theory and dynamical systems. In the information theory side, we consider rate distortion theory, which studies lossy data compression of stochastic processes under distortion constraints. In the dynamical systems side, we consider mean dimension theory, which studies how many parameters per second we need to describe a dynamical system. The main results are new variational principles connecting rate distortion function to metric mean dimension.

Shifts of finite type as fundamental objects in the theory of shadowing

Abstract: Shifts of finite type and the notion of shadowing, or pseudo-orbit tracing, are powerful tools in the study of dynamical systems. In this paper we prove that there is a deep and fundamental relationship between these two concepts. Let $X$ be a compact totally disconnected space and $f:X\to X$ a continuous map. We demonstrate that $f$ has shadowing if and only if the system $({f},{X})$ is (conjugate to) the inverse limit of a directed system of shifts of finite type. In particular, this implies that, in the case that $X$ is the Cantor set, $f$ has shadowing if and only if $(f,X)$ is the inverse limit of a sequence of shifts of finite type. Moreover, in the general compact metric case, where $X$ is not necessarily totally disconnected, we prove that $f$ has shadowing if and only if $({f},{X})$ is a factor of (i.e. semi-conjugate to) the inverse limit of a sequence of shifts of finite type by a quotient that almost lifts pseudo-orbits.

Data-driven spectral analysis of the Koopman operator.

Abstract: Starting from measured data, we develop a method to compute the fine structure of the spectrum of the Koopman operator with rigorous convergence guarantees. The method is based on the observation that, in the measure-preserving ergodic setting, the moments of the spectral measure associated to a given observable are computable from a single trajectory of this observable. Having finitely many moments available, we use the classical Christoffel-Darboux kernel to separate the atomic and absolutely continuous parts of the spectrum, supported by convergence guarantees as the number of moments tends to infinity. In addition, we propose a technique to detect the singular continuous part of the spectrum as well as two methods to approximate the spectral measure with guaranteed convergence in the weak topology, irrespective of whether the singular continuous part is present or not. The proposed method is simple to implement and readily applicable to large-scale systems since the computational complexity is dominated by inverting an $N\times N$ Hermitian positive-definite Toeplitz matrix, where $N$ is the number of moments, for which efficient and numerically stable algorithms exist; in particular, the complexity of the approach is independent of the dimension of the underlying state-space. We also show how to compute, from measured data, the spectral projection on a given segment of the unit circle, allowing us to obtain a finite-dimensional approximation of the operator that explicitly takes into account the point and continuous parts of the spectrum. Finally, we describe a relationship between the proposed method and the so-called Hankel Dynamic Mode Decomposition, providing new insights into the behavior of the eigenvalues of the Hankel DMD operator. A number of numerical examples illustrate the approach, including a study of the spectrum of the lid-driven two-dimensional cavity flow.

Combining Trajectory Optimization, Supervised Machine Learning, and Model Structure for Mitigating the Curse of Dimensionality in the Control of Bipedal Robots

Abstract: To overcome the obstructions imposed by high-dimensional bipedal models, we embed a stable walking motion in an attractive low-dimensional surface of the system's state space. The process begins with trajectory optimization to design an open-loop periodic walking motion of the high-dimensional model and then adding to this solution, a carefully selected set of additional open-loop trajectories of the model that steer toward the nominal motion. A drawback of trajectories is that they provide little information on how to respond to a disturbance. To address this shortcoming, Supervised Machine Learning is used to extract a low-dimensional state-variable realization of the open-loop trajectories. The periodic orbit is now an attractor of the low-dimensional state-variable model but is not attractive in the full-order system. We then use the special structure of mechanical models associated with bipedal robots to embed the low-dimensional model in the original model in such a manner that the desired walking motions are locally exponentially stable. The design procedure is first developed for ordinary differential equations and illustrated on a simple model. The methods are subsequently extended to a class of hybrid models and then realized experimentally on an Atrias-series 3D bipedal robot.

Measure concentration and the weak Pinsker property

Abstract: Let $(X,\mu)$ be a standard probability space. An automorphism $T$ of $(X,\mu)$ has the weak Pinsker property if for every $\varepsilon > 0$ it has a splitting into a direct product of a Bernoulli shift and an automorphism of entropy less than $\varepsilon$. This property was introduced by Thouvenot, who asked whether it holds for all ergodic automorphisms. This paper proves that it does. The proof actually gives a more general result. Firstly, it gives a relative version: any factor map from one ergodic automorphism to another can be enlarged by arbitrarily little entropy to become relatively Bernoulli. Secondly, using some facts about relative orbit equivalence, the analogous result holds for all free and ergodic measure-preserving actions of a countable amenable group. The key to this work is a new result about measure concentration. Suppose now that $\mu$ is a probability measure on a finite product space $A^n$, and endow this space with its Hamming metric. We prove that $\mu$ may be represented as a mixture of other measures in which (i) most of the weight in the mixture is on measures that exhibit a strong kind of concentration, and (ii) the number of summands is bounded in terms of the difference between the Shannon entropy of $\mu$ and the combined Shannon entropies of its marginals.

Rate-dependent tipping from periodic attractors: partial tipping and connecting orbits

Abstract: We consider how breakdown of the quasistatic approximation for attractors can lead to rate-induced tipping, where a qualitative change in tracking/tipping behaviour of trajectories can be characterised in terms of a critical rate. Associated with rate-induced tipping (where tracking of a branch of quasistatic attractors breaks down) we find a new phenomenon for attractors that are not simply equilibria: partial tipping of the pullback attractor where certain phases of the periodic attractor tip and others track the quasistatic attractor. For a specific model system with a parameter shift between two asymptotically autonomous systems with periodic attractors we characterise thresholds of rate-induced tipping to partial and total tipping. We show these thresholds can be found in terms of certain periodic-to-periodic (PtoP) and periodic-to-equilibrium (PtoE) connections that we determine using Lin's method for an augmented system.

Inheritance of oscillation in chemical reaction networks

Abstract: Some results are presented on how oscillation is inherited by chemical reaction networks (CRNs) when they are built in natural ways from smaller oscillatory networks. The main results describe four important ways in which a CRN can be enlarged while preserving its capacity for oscillation. The results are for general CRNs, not necessarily fully open, but lead to an important corollary for fully open networks: if a fully open CRN R with mass action kinetics admits a nondegenerate (resp., linearly stable) periodic orbit, then so do all such CRNs which include R as an induced subnetwork. This claim holds for other classes of kinetics, but fails, in general, for CRNs which are not fully open. Where analogous results for multistationarity can be proved using the implicit function theorem alone, the results here call on regular and singular perturbation theory. Equipped with these results and with the help of some analysis and numerical simulation, lower bounds are put on the proportion of small fully open CRNs capable of stable oscillation under various assumptions on the kinetics. This exploration suggests that small oscillatory motifs are an important source of oscillation in CRNs.

A proof of Wright's conjecture

Abstract: Wright's conjecture states that the origin is the global attractor for the delay differential equation $y'(t) = - \alpha y(t-1) [ 1 + y(t) ] $ for all $\alpha \in (0,\tfrac{\pi}{2}]$. This has been proven to be true for a subset of parameter values $\alpha$. We extend the result to the full parameter range $\alpha \in (0,\tfrac{\pi}{2}]$, and thus prove Wright's conjecture to be true. Our approach relies on a careful investigation of the neighborhood of the Hopf bifurcation occurring at $\alpha =\tfrac{\pi}{2}$. This analysis fills the gap left by complementary work on Wright's conjecture, which covers parameter values further away from the bifurcation point. Furthermore, we show that the branch of (slowly oscillating) periodic orbits originating from this Hopf bifurcation does not have any subsequent bifurcations (and in particular no folds) for $\alpha\in(\tfrac{\pi}{2} , \tfrac{\pi}{2} + 6.830 \times 10^{-3}]$. When combined with other results, this proves that the branch of slowly oscillating solutions that originates from the Hopf bifurcation at $\alpha=\tfrac{\pi}{2}$ is globally parametrized by $\alpha > \tfrac{\pi}{2}$.