Wavelet Monte Carlo dynamics: a new algorithm for simulating the hydrodynamics of interacting Brownian particles

Abstract

We develop a new algorithm for the Brownian dynamics of soft matter systems that evolves time by spatially correlated Monte Carlo moves. The algorithm uses vector wavelets as its basic moves and produces hydrodynamics in the low Reynolds number regime propagated according to the Oseen tensor. When small moves are removed the correlations closely approximate the Rotne-Prager tensor, itself widely used to correct for deficiencies in Oseen. We also include plane wave moves to provide the longest range correlations, which we detail for both infinite and periodic systems. The computational cost of the algorithm scales competitively with the number of particles simulated, $N$, scaling as $N\ln N$ in homogeneous systems and as $N$ in dilute systems. In comparisons to established lattice Boltzmann and Brownian dynamics algorithms the wavelet method was found to be only a factor of order 1 times more expensive than the cheaper lattice Boltzmann algorithm in marginally semi-dilute simulations, while it is significantly faster than both algorithms at large $N$ in dilute simulations. We also validate the algorithm by checking it reproduces the correct dynamics and equilibrium properties of simple single polymer systems, as well as verifying the effect of periodicity on the mobility tensor.

arXiv comment16 pages, 12 figures, submitted to the Journal of Chemical Physics

References in corpus (9)

• Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part I. Theoretical Foundation
• Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part II. Numerical Results
• Event-chain Monte Carlo algorithms for hard-sphere systems
• Corrections to Scaling in the Hydrodynamic Properties of Dilute Polymer Solutions
• Translational Diffusion of Polymer Chains with Excluded Volume and Hydrodynamic Interactions by Brownian Dynamics Simulation
• Optimisation of a Brownian dynamics algorithm for semidilute polymer solutions
• Implicit and explicit solvent models for the simulation of a single polymer chain in solution: Lattice Boltzmann vs Brownian dynamics
• Monte Carlo Simulation of Dense Polymer Melts Using Event Chain Algorithms
• Multiscale Monte Carlo for simple fluids

Cited by in corpus (3)

• Universal scaling and characterisation of gelation in associative polymer solutions
• Influence of thermal fluctuations on active diffusion at large Péclet numbers
• Multi-scale time-stepping in molecular dynamics