Algorithm for rigid-body Brownian dynamics.

We present an algorithm for performing rigid-body Brownian dynamics that can take into account the hydrodynamic properties (translational and rotational friction tensors and the coupling between them) of each rigid body. In the zero temperature limit, the error term scales as Delta;{4} for time step Delta , while at nonzero temperatures the error scaling is Delta;{5/2} . We test the algorithm by applying it to a molecule of four-aminopyridine in water. We intend to use the algorithm to model the interaction between biological ion channels and other channel blocker molecules, but it may also have applicability to modeling other small particles such as colloids or nanoparticles.