Communication: The highest frequency hydrogen bond vibration and an experimental value for the dissociation energy of formic acid dimer

Spectroscopic Evidence for Doubly Hydrogen-Bonded Cationic Dimers in the Solid, Liquid, and Gaseous Phases of Carboxyl-Functionalized Ionic Liquids

Intermolecular interactions determine whether matter sticks together, gases condense into liquids, or liquids freeze into solids. The most prominent example is hydrogen bonding in water, responsible for the anomalous properties in the liquid phase and polymorphism in ice. The physical properties are also exceptional for ionic liquids (ILs), wherein a delicate balance of Coulomb interactions, hydrogen bonds, and dispersion interactions results in a broad liquid range and the vaporization of ILs as ion pairs. In this study, we show that strong, local, and directional hydrogen bonds govern the structures and arrangements in the solid, liquid, and gaseous phases of carboxyl-functionalized ILs. For that purpose, we explored the H-bonded motifs by X-ray diffraction and attenuated total reflection (ATR) infrared (IR) spectroscopy in the solid state, by ATR and transmission IR spectroscopy in the liquid phase, and by cryogenic ion vibrational predissociation spectroscopy (CIVPS) in the gaseous phase at low temperature. The analysis of the CO stretching bands reveals doubly hydrogen-bonded cationic dimers (c═c), resembling the archetype H-bond motif known for carboxylic acids. The like-charge doubly hydrogen-bonded ion pairs are present in the crystal structure of the IL, survive phase transition into the liquid state, and are still present in the gaseous phase even in (2,1) complexes wherein one counterion is removed and repulsive Coulomb interaction increased. The interpretation of the vibrational spectra is supported by quantum chemical methods. These observations have implications for the fundamental nature of the hydrogen bond between ions of like charge.

Equilibrium structure of a dense trimesic acid monolayer on a homogeneous solid surface: from atomistic simulation to thermodynamics

A general methodology for determining the thermodynamic characteristics of rigid organic crystals on the atomistic level is presented. The proposed approach is based on a combination of grid interpolation of the precalculated intermolecular potential and kinetic Monte Carlo simulation of the gas-crystal system with an explicit interphase. The two-phase system is stabilized in a wide range of external parameters with an imposed external potential and damping field. The damping field reduces the intermolecular potential at the edges of the crystals and turns it off in the gas phase. To determine the thermodynamic characteristics of a crystal the conditions of equality of chemical potentials in coexisting phases are used. The intermolecular pairwise potential can be calculated on the atomistic or quantum level. In the kinetic Monte Carlo simulations, a grid interpolation of the precalculated potential is performed on each iteration of the algorithm. We have applied the approach to the thermodynamic analysis of a dense monolayer of trimesic acid on a homogeneous surface. The calculated free energy and entropy for the dense "superflower" and filled chicken-wire phases obey the Gibbs-Duhem equation, which confirms the thermodynamic consistency of our approach. Using the proposed approach, we have revealed that the dense "superflower" phase becomes metastable at zero pressure and 470-500 K. Under these conditions, the filled chicken-wire structure with partially released hexagonal cages is thermodynamically favourable. The proposed approach is a potentially universal tool for the thermodynamic analysis of crystals formed by "rigid" organic molecules of any complexity on the atomistic level.

Formic acid-methanol complexation vs. esterification: elusive pre-reactive species identified by vibrational spectroscopy

Vibrational spectra of the mixed dimer and the two mixed trimers of methanol and formic acid as well as some of their isotopologues are presented. Out of the eight expected OH stretching fundamentals of these three pre-reactive hydrogen-bonded complexes, the three modes mainly involving an acid OH group bound to the alcohol appear to be missing in the jet-cooled spectra despite a combination of infrared and Raman probing. A possibility of spectral overlap is discussed in the mixed dimer case, but largely discarded. The missing modes correspond to (fractional) concerted elongation of all engaged OH bonds, promoting synchronous degenerate proton transfer between the molecules. One other trimer mode is very tentatively attributed to a broad spectral feature, whereas all OH bonds contacting carbonyl groups can be unambiguously identified by four relatively narrow infrared absorptions. The spectral features are confirmed by vibrational perturbation theory and deviate in a subtle but systematic way from scaled harmonic predictions which were previously validated for the formic acid complex with a more acidic alcohol. Despite being exothermic and exergonic, ester formation can only be detected in the rarefied gas expansions after extended pre-mixture of the gases, which somewhat contrasts the recent microwave spectroscopic evidence of in situ ester formation and in particular the lack of pre-reactive complex signals.

