BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Christian Kuehn

Christian Kuehn contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.DSmath.PRnlin.CDnlin.PSmath.CAnlin.AOmath.NAmath.APNeurons and CognitionBiological PhysicsPopulations and EvolutionNumerical AnalysisQuantitative Methodscond-mat.stat-mechmath-phmath.GNmath.MPMolecular Networksphysics.data-anphysics.med-phphysics.soc-phSocial and Information Networks

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
38works
0followers
22topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2022arXiv

Community Integration Algorithms (CIAs) for Dynamical Systems on Networks

Dynamics of large-scale network processes underlies crucial phenomena ranging across all sciences. Forward simulation of large network models is often computationally prohibitive. Yet, most networks have intrinsic community structure. We exploit these communities and propose a fast simulation algorithm for network dynamics. In particular, aggregating the inputs a node receives constitutes the limiting factor in numerically simulating large-scale network dynamics. We develop community integration algorithms (CIAs) significantly reducing function-evaluations. We obtain a substantial reduction from polynomial to linear computational complexity. We illustrate our results in multiple applications including classical and higher-order Kuramoto-type systems for synchronisation and Cucker--Smale systems exhibiting flocking behaviour on synthetic as well as real-world networks. Numerical comparison and theoretical analysis confirm the robustness and efficiency of CIAs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Estimating rate-induced tipping via asymptotic series and a Melnikov-like method

The paper deals with the study of rate-induced tipping in asymptotically autonomous scalar ordinary differential equations. We prove that, in such a tipping scenario, a solution which limits at a hyperbolic stable equilibrium of the past limit-problem loses uniform asymptotic stability and coincides with a solution which limits at a hyperbolic unstable equilibrium of the future limit-problem. We use asymptotic series to approximate such pairs of solutions and characterize the occurrence of a rate-induced tipping by using only solutions calculable on finite time intervals. Moreover, we show that a Melnikov-inspired method employing the asymptotic series allows to asymptotically approximate the tipping point.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Geometric blow-up of a dynamic Turing instability in the Swift-Hohenberg equation

We present a rigorous analysis of the slow passage through a Turing bifurcation in the Swift-Hohenberg equation using a novel approach based on geometric blow-up. We show that the formally derived multiple scales ansatz which is known from classical modulation theory can be adapted for use in the fast-slow setting, by reformulating it as a blow-up transformation. This leads to dynamically simpler modulation equations posed in the blown-up space, via a formal procedure which directly extends the established approach to the time-dependent setting. The modulation equations take the form of non-autonomous Ginzburg-Landau equations, which can be analysed within the blow-up. The asymptotics of solutions in weighted Sobelev spaces are given in two different cases: (i) A symmetric case featuring a delayed loss of stability, and (ii) A second case in which the symmetry is broken by a source term. In order to characterise the dynamics of the Swift-Hohenberg equation itself we derive rigorous estimates on the error of the dynamic modulation approximation. These estimates are obtained by bounding weak solutions to an evolution equation for the error which is also posed in the blown-up space. Using the error estimates obtained, we are able to infer the asymptotics of a large class of solutions to the dynamic Swift-Hohenberg equation. We provide rigorous asymptotics for solutions in both cases (i) and (ii). We also prove the existence of the delayed loss of stability in the symmetric case (i), and provide a lower bound for the delay time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Self-adapting infectious dynamics on random networks

Self-adaptive dynamics occurs in many physical systems such as socio-economics, neuroscience, or biophysics. We formalize a self-adaptive modeling approach, where adaptation takes place within a set of strategies based on the history of the state of the system. This leads to piecewise deterministic Markovian dynamics coupled to a non-Markovian adaptive mechanism. We apply this framework to basic epidemic models (SIS, SIR) on random networks. We consider a co-evolutionary dynamical network where node-states change through the epidemics and network topology changes through creation and deletion of edges. For a simple threshold base application of lockdown measures we observe large regions in parameter space with oscillatory behavior. For the SIS epidemic model, we derive analytic expressions for the oscillation period from a pairwise closed model. Furthermore, we show that there is a link to self-organized criticality as the basic reproduction number fluctuates around one. We also study the dependence of our results on the underlying network structure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Warning signs for non-Markovian bifurcations: colour blindness and scaling laws

