Molecular dynamics with constrained nuclear electronic orbital density functional theory: Accurate vibrational spectra from efficient incorporation of nuclear quantum effects

Nuclear quantum effects play a crucial role in many chemical and biological systems involving hydrogen atoms yet are difficult to include in practical molecular simulations. In this Letter, we combine our recently developed methods of constrained nuclear-electronic orbital density functional theory (cNEO-DFT) and constrained minimized energy surface molecular dynamics (CMES-MD) to create a new method for accurately and efficiently describing nuclear quantum effects in molecular simulations. Using this new method, dubbed cNEO-MD, the vibrational spectra of a set of small molecules are calculated and compared with those from conventional ab initio molecular dynamics (AIMD) as well as from experiments. With the same formal scaling, cNEO-MD greatly outperforms AIMD in describing the vibrational modes with significant hydrogen motion characters, demonstrating the promise of cNEO-MD for simulating chemical and biological systems with significant nuclear quantum effects.

Liquid-Liquid Phase Transition in Supercooled H2O and D2O: A Path-Integral Computer Simulation Study

We perform path-integral molecular dynamics (PIMD) and classical MD simulations of H2O and D2O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density ρ(T), isothermal compressibility κT(T), and self-diffusion coefficients D(T) of H2O and D2O are in excellent agreement with available experimental data; the isobaric heat capacity CP(T) obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H2O and D2O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H2O and D2O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H2O, from PIMD simulations, is located at Pc = 167±9 MPa, Tc = 159±6 K, and ρc = 1.02±0.01 g/cm3. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D2O is estimated to be Pc = 176 ± 4 MPa, Tc = 177 ± 2 K, and ρc = 1.13±0.01 g/cm3. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, Pc = 203 ± 4 MPa, Tc = 175 ± 2 K, and ρc = 1.03 ± 0.01 g/cm3). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of Tc for D2O and, particularly, H2O suggest that improved water models are needed for the study of supercooled water.

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

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.

Deterministic Quantum Mechanics: Classical nuclear motion explains chemical reactions and spectra

Is a classical description of nuclear motion sufficient when describing chemical reactions? The present paper investigates some phenomena that were previously attributed to nuclear quantum effects. The aim is to show that these phenomena can be modelled with traditional Car-Parrinello molecular dynamics, that is, with a method which treats nuclear motion classically. We find that no additional paradigm is needed for describing chemical reactions. The special reactivity observed for carbenes can be attributed to the special environment represented by a noble gas matrix. Also the infrared spectrum of porphycene is perfectly modelled by traditional Car-Parrinello molecular dynamics. If no more convincing examples are produced, one will stick to deterministic quantum mechanics, as it is the simpler theory which, in addition, is free of paradoxa.

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

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.

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller–Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers’ stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm−1–1700 cm−1 and 2300 cm−1–3400 cm−1 in the gas phase and 600 cm−1–1800 cm−1 and 2200 cm−1–3400 cm−1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm−1–1700 cm−1 and 2300 cm−1–3300 cm−1 for the gas phase and one broad absorption region in the solid state between 700 cm−1 and 3100 cm−1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.