Mid-infrared Laser Spectroscopy of Jet-Cooled Formic Acid Trimer: Mode-Dependent Line Broadening in the C-O Stretching Region

Building on recent progress in the vibrational spectroscopy of the formic acid trimer, we present the first high-resolution measurements of the jet-cooled laser absorption spectrum of (HCOOH)3. The spectra of the lowest- and highest-frequency C-O stretching fundamentals are analyzed whereas the third band is not observed, complicated by monomer and dimer absorptions at 1219 cm-1 (8.2 μm). Vibration-rotation parameters are obtained for the band at 1172.31512(68) cm-1 whereas the C-O stretch at 1246.33(5) cm-1 exhibits a significantly larger breadth, allowing only resolution of the coarse PQR structure. Vibrational predissociation can be ruled out, and intramolecular vibrational redistribution mechanisms are discussed, particularly coupling to the concerted proton exchange within the cyclic dimer subunit. Ultimately, the question remains open. The prospects of high-resolution measurements of other trimer bands or isotope substitution experiments, which might assist in revealing the mode-specificity of the underlying broadening mechanisms, are discussed.

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.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Double proton transfer in hydrated formic acid dimer: Interplay of spatial symmetry and solvent-generated force on reactivity

The double proton transfer (DPT) reaction in the hydrated formic acid dimer (FAD) is investigated at molecular-level detail. For this, a global and reactive machine learned (ML) potential energy surface (PES) is developed to run extensive (more than 100 ns) mixed ML/MM molecular dynamics (MD) simulations in explicit molecular mechanics (MM) solvent at MP2-quality for the solute. Simulations with fixed – as in a conventional empirical force field – and conformationally fluctuating – as available from the ML-based PES – charge models for FAD show a significant impact on the competition between DPT and dissociation of FAD into two formic acid monomers. With increasing temperature the barrier height for DPT in solution changes by about 10% (∼1 kcal mol⁻¹) between 300 K and 600 K. The rate for DPT is largest, ∼1 ns⁻¹, at 350 K and decreases for higher temperatures due to destabilisation and increased probability for dissociation of FAD. The water solvent is found to promote the first proton transfer by exerting a favourable solvent-induced Coulomb force along the O–H⋯O hydrogen bond whereas the second proton transfer is significantly controlled by the O–O separation and other conformational degrees of freedom. Double proton transfer in hydrated FAD is found to involve a subtle interplay and balance between structural and electrostatic factors.

Simulation of double proton transfer in formic acid dimer by reactive ML potential in explicit molecular mechanics water solvent.

Thermodynamically Stable Cationic Dimers in Carboxyl-Functionalized Ionic Liquids: The Paradoxical Case of "Anti-Electrostatic" Hydrogen Bonding

