Comparison of Brownian dynamics algorithms with hydrodynamic interaction

The hydrodynamic interaction is an essential effect to consider in Brownian dynamics simulations of polymer and nanoparticle dilute solutions. Several mathematical approaches can be used to build Brownian dynamics algorithms with hydrodynamic interaction, the most common of them being the exact but time demanding Cholesky decomposition and the Chebyshev polynomial expansion. Recently, Geyer and Winter J. Chem. Phys. 130, 1149051 (2009)] have proposed a new approximation to treat the hydrodynamic interaction that seems quite efficient and is increasingly used. So far, a systematic comparison among those approaches has not been clearly made. In this paper, several features and the efficiency of typical implementations of those approaches are evaluated by using bead-and-spring chain models. The different sensitivity to the bead overlap detected for the different implementations may be of interest to select the suitable algorithm for a given simulation.