Anharmonic quantum nuclear densities from full dimensional vibrational eigenfunctions with application to protonated glycine

Quantum Nature of Ubiquitous Vibrational Features Revealed for Ethylene Glycol

Vibrational properties of molecules are of widespread interest and importance in chemistry and biochemistry. The reliability of widely employed approximate computational methods is questioned here against the complex experimental spectrum of ethylene glycol. Comparisons between quantum vibrational self-consistent field and virtual-state configuration interaction (VSCF/VCI), adiabatically switched semiclassical initial value representation (AS-SCIVR), and thermostatted ring polymer molecular dynamics (TRPMD) calculations are made using a full-dimensional, machine-learned potential energy surface. Calculations are done for five low-lying conformers and compared with the experiment, with a focus on the high-frequency, OH-stretches, and CH-stretches, part of the spectrum. Fermi resonances are found in the analysis of VSCF/VCI eigenstates belonging to the CH-stretching band. Results of comparable accuracy, quality, and level of detail are obtained by means of AS SCIVR. The current VSCF/VCI and AS-SCIVR power spectra largely close the gaps between the experiment and TRPMD and classical MD calculations. Analysis of these results provides guidance on what level of accuracy to expect from TRPMD and classical MD calculations of the vibrational spectra for ubiquitous CH- and OH-stretching bands. This work shows that even general vibrational features require a proper quantum treatment, usually not achievable by the most popular theoretical approaches.

Theory of Nonadiabatic Tunneling Splitting

Estimating tunneling splittings is a long-standing quantum mechanical challenge for theoretical methods. Sometimes splittings are so small, i.e., within a fraction of a wavenumber, pushing the limits of experimental detection and computational precision. Currently, most computational methods are able, at best, to obtain only ground-state tunneling splittings, either for symmetric or asymmetric potentials. In this Letter, we introduce a unified theoretical approach, based on a two-state approximation that can be equally applied to symmetric and asymmetric diabatic potential crossing and for excited states, providing reliable estimates even for states near the energy crossing. The method opens the door to analytic approximations for the tunneling splitting of model potential systems. It provides a framework for the introduction of vibrational perturbation theory to the estimation of nonadiabatic tunneling splittings. It also provides new insight into the semiclassical theory, leading to an instanton based steepest descent expression applicable also to excited states. Numerical tests on model systems are promising, providing the groundwork for implementation to future multidimensional applications.

Semiclassical description of nuclear quantum effects in solvated and condensed phase molecular systems

In this perspective we deal with the challenge of investigating nuclear quantum effects in solvated and condensed phase molecular systems in a computationally affordable way. To this end, semiclassical methods are promising theoretical approaches, as we demonstrate through vibrational spectroscopy and reaction kinetics. We show that quantum vibrational features can be found in hydrates of carbonyl compounds and microsolvated amino acids, and we report quantum estimates of the low-temperature reaction rate constant of a unimolecular reaction taking place in a noble-gas matrix. The hallmark of semiclassical methods is their ability to include nuclear quantum effects into classical molecular dynamics simulations. For this reason, unlike other popular methods, semiclassical approaches are able to account also for real-time quantum contributions and are expected to point out the importance of nuclear quantum effects in complex systems for a wider range of chemical properties.

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.

A Perspective on the Investigation of Spectroscopy and Kinetics of Complex Molecular Systems with Semiclassical Approaches

In this Perspective we show that semiclassical methods provide a rigorous hierarchical way to study the vibrational spectroscopy and kinetics of complex molecular systems. The time averaged approach to spectroscopy and the semiclassical transition state theory for kinetics, which have been first adopted and then further developed in our group, provide accurate quantum results on rigorous physical grounds and can be applied even when dealing with a large number of degrees of freedom. In spectroscopy, the multiple coherent, divide-and-conquer, and adiabatically switched semiclassical approaches have practically permitted overcoming issues related to the convergence of results. In this Perspective we demonstrate the possibility of studying the semiclassical vibrational spectroscopy of a molecule adsorbed on an anatase (101) surface, a system made of 51 atoms. In kinetics, the semiclassical transition state theory is able to account for anharmonicity and the coupling between the reactive and bound modes. Our group has developed this technique for practical applications involving the study of phenomena like kinetic isotope effect, heavy atom tunneling, and elusive conformer lifetimes. Here, we show that our multidimensional anharmonic quantum approach is able to tackle on-the-fly the thermal kinetic rate constant of a 135 degree-of-freedom system. Overall, semiclassical methods open up the possibility to describe at the quantum mechanical level systems characterized by hundreds of degrees of freedom leading to the accurate spectroscopic and kinetic description of biomolecules and complex molecular systems.

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.

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.

Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces

Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.

Coupled Cluster Semiclassical Estimates of Experimental Reaction Rates: The Interconversion of Glycine Conformer VIp to Ip

We apply the full-dimensional Semiclassical Transition State Theory (SCTST) to estimate the rate constant of glycine molecule interconversion between the VIp and Ip conformers. We have reached an electronic structure accuracy up to the explicitly correlated Coupled Cluster method (CCSD(T)-F12b/cc-pVDZ-F12) thanks to our parallel implementation. The reaction has been experimentally investigated in the literature and is known to proceed by quantum mechanical tunneling. The SCTST rates improve over other theoretical methods, and our results align with the experimental measurements, thus confirming the accuracy of the fully coupled anharmonic semiclassical tunneling treatment, providing that the level of electronic structure theory gives a reliable estimate of the reaction barrier height and shape. The comparison with experimental half-life times supports the validity of SCTST for glycine VIp-Ip conformer conversion in the cryogenic temperature range, where this theory is usually not considered applicable due to the onset of the deep tunneling regime.

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.

Vibrational spectroscopy simulation of solvation effects on a G-quadruplex

It is commonly believed that solvation effects on the vibrational properties of a solute are easily accounted for by simple rules of thumbs, that is, solvating a polar molecule in a polar medium has the only effect of red shifting all its spectroscopical features and, similarly, solvating a polar molecule in a nonpolar medium has the opposite effect. In this work, we use theoretical vibrational spectroscopy at quasi-classical and quantum approximate semiclassical level to gain atomistic insights about solvent-solute interactions for 2'-deoxyguanosine and the G-quadruplex. We employ the quasi-classical trajectory method to include full anharmonicity into our calculated spectra, and then introduce quantum nuclear effects by means of divide-and-conquer semiclassical spectroscopy calculations. Solvation is treated explicitly leading to a good reproducibility of the available experimental data and reliable predictions when an experimental reference is missing.Communicated by Ramaswamy H. Sarma.

Elucidating NO x Surface Chemistry at the Anatase (101) Surface in TiO2 Nanoparticles

Understanding NO _x chemistry at titania nanoparticle surfaces is important for photocatalytic environmental remediation processes. We focus on this problem and put forward an experimental-computational approach based on vibrational spectroscopy grounds. Temperature-dependent IR experiments of NO _x adsorption on shape-engineered nanoparticle (101) anatase surfaces are paired with power spectra obtained from Born-Oppenheimer trajectories. Then, the harmonic versus anharmonic vibrational frequencies of several adsorption scenarios are directly compared with the IR experiments. We conclude that molecules are adsorbed mainly by the N-end side and both the intermolecular interactions between adsorbed molecules and (NO)₂ dimer formation are responsible for the main NO adsorption spectroscopic features. We also investigate the spectroscopy and the mechanism of formation on defective anatase surfaces of the long-lived greenhouse gas N₂O.

Semiclassical and VSCF/VCI Calculations of the Vibrational Energies of trans- and gauche-Ethanol Using a CCSD(T) Potential Energy Surface

A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.

A Semiclassical Framework for Mixed Quantum Classical Dynamics

Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.

Quantum Calculations on a New CCSD(T) Machine-Learned Potential Energy Surface Reveal the Leaky Nature of Gas-Phase Trans and Gauche Ethanol Conformers

