An ordered line integral method for computing the quasi-potential in the case of variable anisotropic diffusion
From MaRDI portal
Publication:2000680
Abstract: Nongradient stochastic differential equations (SDEs) with position-dependent and anisotropic diffusion are often used in biological modeling. The quasi-potential is a crucial function in the Large Deviation Theory that allows one to estimate transition rates between attractors of the corresponding ordinary differential equation and find the maximum likelihood transition paths. Unfortunately, the quasi-potential can rarely be found analytically. It is defined as the solution to a certain action minimization problem. In this work, the recently introduced Ordered Line Integral Method (OLIM) is extended for computing the quasi-potential for 2D SDEs with anisotropic and position-dependent diffusion scaled by a small parameter on a regular rectangular mesh. The presented solver employs the dynamical programming principle. At each step, a local action minimization problem is solved using straight line path segments and the midpoint quadrature rule. The solver is tested on two examples where analytic formulas for the quasi-potential are available. The dependence of the computational error on the mesh size, the update factor K (a key parameter of OLIMs), as well as the degree and the orientation of anisotropy is established. The effect of anisotropy on the quasi-potential and the maximum likelihood paths is demonstrated on the Maier-Stein model. The proposed solver is applied to find the quasi-potential and the maximum likelihood transition paths in a model of the genetic switch in Lambda Phage between the lysogenic state where the phage reproduces inside the infected cell without killing it, and the lytic state where the phage destroys the infected cell.
Recommendations
- Ordered line integral methods for computing the quasi-potential
- Finding the quasipotential for nongradient SDEs
- Computing the quasipotential for highly dissipative and chaotic SDEs an application to stochastic Lorenz'63
- An efficient jet marcher for computing the quasipotential for 2D SDEs. Enhancing accuracy and efficiency of quasipotential solvers
- A machine learning method for computing quasi-potential of stochastic dynamical systems
Cites work
- scientific article; zbMATH DE number 3126094 (Why is no real title available?)
- scientific article; zbMATH DE number 852530 (Why is no real title available?)
- A scaling theory of bifurcations in the symmetric weak-noise escape problem.
- Crossing the quasi-threshold manifold of a noise-driven excitable system
- Fast two-scale methods for eikonal equations
- Finding the quasipotential for nongradient SDEs
- Generalisation of the Eyring-Kramers transition rate formula to irreversible diffusion processes
- Limiting Exit Location Distributions in the Stochastic Exit Problem
- Ordered Upwind Methods for Static Hamilton--Jacobi Equations: Theory and Algorithms
- Ordered line integral methods for computing the quasi-potential
- Ordered upwind methods for static Hamilton–Jacobi equations
- Random Perturbations of Dynamical Systems
- SMOLUCHOWSKI–KRAMERS APPROXIMATION AND EXIT PROBLEMS
- Study of noise-induced transitions in the Lorenz system using the minimum action method
- The geometric minimum action method: A least action principle on the space of curves
Cited in
(15)- Ordered line integral methods for solving the eikonal equation
- A machine learning method for computing quasi-potential of stochastic dynamical systems
- Dispersion of the prehistory distribution for non-gradient systems
- Computing the quasipotential for nongradient SDEs in 3D
- Controlling mean exit time of stochastic dynamical systems based on quasipotential and machine learning
- Parametric resonance for enhancing the rate of metastable transition
- Computing the quasipotential for highly dissipative and chaotic SDEs an application to stochastic Lorenz'63
- Unified mechanism of inverse stochastic resonance for monostability and bistability in Hindmarsh-rose neuron
- Ordered line integral methods for computing the quasi-potential
- Optimal control for sampling the transition path process and estimating rates
- Noise-induced transition in the Zeldovich–Semenov model with local and global bifurcations
- Finding the quasipotential for nongradient SDEs
- Noise induced escape from stable invariant tori
- An efficient jet marcher for computing the quasipotential for 2D SDEs. Enhancing accuracy and efficiency of quasipotential solvers
- Sex, ducks, and rock ``n roll: mathematical model of sexual response
This page was built for publication: An ordered line integral method for computing the quasi-potential in the case of variable anisotropic diffusion
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2000680)