Warning signs for tipping points (or critical transitions) have been very actively studied. Although the theory has been applied successfully in models and in experiments for many complex systems such as for tipping in climate systems, there are ongoing debates, when warning signs can be extracted from data. In this work, we shed light on this debate by considering different types of underlying noise. Thereby, we significantly advance the general theory of warning signs for nonlinear stochastic dynamics. A key scenario deals with stochastic systems approaching a bifurcation point dynamically upon slow parameter variation. The stochastic fluctuations are generically able to probe the dynamics near a deterministic attractor to reveal critical slowing down. Using scaling laws near bifurcations, one can then anticipate the distance to a bifurcation. Previous warning signs results assume that the noise is Markovian, most often even white. Here, we study warning signs for non-Markovian systems including coloured noise and $α$-regular Volterra processes (of which fractional Brownian motion and the Rosenblatt process are special cases). We prove that early warning scaling laws can disappear completely or drastically change their exponent based upon the parameters controlling the noise process. This provides a clear explanation, why applying standard warning signs results to reduced models of complex systems may not agree with data-driven studies. We demonstrate our results numerically in the context of a box model of the Atlantic Meridional Overturning Circulation (AMOC).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A qualitative mathematical model of the immune response under the effect of stress

In the last decades, the interest to understand the connection between brain and body has grown notably. For example, in psychoneuroimmunology many studies associate stress, arising from many different sources and situations, to changes in the immune system from the medical or immunological point of view as well as from the biochemical one. In this paper we identify important behaviours of this interplay between the immune system and stress from medical studies and seek to represent them qualitatively in a paradigmatic, yet simple, mathematical model. To that end we develop a differential equation model with two equations for infection level and immune system, which integrates the effects of stress as an additional parameter. We are able to reproduce a stable healthy state for little stress, an oscillatory state between healthy and infected states for high stress, and a "burn-out" or stable sick state for extremely high stress. The mechanism between the different dynamics is controlled by two saddle-node in cycle (SNIC) bifurcations. Furthermore, our model is able to capture an induced infection upon dropping from moderate to low stress, and it predicts increasing infection periods upon increasing before eventually reaching a burn-out state.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical continuation for fractional PDEs: sharp teeth and bloated snakes

Partial differential equations (PDEs) involving fractional Laplace operators have been increasingly used to model non-local diffusion processes and are actively investigated using both analytical and numerical approaches. The purpose of this work is to study the effects of the spectral fractional Laplacian on the bifurcation structure of reaction-diffusion systems on bounded domains. In order to do this we use advanced numerical continuation techniques to compute the solution branches. Since current available continuation packages only support systems involving the standard Laplacian, we first extend the pde2path software to treat fractional PDEs. The new capabilities are then applied to the study of the Allen-Cahn equation, the Swift-Hohenberg equation and the Schnakenberg system (in which the standard Laplacian is each replaced by the spectral fractional Laplacian). Our study reveals some common effects, which contributes to a better understanding of fractional diffusion in generic reaction-diffusion systems. In particular, we investigate the changes in snaking bifurcation diagrams and also study the spatial structure of non-trivial steady states upon variation of the order of the fractional Laplacian. Our results show that the fractional order can induce very significant qualitative and quantitative changes in global bifurcation structures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Random attractors via pathwise mild solutions for stochastic parabolic evolution equations

We investigate the longtime behavior of stochastic partial differential equations (SPDEs) with differential operators that depend on time and the underlying probability space. In particular, we consider stochastic parabolic evolution problems in Banach spaces with additive noise and prove the existence of random exponential attractors. These are compact random sets of finite fractal dimension that contain the global random attractor and are attracting at an exponential rate. In order to apply the framework of random dynamical systems, we use the concept of pathwise mild solutions. This approach is essential for our setting since the stochastic evolution equation cannot be transformed into a family of PDEs with random coefficients via the stationary Ornstein-Uhlenbeck process.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stability analysis of multiplayer games on adaptive simplicial complexes

We analyze the influence of multiplayer interactions and network adaptation on the stability of equilibrium points in evolutionary games. We consider the Snowdrift game with both two-player and three-player interactions on simplicial complexes. The state of the system and the topology of the interactions are both adaptive through best response strategies of nodes and rewiring strategies of edges, respectively. We derive a closed set of low-dimensional differential equations using pairwise moment closure, which yields an approximation of the lower moments of the system. We numerically confirm the validity of these moment equations. Moreover, we demonstrate that the stability of the fixed points remains unchanged for the considered adaption process. This stability result indicates that rational best response strategies in games are very difficult to destabilize, even if higher-order interactions are taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Coupled Dynamics on Hypergraphs: Master Stability of Steady States and Synchronization

