Fingerprint region of the formic acid dimer: variational vibrational computations in curvilinear coordinates

Methane dimer rovibrational states and Raman transition moments

Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible ΔK ≠ 0 transition moments arise in this symmetric top system due to strong rovibrational couplings.

Exact quantum dynamics developments for floppy molecular systems and complexes

Molecular rotation, vibration, internal rotation, isomerization, tunneling, intermolecular dynamics of weakly and strongly interacting systems, intra-to-inter-molecular energy transfer, hindered rotation and hindered translation over surfaces are important types of molecular motions. Their fundamentally correct and detailed description can be obtained by solving the nuclear Schrödinger equation on a potential energy surface. Many of the chemically interesting processes involve quantum nuclear motions which are 'delocalized' over multiple potential energy wells. These 'large-amplitude' motions in addition to the high dimensionality of the vibrational problem represent challenges to the current (ro)vibrational methodology. A review of the quantum nuclear motion methodology is provided, current bottlenecks of solving the nuclear Schrödinger equation are identified, and solution strategies are reviewed. Technical details, computational results, and analysis of these results in terms of limiting models and spectroscopically relevant concepts are highlighted for selected numerical examples.

High-dimensional neural network potentials for accurate vibrational frequencies: the formic acid dimer benchmark

In recent years, machine learning potentials (MLP) for atomistic simulations have attracted a lot of attention in chemistry and materials science. Many new approaches have been developed with the primary aim to transfer the accuracy of electronic structure calculations to large condensed systems containing thousands of atoms. In spite of these advances, the reliability of modern MLPs in reproducing the subtle details of the multi-dimensional potential-energy surface is still difficult to assess for such systems. On the other hand, moderately sized systems enabling the application of tools for thorough and systematic quality-control are nowadays rarely investigated. In this work we use benchmark-quality harmonic and anharmonic vibrational frequencies as a sensitive probe for the validation of high-dimensional neural network potentials. For the case of the formic acid dimer, a frequently studied model system for which stringent spectroscopic data became recently available, we show that high-quality frequencies can be obtained from state-of-the-art calculations in excellent agreement with coupled cluster theory and experimental data.

Slow monomer vibrations in formic acid dimer: Stepping up the ladder with FTIR and Raman jet spectroscopy

In an effort to extend the cold gas phase spectroscopic database of the cyclic formic acid dimer (FAD), we present and analyze the jet-cooled vibrational infrared and Raman spectrum of (HCOOH)₂ in the monomer fingerprint region between 600 and 1500 cm^-1. The present study bridges the gap between the intermolecular dimerization-induced and the carbonyl stretching fundamentals that have already been reexamined using jet-cooled or high-resolution spectroscopy. This completes the characterization of the jet-cooled vibrational (HCOOH)₂ spectrum below the complex OH (CH) stretching fundamentals, and we report resonance-induced FAD combination/overtone transitions that will serve as a valuable reference for a theoretical modeling of its vibrational dynamics. As a by-product, several new formic acid trimer fundamentals are identified in the jet spectra and assigned with the help of second-order vibrational perturbation theory (VPT2). The polar formic acid dimer still eludes detection in a supersonic jet, but we are able to estimate an experimental upper-bound of the polar dimer-to-trimer-to-cyclic dimer intensity ratio to about 1:10:100 under typical expansion conditions. Using VPT2 with resonance treatment (VPT2+K), we reinvestigate the notorious ν₂₂ resonance triad. Generally, we find that VPT2, which is, of course, inadequate for modeling the resonance-rich OH stretching spectrum of FAD, is performing very satisfactorily in predicting fundamental and two-quantum state term values for the slower modes below 1500 cm^-1. As these modes are the building blocks for the ultrafast energy dissipation in the OH stretching region, the present work opens the door for its quantitative understanding.

Transfer learned potential energy surfaces: accurate anharmonic vibrational dynamics and dissociation energies for the formic acid monomer and dimer

The vibrational dynamics of the formic acid monomer (FAM) and dimer (FAD) is investigated from machine-learned potential energy surfaces at the MP2 (PES_MP2) and transfer-learned (PES_TL) to the CCSD(T) levels of theory. The normal mode (MAEs of 17.6 and 25.1 cm⁻¹) and second order vibrational perturbation theory (VPT2, MAEs of 6.7 and 17.1 cm⁻¹) frequencies from PES_TL for all modes below 2000 cm⁻¹ for FAM and FAD agree favourably with experiment. For the OH stretch mode the experimental frequencies are overestimated by more than 150 cm⁻¹ for both FAM and FAD from normal mode calculations. Conversely, VPT2 calculations on PES_TL for FAM reproduce the experimental OH frequency to within 22 cm⁻¹. For FAD the VPT2 calculations find the high-frequency OH stretch at 3011 cm⁻¹, compared with an experimentally reported, broad (∼100 cm⁻¹) absorption band with center frequency estimated at ∼3050 cm⁻¹. In agreement with earlier reports, MD simulations at higher temperature shift the position of the OH-stretch in FAM to the red, consistent with improved sampling of the anharmonic regions of the PES. However, for FAD the OH-stretch shifts to the blue and for temperatures higher than 1000 K the dimer partly or fully dissociates using PES_TL. Including zero-point energy corrections from diffusion Monte Carlo simulations for FAM and FAD and corrections due to basis set superposition and completeness errors yields a dissociation energy of D₀ = −14.23 ± 0.08 kcal mol⁻¹ compared with an experimentally determined value of −14.22 ± 0.12 kcal mol⁻¹.

Neural network based PESs are constructed for formic acid monomer and dimer at the MP2 and transfer learned to the CCSD(T) level of theory. The PESs are used to study the vibrational dynamics and dissociation energy of the molecules.

Performance of a black-box-type rovibrational method in comparison with a tailor-made approach: Case study for the methane-water dimer

The present work intends to join and respond to the excellent and thoroughly documented rovibrational study of X. G. Wang and T. Carrington, Jr. [J. Chem. Phys. 154, 124112 (2021)] that used an approach tailored for floppy dimers with an analytic dimer Hamiltonian and a non-product basis set including Wigner D functions. It is shown in the present work that the GENIUSH black-box-type rovibrational method can approach the performance of the tailor-made computation for the example of the floppy methane-water dimer. Rovibrational transition energies and intensities are obtained in the black-box-type computation with a twice as large basis set and in excellent numerical agreement in comparison with the more efficient tailor-made approach.