The Use of Processor Groups in Molecular Dynamics Simulations to Sample Free-Energy States
Abstract
Molecular dynamics calculations composed of many independent simulations are frequently encountered in free energy simulations, as well as many other simulation approaches. In principal, the availability of a large number of independent tasks should make possible the development of highly scalable parallel code that executes these tasks concurrently. This paper discusses the use of processor groups to write simulation codes of this type and describes results a code that evaluates the volume dependence of the Helmholtz free energy for clusters of an immiscible fluid in a solvent. The results show that very high levels of scalability can be achieved using processor groups with corresponding reductions in the time to completion. The main limitation to scaling appears to be load imbalance due to variations in the execution times of the individual tasks.
- Authors:
- Publication Date:
- Research Org.:
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 908948
- Report Number(s):
- PNNL-SA-51253
KP1303000; TRN: US200722%%828
- DOE Contract Number:
- AC05-76RL01830
- Resource Type:
- Journal Article
- Resource Relation:
- Journal Name: Journal of Chemical Theory and Computation, 3(2):583-592; Journal Volume: 3; Journal Issue: 2
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; 99 GENERAL AND MISCELLANEOUS//MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; COMPUTER CODES; AVAILABILITY; FREE ENERGY; COMPUTERIZED SIMULATION; MOLECULAR DYNAMICS METHOD
Citation Formats
Palmer, Bruce J., Kathmann, Shawn M., Krishnan, Manoj Kumar, Tipparaju, Vinod, and Nieplocha, Jarek. The Use of Processor Groups in Molecular Dynamics Simulations to Sample Free-Energy States. United States: N. p., 2007.
Web. doi:10.1021/ct600260u.
Palmer, Bruce J., Kathmann, Shawn M., Krishnan, Manoj Kumar, Tipparaju, Vinod, & Nieplocha, Jarek. The Use of Processor Groups in Molecular Dynamics Simulations to Sample Free-Energy States. United States. doi:10.1021/ct600260u.
Palmer, Bruce J., Kathmann, Shawn M., Krishnan, Manoj Kumar, Tipparaju, Vinod, and Nieplocha, Jarek. Sun .
"The Use of Processor Groups in Molecular Dynamics Simulations to Sample Free-Energy States". United States.
doi:10.1021/ct600260u.
@article{osti_908948,
title = {The Use of Processor Groups in Molecular Dynamics Simulations to Sample Free-Energy States},
author = {Palmer, Bruce J. and Kathmann, Shawn M. and Krishnan, Manoj Kumar and Tipparaju, Vinod and Nieplocha, Jarek},
abstractNote = {Molecular dynamics calculations composed of many independent simulations are frequently encountered in free energy simulations, as well as many other simulation approaches. In principal, the availability of a large number of independent tasks should make possible the development of highly scalable parallel code that executes these tasks concurrently. This paper discusses the use of processor groups to write simulation codes of this type and describes results a code that evaluates the volume dependence of the Helmholtz free energy for clusters of an immiscible fluid in a solvent. The results show that very high levels of scalability can be achieved using processor groups with corresponding reductions in the time to completion. The main limitation to scaling appears to be load imbalance due to variations in the execution times of the individual tasks.},
doi = {10.1021/ct600260u},
journal = {Journal of Chemical Theory and Computation, 3(2):583-592},
number = 2,
volume = 3,
place = {United States},
year = {Sun Apr 01 00:00:00 EDT 2007},
month = {Sun Apr 01 00:00:00 EDT 2007}
}
-
Molecular dynamics simulations of united atom Lennard-Jones model for methane, ethane, and propane in TIP3P water have been used to estimate the chemical potentials of aqueous methane, ethane, and propane from 600 to 1200 {degree}C and densities from 0 to 1 g cm{sup -3}. Estimates of the errors in the predictions due to the inadequacy of the models show that this method of prediction gives reasonable accuracy. The diameter ({sigma}) of the water-methane interaction is the most important parameter. The present predictions for methane are compared with a variety of other methods from the literature. An equation with seven adjustablemore »
-
Insights into the structural function of the complex of HIV-1 protease with TMC-126: molecular dynamics simulations and free-energy calculations
The binding properties of the protein-inhibitor complex of human immunodeficiency virus type 1 (HIV-1) protease with the inhibitor TMC-126 are investigated by combining computational alanine scanning (CAS) mutagenesis with binding free-energy decomposition (BFED). The calculated results demonstrate that the flap region (residues 38-58) and the active site region (residues 23-32) in HIV-1 protease contribute 63.72% of the protease to the binding of the inhibitor. In particular, the mechanisms for the interactions of key residues of these species are fully explored and analyzed. Interestingly, the regression analyses show that both CAS and BFED based on the generalized Born model yield similarmore » -
Coordination and Hydrolysis of Plutonium Ions in Aqueous Solution using Car-Parrinello Molecular Dynamics Free Energy Simulations
Car-Parrinello molecular dynamics (CPMD) simulations have been used to examine the hydration structures, coordination energetics and the first hydrolysis constants of Pu3+, Pu4+, PuO2+ and PuO22+ ions in aqueous solution at 300 K. The coordination numbers and structural properties of the first shell of these ions are in good agreement with available experimental estimates. The hexavalent PuO22+ species is coordinated to 5 aquo ligands while the pentavalent PuO2+ complex is coordinated to 4 aquo ligands. The Pu3+ and Pu4+ ions are both coordinated to 8 water molecules. The first hydrolysis constants obtained for Pu3+ and PuO22+ are 6.65 and 5.70more » -
Excess free energy of liquids from molecular dynamics simulations. Application to water models
Thermodynamic integration and perturbation methods have been applied to calculate the excess free energy, /Delta//Alpha//sub e/, of several models of liquid water. The results for the SPC and TIP models agree well with the experimental excess free energy of liquid water computed from the vapor pressure. The precision of the calculation, in which the nonbonded interaction of all molecules is coupled to a forcing parameter, is demonstrated by the agreement of the computed temperature dependence of /Delta//Alpha//sub e/ with that expected. Some results are presented that show the apparent equivalence of the stepwise perturbation and continuous integration algorithms for thismore » -
Understanding the catalytic mechanism of xanthosine methyltransferase in caffeine biosynthesis from QM/MM molecular dynamics and free energy simulations
S-Adenosyl-l-methionine (SAM) dependent xanthosine methyltransferase (XMT) is the key enzyme that catalyzes the first methyl transfer in the caffeine biosynthesis pathway to produce the intermediate 7-methylxanthosine (7mXR). Although XMT has been a subject of extensive discussions, the catalytic mechanism and nature of the substrate involved in the catalysis are still unclear. Here in this paper, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy (potential of mean force or PMF) simulations are undertaken to determine the catalytic mechanism of the XMT-catalyzed reaction. Both xanthosine and its monoanionic form with N3 deprotonated are used as the substrates for the methylation.more »Cited by 2