6533b85dfe1ef96bd12bdc91

RESEARCH PRODUCT

BROWNIAN DYNAMICS SIMULATIONS WITHOUT GAUSSIAN RANDOM NUMBERS

Wolfgang PaulBurkhard Dünweg

subject

Stochastic processMathematical analysisGeneral Physics and AstronomyStatistical and Nonlinear PhysicsOrnstein–Uhlenbeck processBrownian excursionBrownian bridgeComputer Science Applicationssymbols.namesakeComputational Theory and MathematicsWiener processReflected Brownian motionStochastic simulationsymbolsStatistical physicsGaussian processMathematical PhysicsMathematics

description

We point out that in a Brownian dynamics simulation it is justified to use arbitrary distribution functions of random numbers if the moments exhibit the correct limiting behavior prescribed by the Fokker-Planck equation. Our argument is supported by a simple analytical consideration and some numerical examples: We simulate the Wiener process, the Ornstein-Uhlenbeck process and the diffusion in a Φ4 potential, using both Gaussian and uniform random numbers. In these examples, the rate of convergence of the mean first exit time is found to be nearly identical for both types of random numbers.

https://doi.org/10.1142/s0129183191001037