## Abstract An extended system Hamiltonian is proposed to perform molecular dynamics (MD) simulation in the grand canonical ensemble. The Hamiltonian is similar to the one proposed by Lynch and Pettitt (Lynch and Pettitt, J Chem Phys 1997, 107, 8594), which consists of the kinetic and potential ene
Canonical numerical methods for molecular dynamics simulations
β Scribed by Daniel I. Okunbor; Robert D. Skeel
- Publisher
- John Wiley and Sons
- Year
- 1994
- Tongue
- English
- Weight
- 651 KB
- Volume
- 15
- Category
- Article
- ISSN
- 0192-8651
No coin nor oath required. For personal study only.
β¦ Synopsis
We consider the application of canonical numerical integrators to molecular dynamics simulations. Computer experiments are done to analyze the relative merits of using canonical integrators against their noncanonical counterparts. These experiments involve the study of the properties of liquid argon utilizing the Lennard-Jones interaction potential. To accomplish this comparative study of canonical and noncanonical integrators, we computed thermodynamic and structural quantities. Our results indicate that noncanonical methods, with one exception, fail miserably to conserve energy and as a consequence give poor estimates of the other quantities. Also, higher-order canonical methods may offer an advantage over the Stormeri Verlet method.
π SIMILAR VOLUMES
## Abstract The structure of the general classical trajectory FORTRAN code is discussed for use in molecular dynamics simulations of polymer processes. Substantial reductions in the number of mathematical operations in the code are achieved by using statement functions, modifying the ususal chain r
We describe an adaptive method for achieving load balance in parallel computations simulating phenomena which are distributed over a spatially extended region, but are local in nature. We have tested the method on standard short-ranged parallel molecular dynamics calculations. The performance gain w
Molecular dynamics simulations were performed for an actin monomer without \(\mathrm{Ca}^{2+}\) and nucleotide for a period of \(100 \mathrm{psec}\) in order to study the stable structure and fluctuations of actin monomers in solution. For the initial structure for the molecular dynamics simulations