MILCH SHAKE: an efficient method for constraint dynamics applied to alkanes.
We describe a method to impose constraints in a molecular dynamics simulation. A technique developed to solve the special case of a linear topology (MILC SHAKE) is hybridized with the SHAKE algorithm. The methodology, which we term MILC-hybridized SHAKE (or MILCH SHAKE), applies to more complex topologies. Here we consider the important case of all atom models of alkanes. Exploiting the mass difference between carbon and hydrogen we show that for higher alkanes MILCH SHAKE can be an order of magnitude faster than SHAKE.
Bailey
A G
AG
Department of Physics, Imperial College London, South Kensington Campus, Prince Consort Road, London SW7 2AZ, United Kingdom. aimeebailey@gmail.com
Lowe
C P
CP