Escape driven by alpha-stable white noises

Exponential increase of transition rates in metastable systems driven by non-Gaussian noise

Noise-induced escape from metastable states governs a plethora of transition phenomena in physics, chemistry, and biology. While the escape problem in the presence of thermal Gaussian noise has been well understood since the seminal works of Arrhenius and Kramers, many systems, in particular living ones, are effectively driven by non-Gaussian noise for which the conventional theory does not apply. Here we present a theoretical framework based on path integrals that allows the calculation of both escape rates and optimal escape paths for a generic class of non-Gaussian noises. We find that non-Gaussian noise always leads to more efficient escape and can enhance escape rates by many orders of magnitude compared with thermal noise, highlighting that away from equilibrium escape rates cannot be reliably modelled based on the traditional Arrhenius-Kramers result. Our analysis also identifies a new universality class of non-Gaussian noises, for which escape paths are dominated by large jumps.

Lévy on-off intermittency

We present an alternative form of intermittency, Lévy on-off intermittency, which arises from multiplicative α-stable white noise close to an instability threshold. We study this problem in the linear and nonlinear regimes, both theoretically and numerically, for the case of a pitchfork bifurcation with fluctuating growth rate. We compute the stationary distribution analytically and numerically from the associated fractional Fokker-Planck equation in the Stratonovich interpretation. We characterize the system in the parameter space (α,β) of the noise, with stability parameter α∈(0,2) and skewness parameter β∈[-1,1]. Five regimes are identified in this parameter space, in addition to the well-studied Gaussian case α=2. Three regimes are located at 1<α<2, where the noise has finite mean but infinite variance. They are differentiated by β and all display a critical transition at the deterministic instability threshold, with on-off intermittency close to onset. Critical exponents are computed from the stationary distribution. Each regime is characterized by a specific form of the density and specific critical exponents, which differ starkly from the Gaussian case. A finite or infinite number of integer-order moments may converge, depending on parameters. Two more regimes are found at 0<α≤1. There, the mean of the noise diverges, and no critical transition occurs. In one case, the origin is always unstable, independently of the distance μ from the deterministic threshold. In the other case, the origin is conversely always stable, independently of μ. We thus demonstrate that an instability subject to nonequilibrium, power-law-distributed fluctuations can display substantially different properties than for Gaussian thermal fluctuations, in terms of statistics and critical behavior.

Statistics of residence time for Lévy flights in unstable parabolic potentials

We analyze the residence time problem for an arbitrary Markovian process describing nonlinear systems without a steady state. We obtain exact analytical results for the statistical characteristics of the residence time. For diffusion in a fully unstable potential profile in the presence of Lévy noise we get the conditional probability density of the particle position and the average residence time. The noise-enhanced stability phenomenon is observed in the system investigated. Results from numerical simulations are in very good agreement with analytical ones.

Peculiarities of escape kinetics in the presence of athermal noises

Stochastic evolution of various dynamic systems and reaction networks is commonly described in terms of noise assisted escape of an overdamped particle from a potential well, as devised by the paradigmatic Langevin equation in which additive Gaussian stochastic force reproduces effects of thermal fluctuations from the reservoir. When implemented for systems close to equilibrium, the approach correctly explains the emergence of the Boltzmann distribution for the ensemble of trajectories generated by the Langevin equation and relates the intensity of the noise strength to the mobility. This scenario can be further generalized to include effects of non-Gaussian, burstlike forcing modeled by Lévy noise. In this case, however, the pulsatile additive noise cannot be treated as the internal (thermal) since the relation between the strength of the friction and variance of the noise is violated. Heavy tails of Lévy noise distributions not only facilitate escape kinetics, but also, more importantly, change the escape protocol by altering the final stationary state to a non-Boltzmann, nonequilibrium form. As a result, contrary to the kinetics induced by a Gaussian white noise, escape rates in environments with Lévy noise are determined not by the barrier height, but instead by the barrier width. We further discuss consequences of simultaneous action of thermal and Lévy noises on statistics of passage times and population of reactants in double-well potentials.

Multimodal stationary states under Cauchy noise

A Lévy noise is an efficient description of out-of-equilibrium systems. The presence of Lévy flights results in a plenitude of noise-induced phenomena. Among others, Lévy flights can produce stationary states with more than one modal value in single-well potentials. Here we explore stationary states in special double-well potentials demonstrating that a sufficiently high potential barrier separating potential wells can produce bimodal stationary states in each potential well. Furthermore, we explore how the decrease in the barrier height affects the multimodality of stationary states. Finally, we explore the role of multimodality of stationary states on noise-induced escape over the static potential barrier.

