Brownian Dynamics Simulation of Polydisperse Hard Spheres
  Standard algorithms for the numerical integration of the Langevin equation require that interactions are slowly varying during to the integration timestep. This in not the case for hard-body systems, where there is no clearcut between the correlation time of the noise and the timescale of the interactions. Starting from a short time approximation of the Smoluchowsky equation, we introduce an algorithm for the simulation of the overdamped Brownian dynamics of polydisperse hard-spheres in absence of hydrodynamics interactions and briefly discuss the extension to the case of external drifts. 
Antonio Scala
cond-mat.soft cond-mat.stat-mech physics.comp-ph