Random walk approximation for irreversible drift-diffusion process on manifold: ergodicity, unconditional stability and convergence

by   Yuan Gao, et al.

Irreversible drift-diffusion processes are very common in biochemical reactions. They have a non-equilibrium stationary state (invariant measure) which does not satisfy detailed balance. For the corresponding Fokker-Planck equation on a closed manifold, via Voronoi tessellation, we propose two upwind finite volume schemes with or without the information of the invariant measure. Both two schemes enjoy stochastic Q-matrix structures and can be decomposed as a gradient flow part and a Hamiltonian flow part, which enable us to prove unconditional stability, ergodicity and error estimates. Based on two upwind schemes, several numerical examples - including sampling accelerated by a mixture flow, image transformations and simulations for stochastic model of chaotic system - are conducted. These two structure-preserving schemes also give a natural random walk approximation for a generic irreversible drift-diffusion process on a manifold. Thus they can be adapted to manifold-related computations induced from high dimensional molecular dynamics.



There are no comments yet.


page 18

page 19

page 20

page 21

page 23

page 24


Global random walk solvers for fully coupled flow and transport in saturated/unsaturated porous media (extended version)

In this article, we present new random walk methods to solve flow and tr...

Inbetweening auto-animation via Fokker-Planck dynamics and thresholding

We propose an equilibrium-driven deformation algorithm (EDDA) to simulat...

Global random walk solvers for reactive transport and biodegradation processes in heterogeneous porous media

Flow and multicomponent reactive transport in saturated/unsaturated poro...

An Optimal Mass Transport Method for Random Genetic Drift

We propose and analyze an optimal mass transport method for a random gen...

Numerical analysis of several schemes for a ionic liquid model

In this paper, we consider an unipolar degenerated drift-diffusion syste...

Disentangling by Subspace Diffusion

We present a novel nonparametric algorithm for symmetry-based disentangl...

Branching random walk solutions to the Wigner equation

The stochastic solutions to the Wigner equation, which explain the nonlo...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.