Anomalous transport effects on switching currents of graphene-based Josephson junctions

We explore the effect of noise on the ballistic graphene-based small Josephson junctions in the framework of the resistively and capacitively shunted model. We use the non-sinusoidal current-phase relation specific for graphene layers partially covered by superconducting electrodes. The noise induced escapes from the metastable states, when the external bias current is ramped, given the switching current distribution, i.e. the probability distribution of the passages to finite voltage from the superconducting state as a function of the bias current, that is the information more promptly available in the experiments. We consider a noise source that is a mixture of two different types of processes: a Gaussian contribution to simulate an uncorrelated ordinary thermal bath, and non-Gaussian, α-stable (or Lévy) term, generally associated to non-equilibrium transport phenomena. We find that the analysis of the switching current distribution makes it possible to efficiently detect a non-Gaussian noise component in a Gaussian background.

The Switch in a Genetic Toggle System with Lévy Noise

A bistable toggle switch is a paradigmatic model in the field of biology. The dynamics of the system induced by Gaussian noise has been intensively investigated, but Gaussian noise cannot incorporate large bursts typically occurring in real experiments. This paper aims to examine effects of variations from one protein imposed by a non-Gaussian Lévy noise, which is able to describe even large jumps, on the coherent switch and the on/off switch via the steady-state probability density, the joint steady-state probability density, and the mean first passage time. We find that a large burst of one protein due to the Lévy noises can induce coherent switches even with small noise intensities in contrast to the Gaussian case which requires large intensities for this. The influences of the stability index, skewness parameter and noise intensity on the on/off switch are analyzed, leading to an adjustment of the concentrations of both proteins and a decision which stable point to stay most. The mean first passage times show complex effects under Lévy noise, especially the stability index and skewness parameter. Our results also imply that the presence of non-Gaussian Lévy noises has fundamentally changed the escape mechanism in such a system compared with Gaussian noise.

Escape process in systems characterized by stable noises and position-dependent resting times

Stochastic systems characterized by a random driving in a form of the general stable noise are considered. The particle experiences long rests due to the traps the density of which is position dependent and obeys a power-law form attributed to the underlying self-similar structure. Both the one- and two-dimensional cases are analyzed. The random walk description involves a position-dependent waiting time distribution. On the other hand, the stochastic dynamics is formulated in terms of the subordination technique where the random time generator is position dependent. The first passage time problem is addressed by evaluating a first passage time density distribution and an escape rate. The influence of the medium nonhomogeneity on those quantities is demonstrated; moreover, the dependence of the escape rate on the stability index and the memory parameter is evaluated. Results indicate essential differences between the Gaussian case and the case involving Lévy flights.

Infinite product expansion of the Fokker–Planck equation with steady-state solution

We present an analytical technique for solving Fokker–Planck equations that have a steady-state solution by representing the solution as an infinite product rather than, as usual, an infinite sum. This method has many advantages: automatically ensuring positivity of the resulting approximation, and by design exactly matching both the short- and long-term behaviour. The efficacy of the technique is demonstrated via comparisons with computations of typical examples.

Stepping molecular motor amid Lévy white noise

We consider a model of a stepping molecular motor consisting of two connected heads. Directional motion of the stepper takes place along a one-dimensional track. Each head is subject to a periodic potential without spatial reflection symmetry. When the potential for one head is switched on, it is switched off for the other head. Additionally, the system is subject to the influence of symmetric, white Lévy noise that mimics the action of external random forcing. The stepper exhibits motion with a preferred direction which is examined by analyzing the median of the displacement of a midpoint between the positions of the two heads. We study the modified dynamics of the stepper by numerical simulations. We find flux reversals as noise parameters are changed. Speed and direction appear to very sensitively depend on characteristics of the noise.

Quantifying a resonant-activation-like phenomenon in non-Markovian systems

Resonant activation is an effect of a noise-induced escape over a modulated potential barrier. The modulation of an energy landscape facilitates the escape kinetics and makes it optimal as measured by the mean first-passage time. A canonical example of resonant activation is a Brownian particle moving in a time-dependent potential under action of Gaussian white noise. Resonant activation is observed not only in typical Markovian-Gaussian systems but also in far-from-equilibrium and far-from-Markovianity regimes. We demonstrate that using an alternative to the mean first-passage time, robust measures of resonant activation, the signature of this effect can be observed in general continuous-time random walks in modulated potentials, even in situations when the mean first-passage time diverges.

