Vibrational spectra of highly anharmonic water clusters: Molecular dynamics and harmonic analyses revisited with constrained nuclear-electronic orbital methods

Yang Yang,Yuzhe Zhang,Yiwen Wang,Zehua Chen,Xi Xu

DOI: https://doi.org/10.26434/chemrxiv-2023-0lb60-v3

2023-11-15

Abstract:Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological and material systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.

Chemistry

What problem does this paper attempt to address?

The paper attempts to address the problem of accurately calculating vibrational spectra in highly anharmonic water clusters, such as neutral water dimers, protonated water trimers, and protonated water tetramers. Existing methods, including traditional ab initio molecular dynamics (AIMD) and harmonic analysis, have limitations when dealing with these systems, particularly in describing quantum nuclear delocalization effects and vibrational mode coupling effects. The authors propose a microcanonical molecular dynamics framework based on constrained nuclear-electronic orbital density functional theory (CNEO-DFT), combined with a direct NVE scheme and its improved version, the NVE-eqp method, to calculate the vibrational spectra of these complex systems more efficiently and accurately. Specifically, the paper addresses the following issues: 1. **Improving computational efficiency**: By using the direct NVE scheme and its improved version, the NVE-eqp method, the simulation time required to reach convergence is reduced, thereby improving computational efficiency. 2. **Enhancing spectral accuracy**: The CNEO-MD method better describes quantum nuclear delocalization effects and vibrational mode coupling effects, thereby improving the accuracy of vibrational spectrum calculations. 3. **Validating the effectiveness of the new method**: The effectiveness and superiority of the CNEO-MD combined with the direct NVE-eqp method are validated in three challenging water cluster systems. In summary, this paper aims to develop an efficient and accurate method for calculating the vibrational spectra of highly anharmonic water clusters to overcome the shortcomings of existing methods.

Vibrational spectra of highly anharmonic water clusters: Molecular dynamics and harmonic analyses revisited with constrained nuclear-electronic orbital methods

Ab Initio Investigation Of Structure, Stability, Thermal Behavior, Bonding, And Infrared Spectra Of Ionized Water Cluster (H2o)(6)(+) (Vol 145, 154307, 2016)

Ab Initio Investigation of the Lower Energy Candidate Structures for (h2o)(5)+ Water Cluster

Ab Initio Investigation of Possible Candidate Structures and Properties of Water Cluster (H2O)7+ Via Particle Swarm Optimization Method

Reduced Rovibrational Coupling Cartesian Dynamics for Semiclassical Calculations: Application to the Spectrum of the Zundel Cation

Understanding the Origins of Dipolar Couplings and Correlated Motion in the Vibrational Spectrum of Water

H2O Trimer: Rigorous 12D Quantum Calculations of Intermolecular Vibrational States, Tunneling Splittings, and Low-Frequency Spectrum

Dynamics and structural changes of small water clusters on ionization

Theoretical Study of Binding Interactions and Vibrational Raman Spectra of Water in Hydrogen-Bonded Anionic Complexes: (H2o)n- (n = 2 and 3), H2O···X-(X = F, Cl, Br, and I), and H2O···M- (M = Cu, Ag, and Au)

Disentangling the Complex Vibrational Spectrum of the Protonated Water Trimer, H+(H2O)3, with Two-Color IR-IR Photodissociation of the Bare Ion and Anharmonic VSCF/VCI Theory

Cluster-model DFT Simulations of the Infrared Spectra of Triazine-Based Molecular Crystals.

Full dimensional (15D) quantum-dynamical simulation of the protonated water-dimer II: infrared spectrum and vibrational dynamics

Vibrational Spectra of the OH Radical in Water: Ab Initio Molecular Dynamics Simulations and Quantum Chemical Calculations Using Hybrid Functionals

Impact of Nuclear Quantum Effects on the Structural Properties of Protonated Water Clusters

Accurate Calculation of Tunneling Splittings in Water Clusters Using Path-Integral Based Methods.

Molecular Vibrational Spectral Simulation Connects Theoretical Cluster Structure Identification and Vibrational Spectral Evidence

Digital Quantum Simulation Of Molecular Vibrations

Molecular-Dynamics and First-Principles Calculations of Raman Spectra and Molecular and Electronic Structure of Hydrogen Clusters in Hydrogen Clathrate Hydrate

Full-dimensional Quantum Studies of Vibrational Energy Transfer Dynamics Between H2O and Ar: Theory Assessing Experiment.

Intermolecular Bending States and Tunneling Splittings of Water Trimer from Rigorous 9D Quantum Calculations: I. Methodology, Energy Levels, and Low-Frequency Spectrum

Vibrational strong coupling in liquid water from cavity molecular dynamics