In the study of dynamical systems on networks/graphs, a key theme is how the network topology influences stability for steady states or synchronized states. Ideally, one would like to derive conditions for stability or instability that instead of microscopic details of the individual nodes/vertices rather make the influence of the network coupling topology visible. The master stability function is an important such tool to achieve this goal. Here we generalize the master stability approach to hypergraphs. A hypergraph coupling structure is important as it allows us to take into account arbitrary higher-order interactions between nodes. As for instance in the theory of coupled map lattices, we study Laplace type interaction structures in detail. Since the spectral theory of Laplacians on hypergraphs is richer than on graphs, we see the possibility of new dynamical phenomena. More generally, our arguments provide a blueprint for how to generalize dynamical structures and results from graphs to hypergraphs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Network Dynamics on Graphops

In this brief note, we report a formal mathematical observation: we are about to breach a major century-old barrier in the analysis of interacting particle systems. More precisely, it is well-known that in well-mixed/homogeneous/all-to-all-coupled systems, one may derive mean-field limit equations such as Vlasov-Fokker-Planck equations (VFPEs). A mesoscopic VFPE describes the probability of finding a single vertex/particle in a certain state, forming a bridge between microscopic statistical physics and macroscopic fluid-type approximations. One major obstacle in this framework is to incorporate complex network structures into limiting equations. In many cases, only heuristic approximations exist, or the limits rely on particular classes of integral operators. In this paper, we notice that there is a much more elegant, and profoundly more general, way available due to recent progress in the theory of graph limits. In particular, we show how one may easily enter complex network dynamics via graphops (graph operators) into VFPEs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Random Attractors for Stochastic Partly Dissipative Systems

We prove the existence of a global random attractor for a certain class of stochastic partly dissipative systems. These systems consist of a partial (PDE) and an ordinary differential equation (ODE), where both equations are coupled and perturbed by additive white noise. The deterministic counterpart of such systems and their long-time behaviour have already been considered but there is no theory that deals with the stochastic version of partly dissipative systems in their full generality. We also provide several examples for the application of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Slow Manifolds for Infinite-Dimensional Evolution Equations

We extend classical finite-dimensional Fenichel theory in two directions to infinite dimensions. Under comparably weak assumptions we show that the solution of an infinite-dimensional fast-slow system is approximated well by the corresponding slow flow. After that we construct a two-parameter family of slow manifolds $S_{ε,ζ}$ under more restrictive assumptions on the linear part of the slow equation. The second parameter $ζ$ does not appear in the finite-dimensional setting and describes a certain splitting of the slow variable space in a fast decaying part and its complement. The finite-dimensional setting is contained as a special case in which $S_{ε,ζ}$ does not depend on $ζ$. Finally, we apply our new techniques to three examples of fast-slow systems of partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Travelling waves for discrete stochastic bistable equations

Many physical, chemical and biological systems have an inherent discrete spatial structure that strongly influences their dynamical behaviour. Similar remarks apply to internal or external noise, as well as to nonlocal coupling. In this paper we study the combined effect of nonlocal spatial discretization and stochastic perturbations on travelling waves in the Nagumo equation, which is a prototypical model for bistable reaction-diffusion partial differential equations (PDEs). We prove that under suitable parameter conditions, various discrete-stochastic variants of the Nagumo equation have solutions, which stay close on long time scales to the classical monotone Nagumo front with high probability if the noise level and spatial discretization are sufficiently small.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An adaptive voter model on simplicial complexes

