An ordered line integral method for computing the quasi-potential in the case of variable anisotropic diffusion

DOI10.1016/J.PHYSD.2018.07.002zbMATH Open1415.65017arXiv1806.05321OpenAlexW2808097552WikidataQ129565169 ScholiaQ129565169MaRDI QIDQ2000680FDOQ2000680

Maria Kourkina Cameron, D. Dahiya

Publication date: 28 June 2019

Published in: Physica D (Search for Journal in Brave)

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.

Full work available at URL: https://arxiv.org/abs/1806.05321

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

zbMATH Keywords

quasipotentiallambda phageminimum action pathMaier-Stein modelordered line integral methodvariable and anisotropic diffusion

Mathematics Subject Classification ID

Stochastic partial differential equations (aspects of stochastic analysis) (60H15) Biochemistry, molecular biology (92C40) Numerical solutions to stochastic differential and integral equations (65C30)

Cites Work

• Title not available (Why is that?)
• Ordered Upwind Methods for Static Hamilton--Jacobi Equations: Theory and Algorithms
• Random Perturbations of Dynamical Systems
• Generalisation of the Eyring-Kramers transition rate formula to irreversible diffusion processes
• The geometric minimum action method: A least action principle on the space of curves
• Study of noise-induced transitions in the Lorenz system using the minimum action method
• A scaling theory of bifurcations in the symmetric weak-noise escape problem.
• Finding the quasipotential for nongradient SDEs
• SMOLUCHOWSKI–KRAMERS APPROXIMATION AND EXIT PROBLEMS
• Ordered upwind methods for static Hamilton–Jacobi equations
• Fast two-scale methods for eikonal equations
• Limiting Exit Location Distributions in the Stochastic Exit Problem
• Title not available (Why is that?)
• Ordered line integral methods for computing the quasi-potential
• Crossing the quasi-threshold manifold of a noise-driven excitable system

Cited In (13)

• Parametric Resonance for Enhancing the Rate of Metastable Transition
• 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
• Computing the quasipotential for highly dissipative and chaotic SDEs an application to stochastic Lorenz'63
• Optimal control for sampling the transition path process and estimating rates
• Noise-induced transition in the Zeldovich–Semenov model with local and global bifurcations
• Unified mechanism of inverse stochastic resonance for monostability and bistability in Hindmarsh–Rose neuron
• 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

Uses Software

• QPot

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)