Path-integral formulation for Lévy flights: evaluation of the propagator for free, linear, and harmonic potentials in the over- and underdamped limits

Lévy flights can be described using a Fokker-Planck equation, which involves a fractional derivative operator in the position coordinate. Such an operator has its natural expression in the Fourier domain. Starting with this, we show that the solution of the equation can be written as a Hamiltonian path integral. Though this has been realized in the literature, the method has not found applications as the path integral appears difficult to evaluate. We show that a method in which one integrates over the position coordinates first, after which integration is performed over the momentum coordinates, can be used to evaluate several path integrals that are of interest. Using this, we evaluate the propagators for (a) free particle, (b) particle subjected to a linear potential, and (c) harmonic potential. In all the three cases, we have obtained results for both overdamped and underdamped cases.

Anomalous diffusion in systems driven by the stable Lévy noise with a finite noise relaxation time and inertia

Dynamical systems driven by a general Lévy stable noise are considered. The inertia is included and the noise, represented by a generalized Ornstein-Uhlenbeck process, has a finite relaxation time. A general linear problem (the additive noise) is solved: the resulting distribution converges with time to the distribution for the white-noise, massless case. Moreover, a multiplicative noise is discussed. It can make the distribution steeper and the variance, which is finite, depends sublinearly on time (subdiffusion). For a small mass, a white-noise limit corresponds to the Stratonovich interpretation. On the other hand, the distribution tails agree with the Itô interpretation if the inertia is very large. An escape time from the potential well is calculated.

FIRST EXIT TIMES OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS DRIVEN BY MULTIPLICATIVE LÉVY NOISE WITH HEAVY TAILS

In this paper, we study first exit times from a bounded domain of a gradient dynamical system Ẏt = -∇U(Yt) perturbed by a small multiplicative Lévy noise with heavy tails. A special attention is paid to the way the multiplicative noise is introduced. In particular, we determine the asymptotics of the first exit time of solutions of Itô, Stratonovich and Marcus canonical SDEs.

Hopping over a heat barrier

We analyze diffusion in a finite domain with a position-dependent diffusion coefficient in terms of a stochastic hopping process. Via a coordinate transformation, we map the original system onto a problem with constant diffusion but nontrivial potential. In this way we show that a regime with enhanced diffusion acts as a potential barrier. We compute first-passage time distributions, hopping rates, and eigenvalues of the Fokker-Planck operator, and thereby verify that diffusion with a heat barrier is equivalent to a hopping process between metastable states.

Escape from the potential well: competition between long jumps and long waiting times

Within a concept of the fractional diffusion equation and subordination, the paper examines the influence of a competition between long waiting times and long jumps on the escape from the potential well. Applying analytical arguments and numerical methods, we demonstrate that the presence of long waiting times distributed according to a power-law distribution with a diverging mean leads to very general asymptotic properties of the survival probability. The observed survival probability asymptotically decays like a power law whose form is not affected by the value of the exponent characterizing the power law jump length distribution. It is demonstrated that this behavior is typical of and generic for systems exhibiting long waiting times. We also show that the survival probability has a universal character not only asymptotically, but also at small times. Finally, it is indicated which properties of the first passage time density are sensitive to the exact value of the exponent characterizing the jump length distribution.

Nonlinear stochastic equations with multiplicative Lévy noise

The Langevin equation with a multiplicative Lévy white noise is solved. The noise amplitude and the drift coefficient have a power-law form. A validity of ordinary rules of the calculus for the Stratonovich interpretation is discussed. The solution has the algebraic asymptotic form and the variance may assume a finite value for the case of the Stratonovich interpretation. The problem of escaping from a potential well is analyzed numerically; predictions of different interpretations of the stochastic integral are compared.

Multiplicative Lévy processes: Itô versus Stratonovich interpretation

Langevin equation with a multiplicative stochastic force is considered. That force is uncorrelated, it has the Lévy distribution and the power-law intensity. The Fokker-Planck equations, which correspond both to the Itô and Stratonovich interpretation, are presented. They are solved for the case without drift and for the harmonic oscillator potential. The variance is evaluated; it is always infinite for the Itô case whereas for the Stratonovich one it can be finite and rise with time slower that linearly, which indicates subdiffusion. Analytical results are compared with numerical simulations.

Lévy noises: double stochastic resonance in a single-well potential

We study properties of a single-well fourth-order potential perturbed by a periodically modulated stable noise. Periodic modulation of the stable noise asymmetry results in an occurrence of the dynamical hysteresis which is the manifestation of the stochastic resonance in the system at hand. We show that the single-well potential with time modulated stable driving is a minimalistic setup, allowing the occurrence of the stochastic resonance (as measured by the hysteresis loop area). Finally, we demonstrate that the observed stochastic resonance is of the double type, i.e., the system efficiency measured by the hysteresis loop area depends in a nonmonotonous way both on the scale parameter (noise intensity) and on the stability exponent characterizing tails asymptotic of noise pulses.

Continuous-time multidimensional Markovian description of Lévy walks

The paper presents a multidimensional model for nonlinear Markovian random walks that generalizes the one we developed previously [I. Lubashevsky, R. Friedrich, and A. Heuer, Phys. Rev. E 79, 011110 (2009)] in order to describe the Lévy-type stochastic processes in terms of continuous trajectories of walker motion. This approach may open a way to treat Lévy flights or Lévy random walks in inhomogeneous media or systems with boundaries in the future. The proposed model assumes the velocity of a wandering particle to be affected by a linear friction and a nonlinear Langevin force whose intensity is proportional to the magnitude of the velocity for its large values. Based on the singular perturbation technique, the corresponding Fokker-Planck equation is analyzed and the relationship between the system parameters and the Lévy exponent is found. Following actually the previous paper we demonstrate also that anomalously long displacements of the wandering particle are caused by extremely large fluctuations in the particle velocity whose duration is determined by the system parameters rather than the duration of the observation interval. In this way we overcome the problem of ascribing to Lévy random-walk non-Markov properties.

Numerical solution of stochastic differential equations with Poisson and Lévy white noise

