Oscillatory survival probability: analytical and numerical study of a non-Poissonian exit time

DOI10.1137/151004100zbMATH Open1343.60121arXiv1411.6106OpenAlexW2351176355WikidataQ58051978 ScholiaQ58051978MaRDI QIDQ2806426FDOQ2806426

Zeev Schuss, D. Holcman, K. Dao Duc

Publication date: 18 May 2016

Published in: Multiscale Modeling \& Simulation (Search for Journal in Brave)

Abstract: We study the escape of Brownian motion from the domain of attraction $Omega$ of a stable focus with a strong drift. The boundary $partialOmega$ of $Omega$ is an unstable limit cycle of the drift and the focus is very close to the limit cycle. We find a new phenomenon of oscillatory decay of the peaks of the survival probability of the Brownian motion in $Omega$. We compute explicitly the complex-valued second eigenvalue $lambd{a}_2(Omega$) of the Fokker-Planck operator with Dirichlet boundary conditions and show that it is responsible for the peaks. Specifically, we demonstrate that the dominant oscillation frequency equals $1/mathfrakImlambd{a}_2(Omega)$ and is independent of the relative noise strength. We apply the analysis to a canonical system and compare the density of exit points on $partialOmega$ to that obtained from stochastic simulations. We find that this density is concentrated in a small portion of the boundary, thus rendering the exit a narrow escape problem. Unlike the case in the classical activated escape problem, the principal eigenvalue does not necessarily decay exponentially as the relative noise strength decays. The oscillatory narrow escape problem arises in a mathematical model of neural networks with synaptic depression. We identify oscillation peaks in the density of the time the network spends in a specific state. This observation explains the oscillations of stochastic trajectories around a focus prior to escape and also the non-Poissonian escape times. This phenomenon has been observed and reported in the neural network literature.

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

zbMATH Keywords

neuroscienceeigenvalueFokker-Planck equationlimit cycleoscillationmean first passage timeWKB asymptoticsexit pointsplanar diffusion process

Mathematics Subject Classification ID

Diffusion processes (60J60) Neural biology (92C20) Stochastic ordinary differential equations (aspects of stochastic analysis) (60H10) Ordinary differential equations and systems with randomness (34F05) Boundary value problems for second-order elliptic systems (35J57) Fokker-Planck equations (35Q84) Numerical solution of boundary value problems involving ordinary differential equations (65L10) Numerical solution of eigenvalue problems involving ordinary differential equations (65L15)

Cites Work

• Title not available (Why is that?)
• Elements of applied bifurcation theory
• Singular Perturbation Methods in Stochastic Differential Equations of Mathematical Physics
• Theory and applications of stochastic processes. An analytical approach
• Diffusion Across Characteristic Boundaries with Critical Points
• Diffusion Across Characteristic Boundaries
• The Narrow Escape Problem
• Oscillatory Behavior of the Rate of Escape through an Unstable Limit Cycle
• Strong Localized Perturbations of Eigenvalue Problems
• Summing Logarithmic Expansions for Singularly Perturbed Eigenvalue Problems
• A Singular Perturbation Method for the Computation of the Mean First Passage Time in a Nonlinear Filter
• Nonlinear filtering and optimal phase tracking
• Stochastic Narrow Escape in Molecular and Cellular Biology
• The Exit Problem for Randomly Perturbed Dynamical Systems
• Cycling and skewing of exit measures for planar systems
• Title not available (Why is that?)

Cited In (3)

• Exit versus escape for stochastic dynamical systems and application to the computation of the bursting time duration in neuronal networks
• Tunnel effect and analysis of the survival amplitude in the nonlinear Winter's model
• Dynamic Looping of a Free-Draining Polymer

This page was built for publication: Oscillatory survival probability: analytical and numerical study of a non-Poissonian exit time

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2806426)