Computer simulation of solid C60 using multiple time‐step algorithms

The reversible reference system propagator algorithm (r‐RESPA), based on a Trotter factorization of the classical propagator, is tested in a molecular dynamics simulation of solid C60. We show how, with an appropriate subdivision of the interaction potential and with a careful balancing of the integration parameters, one can adopt large time steps and impressively efficient r‐RESPA integrators which yield the same dynamics obtained by means of the small time‐step Verlet algorithm. The results presented here show that the use of r‐RESPA integrators speeds up the simulation by a factor of between 20–40 with respect to the standard Verlet algorithms.

[1]  Watanabe,et al.  Direct dynamical calculation of entropy and free energy by adiabatic switching. , 1990, Physical review letters.

[2]  Meijer,et al.  New orientationally ordered low-temperature superstructure in high-purity C60. , 1992, Physical review letters.

[3]  B. Chase,et al.  Vibrational spectroscopy of fullerenes (C60 and C70). Temperature dependant studies , 1992 .

[4]  Sven J. Cyvin,et al.  Molecular vibrations of footballene , 1988 .

[5]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[6]  Mark E. Tuckerman,et al.  Molecular dynamics in systems with multiple time scales: Systems with stiff and soft degrees of freedom and with short and long range forces , 1991 .

[7]  S. C. O'brien,et al.  C60: Buckminsterfullerene , 1985, Nature.

[8]  Hans De Raedt,et al.  Applications of the generalized Trotter formula , 1983 .

[9]  Michiel Sprik,et al.  Modeling the orientational ordering transition in solid fullerene (C60) , 1992 .

[10]  W. Krätschmer,et al.  The infrared and ultraviolet absorption spectra of laboratory-produced carbon dust: evidence for the presence of the C60 molecule , 1990 .

[11]  M. Kappes,et al.  The Structure of the C60 Molecule: X-Ray Crystal Structure Determination of a Twin at 110 K , 1991, Science.

[12]  H. Trotter.  On the product of semi-groups of operators , 1959 .

[13]  E. Carter,et al.  Time-Reversible Multiple Time Scale ab Initio Molecular Dynamics , 1993 .

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

[15]  G. Marconi,et al.  A Vibrational Potential Function for Molecular Dynamics Simulation of Buckminsterfullerene , 1993 .

[16]  Masakatsu Watanabe,et al.  Dynamics of molecules with internal degrees of freedom by multiple time-step methods , 1993 .

[17]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[18]  P. Procacci,et al.  Vibrational frequencies of C70 , 1992 .

[19]  P. Procacci,et al.  The far-infrared spectrum of crystalline C[sub 60] , 1993 .

[20]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with disparate masses , 1991 .

[21]  Smith,et al.  Orientational ordering transition in solid C60. , 1991, Physical review letters.

[22]  Z. Gamba.  A first approach to a simple intermolecular potential model for the condensed phases of C60 , 1992 .

[23]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with long range forces , 1991 .

[24]  Ailan Cheng,et al.  Molecular dynamics simulations of solid buckminsterfullerenes , 1991 .