Sampling rare trajectories using stochastic bridges

Abstract

The numerical quantification of the statistics of rare events in stochastic processes is a challenging computational problem. We present a sampling method that constructs an ensemble of stochastic trajectories that are constrained to have fixed start and end points (so-called stochastic bridges). We then show that by carefully choosing a set of such bridges and assigning an appropriate statistical weight to each bridge, one can focus more processing power on the rare events of a target stochastic process while faithfully preserving the statistics of these rate trajectories. Further, we also compare the stochastic bridges produced using our method to the Wentzel-Kramers-Brillouin (WKB) optimal paths of the target process, derived in the limit of low noise. We see that the paths produced using our method, encoding the full statistics of the process, collapse onto the WKB optimal path as the level of noise is reduced. We propose that our method can be used to judge the accuracy of the WKB approximation at finite levels of noise.

References in corpus (7)

• The large deviation approach to statistical mechanics
• Extinction Rates for Fluctuation-Induced Metastabilities : A Real-Space WKB Approach
• Convex Hulls of Random Walks: Large-Deviation Properties
• Non-intersecting Brownian bridges in the flat-to-flat geometry
• Generating discrete-time constrained random walks and Lévy flights
• Joint distributions of partial and global maxima of a Brownian Bridge
• Generating constrained run-and-tumble trajectories