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

Formic acid stability in different solvents by DFT calculations

CONTEXT

New technologies have been developed toward the use of green energies. The production of formic acid (FA) from carbon dioxide (CO[Formula: see text]) hydrogenation with H[Formula: see text] is a sustainable process for H[Formula: see text] storage. However, the FA adduct stabilization is thermodynamically dependent on the type of solvent and thermodynamic conditions. The results suggest a wide range of dielectric permittivity values between the dimethyl sulfoxide (DMSO) and water solvents to stabilize the FA in the absence of base. The thermodynamics analysis and the infrared and charge density difference results show that the formation of the FA complex with H[Formula: see text]O is temperature dependent and has a major influence on aqueous solvents compared to the FA adduct with amine, in good agreement with the experiment. In these conditions, the stability thermodynamic of the FA molecule may be favorable at non-organic solvents and dielectric permittivity values closer to water. Therefore, a mixture of aqueous solvents with possible ionic composition could be used to increase the thermodynamic stability of H[Formula: see text] storage in CO[Formula: see text] conversion processes.

METHODS

Using the Quantum ESPRESSO package, density functional theory (DFT) calculations were performed with periodic boundary conditions, and the electronic wave functions were expanded in plane waves. For the exchange-correlation functional, we use the vdW-DF functional with the inclusion of van der Waals (vdW) forces. Electron-ion interactions are treated by the projector augmented wave (PAW) method with pseudopotentials available in the PSlibrary repository. The wave functions and the electronic densities were expanded employing accurate cut-off energies of 6.80[Formula: see text]10[Formula: see text] and 5.44[Formula: see text]10[Formula: see text] eV, respectively. The electronic density was computed from the wave functions calculated at the [Formula: see text]-point in the first Brillouin-zone. Each structural optimization was minimized according to the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm, with force and energy convergence criteria of 25 meV[Formula: see text]Å[Formula: see text] and 1.36 meV, respectively. The electrostatic solvation effects were performed by the [Formula: see text] package with the Self-Consistent Continuum Solvation (SCCS) approach.

Vibrationally correlated calculations in polyspherical coordinates: Taylor expansion-based kinetic energy operators

The efficiency of quantum chemical simulations of nuclear motion can in many cases greatly benefit from the application of curvilinear coordinate systems. This is rooted in the fact that a set of smartly selected curvilinear coordinates may represent the motion naturally well, thus decreasing the couplings between motions in these coordinates. In this study, we assess the validity of different Taylor expansion-based approximations of kinetic energy operators in a (curvilinear) polyspherical parametrization. To this end, we investigate the accuracy as well as the numerical performance of the approximations in time-independent vibrational coupled cluster and full vibrational interaction calculations for several test cases ranging from tri- to penta-atomic molecules. We find that several of the proposed schemes reproduce the vibrational ground state and excitation energies to a decent accuracy, justifying their application in future investigations. Furthermore, due to the restricted mode coupling and their inherent sum-of-products form, the new approximations open up the possibility of treating large molecular systems with efficient vibrational coupled cluster schemes in general coordinates.

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.

On the vibrations of formic acid predicted from first principles

In this article, we review recent first principles, anharmonic studies on the molecular vibrations of gaseous formic acid in its monomer form. Transitions identified as fundamentals for both cis- and trans form reported in these studies are collected and supported by results from high-resolution experiments. Attention is given to the effect of coordinate coupling on the convergence of the computed vibrational states.

Intramolecular Vibrational Redistribution in Formic Acid and its Deuterated Forms

The intramolecular vibrational relaxation dynamics of formic acid and its deuterated isotopologues is simulated on the full-dimensional potential energy surface of Richter and Carbonnière [F. Richter and P. Carbonnière, J. Chem. Phys. 148, 064303 (2018)] using the Heidelberg MCTDH package. We focus on couplings with the torsion vibrational modes close to the trans- cis isomerisation coordinate from the dynamics of artificially excited vibrational mode overtones. The C-O stretch bright vibrational mode is coupled to the out-of-the plane torsion mode in HCOOH, where this coupling could be exploited for laser-induced trans-to- cis isomerisation. Strong isotopic effects are observed: deuteration of the hydroxyl group, i.e., in HCOOD and DCOOD, destroys the C-O stretch to torsion mode coupling whereas in DCOOH, little to no effect is observed.

Adaptive Fitting of Potential Energy Surfaces of Small to Medium-Sized Molecules in Sum-of-Product Form: Application to Vibrational Spectroscopy

A semi-automatic sampling and fitting procedure for generating sum-of-product (Born-Oppenheimer) potential energy surfaces based on a high-dimensional model representation is presented. The adaptive sampling procedure and subsequent fitting relies on energies only and can be used for re-fitting existing analytic potential energy surfaces in sum-of-product form or for direct fits from ab initio computa- tions. The method is tested by fitting ground electronic state potential energy surfaces for small to medium sized semi-rigid molecules, i.e., HFCO, HONO, and HCOOH, based upon ab initio computations at the CCSD(T)-F12/cc-pVTZ-F12 or MP2/aug-cc-pVTZ levels of theory. Vibrational eigenstates are computed using block improved relaxation in the Heidelberg MCTDH package and compared to available experimental and theoretical data. The new potential energy surfaces are compared to the best ones currently available for these molecules, in terms of accuracy, including of resulting vibrational states, required numbers of sampling points, and number of fitting parameters. The present procedure leads to compact expansions and scales well with the number of dimensions for simple potentials such as single or double wells.

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.

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.

Quantifying Conformational Isomerism in Chain Molecules by Linear Raman Spectroscopy: The Case of Methyl Esters

The conformational preferences of the ester group have the potential to facilitate the large amplitude folding of long alkyl chains in the gas phase. They are monitored by Raman spectroscopy in supersonic jet expansions for the model system methyl butanoate, after establishing a quantitative relationship with quantum-chemical predictions for methyl methanoate. This requires a careful analysis of experimental details, and a simulation of the rovibrational contours for near-symmetric top molecules. The technique is shown to be complementary to microwave spectroscopy in quantifying coexisting conformations. It confirms that a C-O-C(=O)-C-C chain segment can be collapsed into a single all-trans conformation by collisional cooling, whereas alkyl chain isomerism beyond this five-membered chain largely survives the jet expansion. This sets the stage for the investigation of linear alkyl alkanoates in terms of dispersion-induced stretched-chain to hairpin transitions by Raman spectroscopy.