Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine.

We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.

A time averaged semiclassical approach to IR spectroscopy.

We propose a new semiclassical approach to the calculation of molecular IR spectra. The method employs the time averaging technique of Kaledin and Miller upon symmetrization of the quantum dipole-dipole autocorrelation function. Spectra at high and low temperatures are investigated. In the first case, we are able to point out the possible presence of hot bands in the molecular absorption line shape. In the second case, we are able to reproduce accurate IR spectra as demonstrated by a calculation of the IR spectrum of the water molecule, which is within 4% of the exact intensity. Our time averaged IR spectra can be directly compared to time averaged semiclassical power spectra as shown in an application to the CO2 molecule, which points out the differences between IR and power spectra and demonstrates that our new approach can identify active IR transitions correctly. Overall, the method features excellent accuracy in calculating absorption intensities and provides estimates for the frequencies of vibrations in agreement with the corresponding power spectra. In perspective, this work opens up the possibility to interface the new method with the semiclassical techniques developed for power spectra, such as the divide-and-conquer one, to get accurate IR spectra of complex and high-dimensional molecular systems.

The Fate of the Formic Acid Proton on the Anatase TiO2(101) Surface.

As a prototype adsorption reaction of gas Brønsted acid on oxides, we study the adsorption of formic acid on anatase. We perform infrared spectroscopy measurements of adsorbed HCOOH and HCOOD on TiO2 nanopowders, from 13 K up to room temperature in an ultra-high vacuum chamber. We assign the IR signals via computed spectra from nuclear quantum dynamics simulations using our divide-and-conquer semiclassical ab initio molecular dynamics method. The acid proton forms an extraordinarily short and strong hydrogen bond with the surface oxygen. The strength of this hydrogen bond, that compares to H bonds in ice at high pressures, is at the root of a substantial redshift with respect to the typical free OH stretching frequency, which eludes its straightforward detection.

The first HyDRA challenge for computational vibrational spectroscopy.

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge. In the latter case, the final result is unknown to modellers, providing an unbiased testing opportunity for quantum chemical models. In this work, we present the spectroscopic and computational results for the first HyDRA blind challenge. The latter deals with the prediction of water donor stretching vibrations in monohydrates of organic molecules. This edition features a test set of 10 systems. Experimental water donor OH vibrational wavenumbers for the vacuum-isolated monohydrates of formaldehyde, tetrahydrofuran, pyridine, tetrahydrothiophene, trifluoroethanol, methyl lactate, dimethylimidazolidinone, cyclooctanone, trifluoroacetophenone and 1-phenylcyclohexane-cis-1,2-diol are provided. The results of the challenge show promising predictive properties in both purely quantum mechanical approaches as well as regression and other machine learning strategies.

Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region.

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.

Zwitter Ionization of Glycine at Outer Space Conditions due to Microhydration by Six Water Molecules.

We investigate glycine microsolvation with water molecules, mimicking astrophysical conditions, in our laboratory by embedding these clusters in helium nanodroplets at 0.37 K. We recorded mass selective infrared spectra in the frequency range 1500-1800 cm^{-1} where two bands centered at 1630 and 1724 cm^{-1} were observed. By comparison with the extensive accompanying calculations, the band at 1630 cm^{-1} was assigned to the COO^{-} asymmetric stretching mode of the zwitter ion and the band at 1724 cm^{-1} was assigned to redshifted C=O stretch within neutral clusters. We show that zwitter ion formation of amino acids readily occurs with only few water molecules available even under extreme conditions.

Caldeira-Leggett model vs ab initio potential: A vibrational spectroscopy test of water solvation.

We present a semiclassically approximate quantum treatment of solvation with the purpose of investigating the accuracy of the Caldeira-Leggett model. We do that by simulating the vibrational features of water solvation by means of two different approaches. One is entirely based on the adoption of an accurate ab initio potential to describe water clusters of increasing dimensionality. The other one consists of a model made of a central water molecule coupled to a high-dimensional Caldeira-Leggett harmonic bath. We demonstrate the role of quantum effects in the detection of water solvation and show that the computationally cheap approach based on the Caldeira-Leggett bath is only partially effective. The main conclusion of the study is that quantum methods associated with high-level potential energy surfaces are necessary to correctly study solvation features, while simplified models, even if attractive owing to their reduced computational cost, can provide some useful insights but are not able to come up with a comprehensive description of the solvation phenomenon.

On-the-fly adiabatically switched semiclassical initial value representation molecular dynamics for vibrational spectroscopy of biomolecules.

Semiclassical (SC) vibrational spectroscopy is a technique capable of reproducing quantum effects (such as zero-point energies, quantum resonances, and anharmonic overtones) from classical dynamics runs even in the case of very large dimensional systems. In a previous study [Conte et al. J. Chem. Phys. 151, 214107 (2019)], a preliminary sampling based on adiabatic switching has been shown to be able to improve the precision and accuracy of semiclassical results for challenging model potentials and small molecular systems. In this paper, we investigate the possibility to extend the technique to larger (bio)molecular systems whose dynamics must be integrated by means of ab initio "on-the-fly" calculations. After some preliminary tests on small molecules, we obtain the vibrational frequencies of glycine improving on pre-existing SC calculations. Finally, the new approach is applied to 17-atom proline, an amino acid characterized by a strong intramolecular hydrogen bond.

Machine learning for vibrational spectroscopy via divide-and-conquer semiclassical initial value representation molecular dynamics with application to N-methylacetamide.

A machine learning algorithm for partitioning the nuclear vibrational space into subspaces is introduced. The subdivision criterion is based on Liouville's theorem, i.e., the best preservation of the unitary of the reduced dimensionality Jacobian determinant within each subspace along a probe full-dimensional classical trajectory. The algorithm is based on the idea of evolutionary selection, and it is implemented through a probability graph representation of the vibrational space partitioning. We interface this customized version of genetic algorithms with our divide-and-conquer semiclassical initial value representation method for the calculation of molecular power spectra. First, we benchmark the algorithm by calculating the vibrational power spectra of two model systems, for which the exact subspace division is known. Then, we apply it to the calculation of the power spectrum of methane. Exact calculations and full-dimensional semiclassical spectra of this small molecule are available and provide an additional test of the accuracy of the new approach. Finally, the algorithm is applied to the divide-and-conquer semiclassical calculation of the power spectrum of 12-atom trans-N-methylacetamide.

Representing molecular ground and excited vibrational eigenstates with nuclear densities obtained from semiclassical initial value representation molecular dynamics.

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 the protonated glycine molecule [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 the 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, also 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.

Anharmonic calculations of vibrational spectra for molecular adsorbates: A divide-and-conquer semiclassical molecular dynamics approach.

The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born-Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H2O on the anatase-TiO2 (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H2O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach.

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

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.

Semiclassical vibrational spectroscopy with Hessian databases.

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.

An effective semiclassical approach to IR spectroscopy.

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.

Improved semiclassical dynamics through adiabatic switching trajectory sampling.

We introduce an improved semiclassical dynamics approach to quantum vibrational spectroscopy. In this method, a harmonic-based phase space sampling is preliminarily driven toward non-harmonic quantization by slowly switching on the actual potential. The new coordinates and momenta serve as initial conditions for the semiclassical dynamics calculation, leading to a substantial decrease in the number of chaotic trajectories to deal with. Applications are presented for model and molecular systems of increasing dimensionality characterized by moderate or high chaoticity. They include a bidimensional Henon-Heiles potential, water, formaldehyde, and methane. The method improves accuracy and precision of semiclassical results and it can be easily interfaced with all pre-existing semiclassical theories.

Simplified approach to the mixed time-averaging semiclassical initial value representation for the calculation of dense vibrational spectra.

We present and test an approximate method for the semiclassical calculation of vibrational spectra. The approach is based on the mixed time-averaging semiclassical initial value representation method, which is simplified to a form that contains a filter to remove contributions from approximately harmonic environmental degrees of freedom. This filter comes at no additional numerical cost, and it has no negative effect on the accuracy of peaks from the anharmonic system of interest. The method is successfully tested for a model Hamiltonian and then applied to the study of the frequency shift of iodine in a krypton matrix. Using a hierarchic model with up to 108 normal modes included in the calculation, we show how the dynamical interaction between iodine and krypton yields results for the lowest excited iodine peaks that reproduce experimental findings to a high degree of accuracy.

"Divide-and-conquer" semiclassical molecular dynamics: An application to water clusters.

We present an investigation of vibrational features in water clusters performed by means of our recently established divide-and-conquer semiclassical approach [M. Ceotto, G. Di Liberto, and R. Conte, Phys. Rev. Lett. 119, 010401 (2017)]. This technique allows us to simulate quantum vibrational spectra of high-dimensional systems starting from full-dimensional classical trajectories and projection of the semiclassical propagator onto a set of lower dimensional subspaces. The potential energy surface employed is a many-body representation up to three-body terms, in which monomers and two-body interactions are described by the high level Wang-Huang-Braams-Bowman (WHBB) water potential, while, for three-body interactions, calculations adopt a fast permutationally invariant ab initio surface at the same level of theory of the WHBB 3-body potential. Applications range from the water dimer up to the water decamer, a system made of 84 vibrational degrees of freedom. Results are generally in agreement with previous variational estimates in the literature. This is particularly true for the bending and the high-frequency stretching motions, while estimates of modes strongly influenced by hydrogen bonding are red shifted, in a few instances even substantially, as a consequence of the dynamical and global picture provided by the semiclassical approach.

"Divide and conquer" semiclassical molecular dynamics: A practical method for spectroscopic calculations of high dimensional molecular systems.

We extensively describe our recently established "divide-and-conquer" semiclassical method [M. Ceotto, G. Di Liberto, and R. Conte, Phys. Rev. Lett. 119, 010401 (2017)] and propose a new implementation of it to increase the accuracy of results. The technique permits us to perform spectroscopic calculations of high-dimensional systems by dividing the full-dimensional problem into a set of smaller dimensional ones. The partition procedure, originally based on a dynamical analysis of the Hessian matrix, is here more rigorously achieved through a hierarchical subspace-separation criterion based on Liouville's theorem. Comparisons of calculated vibrational frequencies to exact quantum ones for a set of molecules including benzene show that the new implementation performs better than the original one and that, on average, the loss in accuracy with respect to full-dimensional semiclassical calculations is reduced to only 10 wavenumbers. Furthermore, by investigating the challenging Zundel cation, we also demonstrate that the "divide-and-conquer" approach allows us to deal with complex strongly anharmonic molecular systems. Overall the method very much helps the assignment and physical interpretation of experimental IR spectra by providing accurate vibrational fundamentals and overtones decomposed into reduced dimensionality spectra.

Anharmonic vibrational eigenfunctions and infrared spectra from semiclassical molecular dynamics.

We describe a new approach based on semiclassical molecular dynamics that allows simulating infrared absorption or emission spectra of molecular systems with inclusion of anharmonic intensities. This is achieved from semiclassical power spectra by computing first the vibrational eigenfunctions as a linear combination of harmonic states, and then the oscillator strengths associated with the vibrational transitions. We test the approach against a 1D Morse potential and apply it to the water molecule with results in excellent agreement with discrete variable representation quantum benchmarks. The method does not require any grid calculations, and it is directly extendable to high dimensional systems. The usual exponential scaling of the basis set size with the dimensionality of the system can be avoided by means of an appropriate truncation scheme. Furthermore, the approach has the advantage to provide IR spectra beyond the harmonic approximation without losing the possibility of an intuitive assignment of absorption peaks in terms of normal modes of vibration.

A quantum mechanical insight into SN2 reactions: Semiclassical initial value representation calculations of vibrational features of the Cl-⋯CH3Cl pre-reaction complex with the VENUS suite of codes.

The role of vibrational excitation of reactants in driving reactions involving polyatomic species has been often studied by means of classical or quasi-classical trajectory simulations. We propose a different approach based on investigation of vibrational features of the Cl-⋯CH3Cl pre-reaction complex for the Cl- + CH3Cl SN2 reaction. We present vibrational power spectra and frequency estimates for the title pre-reaction complex calculated at the level of classical, semiclassical, and second-order vibrational perturbation theory on a pre-existing analytical potential energy surface. The main goals of the paper are the study of anharmonic effects and understanding of vibrational couplings that permit energy transfer between the collisional kinetic energy and the internal vibrations of the reactants. We provide both classical and quantum pictures of intermode couplings and show that the SN2 mechanism is favored by the coupling of a C-Cl bend involving the Cl- projectile with the CH3 rocking motion of the target molecule. We also illustrate how the routines needed for semiclassical vibrational spectroscopy simulations can be interfaced in a user-friendly way to pre-existing molecular dynamics software. In particular, we present an implementation of semiclassical spectroscopy into the VENUS suite of codes, thus providing a useful computational tool for users who are not experts of semiclassical dynamics.