We show that carboxyl-functionalized ionic liquids (ILs) form doubly hydrogen-bonded cationic dimers (c⁺=c⁺) despite the repulsive forces between ions of like charge and competing hydrogen bonds between cation and anion (c⁺–a⁻). This structural motif as known for formic acid, the archetype of double hydrogen bridges, is present in the solid state of the IL 1−(carboxymethyl)pyridinium bis(trifluoromethylsulfonyl)imide [HOOC−CH₂−py][NTf₂]. By means of quantum chemical calculations, we explored different hydrogen-bonded isomers of neutral (HOOC–(CH₂)_n–py⁺)₂(NTf₂⁻)₂, single-charged (HOOC–(CH₂)_n–py⁺)₂(NTf₂⁻), and double-charged (HOOC– (CH₂)_n−py⁺)₂ complexes for demonstrating the paradoxical case of “anti-electrostatic” hydrogen bonding (AEHB) between ions of like charge. For the pure doubly hydrogen-bonded cationic dimers (HOOC– (CH₂)_n−py⁺)₂, we report robust kinetic stability for n = 1–4. At n = 5, hydrogen bonding and dispersion fully compensate for the repulsive Coulomb forces between the cations, allowing for the quantification of the two equivalent hydrogen bonds and dispersion interaction in the order of 58.5 and 11 kJmol⁻¹, respectively. For n = 6–8, we calculated negative free energies for temperatures below 47, 80, and 114 K, respectively. Quantum cluster equilibrium (QCE) theory predicts the equilibria between cationic monomers and dimers by considering the intermolecular interaction between the species, leading to thermodynamic stability at even higher temperatures. We rationalize the H-bond characteristics of the cationic dimers by the natural bond orbital (NBO) approach, emphasizing the strong correlation between NBO-based and spectroscopic descriptors, such as NMR chemical shifts and vibrational frequencies.

Thermodynamics of self-assembled molecular layers of trimesic acid from fields-supported kinetic Monte Carlo simulation

A technique has been developed for calculating the thermodynamic characteristics of rigid self-assembled organic adsorption layers and the parameters of polymorphic transitions using two types of external fields and the kinetic Monte Carlo method.

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.

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol^-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.

Double Proton Transfer Across a Table: The Formic Acid Dimer-Fluorobenzene Complex

Proton transfer via tunneling is a fundamental quantum‐mechanical phenomenon. We report rotational spectroscopy measurements of this process in the complex of the formic acid dimer with fluorobenzene. The assignment of the spectrum indicates that this complex exists in the form of a π–π stacked structure. Each rotational transition of the parent isotopologue exhibits splitting. Isotopic substitution experiments show that the spectral splitting results from double‐proton transfer tunneling in the formic acid dimer. Presence of fluorobenzene as a neighboring molecule does not quench the double proton transfer in the formic acid dimer but decreases its tunneling splitting from 341(3) MHz to 267.608(1) MHz. Calculations suggest that the presence of the weakly bounded fluorobenzene does not influence the activation energy of the proton transfer. The fluorobenzene is reoriented with respect to the formic acid dimer during the course of the reaction, slowing down the proton transfer motion.

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.

CC-stretched formic acid: isomerisation, dimerisation, and carboxylic acid complexation

The cis-trans-isomerism of the propiolic acid monomer (HC[triple bond, length as m-dash]C-COOH) is examined with linear Raman jet spectroscopy, yielding the first environment-free vibrational band centres of a higher-energy cis-rotamer beyond formic acid (HCOOH) in addition to all fundamentals and a large number of hot and combination/overtone bands of the trans-conformer. Two near-isoenergetic trans-fundamentals of different symmetry (CC[double bond, length as m-dash]O bend and OH torsion) prove to be a sensitive benchmarking target, as their energetic order is susceptible to the choice of electronic structure method, basis set size, and inclusion of vibrational anharmonicity. For the infrared- and Raman-active C[double bond, length as m-dash]O stretching fundamentals of the cyclic (C_2h) trans-propiolic acid dimer, resonance couplings are found that in part extend to the C_s-symmetric heterodimer of trans-propiolic and trans-formic acid. Exploratory vibrational perturbation theory (VPT2) calculations show that all perturbing states involve displacements of the OH moieties located on the doubly hydrogen bonded ring. The comparison of the infrared spectra of the propiolic acid dimer and its heterodimer with formic acid to that of several other carboxylic acid dimers from the literature reveals a notable similarity regarding a non-fundamental dimer band around 1800 cm^-1, which in most cases is so far unassigned. VPT2 calculations and a simple harmonic model suggest an assignment to a combination vibration of the symmetric and antisymmetric OH torsion.