Collective decision making processes lie at the heart of many social, political and economic challenges. The classical voter model is a well-established conceptual model to study such processes. In this work, we define a new form of adaptive (or co-evolutionary) voter model posed on a simplicial complex, i.e., on a certain class of hypernetworks/hypergraphs. We use the persuasion rule along edges of the classical voter model and the recently studied re-wiring rule of edges towards like-minded nodes, and introduce a new peer pressure rule applied to three nodes connected via a 2-simplex. This simplicial adaptive voter model is studied via numerical simulation. We show that adding the effect of peer pressure to an adaptive voter model leaves its fragmentation transition, i.e., the transition upon varying the re-wiring rate from a single majority state into to a fragmented state of two different opinion subgraphs, intact. Yet, above and below the fragmentation transition, we observe that the peer pressure has substantial quantitative effects. It accelerates the transition to a single-opinion state below the transition and also speeds up the system dynamics towards fragmentation above the transition. Furthermore, we quantify that there is a multiscale hierarchy in the model leading to the depletion of 2-simplices, before the depletion of active edges. This leads to the conjecture that many other dynamic network models on simplicial complexes may show a similar behaviour with respect to the sequential evolution of simplicies of different dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On fast-slow consensus networks with a dynamic weight

We study dynamic networks under an undirected consensus communication protocol and with one state-dependent weighted edge. We assume that the aforementioned dynamic edge can take values over the whole real numbers, and that its behaviour depends on the nodes it connects and on an extrinsic slow variable. We show that, under mild conditions on the weight, there exists a reduction such that the dynamics of the network are organized by a transcritical singularity. As such, we detail a slow passage through a transcritical singularity for a simple network, and we observe that an exchange between consensus and clustering of the nodes is possible. In contrast to the classical planar fast-slow transcritical singularity, the network structure of the system under consideration induces the presence of a maximal canard. Our main tool of analysis is the blow-up method. Thus, we also focus on tracking the effects of the blow-up transformation on the network's structure. We show that on each blow-up chart one recovers a particular dynamic network related to the original one. We further indicate a numerical issue produced by the slow passage through the transcritical singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Sample Paths Estimates for Stochastic Fast-Slow Systems driven by Fractional Brownian Motion

We analyze the effect of additive fractional noise with Hurst parameter $H > \frac{1}{2}$ on fast-slow systems. Our strategy is based on sample paths estimates, similar to the approach by Berglund and Gentz in the Brownian motion case. Yet, the setting of fractional Brownian motion does not allow us to use the martingale methods from fast-slow systems with Brownian motion. We thoroughly investigate the case where the deterministic system permits a uniformly hyperbolic stable slow manifold. In this setting, we provide a neighborhood, tailored to the fast-slow structure of the system, that contains the process with high probability. We prove this assertion by providing exponential error estimates on the probability that the system leaves this neighborhood. We also illustrate our results in an example arising in climate modeling, where time-correlated noise processes have become of greater relevance recently.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Early-Warning Signs for Pattern-Formation in Stochastic Partial Differential Equations

There have been significant recent advances in our understanding of the potential use and limitations of early-warning signs for predicting drastic changes, so called critical transitions or tipping points, in dynamical systems. A focus of mathematical modeling and analysis has been on stochastic ordinary differential equations, where generic statistical early-warning signs can be identified near bifurcation-induced tipping points. In this paper, we outline some basic steps to extend this theory to stochastic partial differential equations with a focus on analytically characterizing basic scaling laws for linear SPDEs and comparing the results to numerical simulations of fully nonlinear problems. In particular, we study stochastic versions of the Swift-Hohenberg and Ginzburg-Landau equations. We derive a scaling law of the covariance operator in a regime where linearization is expected to be a good approximation for the local fluctuations around deterministic steady states. We compare these results to direct numerical simulation, and study the influence of noise level, noise color, distance to bifurcation and domain size on early-warning signs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From random Poincaré maps to stochastic mixed-mode-oscillation patterns

We quantify the effect of Gaussian white noise on fast--slow dynamical systems with one fast and two slow variables, which display mixed-mode oscillations owing to the presence of a folded-node singularity. The stochastic system can be described by a continuous-space, discrete-time Markov chain, recording the returns of sample paths to a Poincaré section. We provide estimates on the kernel of this Markov chain, depending on the system parameters and the noise intensity. These results yield predictions on the observed random mixed-mode oscillation patterns. Our analysis shows that there is an intricate interplay between the number of small-amplitude oscillations and the global return mechanism. In combination with a local saturation phenomenon near the folded node, this interplay can modify the number of small-amplitude oscillations after a large-amplitude oscillation. Finally, sufficient conditions are derived which determine when the noise increases the number of small-amplitude oscillations and when it decreases this number.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Bounded Positive Stationary Solutions for a Nonlocal Fisher-KPP Equation

We study the existence of stationary solutions for a nonlocal version of the Fisher-Kolmogorov-Petrovskii-Piscounov (FKPP) equation. The main motivation is a recent study by Berestycki et {al.} [Nonlinearity 22 (2009), {pp.}~2813--2844] where the nonlocal FKPP equation has been studied and it was shown for the spatial domain $\mathbb{R}$ andsufficiently small nonlocality that there are only two bounded non-negative stationary solutions. Here we provide a similar result for $\mathbb{R}^d$ using a completely different approach. In particular, an abstract perturbation argument is used in suitable weighted Sobolev spaces. One aim of the alternative strategy is that it can eventually be generalized to obtain persistence results for hyperbolic invariant sets for other nonlocal evolution equations on unbounded domains with small nonlocality, {i.e.}, to improve our understanding in applications when a small nonlocal influence alters the dynamics and when it does not.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Critical Transitions in Social Network Activity

A large variety of complex systems in ecology, climate science, biomedicine and engineering have been observed to exhibit tipping points, where the internal dynamical state of the system abruptly changes. For example, such critical transitions may result in the sudden change of ecological environments and climate conditions. Data and models suggest that detectable warning signs may precede some of these drastic events. This view is also corroborated by abstract mathematical theory for generic bifurcations in stochastic multi-scale systems. Whether the stochastic scaling laws used as warning signs are also present in social networks that anticipate a-priori {\it unknown} events in society is an exciting open problem, to which at present only highly speculative answers can be given. Here, we instead provide a first step towards tackling this formidable question by focusing on a-priori {\it known} events and analyzing a social network data set with a focus on classical variance and autocorrelation warning signs. Our results thus pertain to one absolutely fundamental question: Can the stochastic warning signs known from other areas also be detected in large-scale social network data? We answer this question affirmatively as we find that several a-priori known events are preceded by variance and autocorrelation growth. Our findings thus clearly establish the necessary starting point to further investigate the relation between abstract mathematical theory and various classes of critical transitions in social networks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Normal Hyperbolicity and Unbounded Critical Manifolds

This work is motivated by mathematical questions arising in differential equation models for autocatalytic reactions. In particular, this paper answers an open question posed by Guckenheimer and Scheper [SIAM J. Appl. Dyn. Syst. 10-1 (2011), pp. 92-128] and provides a more general theoretical approach to parts of the work by Gucwa and Szmolyan [Discr. Cont. Dyn. Sys.-S. 2-4 (2009), pp. 783-806]. We extend the local theory of singularities in fast-slow polynomial vector fields to classes of unbounded manifolds which lose normal hyperbolicity due to an alignment of the tangent and normal bundles. A projective transformation is used to localize the unbounded problem. Then the blow-up method is employed to characterize the loss of normal hyperbolicity for the transformed slow manifolds. Our analysis yields a rigorous scaling law for all unbounded manifolds which exhibit a power-law decay for the alignment with a fast subsystem domain. Furthermore, the proof also provides a technical extension of the blow-up method itself by augmenting the analysis with an optimality criterion for the blow-up exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A mathematical framework for critical transitions: normal forms, variance and applications

Critical transitions occur in a wide variety of applications including mathematical biology, climate change, human physiology and economics. Therefore it is highly desirable to find early-warning signs. We show that it is possible to classify critical transitions by using bifurcation theory and normal forms in the singular limit. Based on this elementary classification, we analyze stochastic fluctuations and calculate scaling laws of the variance of stochastic sample paths near critical transitions for fast subsystem bifurcations up to codimension two. The theory is applied to several models: the Stommel-Cessi box model for the thermohaline circulation from geoscience, an epidemic-spreading model on an adaptive network, an activator-inhibitor switch from systems biology, a predator-prey system from ecology and to the Euler buckling problem from classical mechanics. For the Stommel-Cessi model we compare different detrending techniques to calculate early-warning signs. In the epidemics model we show that link densities could be better variables for prediction than population densities. The activator-inhibitor switch demonstrates effects in three time-scale systems and points out that excitable cells and molecular units have information for subthreshold prediction. In the predator-prey model explosive population growth near a codimension two bifurcation is investigated and we show that early-warnings from normal forms can be misleading in this context. In the biomechanical model we demonstrate that early-warning signs for buckling depend crucially on the control strategy near the instability which illustrates the effect of multiplicative noise.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Characterizing Slow Exit and Entrance Points

Geometric Singular Perturbation Theory (GSPT) and Conley Index Theory are two powerful techniques to analyze dynamical systems. Conley already realized that using his index is easier for singular perturbation problems. In this paper, we will revisit Conley's results and prove that the GSPT technique of Fenichel Normal Form can be used to simplify the application of Conley index techniques even further. We also hope that our results provide a better bridge between the different fields. Furthermore we show how to interpret Conley's conditions in terms of averaging. The result are illustrated by the two-dimensional van der Pol equation and by a three-dimensional Morris-Lecar model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computing Slow Manifolds of Saddle Type

Slow manifolds are important geometric structures in the state spaces of dynamical systems with multiple time scales. This paper introduces an algorithm for computing trajectories on slow manifolds that are normally hyperbolic with both stable and unstable fast manifolds. We present two examples of bifurcation problems where these manifolds play a key role and a third example in which saddle-type slow manifolds are part of a traveling wave profile of a partial differential equation. Initial value solvers are incapable of computing trajectories on saddle-type slow manifolds, so the slow manifold of saddle type (SMST) algorithm presented here is formulated as a boundary value method. We take an empirical approach here to assessing the accuracy and effectiveness of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Deterministic continutation of stochastic metastable equilibria via Lyapunov equations and ellipsoids

Numerical continuation methods for deterministic dynamical systems have been one of the most successful tools in applied dynamical systems theory. Continuation techniques have been employed in all branches of the natural sciences as well as in engineering to analyze ordinary, partial and delay differential equations. Here we show that the deterministic continuation algorithm for equilibrium points can be extended to track information about metastable equilibrium points of stochastic differential equations (SDEs). We stress that we do not develop a new technical tool but that we combine results and methods from probability theory, dynamical systems, numerical analysis, optimization and control theory into an algorithm that augments classical equilibrium continuation methods. In particular, we use ellipsoids defining regions of high concentration of sample paths. It is shown that these ellipsoids and the distances between them can be efficiently calculated using iterative methods that take advantage of the numerical continuation framework. We apply our method to a bistable neural competition model and a classical predator-prey system. Furthermore, we show how global assumptions on the flow can be incorporated - if they are available - by relating numerical continuation, Kramers' formula and Rayleigh iteration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From First Lyapunov Coefficients to Maximal Canards

Hopf bifurcations in fast-slow systems of ordinary differential equations can be associated with surprising rapid growth of periodic orbits. This process is referred to as canard explosion. The key step in locating a canard explosion is to calculate the location of a special trajectory, called a maximal canard, in parameter space. A first-order asymptotic expansion of this location was found by Krupa and Szmolyan in the framework of a "canard point"-normal-form for systems with one fast and one slow variable. We show how to compute the coefficient in this expansion using the first Lyapunov coefficient at the Hopf bifurcation thereby avoiding use of this normal form. Our results connect the theory of canard explosions with existing numerical software, enabling easier calculations of where canard explosions occur.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Homoclinic Orbits of the FitzHugh-Nagumo Equation: Bifurcations in the Full System

This paper investigates travelling wave solutions of the FitzHugh-Nagumo equation from the viewpoint of fast-slow dynamical systems. These solutions are homoclinic orbits of a three dimensional vector field depending upon system parameters of the FitzHugh-Nagumo model and the wave speed. Champneys et al. [A.R. Champneys, V. Kirk, E. Knobloch, B.E. Oldeman, and J. Sneyd, When Shilnikov meets Hopf in excitable systems, SIAM Journal of Applied Dynamical Systems, 6(4), 2007] observed sharp turns in the curves of homoclinic bifurcations in a two dimensional parameter space. This paper demonstrates numerically that these turns are located close to the intersection of two curves in the parameter space that locate non-transversal intersections of invariant manifolds of the three dimensional vector field. The relevant invariant manifolds in phase space are visualized. A geometrical model inspired by the numerical studies displays the sharp turns of the homoclinic bifurcations curves and yields quantitative predictions about multi-pulse and homoclinic orbits and periodic orbits that have not been resolved in the FitzHugh-Nagumo model. Further observations address the existence of canard explosions and mixed-mode oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Homoclinic Orbits of the FitzHugh-Nagumo Equation: The Singular-Limit

