Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userInstantons for rare events in heavy-tailed distributions
129
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Large deviation theory and instanton calculus for stochastic systems are widely used to gain insight into the evolution and probability of rare events. At its core lies the realization that rare events are, under the right circumstances, dominated by their least unlikely realization. Their computation through a saddle-point approximation of the path integral for the corresponding stochastic field theory then reduces an inefficient stochastic sampling problem into a deterministic optimization problem: finding the path of smallest action, the instanton. In the presence of heavy tails, though, standard algorithms to compute the instanton critically fail to converge. The reason for this failure is the divergence of the scaled cumulant generating function (CGF) due to a non-convex large deviation rate function. We propose a solution to this problem by convexifying the rate function through nonlinear reparametrization of the observable, which allows us to compute instantons even in the presence of super-exponential or algebraic tail decay. The approach is generalizable to other situations where the existence of the CGF is required, such as exponential tilting in importance sampling for Monte-Carlo algorithms. We demonstrate the proposed formalism by applying it to rare events in several stochastic systems with heavy tails, including extreme power spikes in fiber optics induced by soliton formation.
rate research
Read More
We study the time until first occurrence, the first-passage time, of rare density fluctuations in diffusive systems. We approach the problem using a model consisting of many independent random walkers on a lattice. The existence of spatial correlations makes this problem analytically intractable. However, for a mean-field approximation in which the walkers can jump anywhere in the system, we obtain a simple asymptotic form for the mean first-passage time to have a given number k of particles at a distinguished site. We show numerically, and argue heuristically, that for large enough k, the mean-field results give a good approximation for first-passage times for systems with nearest-neighbour dynamics, especially for two and higher spatial dimensions. Finally, we show how the results change when density fluctuations anywhere in the system, rather than at a specific distinguished site, are considered.
Rare events in stochastic processes with heavy-tailed distributions are controlled by the big jump principle, which states that a rare large fluctuation is produced by a single event and not by an accumulation of coherent small deviations. The principle has been rigorously proved for sums of independent and identically distributed random variables and it has recently been extended to more complex stochastic processes involving Levy distributions, such as Levy walks and the Levy-Lorentz gas, using an effective rate approach. We review the general rate formalism and we extend its applicability to continuous time random walks and to the Lorentz gas, both with stretched exponential distributions, further enlarging its applicability. We derive an analytic form for the probability density functions for rare events in the two models, which clarify specific properties of stretched exponentials.
Starting from a master equation, we derive the evolution equation for the size distribution of elements in an evolving system, where each element can grow, divide into two, and produce new elements. We then probe general solutions of the evolution quation, to obtain such skew distributions as power-law, log-normal, and Weibull distributions, depending on the growth or division and production. Specifically, repeated production of elements of uniform size leads to power-law distributions, whereas production of elements with the size distributed according to the current distribution as well as no production of new elements results in log-normal distributions. Finally, division into two, or binary fission, bears Weibull distributions. Numerical simulations are also carried out, confirming the validity of the obtained solutions.
In a microcanonical ensemble (constant $NVE$, hard reflecting walls) and in a molecular dynamics ensemble (constant $NVEmathbf{PG}$, periodic boundary conditions) with a number $N$ of smooth elastic hard spheres in a $d$-dimensional volume $V$ having a total energy $E$, a total momentum $mathbf{P}$, and an overall center of mass position $mathbf{G}$, the individual velocity components, velocity moduli, and energies have transformed beta distributions with different arguments and shape parameters depending on $d$, $N$, $E$, the boundary conditions, and possible symmetries in the initial conditions. This can be shown marginalizing the joint distribution of individual energies, which is a symmetric Dirichlet distribution. In the thermodynamic limit the beta distributions converge to gamma distributions with different arguments and shape or scale parameters, corresponding respectively to the Gaussian, i.e., Maxwell-Boltzmann, Maxwell, and Boltzmann or Boltzmann-Gibbs distribution. These analytical results agree with molecular dynamics and Monte Carlo simulations with different numbers of hard disks or spheres and hard reflecting walls or periodic boundary conditions. The agreement is perfect with our Monte Carlo algorithm, which acts only on velocities independently of positions with the collision versor sampled uniformly on a unit half sphere in $d$ dimensions, while slight deviations appear with our molecular dynamics simulations for the smallest values of $N$.
The granular Leidenfrost effect (B. Meerson et al, Phys. Rev. Lett. {bf 91}, 024301 (2003), P. Eshuis et al, Phys. Rev. Lett. {bf 95}, 258001 (2005)) is the levitation of a mass of granular matter when a wall below the grains is vibrated giving rise to a hot granular gas below the cluster. We find by simulation that for a range of parameters the system is bistable: the levitated cluster can occasionally break and give rise to two clusters and a hot granular gas above and below. We use techniques from the theory of rare events to compute the mean transition time for breaking to occur. This requires the introduction of a two-component reaction coordinate.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
