The effect of molecular dynamics sampling on the calculated observable gas-phase structures

2016 ◽  
Vol 18 (27) ◽  
pp. 18237-18245 ◽  
Author(s):  
Denis S. Tikhonov ◽  
Arseniy A. Otlyotov ◽  
Vladimir V. Rybkin
Keyword(s):  
Molecular Dynamics ◽  
Ab Initio ◽  
Gas Phase ◽  
Sampling Methods ◽  
Quantum Effects ◽  
Phase Structures ◽  
Nuclear Quantum Effects

We evaluate the performance of various ab initio molecular dynamics sampling methods for the calculation of observable gas-phase structures and probe the nuclear quantum effects.

scholarly journals Classical Ab Initio Molecular Dynamics Run at an Elevated Temperature is Not a Good Model for the Nuclear Quantum Effects in Water at Ambient Temperature

2021 ◽  
Author(s):  
Chenghan Li ◽  
Francesco Paesani ◽  
Gregory A. Voth
Keyword(s):  
Molecular Dynamics ◽  
Elevated Temperature ◽  
Ambient Temperature ◽  
Ab Initio ◽  
Density Functional ◽  
Quantum Effects ◽  
Ab Initio Molecular Dynamics ◽  
Higher Temperature ◽  
Three Body ◽  
Nuclear Quantum Effects

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the over-structuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic ensemble are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. We find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-body correlations as well as dynamics. In many of these cases, the effects of NQEs are exactly the opposite of the effects of carrying out the simulations at an elevated temperature.

scholarly journals Sensitivity of Intra- and Intermolecular Interactions of Benzo[h]quinoline from Car–Parrinello Molecular Dynamics and Electronic Structure Inspection

2021 ◽  
Vol 22 (10) ◽  
pp. 5220
Author(s):  
Jarosław J. Panek ◽  
Joanna Zasada ◽  
Bartłomiej M. Szyja ◽  
Beata Kizior ◽  
Aneta Jezierska
Keyword(s):  
Molecular Dynamics ◽  
Electronic Structure ◽  
Hydrogen Bonds ◽  
Density Functional ◽  
Quantum Effects ◽  
Potential Of Mean Force ◽  
Probability Density Distribution ◽  
Donor And Acceptor ◽  
Nuclear Quantum Effects ◽  
Vibrational Features

The O-H...N and O-H...O hydrogen bonds were investigated in 10-hydroxybenzo[h]quinoline (HBQ) and benzo[h]quinoline-2-methylresorcinol complex in vacuo, solvent and crystalline phases. The chosen systems contain analogous donor and acceptor moieties but differently coupled (intra- versus intermolecularly). Car–Parrinello molecular dynamics (CPMD) was employed to shed light onto principle components of interactions responsible for the self-assembly. It was applied to study the dynamics of the hydrogen bonds and vibrational features as well as to provide initial geometries for incorporation of quantum effects and electronic structure studies. The vibrational features were revealed using Fourier transformation of the autocorrelation function of atomic velocity and by inclusion of nuclear quantum effects on the O-H stretching solving vibrational Schrödinger equation a posteriori. The potential of mean force (Pmf) was computed for the whole trajectory to derive the probability density distribution and for the O-H stretching mode from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling and nuclear quantum effects. The electronic structure changes of the benzo[h]quinoline-2-methylresorcinol dimer and trimers were studied based on Constrained Density Functional Theory (CDFT) whereas the Electron Localization Function (ELF) method was applied for all systems. It was found that the bridged proton is localized on the donor side in both investigated systems in vacuo. The crystalline phase simulations indicated bridged proton-sharing and transfer events in HBQ. These effects are even more pronounced when nuclear quantization is taken into account, and the quantized Pmf allows the proton to sample the acceptor area more efficiently. The CDFT indicated the charge depletion at the bridged proton for the analyzed dimer and trimers in solvent. The ELF analysis showed the presence of the isolated proton (a signature of the strongest hydrogen bonds) only in some parts of the HBQ crystal simulation. The collected data underline the importance of the intramolecular coupling between the donor and acceptor moieties.

scholarly journals On-the-Fly Ring-Polymer Molecular Dynamics Calculations of the Dissociative Photodetachment Process of the Oxalate Anion

Molecules ◽  
2021 ◽  
Vol 26 (23) ◽  
pp. 7250
Author(s):  
Yukinobu Takahashi ◽  
Yu Hashimoto ◽  
Kohei Saito ◽  
Toshiyuki Takayanagi
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Quantum Effects ◽  
Experimental Spectrum ◽  
Computational Time ◽  
Oxalate Anion ◽  
Dissociative Photodetachment ◽  
Ring Polymer ◽  
Ring Polymer Molecular Dynamics ◽  
Nuclear Quantum Effects

The dissociative photodetachment dynamics of the oxalate anion, C2O4H− + hν → CO2 + HOCO + e−, were theoretically studied using the on-the-fly path-integral and ring-polymer molecular dynamics methods, which can account for nuclear quantum effects at the density-functional theory level in order to compare with the recent experimental study using photoelectron–photofragment coincidence spectroscopy. To reduce computational time, the force acting on each bead of ring-polymer was approximately calculated from the first and second derivatives of the potential energy at the centroid position of the nuclei beads. We find that the calculated photoelectron spectrum qualitatively reproduces the experimental spectrum and that nuclear quantum effects are playing a role in determining spectral widths. The calculated coincidence spectrum is found to reasonably reproduce the experimental spectrum, indicating that a relatively large energy is partitioned into the relative kinetic energy between the CO2 and HOCO fragments. This is because photodetachment of the parent anion leads to Franck–Condon transition to the repulsive region of the neutral potential energy surface. We also find that the dissociation dynamics are slightly different between the two isomers of the C2O4H− anion with closed- and open-form structures.

scholarly journals Compaction of RNA Hairpins and Their Kissing Complexes in Native Electrospray Mass Spectrometry

2020 ◽  
Author(s):  
Joséphine Abi-Ghanem ◽  
Clémence Rabin ◽  
Massimiliano Porrini ◽  
Frédéric Rosu ◽  
Valerie Gabelica
Keyword(s):  
Mass Spectrometry ◽  
Molecular Dynamics ◽  
Gas Phase ◽  
Electrospray Mass Spectrometry ◽  
Initial Structure ◽  
Native Solution ◽  
Phase Structures ◽  
Rna Hairpins ◽  
Level Calculation ◽  
Structure In Solution

When electrosprayed from native solution conditions, RNA hairpins and kissing complexes acquire charge states at which they get significantly more compact in the gas phase than their initial structure in solution. Here we show the limits of using force field molecular dynamics to interpret the gas-phase structures of nucleic acid complexes in the gas phase, and we suggest that higher-level calculation levels should be used in the future.<br>

Sign in / Sign up

Export Citation Format

Share Document