Search

The main goal of this work is to provide a description of transitions from uniform to non-uniform snychronization in diffusions based on large deviation estimates for finite time Lyapunov exponents. These can be characterized in terms of moment Lyapunov exponents which are principal eigenvalues of the generator of the tilted (Feynman-Kac) semigroup. Using a computer assisted proof, we demonstrate how to determine these eigenvalues and investigate the rate function which is the Legendre-Fenichel transform of the moment Lyapunov function. We apply our results to two case studies: the pitchfork bifurcation and a two-dimensional toy model, also considering the transition to a positive asymptotic Lyapunov exponent.

For overparameterized optimization tasks, such as the ones found in modern machine learning, global minima are generally not unique. In order to understand generalization in these settings, it is vital to study to which minimum an optimization algorithm converges. The possibility of having minima that are unstable under the dynamics imposed by the optimization algorithm limits the potential minima that the algorithm can find. In this paper, we characterize the global minima that are dynamically stable/unstable for both deterministic and stochastic gradient descent (SGD). In particular, we introduce a characteristic Lyapunov exponent which depends on the local dynamics around a global minimum and rigorously prove that the sign of this Lyapunov exponent determines whether SGD can accumulate at the respective global minimum.

We study the synchronization behavior of discrete-time Markov chains on countable state spaces. Representing a Markov chain in terms of a random dynamical system, which describes the collective dynamics of trajectories driven by the same noise, allows for the characterization of synchronization via random attractors. We establish the existence and uniqueness of a random attractor under mild conditions and show that forward and pullback attraction are equivalent in our setting. Additionally, we provide necessary and sufficient conditions for reaching the random attractor, or synchronization respectively, in a time of finite mean. By introducing insulated and synchronizing sets, we structure the state space with respect to the synchronization behavior and characterize the shape of the random attractor.

We study a prototypical non-polynomial decision-making model for which agents in a population potentially alternate between two consumption strategies, one related to the exploitation of an unlimited but considerably expensive resource and the other a comparably cheaper but restricted and slowly renewable source. In particular, we study a model following a Boltzmann-like exploration policy, enhancing the accuracy at which the exchange rates are captured with respect to classical polynomial approaches by considering sigmoidal functions to represent the cost-profit relation in both exploit strategies. Additionally, given the intrinsic timescale separation between the decision-making process and recovery rates of the renewable resource, we use geometric singular perturbation theory to analyze the model. We further use numerical analysis to determine parameter ranges for which the model undergoes bifurcations. These bifurcations, being related to critical states of the system, are relevant to the fast transitions between strategies. Hence, we design controllers to regulate such rapid transitions by taking advantage of the system's criticality., Comment: 40 pages, 20 figures

Coherent sets are time-dependent regions in the physical space of nonautonomous flows that exhibit little mixing with their neighborhoods, robustly under small random perturbations of the flow. They thus characterize the global long-term transport behavior of the system. We propose a framework to extract such time-dependent families of coherent sets for nonautonomous systems with an ergodic driving dynamics and (small) Brownian noise in physical space. Our construction involves the assembly and analysis of an operator on functions over the augmented space of the associated skew product that, for each fixed state of the driving, propagates distributions on the corresponding physical-space fibre according to the dynamics. This time-dependent operator has the structure of a semigroup (it is called the Mather semigroup), and we show that a spectral analysis of its generator allows for a trajectory-free computation of coherent families, simultaneously for all states of the driving. Additionally, for quasi-periodically driven torus flows, we propose a tailored Fourier discretization scheme for this generator and demonstrate our method by means of three examples of two-dimensional flows.

The combination of functional limit theorems with the pathwise analysis of deterministic and stochastic differential equations has proven to be a powerful approach to the analysis of fast-slow systems. In a multivariate setting, this requires rough path ideas, as already suggested in the seminal work [Melbourne-Stuart, Nonlinearity, 24, 2011]. This initiated a program pursued by numerous authors and which has required substantial results on invariance principles (also known as functional central limit theorems) in rough path topologies. We take a unified point of view and provide simple and exact formulas, of the Green-Kubo type, that characterize the relevant Brownian rough path limits and discuss how they naturally apply in different settings., Comment: 23 pages

In this work we study the Brusselator - a prototypical model for chemical oscillations - under the assumption that the bifurcation parameter is of order $O(1/\epsilon)$ for positive $\epsilon\ll 1$. The dynamics of this mathematical model exhibits a time scale separation visible via fast and slow regimes along its unique attracting limit cycle. Noticeably this limit cycle accumulates at infinity as $\epsilon\rightarrow 0$, so that in polar coordinates $(\theta,r)$, and by doing a further change of variable $r\mapsto r^{-1}$, we analyse the dynamics near the line at infinity, corresponding to the set $\{r=0\}$. This object becomes a nonhyperbolic invariant manifold for which we use a desingularising rescaling, in order to study the closeby dynamics. Further use of geometric singular perturbation techniques allows us to give a decomposition of the Brusselator limit cycle in terms of four different fully quantified time scales.

We consider a stochastic version of the so-called Brusselator - a mathematical model for a two-dimensional chemical reaction network - in which one of its parameters is assumed to vary randomly. It has been suggested via numerical explorations that the system exhibits noise-induced synchronization when time goes to infinity. Complementing this perspective, in this work we explore some of its finite-time features from a random dynamical systems perspective. In particular, we focus on the deviations that orbits of neighboring initial conditions exhibit under the influence of the same noise realization. For this, we explore its local instabilities via finite-time Lyapunov exponents. Furthermore, we present the stochastic Brusselator as a fast-slow system in the case that one of the parameters is much larger than the other one. In this framework, an apparent mechanism for generating the stochastic instabilities is revealed, being associated to the transition between the slow and fast regimes.

We consider linear cocycles taking values in $\textup{SL}_d(\mathbb{R})$ driven by homeomorphic transformations of a smooth manifold, in discrete and continuous time. We show that any discrete-time cocycle can be extended to a continuous-time cocycle, while preserving its characteristic properties. We provide a necessary and sufficient condition under which this extension is natural in the sense that the base is extended to an associated suspension flow and that the dimension of the cocycle does not change. Further, we refine our general result for the case of (quasi-)periodic driving. As an example, we present a discrete-time cocycle due to Michael Herman. The Furstenberg--Kesten limits of this cocycle do not exist everywhere and its Oseledets splitting is discontinuous. Our results on the continuous-time extension of discrete-time cocycles allow us to construct a continuous-time cocycle with analogous properties.

Canards are a well-studied phenomenon in fast-slow ordinary differential equations implying the delayed loss of stability after the slow passage through a singularity. Recent studies have shown that the corresponding maps stemming from explicit Runge-Kutta discretizations, in particular the forward Euler scheme, exhibit significant distinctions to the continuous-time behavior: for folds, the delay in loss of stability is typically shortened whereas, for transcritical singularities, it is arbitrarily prolonged. We employ the method of modified equations, which correspond with the fixed discretization schemes up to higher order, to understand and quantify these effects directly from a fast-slow ODE, yielding consistent results with the discrete-time behavior and opening a new perspective on the wide range of (de-)stabilization phenomena along canards.

We study the behavior of a stochastic variant of replicator dynamics in two-agent zero-sum games. We characterize the statistics of such systems by their invariant measures which can be shown to be entirely supported on the boundary of the space of mixed strategies. Depending on the noise strength we can furthermore characterize these invariant measures by finding accumulation of mass at specific parts of the boundary. In particular, regardless of the magnitude of noise, we show that any invariant probability measure is a convex combination of Dirac measures on pure strategy profiles, which correspond to vertices/corners of the agents' simplices. Thus, in the presence of stochastic perturbations, even in the most classic zero-sum settings, such as Matching Pennies, we observe a stark disagreement between the axiomatic prediction of Nash equilibrium and the evolutionary emergent behavior derived by an assumption of stochastically adaptive, learning agents.

We prove the positivity of Lyapunov exponents for the normal form of a Hopf bifurcation, perturbed by additive white noise, under sufficiently strong shear strength. This completes a series of related results for simplified situations which we can exploit by studying suitable limits of the shear and noise parameters. The crucial technical ingredient for making this approach rigorous is a result on the continuity of Lyapunov exponents via Furstenberg Khasminskii formulas., Comment: 35 pages

We study a singularly perturbed fast-slow system of two partial differential equations (PDEs) of reaction-diffusion type on a bounded domain. We assume that the reaction terms in the fast variable contain a fold singularity, whereas the slow variable assumes the role of a dynamic bifurcation parameter, thus extending the classical analysis of a fast-slow dynamic fold bifurcation to an infinite-dimensional setting. Our approach combines a spectral Galerkin discretisation with techniques from Geometric Singular Perturbation Theory (GSPT) which are applied to the resulting high-dimensional systems of ordinary differential equations (ODEs). In particular, we show the existence of invariant manifolds away from the fold singularity, while the dynamics in a neighbourhood of the singularity is described by geometric desingularisation, via the blow-up technique. Finally, we relate the Galerkin manifolds that are obtained after the discretisation to the invariant manifolds which exist in the phase space of the original system of PDEs., Comment: 41 pages with 9 figures. Corresponding author: Thomas Zacharis

This work explores a synchronization-like phenomenon induced by common noise for continuous-time Markov jump processes given by chemical reaction networks. A corresponding random dynamical system is formulated in a two-step procedure, at first for the states of the embedded discrete-time Markov chain and then for the augmented Markov chain including also random jump times. We uncover a time-shifted synchronization in the sense that -- after some initial waiting time -- one trajectory exactly replicates another one with a certain time delay. Whether or not such a synchronization behaviour occurs depends on the combination of the initial states. We prove this partial time-shifted synchronization for the special setting of a birth-death process by analyzing the corresponding two-point motion of the embedded Markov chain and determine the structure of the associated random attractor. In this context, we also provide general results on existence and form of random attractors for discrete-time, discrete-space random dynamical systems.

We establish the existence of a full spectrum of Lyapunov exponents for memoryless random dynamical systems with absorption. To this end, we crucially embed the process conditioned to never being absorbed, the Q-process, into the framework of random dynamical systems, allowing us to study multiplicative ergodic properties. We show that the finite-time Lyapunov exponents converge in conditioned probability and apply our results to iterated function systems and stochastic differential equations., Comment: 36 pages, 0 figures

We study in detail the role of covariant Lyapunov vectors and their respective angles for detecting transitions between metastable states in dynamical systems, as recently discussed in several atmospheric science applications. The underlying models are built from data by the dynamical clustering method, called FEM-BV-VAR, and the Lyapunov vectors are approximated based on these models. We test this data-based numerical approach at the hand of three well-understood example systems with increasing dynamical complexity, identifying crucial properties that allow for a successful application of the method: in particular, it turns out that the method requires a clear multiple time scale structure with fast transitions between slow subsystems which can be dynamically characterized by invariant neutral directions of the linear approximation model.