Probability flow solution of the Fokker-Planck equation

06/09/2022
by   Nicholas M. Boffi, et al.
33

The method of choice for integrating the time-dependent Fokker-Planck equation in high-dimension is to generate samples from the solution via integration of the associated stochastic differential equation. Here, we introduce an alternative scheme based on integrating an ordinary differential equation that describes the flow of probability. Unlike the stochastic dynamics, this equation deterministically pushes samples from the initial density onto samples from the solution at any later time. The method has the advantage of giving direct access to quantities that are challenging to estimate only given samples from the solution, such as the probability current, the density itself, and its entropy. The probability flow equation depends on the gradient of the logarithm of the solution (its "score"), and so is a-priori unknown. To resolve this dependence, we model the score with a deep neural network that is learned on-the-fly by propagating a set of particles according to the instantaneous probability current. Our approach is based on recent advances in score-based diffusion for generative modeling, with the important difference that the training procedure is self-contained and does not require samples from the target density to be available beforehand. To demonstrate the validity of the approach, we consider several examples from the physics of interacting particle systems; we find that the method scales well to high-dimensional systems, and accurately matches available analytical solutions and moments computed via Monte-Carlo.

READ FULL TEXT

page 12

page 14

page 29

page 30

page 31

page 33

page 34

page 35

research
11/10/2022

Artificial Neural Network Solver for Time-Dependent Fokker-Planck Equations

Stochastic differential equations play an important role in various appl...
research
03/20/2021

Adaptive deep density approximation for Fokker-Planck equations

In this paper we present a novel adaptive deep density approximation str...
research
06/26/2019

Piecewise polynomial approximation of probability density functions with application to uncertainty quantification for stochastic PDEs

The probability density function (PDF) associated with a given set of sa...
research
06/06/2023

Fisher information and shape-morphing modes for solving the Fokker-Planck equation in higher dimensions

The Fokker-Planck equation describes the evolution of the probability de...
research
12/06/2022

Proposal of a Score Based Approach to Sampling Using Monte Carlo Estimation of Score and Oracle Access to Target Density

Score based approaches to sampling have shown much success as a generati...
research
11/07/2022

Posterior samples of source galaxies in strong gravitational lenses with score-based priors

Inferring accurate posteriors for high-dimensional representations of th...
research
09/10/2021

Physics-based machine learning for modeling stochastic IP3-dependent calcium dynamics

We present a machine learning method for model reduction which incorpora...

Please sign up or login with your details

Forgot password? Click here to reset