Dataset Information

Testing for physical validity in molecular simulations.

ABSTRACT: Advances in recent years have made molecular dynamics (MD) and Monte Carlo (MC) simulations powerful tools in molecular-level research, allowing the prediction of experimental observables in the study of systems such as proteins, membranes, and polymeric materials. However, the quality of any prediction based on molecular dynamics results will strongly depend on the validity of underlying physical assumptions. Unphysical behavior of simulations can have significant influence on the results and reproducibility of these simulations, such as folding of proteins and DNA or properties of lipid bilayers determined by cutoff treatment, dynamics of peptides and polymers affected by the choice of thermostat, or liquid properties depending on the simulation time step. Motivated by such examples, we propose a two-fold approach to increase the robustness of molecular simulations. The first part of this approach involves tests which can be performed by the users of MD programs on their respective systems and setups. We present a number of tests of different complexity, ranging from simple post-processing analysis to more involved tests requiring additional simulations. These tests are shown to significantly increase the reliability of MD simulations by catching a number of common simulation errors violating physical assumptions, such as non-conservative integrators, deviations from the Boltzmann ensemble, and lack of ergodicity between degrees of freedom. To make the usage as easy as possible, we have developed an open-source and platform-independent Python library (https://physical-validation.readthedocs.io) implementing these tests. The second part of the approach involves testing for code correctness. While unphysical behavior can be due to poor or incompatible choices of parameters by the user, it can just as well originate in coding errors within the program. We therefore propose to include physical validation tests in the code-checking mechanism of MD software packages. We have implemented such a validation for the GROMACS software package, ensuring that every major release passes a number of physical sanity checks performed on selected representative systems before shipping. It is, to our knowledge, the first major molecular mechanics software package to run such validation routinely. The tests are, as the rest of the package, open source software, and can be adapted for other software packages.

SUBMITTER: Merz PT

PROVIDER: S-EPMC6126824 | biostudies-literature | 2018

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Testing for physical validity in molecular simulations.

Merz Pascal T PT Shirts Michael R MR

PloS one 20180906 9

Advances in recent years have made molecular dynamics (MD) and Monte Carlo (MC) simulations powerful tools in molecular-level research, allowing the prediction of experimental observables in the study of systems such as proteins, membranes, and polymeric materials. However, the quality of any prediction based on molecular dynamics results will strongly depend on the validity of underlying physical assumptions. Unphysical behavior of simulations can have significant influence on the results and r ...[more]

PMID: 30188933

Similar Datasets

Project description:Many different molecular species have been observed in the interstellar medium. These range from simple diatomic species to saturated organic molecules with several carbon atoms. The latter molecules are assumed to be formed predominantely on the surface of interstellar dust grains. All surface reactions that can proceed under the low interstellar temperatures are exothermic. Their exothermicity can be as high as a few electron volts, which is considerable compared to the thermal energy of the molecules at 10 K. It is postulated that this exothermicity can be used for the desorption of reaction products from the grain. In previous studies, we have shown that translational excitation can lead to desorption, whereas vibrational and rotational excitations are much less efficient in the desorption of surface products. Vibrational excitation is however much more likely upon bond formation than translational excitation. The present study follows energy dissipation upon translational, vibrational, or rotational excitation of admolecules on a surface and its conversion, or lack thereof, to different energy contributions. To this end, thousands of molecular dynamics simulations were performed with an admolecule on top of a surface that received a fixed amount of energy, vibrational, rotational, or translational. Three different surface species have been considered, CO2, H2O, and CH4, spanning a range in binding energies, the number of internal degrees of freedom, and molecular weights. A fast exchange of energy between vibrational stretches is observed, but only very limited exchange to rotational or translation excitation has been found. For the dissipation of energy to the surface, excitation of the surface-admolecule bond is critical. Astrochemical models often assume instantaneous equipartition of energy after a reaction process to estimate the amount of available energy for chemical desorption. Based on the present study, we conclude that this assumption is not justified.

Dataset Information

Testing for physical validity in molecular simulations.

Publications

Testing for physical validity in molecular simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets