Construction of quasi-potentials for stochastic dynamical systems: an optimization approach
arXiv:1805.07273 · doi:10.1103/PhysRevE.98.022136
Abstract
The construction of effective and informative landscapes for stochastic dynamical systems has proven a long-standing and complex problem. In many situations, the dynamics may be described by a Langevin equation while constructing a landscape comes down to obtaining the quasi-potential, a scalar function that quantifies the likelihood of reaching each point in the state-space. In this work we provide a novel method for constructing such landscapes by extending a tool from control theory: the Sum-of-Squares method for generating Lyapunov functions. Applicable to any system described by polynomials, this method provides an analytical polynomial expression for the potential landscape, in which the coefficients of the polynomial are obtained via a convex optimization problem. The resulting landscapes are based upon a decomposition of the deterministic dynamics of the original system, formed in terms of the gradient of the potential and a remaining "curl" component. By satisfying the condition that the inner product of the gradient of the potential and the remaining dynamics is everywhere negative, our derived landscapes provide both upper and lower bounds on the true quasi-potential; these bounds becoming tight if the decomposition is orthogonal. The method is demonstrated to correctly compute the quasi-potential for high-dimensional linear systems and also for a number of nonlinear examples.
arXiv comment12 pages, 4 figures
References in corpus (20)
- JuMP: A Modeling Language for Mathematical Optimization
- Predator-prey cycles from resonant amplification of demographic stochasticity
- Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
- Stochastic Dynamical Structure (SDS) of Nonequilibrium Processes in the Absence of Detailed Balance. II: construction of SDS with nonlinear force and multiplicative noise
- Structure of Stochastic Dynamics near Fixed Points
- A Scaling Theory of Bifurcations in the Symmetric Weak-Noise Escape Problem
- Generalisation of the Eyring-Kramers transition rate formula to irreversible diffusion processes
- Stochastic phenotype transition of a single cell in an intermediate region of gene-state switching
- Diffusive dynamics and stochastic models of turbulent axisymmetric wakes
- Realization of Waddington's Metaphor: Potential Landscape, Quasi-potential, A-type Integral and Beyond
- Designing Attractive Models via Automated Identification of Chaotic and Oscillatory Dynamical Regimes
- Stochastic Dynamical Structure (SDS) of Nonequilibrium Processes in the Absence of Detailed Balance. III: potential function in local stochastic dynamics and in steady state of Boltzmann-Gibbs type distribution function
- Topological resilience in non-normal networked systems
- Enhanced stochastic oscillations in autocatalytic reactions
- Constructive Proof of Global Lyapunov Function as Potential Function
- Ordered Line Integral Methods for Computing the Quasi-potential
- Noise driven neuromorphic tuned amplifier
- Minimum Action Path theory reveals the details of stochastic biochemical transitions out of oscillatory cellular states
- Non-normality, reactivity, and intrinsic stochasticity in neural dynamics: a non-equilibrium potential approach
- Construction of the quasi-potential for linear SDEs using false quasi-potentials and a geometric recursion