ALGORITHMS FOR BROWNIAN DYNAMICS

A third order algorithm for brownian dynamics (BD) simulations is proposed, which is identical to the powerful molecular dynamics (MD) algorithm due to Verlet in the limit of infinitely small friction coefficient γ. In contrast to most BD algorithms used up till now, the integration time step Δt is not limited by the condition Δt ≪ γ-1. It is shown how constraints, such as bondlength or bond angle constraints, can be incorporated in the computational scheme. For the molecule considered here the proposed BD algorithm is about ten times more efficient than the second order ones that are generally used. The introduction of bondlength constraints saves about a factor of three in computing time, as in the MD case. The application of bond angle constraints is not recommended.

[1]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[2]  J. H. Weiner,et al.  A generalization of Kramers’ rate formula to include some anharmonic effects , 1978 .

[3]  W. Vangunsteren.  Constrained dynamics of flexible molecules , 1980 .

[4]  Athanasios Papoulis,et al.  Probability, Random Variables and Stochastic Processes , 1965 .

[5]  E. Helfand,et al.  Brownian dynamics study of transitions in a polymer chain of bistable oscillators , 1978 .

[6]  D. Lévesque,et al.  Transport properties and the time evolution of electrolyte solutions in the Brownian dynamics approximation , 1979 .

[7]  J. H. Weiner,et al.  Computer Simulation of Conformational Transitions in an Idealized Polymer Model , 1977 .

[8]  Thomas A. Weber,et al.  Brownian Dynamics Study of Polymer Conformational Transitions , 1979 .

[9]  Thomas A. Weber,et al.  Simulation of n‐butane using a skeletal alkane model , 1978 .

[10]  J. H. Weiner,et al.  Rate theory for solids. IV. Classical Brownian-motion model , 1974 .

[11]  R. Kubo.  The fluctuation-dissipation theorem , 1966 .

[12]  D. Ermak,et al.  Numerical integration of the Langevin equation: Monte Carlo simulation , 1980 .

[13]  E. Helfand,et al.  Numerical integration of stochastic differential equations — ii , 1979, The Bell System Technical Journal.

[14]  T. Schneider,et al.  Molecular-dynamics study of a three-dimensional one-component model for distortive phase transitions , 1978 .

[15]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[16]  G. Torrie,et al.  Stochastic trajectory simulation of iodine recombination in liquids , 1980 .

[17]  Harold L. Friedman,et al.  Brownian dynamics: Its application to ionic solutions , 1977 .

[18]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[19]  H. Berendsen,et al.  STOCHASTIC DYNAMICS FOR MOLECULES WITH CONSTRAINTS BROWNIAN DYNAMICS OF NORMAL-ALKANES , 1981 .

[20]  Abraham Nitzan,et al.  Dynamics of gas–solid interactions: Calculations of energy transfer and sticking , 1977 .

[21]  G. Gilmer,et al.  Molecular dynamics of surface diffusion. I. The motion of adatoms and clusters , 1979 .

[22]  Jean-Paul Ryckaert,et al.  Molecular dynamics of liquid alkanes , 1978 .

[23]  Irwin Oppenheim,et al.  Molecular Theory of Brownian Motion for Several Particles , 1971 .