Transfer Learning to CCSD(T): Accurate Anharmonic Frequencies from Machine Learning Models

The calculation of the anharmonic modes of small- to medium-sized molecules for assigning experimentally measured frequencies to the corresponding type of molecular motions is computationally challenging at sufficiently high levels of quantum chemical theory. Here, a practical and affordable way to calculate coupled-cluster quality anharmonic frequencies using second-order vibrational perturbation theory (VPT2) from machine-learned models is presented. The approach, referenced as "NN + VPT2", uses a high-dimensional neural network (PhysNet) to learn potential energy surfaces (PESs) at different levels of theory from which harmonic and VPT2 frequencies can be efficiently determined. The NN + VPT2 approach is applied to eight small- to medium-sized molecules (H₂CO, trans-HONO, HCOOH, CH₃OH, CH₃CHO, CH₃NO₂, CH₃COOH, and CH₃CONH₂) and frequencies are reported from NN-learned models at the MP2/aug-cc-pVTZ, CCSD(T)/aug-cc-pVTZ, and CCSD(T)-F12/aug-cc-pVTZ-F12 levels of theory. For the largest molecules and at the highest levels of theory, transfer learning (TL) is used to determine the necessary full-dimensional, near-equilibrium PESs. Overall, NN + VPT2 yields anharmonic frequencies to within 20 cm^-1 of experimentally determined frequencies for close to 90% of the modes for the highest quality PES available and to within 10 cm^-1 for more than 60% of the modes. For the MP2 PESs only ∼60% of the NN + VPT2 frequencies were within 20 cm^-1 of the experiment, with outliers up to ∼150 cm^-1, compared to the experiment. It is also demonstrated that the approach allows to provide correct assignments for strongly interacting modes such as the OH bending and the OH torsional modes in formic acid monomer and the CO-stretch and OH-bend mode in acetic acid.

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

Curvilinear kinetic energy models are developed for variational nuclear motion computations including the inter- and the low-frequency intra-molecular degrees of freedom of the formic acid dimer. The coupling of the inter- and intra-molecular modes is studied by solving the vibrational Schrödinger equation for a series of vibrational models, from two up to ten active vibrational degrees of freedom by selecting various combinations of active modes and constrained coordinate values. Vibrational states, nodal assignment, and infrared vibrational intensity information is computed using the full-dimensional potential energy surface (PES) and electric dipole moment surface developed by Qu and Bowman [Phys. Chem. Chem. Phys., 2016, 18, 24835; J. Chem. Phys., 2018, 148, 241713]. Good results are obtained for several fundamental and combination bands in comparison with jet-cooled vibrational spectroscopy experiments, but the description of the ν₈ and ν₉ fundamental vibrations, which are close in energy and have the same symmetry, appears to be problematic. For further progress in comparison with experiment, the potential energy surface, and in particular, its multi-dimensional couplings representation, requires further improvement.

Understanding dispersion interactions in molecular chemistry

This themed collection contains a selection of articles on the topic of Understanding Dispersion Interactions in Molecular Chemistry.

A new instrument for kinetics and branching ratio studies of gas phase collisional processes at very low temperatures

A new instrument dedicated to the kinetic study of low-temperature gas phase neutral-neutral reactions, including clustering processes, is presented. It combines a supersonic flow reactor with vacuum ultra-violet synchrotron photoionization time-of-flight mass spectrometry. A photoion-photoelectron coincidence detection scheme has been adopted to optimize the particle counting efficiency. The characteristics of the instrument are detailed along with its capabilities illustrated through a few results obtained at low temperatures (<100 K) including a photoionization spectrum of n-butane, the detection of formic acid dimer formation, and the observation of diacetylene molecules formed by the reaction between the C₂H radical and C₂H₂.

Increasing the weights in the molecular work-out of cis- and trans-formic acid: extension of the vibrational database via deuteration

The higher-energy cis- as well as the global minimum trans-rotamers of the four H/D isotopologues of the formic acid monomer have been examined with Raman jet spectroscopy extending the vibrational gas phase reference database by eleven new cis-band positions for HCOOD, DCOOH, and DCOOD. With these new additions, all O-H/D, C-H/D, and C[double bond, length as m-dash]O stretching as well as the O-D in-plane bending vibrations of these higher-energy rotamers are known in addition to the previously determined C-O stretch and OH torsion of cis-HCOOH. Further, a comparison of the vibrational spectra of all four H/D isotopologues of the globally stable trans-rotamer of formic acid is shown to be very helpful in revealing similarities and differences in these systems, particularly with regard to Fermi resonances. Amongst the most prominent ones is the ν5/2ν9 resonance doublet of trans-HCOOH, for which we provide more insight into a recently suggested label switch of the resonance partners via the comparison of infrared and Raman jet spectra.

Carboxylation Enhances Fragmentation of Furan upon Resonant Electron Attachment

We report a dissociative electron attachment study to 2-furoic acid (C5H4O3) isolated in a gas phase, which is a model molecule consisting of a carboxylic group and a furan ring. Dissociation of furan by low energy electrons is accessible only via electronic excited Feshbach resonances at energies of incident electrons above 5 eV. On the other hand, carboxylic acids are well-known to dissociate via attachment of electrons at subexcitation energies. Here we elucidate how the electron and proton transfer reactions induced by carboxylation influence stability of the furan ring. Overlap of the furan and carboxyl π orbitals results in transformation of the nondissociative π2 resonance of the furan ring to a dissociative resonance. The interpretation of hydrogen transfer reactions is supported by experimental studies of 3-methyl-2-furoic and 5-methyl-2-furoic acids (C6H6O3) and density functional theory (DFT) calculations.

Stretching force constants as descriptors of energy and geometry of F···HF hydrogen bonds

In this work applicability of proton donor group stretching vibration force constants k_s and intermolecular stretching force constants k_σ for evaluations of hydrogen bond strength and geometry are discussed. For a set of 30 complexes with F···HF hydrogen bonds in a wide range 0.5-48 kcal/mol by means of quantum chemical calculations equilibrium geometries, complexation energies, vibrational frequencies and corresponding force constants were calculated (MP2/aug-cc-pVTZ). It is shown, that properties of a hydrogen bond are more strictly correlated with the values of force constants than with vibrational frequencies. Easy-to-use equations for estimations of hydrogen bond energy ∆E and geometry (r_FH, r_FF) based on k_s and k_σ values are proposed.

Proton Transfer vs Complex Formation Channels in Ionized Formic Acid Dimer: A Direct Ab Initio Molecular Dynamics Study

Photoirradiation to a hydrogen-bonded system plays an important role in the initial DNA and enzyme damage processes. The formic acid (FA) dimer is a model compound of double proton transfer systems, such as DNA base pairs. In the present study, the reactions of the FA dimer cation, formed upon ionization of the neutral dimer, have been investigated by the direct ab initio molecular dynamics method. Two reaction channels were identified for the FA dimer cation: complex formation and proton transfer (PT). In the complex formation channel, the carbonyl oxygen atoms of the two FA monomers were bound symmetrically, and a face-to-face complex was formed. In the PT channel, the proton of FA⁺ was transferred to FA, forming the H⁺(HCOOH)--HCO₂ radical cation as product. At low temperature, the complex channel was dominant, whereas the PT channel increased with increasing temperature. The asymmetric spin distribution on the FA dimer cation exhibited a strong correlation with the PT channel.

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.

High-Resolution Infrared Synchrotron Investigation of (HCN)2 and a Semi-Experimental Determination of the Dissociation Energy D0

