We present a novel approach to calculate molecular infrared (IR) spectra based on semiclassical (SC) molecular dynamics. The main advance from a previous SC method [M. Micciarelli et al. J. Chem. Phys. 149, 064115 (2018)] consists of the possibility to avoid state-to-state calculations making applications to systems characterized by sizable densities of vibrational states feasible. Furthermore, this new method accounts not only for positions and intensities of the several absorption bands which make up the IR spectrum but also for their shapes. We show that accurate SC IR spectra including quantum effects and anharmonicities for both frequencies and intensities can be obtained starting from SC power spectra. The approach is first tested against the water molecule and then applied to the 10-atom glycine amino acid.

An effective semiclassical approach to IR spectroscopy.

A mixed semiclassical initial value representation expression for spectroscopic calculations is derived. The formulation takes advantage of the time-averaging filtering and the hierarchical properties of different trajectory based propagation methods. A separable approximation is then introduced that greatly reduces (about an order of magnitude) the computational cost compared with a full Herman-Kluk time-averaging semiclassical calculation for the same systems. The expression is exact for the harmonic case and it is tested numerically for a Morse potential coupled to one or two additional harmonic degrees of freedom. Results are compared to full Herman-Kluk time-averaging calculations and exact quantum wavepacket propagations. We found the peak positions of the mixed semiclassical approximations to be always in very good agreement with full quantum calculations, while overtone peak intensities are lower with respect to the exact ones. Given the reduced computational effort required by this new mixed semiclassical approximation, we believe the present method to make spectroscopic calculations available for higher dimensional systems than accessible before.

Mixed Semiclassical Initial Value Representation Time-Averaging Propagator for Spectroscopic Calculations

In this work, we report a computational study of the vibrational features of four different nucleobases employing the divide-and-conquer semiclassical initial value representation molecular dynamics method. Calculations are performed on uracil, cytosine, thymine, and adenine. Results show that the overall accuracy with respect to experiments is within 20 wavenumbers, regardless of the dimensionality of the nucleobase. Vibrational estimates are accurate even in the complex case of cytosine, where two relevant conformers are taken into account. These results are promising in the perspective of future studies on more complex systems, such as nucleotides or nucleobase pairs.

/pdf/vibrational-investigation-of-nucleobases-by-means-of-divide-4xyhcozxb6.pdf

Vibrational investigation of nucleobases by means of divide and conquer semiclassical dynamics

We report on a new approach to ease the computational overhead of ab initio “on-the-fly” semiclassical dynamics simulations for vibrational spectroscopy. The well known bottleneck of such computations lies in the necessity to estimate the Hessian matrix for propagating the semiclassical pre-exponential factor at each step along the dynamics. The procedure proposed here is based on the creation of a dynamical database of Hessians and associated molecular geometries able to speed up calculations while preserving the accuracy of results at a satisfactory level. This new approach can be interfaced to both analytical potential energy surfaces and on-the-fly dynamics, allowing one to study even large systems previously not achievable. We present results obtained for semiclassical vibrational power spectra of methane, glycine, and N-acetyl-L-phenylalaninyl-L-methionine-amide, a molecule of biological interest made of 46 atoms.

Semiclassical vibrational spectroscopy with Hessian databases

In this review, we give a comprehensive comparison of the most widely used coherent state (CS) based methods to solve the time-dependent Schrodinger equation (TDSE). Starting from the fully variati...

M. Buchholz

Papers

Herman-Kluk propagator is free from zero-point energy leakage

Interference in semiclassical electron dynamics