Dataset Information

Molecular Dynamics Simulations of Energy Dissipation on Amorphous Solid Water: Testing the Validity of Equipartition.

ABSTRACT: 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, CO₂, H₂O, and CH₄, 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.

SUBMITTER: Fredon A

PROVIDER: S-EPMC8381352 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description: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.

Project description:The Internal Coordinate Molecular Dynamics (ICMD) method is an attractive molecular dynamics (MD) method for studying the dynamics of bonded systems such as proteins and polymers. It offers a simple venue for coarsening the dynamics model of a system at multiple hierarchical levels. For example, large scale protein dynamics can be studied using torsional dynamics, where large domains or helical structures can be treated as rigid bodies and the loops connecting them as flexible torsions. ICMD with such a dynamic model of the protein, combined with enhanced conformational sampling method such as temperature replica exchange, allows the sampling of large scale domain motion involving high energy barrier transitions. Once these large scale conformational transitions are sampled, all-torsion, or even all-atom, MD simulations can be carried out for the low energy conformations sampled via coarse grained ICMD to calculate the energetics of distinct conformations. Such hierarchical MD simulations can be carried out with standard all-atom forcefields without the need for compromising on the accuracy of the forces. Using constraints to treat bond lengths and bond angles as rigid can, however, distort the potential energy landscape of the system and reduce the number of dihedral transitions as well as conformational sampling. We present here a two-part solution to overcome such distortions of the potential energy landscape with ICMD models. To alleviate the intrinsic distortion that stems from the reduced phase space in torsional MD, we use the Fixman compensating potential. To additionally alleviate the extrinsic distortion that arises from the coupling between the dihedral angles and bond angles within a force field, we propose a hybrid ICMD method that allows the selective relaxing of bond angles. This hybrid ICMD method bridges the gap between all-atom MD and torsional MD. We demonstrate with examples that these methods together offer a solution to eliminate the potential energy distortions encountered in constrained ICMD simulations of peptide molecules.

Dataset Information

Molecular Dynamics Simulations of Energy Dissipation on Amorphous Solid Water: Testing the Validity of Equipartition.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets