Jet-cooled infrared spectra of the formic acid dimer by cavity ring-down spectroscopy: observation of the OH stretching region

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.

Coupled proton vibrations between two weak acids: the hinge complex between formic acid and trifluoroethanol

Raman and FTIR spectra of an acid–alcohol complex show complementary signatures from acidic and alcoholic OH stretching, proving its existence.

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.

Dissecting intermolecular interactions in the condensed phase of ibuprofen and related compounds: the specific role and quantification of hydrogen bonding and dispersion forces

Ibuprofen is a well-established non-steroidal anti-inflammatory drug, inhibiting the prostaglandin-endoperoxide synthase. One of the key features defining the ibuprofen structure is the doubly intermolecular O-HO[double bond, length as m-dash]C hydrogen bond in cyclic dimers as know from carboxylic acids and confirmed by X-ray analysis. Until now, there was neither information about the vaporization enthalpy of ibuprofen nor about how this thermal property is determined by the subtle balance between different types of intermolecular interaction. In this study we derive the vaporization enthalpy of ibuprofen from thermochemical experiments to be . We dissected the hydrogen bond energy, EHB = 45.0 kJ mol-1, exclusively from measured vaporization enthalpies of related aliphatic carboxylic acids, their homomorph methyl esters and alkyl acetates, respectively. This contribution from hydrogen bonding could be confirmed almost quantitatively from quantum chemical calculations of ibuprofen clusters, which also suggest dispersion interaction of similar order (Edisp = 47 kJ mol-1). Following the full analysis of the gas-vapor transition enthalpy, we studied the changing structural components from the solid to the liquid phase of ibuprofen by means of Attenuated Total Reflection Infrared (ATR-IR) spectroscopy. The cyclic dimers as observed in the X-ray patterns are essentially preserved in the liquid state just above the melting point. However, with increasing temperature the doubly hydrogen-bonded cyclic dimers are replaced by singly hydrogen-bonded linear dimers in the liquid ibuprofen. The transfer enthalpy from the temperature-dependent equilibria of both dimers as obtained from the IR intensity ratios of the vibrational bands quantifies for the first time the energy of the released, single hydrogen bond to be EHB = 21.0 kJ mol-1. Overall, we show that a combination of thermodynamics, infrared spectroscopy and quantum chemistry provides quantification and detailed understanding of structure and molecular interaction in ibuprofen and related compounds.

Shifting formic acid dimers into perspective: vibrational scrutiny in helium nanodroplets

A metastable dimer of formic acid has been prepared inside superfluid helium nanodroplets and examined using IR spectroscopy and quantum chemical calculations.

Concerted Pair Motion Due to Double Hydrogen Bonding: The Formic Acid Dimer Case

Formic acid dimer as the prototypical doubly hydrogen-bonded gas-phase species is discussed from the perspective of the three translational and the three rotational degrees of freedom which are lost when two formic acid molecules form a stable complex. The experimental characterisation of these strongly hindered translations and rotations is reviewed, as are attempts to describe the associated fundamental vibrations, their combinations, and their thermal shifts by different electronic structure calculations and vibrational models. A remarkable match is confirmed for the combination of a CCSD(T)-level harmonic treatment and an MP2-level anharmonic VPT2 correction. Qualitatively correct thermal shifts of the vibrational spectra can be obtained from classical molecular dynamics in CCSD(T)-quality force fields. A detailed analysis suggests that this agreement between experiment and composite theoretical treatment is not strongly affected by fortuitous error cancellation but fully converged variational treatments of the six pair or intermolecular modes and their overtones and combinations in this model system would be welcome.

Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases?

The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.

Teaching vibrational spectra to assign themselves

A new paradigm for assigning vibrational spectra is described. Instead of proceeding from potential to Hamiltonian to eigenvalues/eigenvectors/intensities to spectrum, the new method goes "backwards" directly from spectrum to eigenvectors. The eigenvectors then "assign" the spectrum, in that they identify the basis states that contribute to each eigenvalue. To start, we demonstrate an algorithm that can obtain useful estimates of the eigenvectors connecting a real, symmetric Hamiltonian to its eigenvalues even if the only available information about the Hamiltonian is its diagonal elements. When this algorithm is augmented with information about transition intensities, it can be used to assign a complex vibrational spectrum using only information about (1) eigenvalues (the peak centers of the spectrum) and (2) a harmonic basis set (taken to be the diagonal elements of the Hamiltonian). Examples will be discussed, including application to the vibrationally complex spectral region of the formic acid dimer.

IR Spectra of (HCOOH)2 and (DCOOH)2: Experiment, VSCF/VCI, and Ab Initio Molecular Dynamics Calculations Using Full-Dimensional Potential and Dipole Moment Surfaces

