Scalability of the second-order reliability method for stochastic differential equations with multiplicative noise

TL;DR

This paper develops a scalable numerical method using automatic differentiation to efficiently compute sharp probability estimates for extreme events in high-dimensional stochastic differential equations with multiplicative noise, based on second-order reliability methods.

Contribution

It introduces a scalable, high-dimensional numerical approach for second-order reliability estimates in SDEs, utilizing automatic differentiation and operator trace computations.

Findings

Efficient computation of tail probabilities in high-dimensional SDEs.

Demonstration of the method on complex SPDE examples.

Public availability of the JAX implementation.

Abstract

We show how to efficiently compute asymptotically sharp estimates of extreme event probabilities in stochastic differential equations (SDEs) with small multiplicative Brownian noise. The underlying approximation is known as sharp large deviation theory or precise Laplace asymptotics in mathematics, the second-order reliability method (SORM) in reliability engineering, and the instanton or optimal fluctuation method with 1-loop corrections in physics. It is based on approximating the tail probability in question with the most probable realization of the stochastic process, and local perturbations around this realization. We first recall and contextualize the relevant classical theoretical result on precise Laplace asymptotics of diffusion processes [Ben Arous (1988), Stochastics, 25(3), 125-153], and then show how to compute the involved infinite-dimensional quantities - operator traces and Carleman-Fredholm determinants - numerically in a way that is scalable with respect to the time discretization and remains feasible in high spatial dimensions. Using tools from automatic differentiation, we achieve a straightforward black-box numerical computation of the SORM estimates in JAX. The method is illustrated in examples of SDEs and stochastic partial differential equations, including a two-dimensional random advection-diffusion model of a passive scalar. We thereby demonstrate that it is possible to obtain efficient and accurate SORM estimates for very high-dimensional problems, as long as the infinite-dimensional structure of the problem is correctly taken into account. Our JAX implementation of the method is made publicly available.