The FitzHugh-Nagumo equation has been investigated with a wide array of different methods in the last three decades. Recently a version of the equations with an applied current was analyzed by Champneys, Kirk, Knobloch, Oldeman and Sneyd using numerical continuation methods. They obtained a complicated bifurcation diagram in parameter space featuring a C-shaped curve of homoclinic bifurcations and a U-shaped curve of Hopf bifurcations. We use techniques from multiple time-scale dynamics to understand the structures of this bifurcation diagram based on geometric singular perturbation analysis of the FitzHugh-Nagumo equation. Numerical and analytical techniques show that if the ratio of the time-scales in the FitzHugh-Nagumo equation tends to zero, then our singular limit analysis correctly represents the observed CU-structure. Geometric insight from the analysis can even be used to compute bifurcation curves which are inaccessible via continuation methods. The results of our analysis are summarized in a singular bifurcation diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scaling of Saddle-Node Bifurcations: Degeneracies and Rapid Quantitative Changes

The scaling of the time delay near a "bottleneck" of a generic saddle-node bifurcation is well-known to be given by an inverse square-root law. We extend the analysis to several non-generic cases for smooth vector fields. We proceed to investigate $C^0$ vector fields. Our main result is a new phenomenon in two-parameter families having a saddle-node bifurcation upon changing the first parameter. We find distinct scalings for different values of the second parameter ranging from power laws with exponents in (0,1) to scalings given by O(1). We illustrate this rapid quantitative change of the scaling law by a an overdamped pendulum with varying length.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Warning signs for wave speed transitions of noisy Fisher-KPP invasion fronts

Invasion waves are a fundamental building block of theoretical ecology. In this study we aim to take the first steps to link propagation failure and fast acceleration of traveling waves to critical transitions (or tipping points). The approach is based upon a detailed numerical study of various versions of the Fisher-Kolmogorov-Petrovskii-Piscounov (FKPP) equation. The main motivation of this work is to contribute to the following question: how much information do statistics, collected by a stationary observer, contain about the speed and bifurcations of traveling waves? We suggest warning signs based upon closeness to carrying capacity, second-order moments and transients of localized initial invasions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A mathematical framework for critical transitions: bifurcations, fast-slow systems and stochastic dynamics

Bifurcations can cause dynamical systems with slowly varying parameters to transition to far-away attractors. The terms ``critical transition'' or ``tipping point'' have been used to describe this situation. Critical transitions have been observed in an astonishingly diverse set of applications from ecosystems and climate change to medicine and finance. The main goal of this paper is to give an overview which standard mathematical theories can be applied to critical transitions. We shall focus on early-warning signs that have been suggested to predict critical transitions and point out what mathematical theory can provide in this context. Starting from classical bifurcation theory and incorporating multiple time scale dynamics one can give a detailed analysis of local bifurcations that induce critical transitions. We suggest that the mathematical theory of fast-slow systems provides a natural definition of critical transitions. Since noise often plays a crucial role near critical transitions the next step is to consider stochastic fast-slow systems. The interplay between sample path techniques, partial differential equations and random dynamical systems is highlighted. Each viewpoint provides potential early-warning signs for critical transitions. Since increasing variance has been suggested as an early-warning sign we examine it in the context of normal forms analytically, numerically and geometrically; we also consider autocorrelation numerically. Hence we demonstrate the applicability of early-warning signs for generic models. We end with suggestions for future directions of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical analysis of evolution equations in generalized models

Generalized models provide a framework for the study of evolution equations without specifying all functional forms. The generalized formulation of problems has been shown to facilitate the analytical investigation of local dynamics and has been used successfully to answer applied questions. Yet their potential to facilitate analytical computations has not been realized in the mathematical literature. In the present paper we introduce the method of generalized modeling in mathematical terms, supporting the key steps of the procedure by rigorous proofs. Further, we point out open questions that are in the scope of present mathematical research and, if answered could greatly increase the predictive power of generalized models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonlocal Generalized Models of Predator-Prey Systems