Ethanol is a molecule of fundamental interest in combustion, astrochemistry, and condensed phase as a solvent. It is characterized by two methyl rotors and trans (anti) and gauche conformers, which are known to be very close in energy. Here we show that based on rigorous quantum calculations of the vibrational zero-point state, using a new ab initio potential energy surface (PES), the ground state resembles the trans conformer, but substantial delocalization to the gauche conformer is present. This explains experimental issues about identification and isolation of the two conformers. This “leak” effect is partially quenched when deuterating the OH group, which further demonstrates the need for a quantum mechanical approach. Diffusion Monte Carlo and full-dimensional semiclassical dynamics calculations are employed. The new PES is obtained by means of a Δ-machine learning approach starting from a pre-existing low level density functional theory surface. This surface is brought to the CCSD(T) level of theory using a relatively small number of ab initio CCSD(T) energies. Agreement between the corrected PES and direct ab initio results for standard tests is excellent. One- and two-dimensional discrete variable representation calculations focusing on the trans–gauche torsional motion are also reported, in reasonable agreement with experiment.

Quantum Anharmonic Calculations of Vibrational Spectra for Water Adsorbed on Titania Anatase(101) Surface: Dissociative versus Molecular Adsorption

The interaction of water molecules and hydroxyl groups with titanium dioxide (TiO₂) surfaces is ubiquitous and very important in anatase nanoparticle photocatalytic processes. Infrared spectroscopy, assisted by ab initio calculations of vibrational frequencies, can be a powerful tool to elucidate the mechanisms behind water adsorption. However, a straightforward comparison between measurements and calculations remains a challenging task because of the complexity of the physical phenomena occurring on nanoparticle surfaces. Consequently, severe computational approximations, such as harmonic vibrational ones, are usually employed. In the present work we partially address this complexity issue by overcoming some of the standard approximations used in theoretical simulations and employ the Divide and Conquer Semiclassical Initial Value Representation (DC-SCIVR) molecular dynamics. This method allows to perform simulations of vibrational spectra of large dimensional systems accounting not only for anharmonicities, but also for nuclear quantum effects. We apply this computational method to water and deuterated water adsorbed on the ideal TiO₂ anatase(101) surface, contemplating both the molecular and the dissociated adsorption processes. The results highlight not only the presence of an anharmonic shift of the frequencies in agreement with the experiments, but also complex quantum mechanical spectral signatures induced by the coupling of molecular vibrational modes with the surface ones, which are different in the hydrogenated case from the deuterated one. These couplings are further analyzed by exploiting the mode subdivision performed during the divide and conquer procedure.

The complex vibrational spectrum of proline explained through the adiabatically switched semiclassical initial value representation

Proline, a 17-atom amino acid with a closed-ring side chain, has a complex potential energy surface characterized by several minima. Its IR experimental spectrum, reported in the literature, is of difficult and controversial assignment. In particular, the experimental signal at 3559 cm−1 associated with the OH stretch is interesting because it is inconsistent with the global minimum, trans-proline conformer. This suggests the possibility that multiple conformers may contribute to the IR spectrum. The same conclusion is obtained by investigating the splitting of the CO stretch at 1766 and 1789 cm−1 and other, more complex spectroscopic features involving CH stretches and COH/CNH bendings. In this work, we perform full-dimensional, on-the-fly adiabatically switched semiclassical initial value representation simulations employing the ab initio dft-d3-B3LYP level of theory with aug-cc-pVDZ basis set. We reconstruct the experimental spectrum of proline in its main features by studying the vibrational features of trans-proline and cis1-proline, and provide a new assignment for the OH stretch of trans-proline.

Quantum Vibrational Spectroscopy of Explicitly Solvated Thymidine in Semiclassical Approximation

In this paper, we demonstrate the possibility to perform spectroscopy simulations of solvated biological species taking into consideration quantum effects and explicit solvation. We achieve this goal by interfacing our recently developed divide-and-conquer approach for semiclassical initial value representation molecular dynamics with the polarizable AMOEBABIO18 force field. The method is applied to the study of solvation of the thymidine nucleoside in two different polar solvents, water and N,N-dimethylformamide. Such systems are made of up to 2476 atoms. Experimental evidence concerning the different behavior of thymidine in the two solvents is well reproduced by our study, even though quantitative estimates are hampered by the limited accuracy of the classical force field employed. Overall, this study shows that semiclassically approximate quantum dynamical studies of explicitly solvated biological systems are both computationally affordable and insightful.

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.