A fixed time step method is developed for integrating stochastic differential equations (SDE's) with Poisson white noise (PWN) and Lévy white noise (LWN). The method for integrating SDE's with PWN has the same structure as that proposed by Kim [Phys. Rev. E 76, 011109 (2007)], but is established by using different arguments. The integration of SDE's with LWN is based on a representation of Lévy processes by sums of scaled Brownian motions and compound Poisson processes. It is shown that the numerical solutions of SDE's with PWN and LWN converge weakly to the exact solutions of these equations, so that they can be used to estimate not only marginal properties but also distributions of functionals of the exact solutions. Numerical examples are used to demonstrate the applications and the accuracy of the proposed integration algorithms.

Realization of Lévy walks as Markovian stochastic processes

Based on multivariate Langevin processes we present a realization of Lévy flights as a continuous process. For the simple case of a particle moving under the influence of friction and a velocity-dependent stochastic force we explicitly derive the generalized Langevin equation and the corresponding generalized Fokker-Planck equation describing Lévy flights. Our procedure is similar to the treatment of the Kramers-Fokker-Planck equation in the Smoluchowski limit. The proposed approach may open a way to treat Lévy flights in inhomogeneous media or systems with boundaries in the future.

Transport in a Lévy ratchet: group velocity and distribution spread

We consider the motion of an overdamped particle in a periodic potential lacking spatial symmetry under the influence of symmetric, white, Lévy noise, being a minimal setup for a "Lévy ratchet." Due to the nonthermal character of the Lévy noise, the particle exhibits a motion with a preferred direction even in the absence of whatever additional time-dependent forces. The examination of the Lévy ratchet has to be based on the characteristics of directionality which are different from typically used measures such as mean current and the dispersion of particle positions, since these become inappropriate when the moments of the noise diverge. To overcome this problem, we discuss robust measures of directionality of transport such as the position of the median of the particle displacement distribution characterizing the group velocity and the interquantile distance giving the measure of the distribution width. Moreover, we analyze the behavior of splitting probabilities for leaving an interval of a given length, unveiling qualitative differences between the noises with Lévy indices below and above unity.

Steady-state Lévy flights in a confined domain

We derive the generalized Fokker-Planck equation associated with a Langevin equation driven by arbitrary additive white noise. We apply our result to study the distribution of symmetric and asymmetric Lévy flights in an infinitely deep potential well. The fractional Fokker-Planck equation for Lévy flights is derived and solved analytically in the steady state. It is shown that Lévy flights are distributed according to the beta distribution, whose probability density becomes singular at the boundaries of the well. The origin of the preferred concentration of flying objects near the boundaries in nonequilibrium systems is clarified.

Leapover lengths and first passage time statistics for Lévy flights

Exact results for the first passage time and leapover statistics of symmetric and one-sided Lévy flights (LFs) are derived. LFs with a stable index alpha are shown to have leapover lengths that are asymptotically power law distributed with an index alpha for one-sided LFs and, surprisingly, with an index alpha/2 for symmetric LFs. The first passage time distribution scales like a power law with an index 1/2 as required by the Sparre-Andersen theorem for symmetric LFs, whereas one-sided LFs have a narrow distribution of first passage times. The exact analytic results are confirmed by extensive simulations.

Stationary states in Langevin dynamics under asymmetric Lévy noises

Properties of systems driven by white non-Gaussian noises can be very different from these of systems driven by the white Gaussian noise. We investigate stationary probability densities for systems driven by alpha-stable Lévy-type noises, which provide natural extension to the Gaussian noise having, however, a new property, namely a possibility of being asymmetric. Stationary probability densities are examined for a particle moving in parabolic, quartic, and in generic double well potential models subjected to the action of alpha-stable noises. Relevant solutions are constructed by methods of stochastic dynamics. In situations where analytical results are known they are compared with numerical results. Furthermore, the problem of estimation of the parameters of stationary densities is investigated.

Dynamic disorder in single-molecule Michaelis-Menten kinetics: The reaction-diffusion formalism in the Wilemski-Fixman approximation

Single-molecule equations for the Michaelis-Menten [Biochem. Z. 49, 333 (1913)] mechanism of enzyme action are analyzed within the Wilemski-Fixman [J. Chem. Phys. 58, 4009 (1973); 60, 866 (1974)] approximation after the effects of dynamic disorder--modeled by the anomalous diffusion of a particle in a harmonic well--are incorporated into the catalytic step of the reaction. The solution of the Michaelis-Menten equations is used to calculate the distribution of waiting times between successive catalytic turnovers in the enzyme beta-galactosidase. The calculated distribution is found to agree qualitatively with experimental results on this enzyme obtained at four different substrate concentrations. The calculations are also consistent with measurements of correlations in the fluctuations of the fluorescent light emitted during the course of catalysis, and with measurements of the concentration dependence of the randomness parameter.

Numerical method for solving stochastic differential equations with Poissonian white shot noise

We propose a numerical integration scheme to solve stochastic differential equations driven by Poissonian white shot noise. Our formula, which is based on an integral equation, which is equivalent to the stochastic differential equation, utilizes a discrete time approximation with fixed integration time step. We show that our integration formula approaches the Euler formula if the Poissonian noise approaches the Gaussian white noise. The accuracy and efficiency of the proposed algorithm are examined by studying the dynamics of an overdamped particle driven by Poissonian white shot noise in a spatially periodic potential. We find that the accuracy of the proposed algorithm only weakly depends on the parameters characterizing the Poissonian white shot noise; this holds true even if the limit of Gaussian white noise is approached.

Barrier crossing driven by Lévy noise: universality and the role of noise intensity

We study the barrier crossing of a particle driven by white symmetric Lévy noise of index alpha and intensity D for three different generic types of potentials: (a) a bistable potential, (b) a metastable potential, and (c) a truncated harmonic potential. For the low noise intensity regime we recover the previously proposed algebraic dependence on D of the characteristic escape time, T_{esc} approximately C(alpha)D;{mu(alpha)} , where C(alpha) is a coefficient. It is shown that the exponent mu(alpha) remains approximately constant, mu approximately 1 for 0<alpha<2 ; at alpha=2 the power-law form of T_{esc} changes into the known exponential dependence on 1D ; it exhibits a divergencelike behavior as alpha approaches 2. In this regime we observe a monotonous increase of the escape time T_{esc} with increasing alpha (keeping the noise intensity D constant). The probability density of the escape time decays exponentially. In addition, for low noise intensities the escape times correspond to barrier crossing by multiple Lévy steps. For high noise intensities, the escape time curves collapse for all values of alpha . At intermediate noise intensities, the escape time exhibits nonmonotonic dependence on the index alpha , while still retaining the exponential form of the escape time density.