The method of generalized modeling has been applied successfully in many different contexts, particularly in ecology and systems biology. It can be used to analyze the stability and bifurcations of steady-state solutions. Although many dynamical systems in mathematical biology exhibit steady-state behaviour one also wants to understand nonlocal dynamics beyond equilibrium points. In this paper we analyze predator-prey dynamical systems and extend the method of generalized models to periodic solutions. First, we adapt the equilibrium generalized modeling approach and compute the unique Floquet multiplier of the periodic solution which depends upon so-called generalized elasticity and scale functions. We prove that these functions also have to satisfy a flow on parameter (or moduli) space. Then we use Fourier analysis to provide computable conditions for stability and the moduli space flow. The final stability analysis reduces to two discrete convolutions which can be interpreted to understand when the predator-prey system is stable and what factors enhance or prohibit stable oscillatory behaviour. Finally, we provide a sampling algorithm for parameter space based on nonlinear optimization and the Fast Fourier Transform which enables us to gain a statistical understanding of the stability properties of periodic predator-prey dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On decomposing mixed-mode oscillations and their return maps

Alternating patterns of small and large amplitude oscillations occur in a wide variety of physical, chemical, biological and engineering systems. These mixed-mode oscillations (MMOs) are often found in systems with multiple time scales. Previous differential equation modeling and analysis of MMOs has mainly focused on local mechanisms to explain the small oscillations. Numerical continuation studies reported different MMO patterns based on parameter variation. This paper aims at improving the link between local analysis and numerical simulation. Our starting point is a numerical study of a singular return map for the Koper model which is a prototypical example for MMOs that also relates to local normal form theory. We demonstrate that many MMO patterns can be understood geometrically by approximating the singular maps with affine and quadratic maps. Motivated by our numerical analysis we use abstract affine and quadratic return map models in combination with two local normal forms that generate small oscillations. Using this decomposition approach we can reproduce many classical MMO patterns and effectively decouple bifurcation parameters for local and global parts of the flow. The overall strategy we employ provides an alternative technique for understanding MMOs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On spatial and temporal multilevel dynamics and scaling effects in epileptic seizures

Epileptic seizures are one of the most well-known dysfunctions of the nervous system. During a seizure, a highly synchronized behavior of neural activity is observed that can cause symptoms ranging from mild sensual malfunctions to the complete loss of body control. In this paper, we aim to contribute towards a better understanding of the dynamical systems phenomena that cause seizures. Based on data analysis and modelling, seizure dynamics can be identified to possess multiple spatial scales and on each spatial scale also multiple time scales. At each scale, we reach several novel insights. On the smallest spatial scale we consider single model neurons and investigate early-warning signs of spiking. This introduces the theory of critical transitions to excitable systems. For clusters of neurons (or neuronal regions) we use patient data and find oscillatory behavior and new scaling laws near the seizure onset. These scalings lead to substantiate the conjecture obtained from mean-field models that a Hopf bifurcation could be involved near seizure onset. On the largest spatial scale we introduce a measure based on phase-locking intervals and wavelets into seizure modelling. It is used to resolve synchronization between different regions in the brain and identifies time-shifted scaling laws at different wavelet scales. We also compare our wavelet-based multiscale approach with maximum linear cross-correlation and mean-phase coherence measures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Time-Scale and Noise Optimality in Self-Organized Critical Adaptive Networks

Recent studies have shown that adaptive networks driven by simple local rules can organize into "critical" global steady states, providing another framework for self-organized criticality (SOC). We focus on the important convergence to criticality and show that noise and time-scale optimality are reached at finite values. This is in sharp contrast to the previously believed optimal zero noise and infinite time scale separation case. Furthermore, we discover a noise induced phase transition for the breakdown of SOC. We also investigate each of the three new effects separately by developing models. These models reveal three generically low-dimensional dynamical behaviors: time-scale resonance (TR), a new simplified version of stochastic resonance - which we call steady state stochastic resonance (SSR) - as well as noise-induced phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hunting French Ducks in a Noisy Environment

We consider the effect of Gaussian white noise on fast-slow dynamical systems with one fast and two slow variables, containing a folded-node singularity. In the absence of noise, these systems are known to display mixed-mode oscillations, consisting of alternating large- and small-amplitude oscillations. We quantify the effect of noise and obtain critical noise intensities above which the small-amplitude oscillations become hidden by fluctuations. Furthermore we prove that the noise can cause sample paths to jump away from so-called canard solutions with high probability before deterministic orbits do. This early-jump mechanism can drastically influence the local and global dynamics of the system by changing the mixed-mode patterns.

0 citations0 reviewsNo code linkedOpen access