Scalable Methods for Computing Sharp Extreme Event Probabilities in Infinite-Dimensional Stochastic Systems

Abstract

We introduce and compare computational techniques for sharp extreme event probability estimates in stochastic differential equations with small additive Gaussian noise. In particular, we focus on strategies that are scalable, i.e. their efficiency does not degrade upon temporal and possibly spatial refinement. For that purpose, we extend algorithms based on the Laplace method for estimating the probability of an extreme event to infinite dimensional path space. The method estimates the limiting exponential scaling using a single realization of the random variable, the large deviation minimizer. Finding this minimizer amounts to solving an optimization problem governed by a differential equation. The probability estimate becomes sharp when it additionally includes prefactor information, which necessitates computing the determinant of a second derivative operator to evaluate a Gaussian integral around the minimizer. We present an approach in infinite dimensions based on Fredholm determinants, and develop numerical algorithms to compute these determinants efficiently for the high-dimensional systems that arise upon discretization. We also give an interpretation of this approach using Gaussian process covariances and transition tubes. An example model problem, for which we provide an open-source python implementation, is used throughout the paper to illustrate all methods discussed. To study the performance of the methods, we consider examples of stochastic differential and stochastic partial differential equations, including the randomly forced incompressible three-dimensional Navier-Stokes equations.

References in corpus (28)

• PyCUDA and PyOpenCL: A Scripting-Based Approach to GPU Run-Time Code Generation
• On the Numerical Evaluation of Fredholm Determinants
• Computation of extreme heat waves in climate models using a large deviation algorithm
• Instantons and Intermittency
• Functional determinants by contour integration methods
• A Scaling Theory of Bifurcations in the Symmetric Weak-Noise Escape Problem
• Rogue Waves and Large Deviations in Deep Sea
• Numerical computation of rare events via large deviation theory
• The instanton method and its numerical implementation in fluid mechanics
• Asymptotics of work distributions: The pre-exponential factor
• Instanton filtering for the stochastic Burgers equation
• Dynamical Landscape and Multistability of a Climate Model
• An Eyring-Kramers law for the stochastic Allen-Cahn equation in dimension two
• A Large Deviation Theory-based Analysis of Heat Waves and Cold Spells in a Simplified Model of the General Circulation of the Atmosphere
• Instanton based importance sampling for rare events in stochastic PDEs
• Instantons and the path to intermittency in turbulent flows
• Spontaneous Symmetry Breaking for Extreme Vorticity and Strain in the 3D Navier-Stokes Equations
• Adaptive high-order splitting schemes for large-scale differential Riccati equations
• Observing symmetry-broken optimal paths of stationary Kardar-Parisi-Zhang interface via a large-deviation sampling of directed polymers in random media
• A Koopman framework for rare event simulation in stochastic differential equations
• Symmetries and zero modes in sample path large deviations
• Gel'fand-Yaglom type equations for calculating fluctuations around Instantons in stochastic systems
• Path integral derivation and numerical computation of large deviation prefactors for non-equilibrium dynamics through matrix Riccati equations
• Instantons for rare events in heavy-tailed distributions
• Efficient Computation of Instantons for Multi-Dimensional Turbulent Flows with Large Scale Forcing
• Scalable Methods for Computing Sharp Extreme Event Probabilities in Infinite-Dimensional Stochastic Systems
• Eddy-Viscous Modeling and the Topology of Extreme Circulation Events in Three-Dimensional Turbulence
• Normalized Gaussian Path Integrals

Cited by in corpus (5)

• Evaluation of transition rates from nonequilibrium instantons
• Scalable Methods for Computing Sharp Extreme Event Probabilities in Infinite-Dimensional Stochastic Systems
• Short-time large deviations of the spatially averaged height of a KPZ interface on a ring
• Instantons, fluctuations and singularities in the supercritical stochastic nonlinear Schroedinger equation
• Synthetic Turbulence via an Instanton Gas Approximation