Physics – Condensed Matter – Soft Condensed Matter
Scientific paper
2008-01-21
J. Chem. Phys. 130 (2009) 114905
Physics
Condensed Matter
Soft Condensed Matter
8 pages, 4 figures; extended version: better approximation + added comparison to Ermak and Fixman HI for bead-spring polymers;
Scientific paper
10.1063/1.3089668
In the Ermak-McCammon algorithm for Brownian Dynamics, the hydrodynamic interactions (HI) between N spherical particles are described by a 3N x 3N diffusion tensor. This tensor has to be factorized at each timestep with a runtime of O(N^3), making the calculation of the correlated random displacements the bottleneck for many-particle simulations. Here we present a faster algorithm for this step, which is based on a truncated expansion of the hydrodynamic multi-particle correlations as two-body contributions. The comparison to the exact algorithm and to the Chebyshev approximation of Fixman verifies that for bead-spring polymers this approximation yields about 95% of the hydrodynamic correlations at an improved runtime scaling of O(N^2) and a reduced memory footprint. The approximation is independent of the actual form of the hydrodynamic tensor and can be applied to arbitrary particle configurations. This now allows to include HI into large many-particle Brownian dynamics simulations, where until now the runtime scaling of the correlated random motion was prohibitive.
Geyer Tihamer
Winter Uwe
No associations
LandOfFree
An O(N^2) Approximation for Hydrodynamic Interactions in Brownian Dynamics Simulations does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with An O(N^2) Approximation for Hydrodynamic Interactions in Brownian Dynamics Simulations, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and An O(N^2) Approximation for Hydrodynamic Interactions in Brownian Dynamics Simulations will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-594334