We report quantum VSCF/VCI and ab initio molecular dynamics (AIMD) calculations of the IR spectra of (HCOOH)₂ and (DCOOH)₂, using full-dimensional, ab initio potential energy and dipole moment surfaces (PES and DMS). These surfaces are fits, using permutationally invariant polynomials, to 13 475 ab initio CCSD(T)-F12a electronic energies and MP2 dipole moments. Here "AIMD" means using these ab initio potential and dipole moment surfaces in the MD calculations. The VSCF/VCI calculations use all (24) normal modes for coupling, with a four-mode representation of the potential. The quantum spectra align well with jet-cooled and room-temperature experimental spectra over the spectral range 600-3600 cm^-1. Analyses of the complex O-H and C-H stretch bands are made based on the mixing of the VSCF/VCI basis functions. The comparisons of the AIMD IR spectra with both experimental and VSCF/VCI ones provide tests of the accuracy of the AIMD approach. These indicate good accuracy for simple bands but not for the complex O-H stretch band, which is upshifted from experimental and VSCF/VCI bands by roughly 300 cm^-1. In addition to testing the AIMD approach, the PES, DMS, and VSCF/VCI calculations for formic acid dimer provide opportunities for testing other methods to represent high-dimensional data and other methods that perform postharmonic vibrational calculations.

Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces

Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)₂, and comparison with classical MD one.

An ab initio potential energy surface for the formic acid dimer: zero-point energy, selected anharmonic fundamental energies, and ground-state tunneling splitting calculated in relaxed 1–4-mode subspaces

We report a full-dimensional, permutationally invariant potential energy surface (PES) for the cyclic formic acid dimer.

Vibrational dynamics of carboxylic acid dimers in gas and dilute solution

Ultrafast mid-IR transient absorption spectroscopy has been used to study the vibrational dynamics of hydrogen-bonded cyclic dimers of trifluoroacetic acid and formic acid in both the gas and solution phases (0.05 M in CCl(4)). Ultrafast excitation of the broad O-H cyclic dimer band leads, in the gas phase, to large-scale structural changes of the dimer creating a species with a distinct free O-H stretching band on 20 ps and 200 ps timescales. These timescales are assigned to ring-opening and dissociation of the dimer, respectively. In the solution phase, no such structural rearrangement occurs and our results are consistent with previous studies. The gas phase dynamics are insensitive to both the specific excitation energy (over a span of 550 cm(-1)) and the chemical identity of the dimer.

Raman jet spectroscopy of formic acid dimers: low frequency vibrational dynamics and beyond

The vibrational dynamics of formic acid dimer is quite regular at low fundamental excitation frequencies, whereas it evolves into a complex and irregular vibrational signature in the OH stretching region. This is evidenced by the first Raman investigation of the jet-cooled formic acid dimer and its three deuterated isotopomers. Subtle isotope effects in the inter-monomer stretching mode, which is directly observed for the first time at 194 cm(-1), find an interpretation based on hydrogen bond weakening due to quantum delocalization of the protons. The reported high-frequency jet spectra should provide essential experimental stepping stones towards a more complete understanding of this planar prototype for strong double hydrogen bonding.

Infrared spectra of the CF3I dimer: A concurrent application of matrix-isolation spectroscopy and cavity ring-down spectroscopy

We have observed infrared spectra of the CF(3)I dimer produced in a supersonic jet by matrix-isolation Fourier transform infrared spectroscopy and infrared cavity ring-down (IR-CRD) spectroscopy. In the matrix-isolation experiments, the dimer was isolated in an Ar matrix by the pulse-deposition method. The recorded spectral range covers the symmetric (nu(1)) and doubly degenerate (nu(4)) C-F stretching regions. From the concentration dependence of the matrix-isolation spectra we have assigned one dimer band for each fundamental region. It was not easy to identify the dimer band for the nu(4) band because of the multiplet feature of the monomeric nu(4) band caused by the site symmetry breaking. The spectra of (CF(3)I)(2) in the nu(4) band region were thus also measured in the gas phase by IR-CRD spectroscopy, where we detected two dimer bands. Comparing the observed band positions with the results of quantum chemical calculations, we have assigned the observed dimer bands to the head-to-head isomer. The structure of (CF(3)I)(2) and its photochemical implications are discussed, in comparison with methyl iodide dimer reported previously [Ito et al., Chem. Phys. Lett. 343, 185 (2001)].

Methyl iodide clusters observed in gas phase by infrared cavity ring-down spectroscopy: The CH3 bending mode at 8μm

