Event-driven Langevin simulations of hard spheres
Abstract
The blossoming of interest in colloids and nanoparticles has given renewed impulse to the study of hard-body systems. In particular, hard spheres have become a real test system for theories and experiments. It is therefore necessary to study the complex dynamics of such systems in presence of a solvent; disregarding hydrodynamic interactions, the simplest model is the Langevin equation. Unfortunately, standard algorithms for the numerical integration of the Langevin equation require that interactions are slowly varying during an integration time step. This is not the case for hard-body systems, where there is no clear-cut distinction between the correlation time of the noise and the time scale of the interactions. Starting first from a splitting of the Fokker-Plank operator associated with the Langevin dynamics, and then from an approximation of the two-body Green's function, we introduce and test two algorithms for the simulation of the Langevin dynamics of hard spheres.
- Publication:
-
Physical Review E
- Pub Date:
- August 2012
- DOI:
- 10.1103/PhysRevE.86.026709
- arXiv:
- arXiv:1111.6197
- Bibcode:
- 2012PhRvE..86b6709S
- Keywords:
-
- 05.10.Gg;
- 05.40.Jc;
- Stochastic analysis methods;
- Brownian motion;
- Physics - Computational Physics;
- Condensed Matter - Soft Condensed Matter;
- Condensed Matter - Statistical Mechanics;
- 82C05;
- 82C80;
- I.6.1;
- G.3
- E-Print:
- Phys. Rev. E vol 86 pag 026709 year 2012