Isomeric Equilibria, Nuclear Quantum Effects, and Vibrational Spectra of M+(H2O)n=1–3 Clusters, with M = Li, Na, K, Rb, and Cs, through Many-Body Representations

2018 ◽  
Vol 122 (27) ◽  
pp. 5811-5821 ◽  
Author(s):  
Marc Riera ◽  
Sandra E. Brown ◽  
Francesco Paesani
2018 ◽  
Author(s):  
Marc Riera ◽  
Sandra E. Brown ◽  
Francesco Paesani

<div> <div> <div> <p>A quantitative characterization of the molecular mechanisms that regulate ion solvation is key to the microscopic understanding of fundamental processes taking place in aqueous environments, with major implications for different fields, from atmospheric chemistry to materials research and biochemistry. This study presents a systematic analysis of isomeric equilibria for small M<sup>+</sup>(H<sub>2</sub>O)<sub>n</sub> clusters, with M = Li, Na, K, Rb, and Cs, from 0 K to 200 K. To determine the relative stability of different isomers of each M<sup>+</sup>(H<sub>2</sub>O)<sub>n </sub>cluster as a function of temperature, replica exchange simulations are carried out at both classical and quantum levels with the recently developed many-body MB-nrg potential energy functions, which have been shown to exhibit chemical accuracy. Anharmonic vibrational spectra are then calculated within the local monomer approximation and found to be in agreement with the available experimental data, providing further support for the accuracy of the MB-nrg potential energy functions. The present analysis indicates that nuclear quantum effects become increasingly im- portant for larger M<sup>+</sup>(H<sub>2</sub>O)<sub>n </sub>clusters containing the heavier alkali metal ions, which is explained in terms of competing ion-water and water-water interactions along with the interplay between energetic and entropic effects. Directly connecting experimental measurements with molecular properties calculated at the quantum mechanical level, this study represents a further step toward the development of a consistent picture of ion hydration from the gas to the condensed phase. </p> </div> </div> </div>


2018 ◽  
Author(s):  
Marc Riera ◽  
Sandra E. Brown ◽  
Francesco Paesani

<div> <div> <div> <p>A quantitative characterization of the molecular mechanisms that regulate ion solvation is key to the microscopic understanding of fundamental processes taking place in aqueous environments, with major implications for different fields, from atmospheric chemistry to materials research and biochemistry. This study presents a systematic analysis of isomeric equilibria for small M<sup>+</sup>(H<sub>2</sub>O)<sub>n</sub> clusters, with M = Li, Na, K, Rb, and Cs, from 0 K to 200 K. To determine the relative stability of different isomers of each M<sup>+</sup>(H<sub>2</sub>O)<sub>n </sub>cluster as a function of temperature, replica exchange simulations are carried out at both classical and quantum levels with the recently developed many-body MB-nrg potential energy functions, which have been shown to exhibit chemical accuracy. Anharmonic vibrational spectra are then calculated within the local monomer approximation and found to be in agreement with the available experimental data, providing further support for the accuracy of the MB-nrg potential energy functions. The present analysis indicates that nuclear quantum effects become increasingly im- portant for larger M<sup>+</sup>(H<sub>2</sub>O)<sub>n </sub>clusters containing the heavier alkali metal ions, which is explained in terms of competing ion-water and water-water interactions along with the interplay between energetic and entropic effects. Directly connecting experimental measurements with molecular properties calculated at the quantum mechanical level, this study represents a further step toward the development of a consistent picture of ion hydration from the gas to the condensed phase. </p> </div> </div> </div>


2021 ◽  
Author(s):  
Xi Xu ◽  
Zehua Chen ◽  
Yang Yang

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.


2019 ◽  
Vol 10 (35) ◽  
pp. 8211-8218 ◽  
Author(s):  
Marc Riera ◽  
Eleftherios Lambros ◽  
Thuong T. Nguyen ◽  
Andreas W. Götz ◽  
Francesco Paesani

Two-body and three-body energies, modulated by higher-body terms and nuclear quantum effects, determine the structure of liquid water and require sub-chemical accuracy that is achieved by the MB-pol model but not by existing DFT functionals.


2021 ◽  
Vol 22 (10) ◽  
pp. 5220
Author(s):  
Jarosław J. Panek ◽  
Joanna Zasada ◽  
Bartłomiej M. Szyja ◽  
Beata Kizior ◽  
Aneta Jezierska

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.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Huziel E. Sauceda ◽  
Valentin Vassilev-Galindo ◽  
Stefan Chmiela ◽  
Klaus-Robert Müller ◽  
Alexandre Tkatchenko

AbstractNuclear quantum effects (NQE) tend to generate delocalized molecular dynamics due to the inclusion of the zero point energy and its coupling with the anharmonicities in interatomic interactions. Here, we present evidence that NQE often enhance electronic interactions and, in turn, can result in dynamical molecular stabilization at finite temperature. The underlying physical mechanism promoted by NQE depends on the particular interaction under consideration. First, the effective reduction of interatomic distances between functional groups within a molecule can enhance the n → π* interaction by increasing the overlap between molecular orbitals or by strengthening electrostatic interactions between neighboring charge densities. Second, NQE can localize methyl rotors by temporarily changing molecular bond orders and leading to the emergence of localized transient rotor states. Third, for noncovalent van der Waals interactions the strengthening comes from the increase of the polarizability given the expanded average interatomic distances induced by NQE. The implications of these boosted interactions include counterintuitive hydroxyl–hydroxyl bonding, hindered methyl rotor dynamics, and molecular stiffening which generates smoother free-energy surfaces. Our findings yield new insights into the versatile role of nuclear quantum fluctuations in molecules and materials.


Sign in / Sign up

Export Citation Format

Share Document