Unsupervised Machine Learning Neural Gas Algorithm for Accurate Evaluations of the Hessian Matrix in Molecular Dynamics

The Hessian matrix of the potential energy of molecular systems is employed not only in geometry optimizations or high-order molecular dynamics integrators but also in many other molecular procedures, such as instantaneous normal mode analysis, force field construction, instanton calculations, and semiclassical initial value representation molecular dynamics, to name a few. Here, we present an algorithm for the calculation of the approximated Hessian in molecular dynamics. The algorithm belongs to the family of unsupervised machine learning methods, and it is based on the neural gas idea, where neurons are molecular configurations whose Hessians are adopted for groups of molecular dynamics configurations with similar geometries. The method is tested on several molecular systems of different dimensionalities both in terms of accuracy and computational time versus calculating the Hessian matrix at each time-step, that is, without any approximation, and other Hessian approximation schemes. Finally, the method is applied to the on-the-fly, full-dimensional simulation of a small synthetic peptide (the 46 atom N-acetyl-l-phenylalaninyl-l-methionine amide) at the level of DFT-B3LYP-D/6-31G* theory, from which the semiclassical vibrational power spectrum is calculated.

Laser-Driven Anharmonic Oscillator: Ground-State Dissociation of the Helium Hydride Molecular Ion by Midinfrared Pulses

The vibrational motion of molecules represents a fundamental example of an anharmonic oscillator. Using a prototype molecular system, HeH^{+}, we demonstrate that appropriate laser pulses make it possible to drive the nuclear motion in the anharmonic potential of the electronic ground state, increasing its energy above the potential barrier and facilitating dissociation by purely vibrational excitation. We find excellent agreement between the frequency-dependent response of the helium hydride molecular cation to both classical and quantum mechanical simulations, thus removing any ambiguities through electronic excitation. Our results provide access to the rich dynamics of anharmonic quantum oscillator systems and pave the way to state-selective control schemes in ground-state chemistry by the adequate choice of the laser parameters.

Laser-Induced Electron Symmetry Restoration in Oriented Molecules Made Simple

Electron symmetry determines many important properties of molecules, from selection rules for photoelectron spectroscopy to symmetry selection rules for chemical reactions. The original electron symmetry is broken if a laser pulse changes the initial state, typically the ground state g, to a superposition of g and an excited state e with different irreducible representations (IRREPs). Quantum dynamics simulations for two examples, the oriented benzene and LiCN molecules, show that the original electron symmetry can be restored by means of a reoptimized π-laser pulse which transfers the component in the excited state e to another state e', or to several others with the same IRREP as the ground state. This method lends itself to much easier experimental applications than all previous ones because it allows the healing of electron symmetry immediately, without any attosecond constraint on the timing of the second pulse.

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.

How many water molecules are needed to solvate one?

Many efforts undertaken to study the solvation process have led to general theories that may describe mean properties, but are unable to provide a detailed understanding at the molecular level. Remarkably, the basic question of how many solvent molecules are necessary to solvate one solute molecule is still open. By exploring several water aggregates of increasing complexity, in this contribution we employ semiclassical spectroscopy to determine on quantum dynamical grounds the minimal network of surrounding water molecules to make the central one display the same vibrational features of liquid water. We find out that double-acceptor double-donor tetrahedral coordination constituting the standard picture is necessary but not sufficient, and that particular care must be reserved for the quantum description of the combination band due to the coupling of the central monomer bending mode with network librations. It is actually our ability to investigate the combination band with a quantum-derived approach that allows us to answer the titular question. The minimal structure eventually responsible for proper solvation is made of a total of 21 water molecules and includes two complete solvation shells, of which the whole first one is tetrahedrally coordinated to the central molecule.

How quantum spectroscopic simulations can explain water solvation by comparison with experimental spectra.

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.