z-logo
Premium
Rigid‐body dynamics in the isothermal‐isobaric ensemble: A test on the accuracy and computational efficiency
Author(s) -
Shinoda Wataru,
Mikami Masuhiro
Publication year - 2003
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.10249
Subject(s) - molecular dynamics , isobaric process , hamiltonian (control theory) , isothermal process , shake , rigid body , statistical physics , algorithm , bilayer , computer science , chemistry , thermodynamics , computational chemistry , physics , mathematics , classical mechanics , mathematical optimization , quantum mechanics , membrane , biochemistry
We have developed a time‐reversible rigid‐body (rRB) molecular dynamics algorithm in the isothermal‐isobaric (NPT) ensemble. The algorithm is an extension of rigid‐body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non‐Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state‐of‐the‐art constraint‐dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time‐step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS‐MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 920–930, 2003

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here