Infrared spectra of methyl iodide clusters produced in a supersonic jet have been observed in the CH3 bending region at 8 mum by cavity ring-down spectroscopy. The dependence of the spectral features on the mixing ratio of CH3I to He and on the stagnation pressure has allowed us to assign the absorption peaks, with the help of the previous results obtained by matrix-isolation technique [Ito et al., Chem. Phys. Lett. 343, 185 (2001)] and infrared cavity ring-down spectroscopy in the C-H stretching region [Ito et al., Chem. Phys. 286, 337 (2003)]. Ab initio calculations at the MP2 level have been carried out up to tetramer to confirm the assignments. It has been found that the frequency shifts upon clustering (relative to monomer) observed in the bending region are not monotonic, in contrast to those in the C-H stretching region. The observed frequency shifts are discussed in terms of dispersion interaction and its variation upon vibrational excitation.

Ground and Excited States of the Monomer and Dimer of Certain Carboxylic Acids

The ground-state properties of the monomer and the dimer of formic acid, acetic acid, and benzoic acid have been investigated using Hartree-Fock (HF) and density functional theory (DFT) methods using the 6-311++G(d,p) basis set. Some of the low-lying excited states have been studied using the time-dependent density functional theory (TDDFT) with LDA and B3LYP functionals and also employing complete-active-space-self-consistent-field (CASSCF) and multireference configuration interaction (MRCI) methodologies. DFT calculations predict the ground-state geometries in quantitative agreement with the available experimental results. The computed binding energies for the three carboxylic acid dimers are also in accord with the known thermodynamic data. The TDDFT predicted wavelengths corresponding to the lowest energy n-pi* transition in formic acid (214 nm) and acetic acid (214 nm) and the pi-pi* transition in benzoic acid (255 nm) are comparable to the experimentally observed absorption maxima. In addition, TDDFT calculations predict qualitatively correctly the blue shift (4-5 nm) in the excitation energy for the pi-pi* transition in going from the monomer to the dimer of formic acid and acetic acid and the red shift (approximately 19 nm) in pi-pi* transition in going from benzoic acid monomer to dimer. This also indicates that the electronic interaction arising from the hydrogen bonds between the monomers is marginal in all three carboxylic acids investigated.

Hydrogen bond-induced vibronic mode mixing in benzoic acid dimer: A laser-induced fluorescence study

Laser-induced dispersed fluorescence spectra of benzoic acid dimer in the cold environment of supersonic jet expansion have been reinvestigated with improved spectral resolution of measurements. The spectra are analyzed with the aid of the normal mode vibrations of the dimer calculated by the ab initio quantum chemistry method at the DFT/B3LYP/6-311+G(*) (*) level of theory. The analysis reveals that the low-frequency intermolecular hydrogen bond modes are mixed extensively with the carboxyl as well as aromatic ring vibrations upon electronic excitation. The mode mixing is manifested as the complete loss of mirror symmetry relation between the fluorescence excitation and dispersed fluorescence spectra of the S(1) origin, and appearance of large number of cross-sequence transitions when the DF spectra are measured by exciting the low-energy vibrations near the S(1) origin. The cross-sequence bands are found in all the cases to be the combinations of two nontotally symmetric fundamentals consisting of one of the intermolecular hydrogen bond modes and the other from the aromatic ring and carboxyl group vibrations. The implications of this mode mixing on the excited state dynamics of the dimer are discussed.

Infrared depletion spectra of 2-aminopyridine⋅2-pyridone, a Watson–Crick mimic of adenine⋅uracil

The 2-aminopyridine2-pyridone (2AP2PY) dimer is linked by N-H...O=C and N-H...N hydrogen bonds, providing a model for the Watson-Crick hydrogen bond configuration of the adenine.thymine and adenine.uracil nucleobase pairs. Mass-specific infrared spectra of 2AP2PY and its seven N-H deuterated isotopomers have been measured between 2550 and 3650 cm(-1) by IR laser depletion combined with UV two-color resonant two-photon ionization. The 2PY amide N-H stretch is a very intense band spread over the range 2700-3000 cm(-1) due to large anharmonic couplings. It is shifted to lower frequency by 710 cm(-1) or approximately 20% upon H bonding to 2AP. On the 2AP moiety, the "bound" amino N-H stretch gives rise to a sharp band at 3140 cm(-1), which is downshifted by 354 cm(-1) or approximately 10% upon H bonding to 2PY. The amino group "free" N-H stretch and the H-N-H bend overtone are sharp bands at approximately 3530 cm(-1) and 3320 cm(-1). Ab initio structures and harmonic vibrations were calculated at the Hartree-Fock level and with the PW91 and B3LYP density functionals. The PW91/6-311++G(d,p) method provides excellent predictions for the frequencies and IR intensities of all the isotopomers.