The high‐resolution infrared absorption spectrum of the donor bending fundamental band ν61 of the homodimer (HCN)₂ has been collected by long‐path static gas‐phase Fourier transform spectroscopy at 207 K employing the highly brilliant 2.75 GeV electron storage ring source at Synchrotron SOLEIL. The rovibrational structure of the ν61 transition has the typical appearance of a perpendicular type band associated with a Σ–Π transition for a linear polyatomic molecule. The total number of 100 assigned transitions are fitted employing a standard semi‐rigid linear molecule Hamiltonian, providing the band origin ν₀ of 779.05182(50) cm⁻¹ together with spectroscopic parameters for the degenerate excited state. This band origin, blue‐shifted by 67.15 cm⁻¹ relative to the HCN monomer, provides the final significant contribution to the change of intra‐molecular vibrational zero‐point energy upon HCN dimerization. The combination with the vibrational zero‐point energy contribution determined recently for the class of large‐amplitude inter‐molecular fundamental transitions then enables a complete determination of the total change of vibrational zero‐point energy of 3.35±0.30 kJ mol⁻¹. The new spectroscopic findings together with previously reported benchmark CCSDT(Q)/CBS electronic energies [Hoobler et al. ChemPhysChem. 19, 3257–3265 (2018)] provide the best semi‐experimental estimate of 16.48±0.30 kJ mol⁻¹ for the dissociation energy D₀ of this prototypical homodimer.

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.

Soft experimental constraints for soft interactions: a spectroscopic benchmark data set for weak and strong hydrogen bonds

An experimental benchmark data base on rotational constants, vibrational properties and energy differences for weakly and more strongly hydrogen-bonded complexes and their constituents from the spectroscopic literature is assembled. It is characterized in detail and finally contracted to a more compact, discriminatory set (ENCH-51, for Experimental Non-Covalent Harmonic with 51 entries). The meeting points between theory and experiment consist of equilibrium rotational constants and harmonic frequencies and energies, which are back-corrected from experimental observables and are very easily accessible by quantum chemical calculations. The relative performance of B3LYP-D3, PBE0-D3 and M06-2X density functional theory predictions with a quadruple-zeta basis set is used to illustrate systematic errors, error compensation and selective performance for structural, vibrational and energetical observables. The current focus is on perspectives and different benchmarking methodologies, rather than on a specific theoretical method or a specific class of compounds. Extension of the data base in chemical, observable and quantum chemical method space is encouraged.

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.

THz spectroscopy of weakly bound cluster molecules in solid para-hydrogen: a sensitive probe of van der Waals interactions

The present work demonstrates that 99.9% enriched solid para-H2 below 3 K provides an excellent inert and transparent medium for the exploration of large-amplitude intermolecular vibrational motion of weakly bound van der Waals cluster molecules in the THz spectral region. THz absorption spectra have been generated for CO2/H2O and CS2/H2O mixtures embedded in enriched solid para-H2 and numerous observed transitions associated with large-amplitude librational motion of the weakly bound binary CO2H2O and CS2H2O van der Waals cluster molecules have been assigned together with tentative assignments for the ternary CS2(H2O)2 system. The interaction strength, directionality and anharmonicity of the weak van der Waals "bonds" between the molecules can be characterized via these THz spectral signatures and yield rigorous benchmarks for high-level ab initio methodologies. It is suggested that even a less stable linear conformation of the ternary CS2(H2O)2 system, where one H2O molecule is linked to each S atom of the CS2 subunit, may be formed due to the kinetics associated with the mobility of free H2O molecules in the soft para-H2 medium. In addition, the spectroscopic observations confirm a linear and planar global intermolecular potential energy minimum for the binary CS2H2O system with C2v symmetry, where the O atom on the H2O molecule is linked to one of the S atoms on the CS2 subunit. A semi-experimental value for the vibrational zero-point energy contribution of 1.93 ± 0.10 kJ mol-1 from the class of large-amplitude intermolecular vibrational modes is proposed. The combination with CCSD(T)/CBS electronic energy predictions provides a semi-experimental estimate of 5.08 ± 0.15 kJ mol-1 for the binding energy D0 of the CS2H2O van der Waals system.

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.

Theory of nanobubble formation and induced force in nanochannels

