Representing Molecular Ground and Excited Vibrational Eigenstates with Nuclear Densities obtained from Semiclassical Initial Value Representation Molecular Dynamics
arXiv:2012.03712 · doi:10.1063/5.0031391
Abstract
We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in protonated glycine molecule [C. Aieta et. al. Nat. Commun. 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, and to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond-lengths, angles and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching the typical harmonic nodal pattern is absent in the anharmonic distribution.
arXiv comment13 pages main text with 11 figures; 7 pages Supplementary material with 9 figures. Accepted by the Journal of Chemical Physics (See link below). Post-print version
References in corpus (21)
- Full dimensional (15D) quantum-dynamical simulation of the protonated water-dimer II: infrared spectrum and vibrational dynamics
- Introduction to the Diffusion Monte Carlo Method
- First-Principles Semiclassical Initial Value Representation Molecular Dynamics
- Semiclassical "Divide-and-Conquer" Method for Spectroscopic Calculations of High Dimensional Molecular Systems
- On-the-fly ab initio semiclassical dynamics: Identifying degrees of freedom essential for emission spectra of oligothiophenes
- 'Divide-and-conquer' semiclassical molecular dynamics: An application to water clusters
- On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature
- "Divide and Conquer" Semiclassical Molecular Dynamics: A practical method for Spectroscopic calculations of High Dimensional Molecular Systems
- Mixed Semiclassical Initial Value Representation Time-Averaging Propagator for Spectroscopic Calculations
- Anharmonic Vibrational Eigenfunctions and Infrared Spectra from Semiclassical Molecular Dynamics
- On-the-fly ab initio semiclassical evaluation of electronic coherences in polyatomic molecules reveals a simple mechanism of decoherence
- Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
- An Effective Semiclassical Approach to IR Spectroscopy
- Simplified Approach to the Mixed Time-averaging Semiclassical Initial Value Representation for the Calculation of Dense Vibrational Spectra
- Validating and Implementing Modified Filinov Phase Filtration in Semiclassical Dynamics
- Graphics processing units accelerated semiclassical initial value representation molecular dynamics
- Application of the Mixed Time-averaging Semiclassical Initial Value Representation method to Complex Molecular Spectra
- Reduced Rovibrational Coupling Cartesian Dynamics for Semiclassical Calculations: Application to the Spectrum of the Zundel Cation
- Vibrational Investigation of Nucleobases by Means of Divide and Conquer Semiclassical Dynamics
- Semiclassical vibrational spectroscopy with Hessian databases
- A quantum mechanical insight into SN2 reactions: Semiclassical initial value representation calculations of vibrational features of the Cl--CHCl pre-reaction complex with the VENUS suite of codes