When evaluating integration schemes used in molecular dynamics (MD) simulations, energy conservation is often cited as the primary criterion by which the integrators should be com- pared. As a result variable stepsize Runge-Kutta methods are often ruled out of consideration due to their characteristic energy drift.
We have shown that by appropriately modifying the stepsize selection strategy in a variable stepsize RK method it is possible for the MD practitioner to obtain substantial control over the energy drift during the course of a simulation. This ability has been previously unreported in the literature, and we present numerical examples to illustrate that it can be achieved without sacrificing computational efficiency under currently obtainable timescales.
Identifer | oai:union.ndltd.org:LACETR/oai:collectionscanada.gc.ca:OTU.1807/25559 |
Date | 31 December 2010 |
Creators | Easley, Kante |
Contributors | Enright, Wayne |
Source Sets | Library and Archives Canada ETDs Repository / Centre d'archives des thèses électroniques de Bibliothèque et Archives Canada |
Language | en_ca |
Detected Language | English |
Type | Thesis |
Page generated in 0.0015 seconds