-
Early warning signs of critical transitions -- The $α$-stable case
Authors:
Lucia S. Layritz,
Ilya Pavlyukevich,
Anja Rammig,
Christian Kuehn
Abstract:
Statistical early warning signs can be used to identify an approaching bifurcation in stochastic dynamical systems and are now regularly employed in applications concerned with the identification of potential rapid, non-linear change or tipping points. However, the reliability of these early warning signs relies on a number of key mathematical assumptions, most notably the presence of Gaussian noi…
▽ More
Statistical early warning signs can be used to identify an approaching bifurcation in stochastic dynamical systems and are now regularly employed in applications concerned with the identification of potential rapid, non-linear change or tipping points. However, the reliability of these early warning signs relies on a number of key mathematical assumptions, most notably the presence of Gaussian noise. We here show that for systems driven by non-Gaussian, $α$-stable noise, the classical early warning signs of rising variance and autocorrelation are not supported by mathematical theory and their use poses the danger of spurious, false-positive results. To address this, we provide a generalized approach by introduce the scaling factor $γ_X$ as an alternative early warning sign. We show that in the case of the Ornstein-Uhlenbeck process, there exists a direct inverse relationship between $γ_{X}$ and the bifurcation parameter, telling us that $γ_{X}$ will increase as we approach the bifurcation. Our numerical simulations confirm theoretical results and show that our findings generalize well to non-linear, non-equilibrium systems. We thus provide a generalized, robust and applicable statistical early warning sign for systems driven by Gaussian and non-Gaussian $α$-stable noise.
△ Less
Submitted 27 November, 2023;
originally announced November 2023.
-
Stochastic energy-balance model with a moving ice line
Authors:
Ilya Pavlyukevich,
Marian Ritsch
Abstract:
In [SIAM J. Appl. Dyn. Sys., 12(4):2068--2092, 2013], Widiasih proposed and analyzed a deterministic one-dimensional Budyko-Sellers energy-balance model with a moving ice-line. In this paper, we extend this model to the stochastic setting and analyze it within the framework of stochastic slow-fast systems. We derive the dynamics for the ice line in the limit of a small parameter as a solution to a…
▽ More
In [SIAM J. Appl. Dyn. Sys., 12(4):2068--2092, 2013], Widiasih proposed and analyzed a deterministic one-dimensional Budyko-Sellers energy-balance model with a moving ice-line. In this paper, we extend this model to the stochastic setting and analyze it within the framework of stochastic slow-fast systems. We derive the dynamics for the ice line in the limit of a small parameter as a solution to a stochastic differential equation. The stochastic approach enables the study of co-existing (metastable) climate states as well as the transition dynamics between them.
△ Less
Submitted 26 November, 2023;
originally announced November 2023.
-
Walsh's Brownian Motion and Donsker Scaling Limits of Perturbed Random Walks
Authors:
Ilya Pavlyukevich,
Andrey Pilipenko
Abstract:
In this paper, we study the Donsker scaling limit of integer-valued random walks perturbed on a finite subset of $\mathbb Z$ called a membrane. Under very mild assumptions about the law of the random walk's increments inside and outside of the membrane we show weak convergence of the scaled processes to a skew Brownian motion and provide the explicit formula for its permeability parameter in terms…
▽ More
In this paper, we study the Donsker scaling limit of integer-valued random walks perturbed on a finite subset of $\mathbb Z$ called a membrane. Under very mild assumptions about the law of the random walk's increments inside and outside of the membrane we show weak convergence of the scaled processes to a skew Brownian motion and provide the explicit formula for its permeability parameter in terms of stationary distributions of certain embedded Markov chains. The proof is based on a representation of the original random walk as a multidimensional coordinate process and its convergence to a Walsh Brownian motion.
△ Less
Submitted 16 October, 2023;
originally announced October 2023.
-
Stochastic selection problem for a Stratonovich SDE with power non-linearity
Authors:
Ilya Pavlyukevich,
Georgiy Shevchenko
Abstract:
In our paper [Bernoulli 26(2), 2020, 1381-1409], we found all strong Markov solutions that spend zero time at $0$ of the Stratonovich stochastic differential equation $d X=|X|^α\circ dB$, $α\in (0,1)$. These solutions have the form $X_t^θ=F(B^θ_t)$, where $F(x)=\frac{1}{1-α}|x|^{1/(1-α)}\text{sign}\, x$ and $B^θ$ is the skew Brownian motion with skewness parameter $θ\in [-1,1]$ starting at…
▽ More
In our paper [Bernoulli 26(2), 2020, 1381-1409], we found all strong Markov solutions that spend zero time at $0$ of the Stratonovich stochastic differential equation $d X=|X|^α\circ dB$, $α\in (0,1)$. These solutions have the form $X_t^θ=F(B^θ_t)$, where $F(x)=\frac{1}{1-α}|x|^{1/(1-α)}\text{sign}\, x$ and $B^θ$ is the skew Brownian motion with skewness parameter $θ\in [-1,1]$ starting at $F^{-1}(X_0)$. In this paper we show how an addition of small external additive noise $\varepsilon W$ restores uniqueness. In the limit as $\varepsilon\to 0$, we recover heterogeneous diffusion corresponding to the physically symmetric case $θ=0$.
△ Less
Submitted 5 May, 2024; v1 submitted 12 August, 2023;
originally announced August 2023.
-
Homogenization of a multivariate diffusion with semipermeable interfaces
Authors:
Olga Aryasova,
Ilya Pavlyukevich,
Andrey Pilipenko
Abstract:
We study the homogenization problem for a system of stochastic differential equation with local time terms that models a multivariate diffusion in presence of semipermeable hyperplane interfaces with oblique penetration. We show that this system has a unique weak solution and determine its weak limit as the distances between the interfaces converge to zero. In the limit, the singular local times t…
▽ More
We study the homogenization problem for a system of stochastic differential equation with local time terms that models a multivariate diffusion in presence of semipermeable hyperplane interfaces with oblique penetration. We show that this system has a unique weak solution and determine its weak limit as the distances between the interfaces converge to zero. In the limit, the singular local times terms vanish and give rise to an additional regular interface-induced drift.
△ Less
Submitted 1 February, 2024; v1 submitted 5 March, 2023;
originally announced March 2023.
-
First Order Linear Marcus SPDEs
Authors:
Lena-Susanne Hartmann,
Ilya Pavlyukevich
Abstract:
In this paper we solve a Lévy driven linear stochastic first order partial differential equation (transport equation) understood in the canonical (Marcus) form. The solution can be obtained with the help of the method of stochastic characteristics. It has the same form as a solution of a deterministic PDE or a solution of a stochastic PDE driven by a Brownian motion studied by Kunita (1984, 1997).
In this paper we solve a Lévy driven linear stochastic first order partial differential equation (transport equation) understood in the canonical (Marcus) form. The solution can be obtained with the help of the method of stochastic characteristics. It has the same form as a solution of a deterministic PDE or a solution of a stochastic PDE driven by a Brownian motion studied by Kunita (1984, 1997).
△ Less
Submitted 1 March, 2023;
originally announced March 2023.
-
Cutoff ergodicity bounds in Wasserstein distance for a viscous energy shell model with Lévy noise
Authors:
Gerardo Barrera,
Michael A. Högele,
Juan Carlos Pardo,
Ilya Pavlyukevich
Abstract:
This article establishes explicit non-asymptotic ergodic bounds in the renormalized Wasserstein-Kantorovich-Rubinstein (WKR) distance for a viscous energy shell lattice model of turbulence with random energy injection. The system under consideration is driven either by a Brownian motion, a symmetric $α$-stable Lévy process, a stationary Gaussian or $α$-stable Ornstein-Uhlenbeck process, or by a ge…
▽ More
This article establishes explicit non-asymptotic ergodic bounds in the renormalized Wasserstein-Kantorovich-Rubinstein (WKR) distance for a viscous energy shell lattice model of turbulence with random energy injection. The system under consideration is driven either by a Brownian motion, a symmetric $α$-stable Lévy process, a stationary Gaussian or $α$-stable Ornstein-Uhlenbeck process, or by a general Lévy process with second moments. The obtained non-asymptotic bounds establish asymptotically abrupt thermalization. The analysis is based on the explicit representation of the solution of the system in terms of convolutions of Bessel functions.
△ Less
Submitted 31 January, 2024; v1 submitted 27 February, 2023;
originally announced February 2023.
-
Limit behaviour of random walks on $\mathbb Z^m$ with two-sided membrane
Authors:
V. Bogdanskii,
I. Pavlyukevich,
A. Pilipenko
Abstract:
We study Markov chains on $\mathbb Z^m$, $m\geq 2$, that behave like a standard symmetric random walk outside of the hyperplane (membrane) $H=\{0\}\times \mathbb Z^{m-1}$. The transition probabilities on the membrane $H$ are periodic and also depend on the incoming direction to $H$, what makes the membrane $H$ two-sided. Moreover, sliding along the membrane is allowed. We show that the natural sca…
▽ More
We study Markov chains on $\mathbb Z^m$, $m\geq 2$, that behave like a standard symmetric random walk outside of the hyperplane (membrane) $H=\{0\}\times \mathbb Z^{m-1}$. The transition probabilities on the membrane $H$ are periodic and also depend on the incoming direction to $H$, what makes the membrane $H$ two-sided. Moreover, sliding along the membrane is allowed. We show that the natural scaling limit of such Markov chains is a $m$-dimensional diffusion whose first coordinate is a skew Brownian motion and the other $m-1$ coordinates is a Brownian motion with a singular drift controlled by the local time of the first coordinate at $0$. In the proof we utilize a martingale characterization of the Walsh Brownian motion and determine the effective permeability and slide direction. Eventually, a similar convergence theorem is established for the one-sided membrane without slides and random iid transition probabilities.
△ Less
Submitted 4 August, 2021;
originally announced August 2021.
-
Geodesic random walks, diffusion processes and Brownian motion on Finsler manifolds
Authors:
Tianyu Ma,
Vladimir S. Matveev,
Ilya Pavlyukevich
Abstract:
We show that geodesic random walks on a complete Finsler manifold of bounded geometry converge to a diffusion process which is, up to a drift, the Brownian motion corresponding to a Riemannian metric.
We show that geodesic random walks on a complete Finsler manifold of bounded geometry converge to a diffusion process which is, up to a drift, the Brownian motion corresponding to a Riemannian metric.
△ Less
Submitted 16 February, 2021;
originally announced February 2021.
-
Moment bounds for dissipative semimartingales with heavy jumps
Authors:
Alexei Kulik,
Ilya Pavlyukevich
Abstract:
In this paper we show that if large jumps of an Itô-semimartingale $X$ have a finite $p$-moment, $p>0$, the radial part of its drift is dominated by $-|X|^κ$ for some $κ\geq -1$, and the balance condition $p+κ>1$ holds true, then under some further natural technical assumptions $\sup_{t\geq 0} \mathbf{E} |X_t|^{p_X}<\infty$ for each $p_X\in(0,p+κ-1)$. The upper bound $p+κ-1$ is generically optimal…
▽ More
In this paper we show that if large jumps of an Itô-semimartingale $X$ have a finite $p$-moment, $p>0$, the radial part of its drift is dominated by $-|X|^κ$ for some $κ\geq -1$, and the balance condition $p+κ>1$ holds true, then under some further natural technical assumptions $\sup_{t\geq 0} \mathbf{E} |X_t|^{p_X}<\infty$ for each $p_X\in(0,p+κ-1)$. The upper bound $p+κ-1$ is generically optimal. The proof is based on the extension of the method of Lyapunov functions to the semimartingale framework. The uniform moment estimates obtained in this paper are indispensable for the analysis of ergodic properties of Lévy driven stochastic differential equations and Lévy driven multi-scale systems.
△ Less
Submitted 13 June, 2021; v1 submitted 26 April, 2020;
originally announced April 2020.
-
Generalized selection problem with Lévy noise
Authors:
Ilya Pavlyukevich,
Andrey Pilipenko
Abstract:
Let $A_\pm>0$, $β\in(0,1)$, and let $Z^{(α)}$ be a strictly $α$-stable Lévy process with the jump measure $ν(\mathrm{d} z)=(C_+\mathbb{I}_{(0,\infty)}(z)+ C_-\mathbb{I}_{(-\infty,0)}(z))|z|^{-1-α}\,\mathrm{d} z$, $α\in (1,2)$, $C_\pm\geq 0$, $C_++C_->0$. The selection problem for the model stochastic differential equation…
▽ More
Let $A_\pm>0$, $β\in(0,1)$, and let $Z^{(α)}$ be a strictly $α$-stable Lévy process with the jump measure $ν(\mathrm{d} z)=(C_+\mathbb{I}_{(0,\infty)}(z)+ C_-\mathbb{I}_{(-\infty,0)}(z))|z|^{-1-α}\,\mathrm{d} z$, $α\in (1,2)$, $C_\pm\geq 0$, $C_++C_->0$. The selection problem for the model stochastic differential equation $\mathrm{d} \bar X^\varepsilon=(A_+\mathbb{I}_{[0,\infty)}(\bar X^\varepsilon) - A_-\mathbb{I}_{(-\infty,0)}(\bar X^\varepsilon))|\bar X^\varepsilon|^β\,\mathrm{d} t +\varepsilon \mathrm{d} Z^{(α)}$ states that in the small noise limit $\varepsilon\to 0$, solutions $\bar X^\varepsilon$ converge weakly to the maximal or minimal solutions of the limiting non-Lipschitzian ordinary differential equation $\mathrm{d} \bar x=(A_+\mathbb{I}_{[0,\infty)}(\bar x)- A_-\mathbb{I}_{(\infty,0)}(\bar x))|\bar x|^β\,\mathrm{d} t$ with probabilities $\bar p_\pm=\bar p_\pm(α,C_+/C_-,β, A_+/A_-)$, see [Pilipenko and Proske, Stat. Probab. Lett., 132:62-73, 2018]. In this paper we solve the generalized selection problem for the stochastic differential equation $\mathrm{d} X^\varepsilon=a(X^\varepsilon)\,\mathrm{d} t+\varepsilon b(X^\varepsilon)\,\mathrm{d} Z$ whose dynamics in the vicinity of the origin in certain sense reminds of dynamics of the model equation. In particular we show that solutions $X^\varepsilon$ also converge to the maximal or minimal solutions of the limiting irregular ordinary differential equation $\mathrm{d} x=a(x) \,\mathrm{d} t$ with the same model selection probabilities $\bar p_\pm$. This means that for a large class of irregular stochastic differential equations, the selection dynamics is completely determined by four local parameters of the drift and the jump measure.
△ Less
Submitted 11 April, 2020;
originally announced April 2020.
-
First order convergence of weak Wong--Zakai approximations of Lévy driven Marcus SDEs
Authors:
Tetyana Kosenkova,
Alexei Kulik,
Ilya Pavlyukevich
Abstract:
For solutions $X=(X_t)_{t\in[0,T]}$ of Lévy-driven Marcus stochastic differential equations we study the Wong--Zakai type time discrete approximations $\bar X=(\bar X_{kh})_{0\leq k\leq T/h}$, $h>0$, and establish the first order convergence $|E f(X_T)-E f(X^h_T)|\leq C h$ for $f\in C_b^4$.
For solutions $X=(X_t)_{t\in[0,T]}$ of Lévy-driven Marcus stochastic differential equations we study the Wong--Zakai type time discrete approximations $\bar X=(\bar X_{kh})_{0\leq k\leq T/h}$, $h>0$, and establish the first order convergence $|E f(X_T)-E f(X^h_T)|\leq C h$ for $f\in C_b^4$.
△ Less
Submitted 7 February, 2020;
originally announced February 2020.
-
Drift Estimation for a Lévy-Driven Ornstein-Uhlenbeck Process with Heavy Tails
Authors:
Alexander Gushchin,
Ilya Pavlyukevich,
Marian Ritsch
Abstract:
We consider the problem of estimation of the drift parameter of an ergodic Ornstein--Uhlenbeck type process driven by a Lévy process with heavy tails. The process is observed continuously on a long time interval $[0,T]$, $T\to\infty$. We prove that the statistical model is locally asymptotic mixed normal and the maximum likelihood estimator is asymptotically efficient.
We consider the problem of estimation of the drift parameter of an ergodic Ornstein--Uhlenbeck type process driven by a Lévy process with heavy tails. The process is observed continuously on a long time interval $[0,T]$, $T\to\infty$. We prove that the statistical model is locally asymptotic mixed normal and the maximum likelihood estimator is asymptotically efficient.
△ Less
Submitted 25 November, 2019;
originally announced November 2019.
-
First passage properties of asymmetric Lévy flights
Authors:
A. Padash,
A. V. Chechkin,
B. Dybiec,
I. Pavlyukevich,
B. Shokri,
R. Metzler
Abstract:
Lévy Flights are paradigmatic generalised random walk processes, in which the independent stationary increments---the "jump lengths"---are drawn from an $α$-stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of Lévy Flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probab…
▽ More
Lévy Flights are paradigmatic generalised random walk processes, in which the independent stationary increments---the "jump lengths"---are drawn from an $α$-stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of Lévy Flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probability to revisit previous points of visitation, rendering Lévy Flights efficient search processes in one and two dimensions. To further quantify their precise property as random search strategies we here study the first-passage time properties of Lévy Flights in one-dimensional semi-infinite and bounded domains for symmetric and asymmetric jump length distributions. To obtain the full probability density function of first-passage times for these cases we employ two complementary methods. One approach is based on the space-fractional diffusion equation for the probability density function, from which the survival probability is obtained for different values of the stable index $α$ and the skewness (asymmetry) parameter $β$. The other approach is based on the stochastic Langevin equation with $α$-stable driving noise. Both methods have their advantages and disadvantages for explicit calculations and numerical evaluation, and the complementary approach involving both methods will be profitable for concrete applications. We also make use of the Skorokhod theorem for processes with independent increments and demonstrate that the numerical results are in good agreement with the analytical expressions for the probability density function of the first-passage times.
△ Less
Submitted 16 October, 2019;
originally announced October 2019.
-
The first passage problem for stable linear delay equations perturbed by power law Lévy noise
Authors:
Michael A. Högele,
Ilya Pavlyukevich
Abstract:
This article studies a linear scalar delay differential equation subject to small multiplicative power tail Lévy noise. We solve the first passage (the Kramers) problem with probabilistic methods and discover an asymptotic loss of memory in this non-Markovian system. Furthermore, the mean exit time increases as the power of the small noise amplitude, whereas the pre-factor accounts for memory effe…
▽ More
This article studies a linear scalar delay differential equation subject to small multiplicative power tail Lévy noise. We solve the first passage (the Kramers) problem with probabilistic methods and discover an asymptotic loss of memory in this non-Markovian system. Furthermore, the mean exit time increases as the power of the small noise amplitude, whereas the pre-factor accounts for memory effects. In particular, we discover a non-linear delay-induced exit acceleration due to a non-normal growth phenomenon. Our results are illustrated in the well-known linear delay oscillator driven by $α$-stable Lévy flights.
△ Less
Submitted 26 March, 2019;
originally announced March 2019.
-
Stratonovich SDE with irregular coefficients: Girsanov's example revisited
Authors:
Ilya Pavlyukevich,
Georgiy Shevchenko
Abstract:
In this paper we study the Stratonovich stochastic differential equation $\mathrm{d} X=|X|^α\circ\mathrm{d} B$, $α\in(-1,1)$, which has been introduced by Cherstvy et al. [New Journal of Physics 15:083039 (2013)] in the context of analysis of anomalous diffusions in heterogeneous media. We determine its weak and strong solutions, which are homogeneous strong Markov processes \chng{spending zero ti…
▽ More
In this paper we study the Stratonovich stochastic differential equation $\mathrm{d} X=|X|^α\circ\mathrm{d} B$, $α\in(-1,1)$, which has been introduced by Cherstvy et al. [New Journal of Physics 15:083039 (2013)] in the context of analysis of anomalous diffusions in heterogeneous media. We determine its weak and strong solutions, which are homogeneous strong Markov processes \chng{spending zero time at $0$: for $α\in (0,1)$, these solutions have the form $$ X_t^θ=\bigl((1-α)B_t^θ\bigr)^{1/(1-α)}, $$ where $B^θ$ is the $θ$-skew Brownian motion driven by $B$ and starting at $\frac{1}{1-α}(X_0)^{1-α}$, $θ\in [-1,1]$,} and $(x)^γ=|x|^γ\operatorname{sign} x$; for $α\in(-1,0]$, only the case $θ=0$ is possible. The central part of the paper consists in the proof of the existence of a quadratic covariation $[f(B^θ),B]$ for a locally square integrable function $f$ and is based on the time-reversion technique for Markovian diffusions.
△ Less
Submitted 30 September, 2019; v1 submitted 13 December, 2018;
originally announced December 2018.
-
Non-Gaussian Limit Theorem for Non-Linear Langevin Equations Driven by Lévy Noise
Authors:
Alexei Kulik,
Ilya Pavlyukevich
Abstract:
In this paper, we study the small noise behaviour of solutions of a non-linear second order Langevin equation $\ddot x^\varepsilon_t +|\dot x^\varepsilon_t|^β=\dot Z^\varepsilon_{\varepsilon t}$, $β\in\mathbb R$, driven by symmetric non-Gaussian Lévy processes $Z^\varepsilon$. This equation describes the dynamics of a one-degree-of-freedom mechanical system subject to non-linear friction and noisy…
▽ More
In this paper, we study the small noise behaviour of solutions of a non-linear second order Langevin equation $\ddot x^\varepsilon_t +|\dot x^\varepsilon_t|^β=\dot Z^\varepsilon_{\varepsilon t}$, $β\in\mathbb R$, driven by symmetric non-Gaussian Lévy processes $Z^\varepsilon$. This equation describes the dynamics of a one-degree-of-freedom mechanical system subject to non-linear friction and noisy vibrations. For a compound Poisson noise, the process $x^\varepsilon$ on the macroscopic time scale $t/\varepsilon$ has a natural interpretation as a non-linear filter which responds to each single jump of the driving process. We prove that a system driven by a general symmetric Lévy noise exhibits essentially the same asymptotic behaviour under the principal condition $α+2β<4$, where $α\in [0,2]$ is the ``uniform'' Blumenthal--Getoor index of the family $\{Z^\varepsilon\}_{\varepsilon>0}$.
△ Less
Submitted 20 July, 2018; v1 submitted 6 July, 2017;
originally announced July 2017.
-
Random Perturbations of a Periodically Driven Nonlinear Oscillator: Escape from a resonance zone
Authors:
Nishanth Lingala,
Navaratnam Sri Namachchivaya,
Ilya Pavlyukevich
Abstract:
The phase space for a periodically driven nonlinear oscillator consists of many resonance zones. Let the strength of periodic excitation and the strength of the damping be indexed by a small parameter $\varepsilon$. It is well known that, as $\varepsilon \to 0$, the measure of the set of initial conditions which lead to 'capture in a resonance zone' goes to zero. In this paper we study the effect…
▽ More
The phase space for a periodically driven nonlinear oscillator consists of many resonance zones. Let the strength of periodic excitation and the strength of the damping be indexed by a small parameter $\varepsilon$. It is well known that, as $\varepsilon \to 0$, the measure of the set of initial conditions which lead to 'capture in a resonance zone' goes to zero. In this paper we study the effect of weak noise on the escape from a resonance zone.
△ Less
Submitted 16 February, 2017; v1 submitted 29 October, 2015;
originally announced October 2015.
-
Spectral Analysis of a Discrete Metastable System Driven by Lévy Flights
Authors:
Toralf Burghoff,
Ilya Pavlyukevich
Abstract:
In this paper we consider a finite state time discrete Markov chain that mimics the behaviour of solutions of the stochastic differential equation $dX=-U'(X)dt+εdL$, where $U$ is a multi-well potential with $n\geq 2$ local minima and L is a symmetric α-stable Lévy process (Lévy flights process). We investigate the spectrum of the generator of this Markov chain in the limit $ε\to 0$ and localize th…
▽ More
In this paper we consider a finite state time discrete Markov chain that mimics the behaviour of solutions of the stochastic differential equation $dX=-U'(X)dt+εdL$, where $U$ is a multi-well potential with $n\geq 2$ local minima and L is a symmetric α-stable Lévy process (Lévy flights process). We investigate the spectrum of the generator of this Markov chain in the limit $ε\to 0$ and localize the top n eigenvalues $λ^ε_1,\dots, λ^ε_n$. These eigenvalues turn out to be of the same algebraic order $O(ε^α)$ and are well separated from the rest of the spectrum by a spectral gap. We also determine the limits $\lim_{ε\to 0}ε^{-α} λ^ε_i$, $1\leq i\leq n$, and show that the corresponding eigenvectors are approximately constant over the domains which correspond to the potential wells of $U$.
△ Less
Submitted 14 January, 2015;
originally announced January 2015.
-
Marcus versus Stratonovich for Systems with Jump Noise
Authors:
Alexei Chechkin,
Ilya Pavlyukevich
Abstract:
The famous Itô-Stratonovich dilemma arises when one examines a dynamical system with a multiplicative white noise. In physics literature, this dilemma is often resolved in favour of the Stratonovich prescription because of its two characteristic properties valid for systems driven by Brownian motion: (i) it allows physicists to treat stochastic integrals in the same way as conventional integrals,…
▽ More
The famous Itô-Stratonovich dilemma arises when one examines a dynamical system with a multiplicative white noise. In physics literature, this dilemma is often resolved in favour of the Stratonovich prescription because of its two characteristic properties valid for systems driven by Brownian motion: (i) it allows physicists to treat stochastic integrals in the same way as conventional integrals, and (ii) it appears naturally as a result of a small correlation time limit procedure. On the other hand, the Marcus prescription [IEEE Trans. Inform. Theory 24, 164 (1978); Stochastics 4, 223 (1981)] should be used to retain (i) and (ii) for systems driven by a Poisson process, Lévy flights or more general jump processes. In present communication we present an in-depth comparison of the Itô, Stratonovich, and Marcus equations for systems with multiplicative jump noise. By the examples of areal-valued linear system and a complex oscillator with noisy frequency (the Kubo-Anderson oscillator) we compare solutions obtained with the three prescriptions.
△ Less
Submitted 31 May, 2014;
originally announced June 2014.
-
Metastability of Morse-Smale dynamical systems perturbed by heavy-tailed Lévy type noise
Authors:
Michael Högele,
Ilya Pavlyukevich
Abstract:
We consider a general class of finite dimensional deterministic dynamical systems with finitely many local attractors $K^i$ each of which supports a unique ergodic probability measure $P^i$, which includes in particular the class of Morse-Smale systems in any finite dimension. The dynamical system is perturbed by a multiplicative non-Gaussian heavy-tailed Lévy type noise of small amplitude…
▽ More
We consider a general class of finite dimensional deterministic dynamical systems with finitely many local attractors $K^i$ each of which supports a unique ergodic probability measure $P^i$, which includes in particular the class of Morse-Smale systems in any finite dimension. The dynamical system is perturbed by a multiplicative non-Gaussian heavy-tailed Lévy type noise of small amplitude $\varepsilon>0$. Specifically we consider perturbations leading to a Itô, Stratonovich and canonical (Marcus) stochastic differential equation. The respective asymptotic first exit time and location problem from each of the domains of attractions $D^i$ in case of inward pointing vector fields in the limit of $\varepsilon \to 0$ was solved by the authors in [J. Stoch. An. Appl. 32(1), 163-190, 2014]. We extend these results to domains with characteristic boundaries and show that the perturbed system exhibits a metastable behavior in the sense that there exits a unique $\varepsilon$-dependent time scale on which the random system converges to a continuous time Markov chain switching between the invariant measures $P^i$. As examples we consider $α$-stable perturbations of the Duffing equation and a chemical system exhibiting a birhythmic behavior.
△ Less
Submitted 21 May, 2014;
originally announced May 2014.
-
Non-standard Skorokhod convergence of Levy-driven convolution integrals in Hilbert spaces
Authors:
Ilya Pavlyukevich,
Markus Riedle
Abstract:
We study the convergence in probability in the non-standard $M_1$ Skorokhod topology of the Hilbert valued stochastic convolution integrals of the type $\int_0^t F_γ(t-s)\,d L(s)$ to a process $\int_0^t F(t-s)\, d L(s)$ driven by a Lévy process $L$. In Banach spaces we introduce strong, weak and product modes of $M_1$-convergence, prove a criterion for the $M_1$-convergence in probability of stoch…
▽ More
We study the convergence in probability in the non-standard $M_1$ Skorokhod topology of the Hilbert valued stochastic convolution integrals of the type $\int_0^t F_γ(t-s)\,d L(s)$ to a process $\int_0^t F(t-s)\, d L(s)$ driven by a Lévy process $L$. In Banach spaces we introduce strong, weak and product modes of $M_1$-convergence, prove a criterion for the $M_1$-convergence in probability of stochastically continuous càdlàg processes in terms of the convergence in probability of the finite dimensional marginals and a good behaviour of the corresponding oscillation functions, and establish criteria for the convergence in probability of Lévy driven stochastic convolutions. The theory is applied to the infinitely dimensional integrated Ornstein--Uhlenbeck processes with diagonalisable generators.
△ Less
Submitted 19 August, 2014; v1 submitted 6 November, 2013;
originally announced November 2013.
-
The exit problem from the neighborhood of a global attractor for heavy-tailed Lévy diffusions
Authors:
Michael Högele,
Ilya Pavlyukevich
Abstract:
We consider a finite dimensional deterministic dynamical system with a global attractor A with a unique ergodic measure P concentrated on it, which is uniformly parametrized by the mean of the trajectories in a bounded set D containing A. We perturbe this dynamical system by a multiplicative heavy tailed Lévy noise of small intensity ε>0 and solve the asymptotic first exit time and location proble…
▽ More
We consider a finite dimensional deterministic dynamical system with a global attractor A with a unique ergodic measure P concentrated on it, which is uniformly parametrized by the mean of the trajectories in a bounded set D containing A. We perturbe this dynamical system by a multiplicative heavy tailed Lévy noise of small intensity ε>0 and solve the asymptotic first exit time and location problem from a bounded domain D around the attractor A in the limit of ε to 0. In contrast to the case of Gaussian perturbations, the exit time has the asymptotically algebraic exit rate as a function of ε, just as in the case when A is a stable fixed point. In the small noise limit, we determine the joint law of the first time and the exit location on the complement of D. As an example, we study the first exit problem from a neighbourhood of a stable limit cycle for the Van der Pol oscillator perturbed by multiplicative α-stable Lévy noise.
△ Less
Submitted 20 March, 2013;
originally announced March 2013.
-
Small noise asymptotics and first passage times of integrated Ornstein-Uhlenbeck processes driven by $α$-stable Lévy processes
Authors:
Robert Hintze,
Ilya Pavlyukevich
Abstract:
In this paper, we study the asymptotic behaviour of one-dimensional integrated Ornstein-Uhlenbeck processes driven by $α$-stable Lévy processes of small amplitude. We prove that the integrated Ornstein-Uhlenbeck process converges weakly to the underlying $α$-stable Lévy process in the Skorokhod $M_1$-topology which secures the weak convergence of first passage times. This result follows from a mor…
▽ More
In this paper, we study the asymptotic behaviour of one-dimensional integrated Ornstein-Uhlenbeck processes driven by $α$-stable Lévy processes of small amplitude. We prove that the integrated Ornstein-Uhlenbeck process converges weakly to the underlying $α$-stable Lévy process in the Skorokhod $M_1$-topology which secures the weak convergence of first passage times. This result follows from a more general result about approximations of an arbitrary Lévy process by continuous integrated Ornstein-Uhlenbeck processes in the $M_1$-topology.
△ Less
Submitted 5 February, 2014; v1 submitted 28 May, 2012;
originally announced May 2012.
-
First exit times of solutions of stochastic differential equations driven by multiplicative Levy noise with heavy tails
Authors:
Ilya Pavlyukevich
Abstract:
In this paper we study first exit times from a bounded domain of a gradient dynamical system $\dot Y_t=-\nabla U(Y_t)$ 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.
In this paper we study first exit times from a bounded domain of a gradient dynamical system $\dot Y_t=-\nabla U(Y_t)$ 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.
△ Less
Submitted 28 May, 2012;
originally announced May 2012.
-
Lévy Ratchet in a Weak Noise Limit: Theory and Simulation
Authors:
Ilya Pavlyukevich,
Bartlomiej Dybiec,
Aleksei V. Chechkin,
Igor M. Sokolov
Abstract:
We study the motion of a particle embedded in a time independent periodic potential with broken mirror symmetry and subjected to a Lévy noise possessing Lévy stable probability law (Lévy ratchet). We develop analytical approach to the problem based on the asymptotic probabilistic method of decomposition proposed by P. Imkeller and I. Pavlyukevich [J. Phys. A {\bf39}, L237 (2006); Stoch. Proc. Appl…
▽ More
We study the motion of a particle embedded in a time independent periodic potential with broken mirror symmetry and subjected to a Lévy noise possessing Lévy stable probability law (Lévy ratchet). We develop analytical approach to the problem based on the asymptotic probabilistic method of decomposition proposed by P. Imkeller and I. Pavlyukevich [J. Phys. A {\bf39}, L237 (2006); Stoch. Proc. Appl. {\bf116}, 611 (2006)]. We derive analytical expressions for the quantities characterizing the particle motion, namely the splitting probabilities of first escape from a single well, the transition probabilities and the particle current. A particular attention is devoted to the interplay between the asymmetry of the ratchet potential and the asymmetry (skewness) of the Lévy noise. Intensive numerical simulations demonstrate a good agreement with the analytical predictions for sufficiently small intensities of the Lévy noise driving the particle.
△ Less
Submitted 25 August, 2010;
originally announced August 2010.
-
Limit theorems for p-variations of solutions of SDEs driven by additive non-Gaussian stable Levy noise
Authors:
C. Hein,
P. Imkeller,
I. Pavlyukevich
Abstract:
In this paper we study the asymptotic properties of the power variations of stochastic processes of the type X=Y+L, where L is an alpha-stable Levy process, and Y a perturbation which satisfies some mild Lipschitz continuity assumptions. We establish local functional limit theorems for the power variation processes of X. In case X is a solution of a stochastic differential equation driven by L,…
▽ More
In this paper we study the asymptotic properties of the power variations of stochastic processes of the type X=Y+L, where L is an alpha-stable Levy process, and Y a perturbation which satisfies some mild Lipschitz continuity assumptions. We establish local functional limit theorems for the power variation processes of X. In case X is a solution of a stochastic differential equation driven by L, these limit theorems provide estimators of the stability index alpha. They are applicable for instance to model fitting problems for paleo-climatic temperature time series taken from the Greenland ice core.
△ Less
Submitted 23 November, 2008;
originally announced November 2008.
-
One-dimensional space-discrete transport subject to Levy perturbations
Authors:
I. Pavlyukevich,
I. M. Sokolov
Abstract:
In this paper we study a one-dimensional space-discrete transport equation subject to additive Levy forcing. The explicit form of the solutions allows their analytic study. In particular we discuss the invariance of the covariance structure of the stationary distribution for Levy perturbations with finite second moment. The situation of more general Levy perturbations lacking the second moment i…
▽ More
In this paper we study a one-dimensional space-discrete transport equation subject to additive Levy forcing. The explicit form of the solutions allows their analytic study. In particular we discuss the invariance of the covariance structure of the stationary distribution for Levy perturbations with finite second moment. The situation of more general Levy perturbations lacking the second moment is considered as well. We moreover show that some of the properties of the solutions are pertinent to a discrete system and are not reproduced by its continuous analogue.
△ Less
Submitted 14 April, 2008;
originally announced April 2008.
-
First exit times for Lévy-driven diffusions with exponentially light jumps
Authors:
Peter Imkeller,
Ilya Pavlyukevich,
Torsten Wetzel
Abstract:
We consider a dynamical system described by the differential equation $\dot{Y}_t=-U'(Y_t)$ with a unique stable point at the origin. We perturb the system by the Lévy noise of intensity $\varepsilon$ to obtain the stochastic differential equation $dX^{\varepsilon}_t=-U'(X^{\varepsilon}_{t-}) dt+\varepsilon dL_t.$ The process $L$ is a symmetric Lévy process whose jump measure $ν$ has exponentiall…
▽ More
We consider a dynamical system described by the differential equation $\dot{Y}_t=-U'(Y_t)$ with a unique stable point at the origin. We perturb the system by the Lévy noise of intensity $\varepsilon$ to obtain the stochastic differential equation $dX^{\varepsilon}_t=-U'(X^{\varepsilon}_{t-}) dt+\varepsilon dL_t.$ The process $L$ is a symmetric Lévy process whose jump measure $ν$ has exponentially light tails, $ν([u,\infty))\sim\exp(-u^α)$, $α>0$, $u\to \infty$. We study the first exit problem for the trajectories of the solutions of the stochastic differential equation from the interval $(-1,1)$. In the small noise limit $\varepsilon\to0$, the law of the first exit time $σ_x$, $x\in(-1,1)$, has exponential tail and the mean value exhibiting an intriguing phase transition at the critical index $α=1$, namely, $\ln\mathbf{E}σ\sim\varepsilon^{-α}$ for $0<α<1$, whereas $\ln\mathbf{E}σ\sim\varepsilon^{- 1}|\ln\varepsilon|^{1-{1}/α}$ for $α>1$.
△ Less
Submitted 10 June, 2009; v1 submitted 6 November, 2007;
originally announced November 2007.
-
Levy Flights, Non-local Search and Simulated Annealing
Authors:
I. Pavlyukevich
Abstract:
We solve a problem of non-convex stochastic optimisation with help of simulated annealing of Levy flights of a variable stability index. The search of the ground state of an unknown potential is non-local due to big jumps of the Levy flights process. The convergence to the ground state is fast due to a polynomial decrease rate of the temperature.
We solve a problem of non-convex stochastic optimisation with help of simulated annealing of Levy flights of a variable stability index. The search of the ground state of an unknown potential is non-local due to big jumps of the Levy flights process. The convergence to the ground state is fast due to a polynomial decrease rate of the temperature.
△ Less
Submitted 26 January, 2007;
originally announced January 2007.
-
Cooling down Levy flights
Authors:
I. Pavlyukevich
Abstract:
Let L(t) be a Levy flights process with a stability index α\in(0,2), and U be an external multi-well potential. A jump-diffusion Z satisfying a stochastic differential equation dZ(t)=-U'(Z(t-))dt+σ(t)dL(t) describes an evolution of a Levy particle of an `instant temperature' σ(t) in an external force field. The temperature is supposed to decrease polynomially fast, i.e. σ(t)\approx t^{-θ} for so…
▽ More
Let L(t) be a Levy flights process with a stability index α\in(0,2), and U be an external multi-well potential. A jump-diffusion Z satisfying a stochastic differential equation dZ(t)=-U'(Z(t-))dt+σ(t)dL(t) describes an evolution of a Levy particle of an `instant temperature' σ(t) in an external force field. The temperature is supposed to decrease polynomially fast, i.e. σ(t)\approx t^{-θ} for some θ>0. We discover two different cooling regimes. If θ<1/α(slow cooling), the jump diffusion Z(t) has a non-trivial limiting distribution as t\to \infty, which is concentrated at the potential's local minima. If θ>1/α(fast cooling) the Levy particle gets trapped in one of the potential wells.
△ Less
Submitted 26 January, 2007;
originally announced January 2007.
-
Metastable Behaviour of Small Noise Levy-Driven Diffusions
Authors:
Peter Imkeller,
Ilya Pavlyukevich
Abstract:
We consider a dynamical system in R driven by a vector field -U', where U is a multi-well potential satisfying some regularity conditions. We perturb this dynamical system by a Levy noise of small intensity and such that the heaviest tail of its Levy measure is regularly varying. We show that the perturbed dynamical system exhibits metastable behaviour i.e. on a proper time scale it reminds of a…
▽ More
We consider a dynamical system in R driven by a vector field -U', where U is a multi-well potential satisfying some regularity conditions. We perturb this dynamical system by a Levy noise of small intensity and such that the heaviest tail of its Levy measure is regularly varying. We show that the perturbed dynamical system exhibits metastable behaviour i.e. on a proper time scale it reminds of a Markov jump process taking values in the local minima of the potential U. Due to the heavy-tail nature of the random perturbation, the results differ strongly from the well studied purely Gaussian case.
△ Less
Submitted 20 July, 2006; v1 submitted 31 January, 2006;
originally announced January 2006.
-
First exit times of solutions of non-linear stochastic differential equations driven by symmetric Levy processes with alpha-stable components
Authors:
Peter Imkeller,
Ilya Pavlyukevich
Abstract:
We study the exit problem of solutions of the stochastic differential equation dX(t)=-U'(X(t))dt+epsilon dL(t) from bounded or unbounded intervals which contain the unique asymptotically stable critical point of the deterministic dynamical system dY=-U'(Y) dt. The process L is composed of a standard Brownian motion and a symmetric alpha-stable Levy process. Using probabilistic estimates we show…
▽ More
We study the exit problem of solutions of the stochastic differential equation dX(t)=-U'(X(t))dt+epsilon dL(t) from bounded or unbounded intervals which contain the unique asymptotically stable critical point of the deterministic dynamical system dY=-U'(Y) dt. The process L is composed of a standard Brownian motion and a symmetric alpha-stable Levy process. Using probabilistic estimates we show that in the small noise limit epsilon->0, the exit time of X from an interval is an exponentially distributed random variable and determine its expected value. Due to the heavy-tail nature of the alpha-stable component of L, the results differ strongly from the well known case in which the deterministic dynamical system undergoes purely Gaussian perturbations.
△ Less
Submitted 31 January, 2006; v1 submitted 15 September, 2004;
originally announced September 2004.
-
Stochastic Resonance in Two-State Markov Chains
Authors:
Peter Imkeller,
Ilya Pavlyukevich
Abstract:
In this paper we introduce a model which provides a new approach to the phenomenon of stochastic resonance. It is based on the study of the properties of the stationary distribution of the underlying stochastic process. We derive the formula for the spectral power amplification coefficient, study its asymptotic properties and dependence on parameters.
In this paper we introduce a model which provides a new approach to the phenomenon of stochastic resonance. It is based on the study of the properties of the stationary distribution of the underlying stochastic process. We derive the formula for the spectral power amplification coefficient, study its asymptotic properties and dependence on parameters.
△ Less
Submitted 28 February, 2003;
originally announced February 2003.
