Algorithms for Brownian dynamics

We derive a class of efficient and stable algorithms of Brownian dynamics using a formula, derived by Suzuki, to express time-ordered operators. These algorithms are simpler than those derived by Helfand from Runge-Kutta algorithms and, like Helfand algorithms, can be combined with SHAKE to describe the Brownian dynamics of constrained systems.