This paper presents a fundamental theory of nanobubble formation and induced force in confined nanochannels. It is shown that nanobubble formation between hydrophobic plates can be predicted from their surface tension and geometry, with estimated values for the surface free energy and the force acting on the plates in good agreement with the results of molecular dynamics simulation and experimentation. When a bubble is formed between two plates, vertical attractive force and horizontal retract force due to the shifted plates are applied to the plates. The net force exerted on the plates is not dependent on the distance between them. The short-range force between hydrophobic surfaces due to hydrophobic interaction appears to correspond to the force estimated by our theory. We compared between experimental and theoretical values for the binding energy of a molecular motor system to validate our theory. The tendency that the binding energy increases as the size of the protein increases is consistent with the theory.

Quantitative probing of subtle interactions among H-bonds in alpha hydroxy carboxylic acid complexes

The alpha OH stretching frequency may be affected upon complexing with water and formic acid.

Dissection of H-bonding interactions in a glycolic acid–water dimer

The binding strength and collective effects of multiple H-bonds in the glycolic acid–water dimer were studied in comparison to the aromatic analog, 9-hydroxy-9-fluorene carboxylic acid (9HFCA).

Spectroscopy and dynamics of double proton transfer in formic acid dimer

We present the isolated gas phase infrared spectra of formic acid dimer, (HCOOH)2, and its deuterated counterpart formic-d acid, (DCOOH)2, at room temperature. The formic acid dimer spectrum was obtained by spectral subtraction of a spectrum of formic acid vapor recorded at low pressure from that recorded at a higher pressure. The spectra of formic acid vapor contain features from both formic acid monomer and formic acid dimer, but at low and high pressures of formic acid, the equilibrium is pushed towards the monomer and dimer, respectively. A similar approach was used for the formic-d acid dimer. Building on the previous development of the Molecular Mechanics with Proton Transfer (MMPT) force field for simulating proton transfer reactions, molecular dynamics (MD) simulations were carried out to interpret the experimental spectra in the OH-stretching region. Within the framework of MMPT, a combination of symmetric single and double minimum potential energy surfaces (PESs) provides a good description of the double proton transfer PES. In a next step, potential morphing together with electronic structure calculations at the B3LYP and MP2 level of theory was used to align the computed and experimentally observed spectral features in the OH-stretching region. From this analysis, a barrier for double proton transfer between 5 and 7 kcal mol(-1) was derived, which compares with a CCSD(T)/aug-cc-pVTZ calculated barrier of 7.9 kcal mol(-1). Such a combination of experimental and computational techniques for estimating barriers for proton transfer in gas phase systems is generic and holds promise for further improved PESs and energetics of these important systems. Additional MD simulations at the semi-empirical DFTB level of theory agree quite well for the center band position but underestimate the width of the OH-stretching band.

Standard free energy of the equilibrium between the trans-monomer and the cyclic-dimer of acetic acid in the gas phase from infrared spectroscopy

Survey jet-cooled spectra of acetic acid have been recorded in the infrared region (200-4000 cm(-1)) over a wide range of expansion conditions. From the variations of the relative intensities of the signals, vibrational transitions have been assigned unambiguously to the trans-monomer and cyclic-dimer. The IR-active fundamental frequencies have been determined at the instrumental accuracy of 0.5 cm(-1). This analysis of the jet-cooled spectra supported by electronic structure calculations permitted us to characterize the trans-monomer/cyclic-dimer equilibrium. From static cell spectra at 298 K, variations of the molar fractions ratio as a function of the total pressure were used to estimate the equilibrium constant and the Gibbs free energy of dimerization at 298 K. The very good agreement with the literature data shows that the present method is able to produce, from a single study, a free energy value as reliable as the one obtained from a large collection of data. In addition, the semi-empirical free energy value was used to estimate the accuracy of electronic structure calculations and in turn the accuracy of the derived useful information such as the dissociation energy of the complex (i.e. the strength of the hydrogen bonds) or the relative energies within the conformational landscape.