Hydrogen Bonding between Formic Acid and Water: Complete Stabilization of the Intrinsically Unstable Conformer

Study on depolymerization kinetics of formic acid dimers in binary mixture

In this study, polarization Raman spectra were collected for binary mixtures of formic acid/methanol and formic acid/acetonitrile with different volume fractions. The broad band of formic acid in the CO vibration region was divided into four vibration peaks, corresponding to CO symmetric and anti-symmetric stretching vibration from cyclic dimer, CO stretching from open dimer, and CO stretching from the free monomer. The experiments showed that as the volume fraction of formic acid in the binary mixture decreased, the cyclic dimer gradually converted to the open dimer, and at a volume fraction of 0.1, fully depolymerized into monomer form (free monomer, solvated monomer, and hydrogen bonding monomer clusters with solvent). The contribution percentage of the total CO stretching intensity of each structure at different concentrations was quantitatively calculated using high resolution infrared spectroscopy, and the results were consistent with the conclusions predicted by polarization Raman spectroscopy. Concentration-triggered 2D-COS synchronous and asynchronous spectra also confirmed the kinetics of formic acid diluted in acetonitrile. This work provides a spectroscopic method for studying the structure of organic compounds in solution and concentration-triggering kinetics in mixtures.

Controlling the Trimerization of the Collagen Triple-Helix by Solvent Switching

Collagen hybridizing peptides (CHPs) are a powerful tool for targeting collagen damage in pathological tissues due to their ability to specifically form a hybrid collagen triple-helix with the denatured collagen chains. However, CHPs have a strong tendency to self-trimerize, requiring preheating or complicated chemical modifications to dissociate their homotrimers into monomers, which hinders their applications. To control the self-assembly of CHP monomers, we evaluated the effects of 22 cosolvents on the triple-helix structure: unlike typical globular proteins, the CHP homotrimers (as well as the hybrid CHP-collagen triple helix) cannot be destabilized by the hydrophobic alcohols and detergents (e.g., SDS) but can be effectively dissociated by the cosolvents that dominate hydrogen bonds (e.g., urea, guanidinium salts, and hexafluoroisopropanol). Our study provided a reference for the solvent effects on natural collagen and a simple effective solvent-switch method, enabling CHP utilization in automated histopathology staining and in vivo imaging and targeting of collagen damage.

Solvation Effects on Quantum Tunneling Reactions

A decisive factor for obtaining high yields and selectivities in organic synthesis is the choice of the proper solvent. Solvent selection is often guided by the intuitive understanding of transition state-solvent interactions. However, quantum-mechanical tunneling can significantly contribute to chemical reactions, circumventing the transition state and thus depriving chemists of their intuitive handle on the reaction kinetics. In this Account, we aim to provide rationales for the effects of solvation on tunneling reactions derived from experiments performed in cryogenic matrices.The tunneling reactions analyzed here cover a broad range of prototypical organic transformations that are subject to strong solvation effects. Examples are the hydrogen tunneling probability for the cis-trans isomerization of formic acid which is strongly reduced upon formation of hydrogen-bonded complexes and the [1,2]H-shift in methylhydroxycarbene where a change in product selectivity is predicted upon interaction with hydrogen bond acceptors.Not only hydrogen but also heavy atom tunneling can exhibit strong solvent effects. The direction of the nearly degenerate valence tautomerization between benzene oxide and oxepin was found to reverse upon formation of a halogen or hydrogen bond with ICF₃ or H₂O. But even in the absence of strong noncovalent interactions such as hydrogen or halogen bonding, solvation can have a decisive effect on tunneling as evidenced by the Cope rearrangement of semibullvalenes via heavy-atom tunneling. Can quantum tunneling be catalyzed? The acceleration of the ring expansion of 1H-bicyclo[3.1.0.]-hexa-3,5-dien-2-one by complexation with Lewis acids provides a proof-of-concept for tunneling catalysis.Two concepts are central for the explanation and prediction of solvation effects on tunneling phenomena: a simple approach expands the Born-Oppenheimer approximation by separating nuclear degrees of freedom into intra- and intermolecular degrees. Intermolecular movements represent the slowest motions within molecular aggregates, thus effectively freezing the position of the solvent in relation to the reactant during the tunneling process. Another useful approach is to treat reactants and products by separate single-well potentials, where the intersection represents the transition state. Thus, stabilization of the reactants via solvation should result in an increase in barrier heights and widths which in turn lowers tunneling probabilities. These simple models can predict trends in tunneling kinetics and provide a rational basis for controlling tunneling reactions via solvation.

Polarization-Dependent Sum-Frequency Generation Spectroscopy for Ångstrom-Scale Depth Profiling of Molecules at Interfaces

The three-dimensional spatial distribution of molecules at soft matter interfaces is crucial for processes ranging from membrane biophysics to atmospheric chemistry. While several techniques can access surface composition, obtaining information on the depth distribution is challenging. We develop a noninvasive, polarization-resolved, surface-specific sum-frequency generation spectroscopy providing quantitative depth information. We demonstrate the technique on formic acid molecules at the air-water interface. With increasing molar fraction from 2.5% to 10%, the formic acid molecules shift, on average, ∼0.9  Å into the bulk. The consistency with the simulation data manifests that the technique allows for probing the Ångstrom-scale depth profile.

Self-association of diphenylpnictoginic acids in solution and solid state: covalent vs. hydrogen bonding

Triphenylpnictogens were oxidized to access diphenylpnictioginic acids Ph₂XOOH (X = P, As, Sb, Bi). It was shown that oxidation with chloramine-T does not lead to the cleavage of a C-pnictogen bond. The preliminary reductive cleavage with sodium in liquid ammonia followed by the oxidation with hydrogen peroxide was successfully utilised for the synthesis of diphenylphosphinic and diphenylarsinic acids. It was shown that in solid state (by means of XRD), all diphenylpnictoginic acids form polymeric chains. Diphenylbismuthinic and diphenylantimonic acids form polymeric covalent adducts, while diphenylphosphinic and diphenylarsinic chains are associated through hydrogen bonding. Unlike diphenylphosphinic acid, diphenilarsinic acid forms two polymorphs of hydrogen-bonded infinite chains. In solution in a polar aprotic solvent diphenylarsinic acid, similarly to dimethylarsinic, forms hydrogen-bonded cyclic dimers together with a small amount of cyclic trimers.

Accurate molecular orientation at interfaces determined by multimode polarization-dependent heterodyne-detected sum-frequency generation spectroscopy via multidimensional orientational distribution function

Many essential processes occur at soft interfaces, from chemical reactions on aqueous aerosols in the atmosphere to biochemical recognition and binding at the surface of cell membranes. The spatial arrangement of molecules specifically at these interfaces is crucial for many of such processes. The accurate determination of the interfacial molecular orientation has been challenging due to the low number of molecules at interfaces and the ambiguity of their orientational distribution. Here, we combine phase- and polarization-resolved sum-frequency generation (SFG) spectroscopy to obtain the molecular orientation at the interface. We extend an exponentially decaying orientational distribution to multiple dimensions, which, in conjunction with multiple SFG datasets obtained from the different vibrational modes, allows us to determine the molecular orientation. We apply this new approach to formic acid molecules at the air–water interface. The inferred orientation of formic acid agrees very well with ab initio molecular dynamics data. The phase-resolved SFG multimode analysis scheme using the multidimensional orientational distribution thus provides a universal approach for obtaining the interfacial molecular orientation.

Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra

The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH₂O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.

Interfacial Water Structure of Binary Liquid Mixtures Reflects Nonideal Behavior

The evaporation of molecules from water-organic solute binary mixtures is key for both atmospheric and industrial processes such as aerosol formation and distillation. Deviations from ideal evaporation energetics can be assigned to intermolecular interactions in solution, yet evaporation occurs from the interface, and the poorly understood interfacial, rather than the bulk, structure of binary mixtures affects evaporation kinetics. Here we determine the interfacial structure of nonideal binary mixtures of water with methanol, ethanol, and formic acid, by combining surface-specific vibrational spectroscopy with molecular dynamics simulations. We find that the free, dangling OH groups at the interfaces of these differently behaving nonideal mixtures are essentially indistinguishable. In contrast, the ordering of hydrogen-bonded interfacial water molecules differs substantially at these three interfaces. Specifically, the interfacial water molecules become more disordered (ordered) in mixtures with methanol and ethanol (formic acid), showing higher (lower) vapor pressure than that predicted by Raoult's law.

An optimized method for sample collection, extraction, and analysis of fentanyl and fentanyl analogs from a non-porous surface

Illicit use of the potent opioid fentanyl and its analogs (fentanyls) are on the rise in the United States. As use increases, drug production tends to also increase, leading to more locations being contaminated with the potentially lethal substance. Because fentanyl-contaminated locations may present a risk to the general public, a method for sampling, identifying, and quantitating these fentanyls from surfaces is in need. This research developed and optimized a surface-wipe collection and extraction method for 17 fentanyls and 10 common fentanyl adulterants from a non-porous surface and quantitated the amount of each compound collected with liquid chromatography tandem mass spectrometry. The final, optimized surface-wipe method resulted in an average collection and extraction efficiency (±SD) of 62.0 (±14.0)%, with a range of 34.1 (±2.6) - 82.5 (±9.6)%. While legislation has yet to be implemented regarding remediation levels for fentanyl-contaminated locations, when such legislation is drafted, this method can be implemented to determine the safety of these locations prior to and after decontamination has occurred.

High level ab initio investigation of the catalytic effect of water on formic acid decomposition and isomerization

Formic acid (FA) is a ubiquitous molecule found in the atmosphere, and is relevant to many important processes. The FA molecule generally exists as the trans isomer, which can decompose into H2O and CO (dehydration). It can also exist in the less favorable cis isomer which can decompose into H2 and CO2 (decarboxylation). Our work examines the complexes formed between each isomer of FA with water. We present geometries and vibrational frequencies obtained at the reliable CCSD(T)/aug-cc-pVTZ level of theory for seven FAwater complexes. We utilize the focal point method to determine CCSDT(Q)/CBS plus corrections binding energies of 7.37, 3.36, and 2.02 kcal mol-1 plus 6.07, 3.79, 2.60, and 2.55 kcal mol-1 for the trans-FAwater and cis-FAwater complexes, respectively. Natural bond orbital analysis is used to further decompose the interactions in each complex and gain insight into their relative strengths. Furthermore, we examine the effect that a single water molecule has on the barrier heights to each decomposition pathway by optimizing the transition states and verifying their connectivity with intrinsic reaction coordinate computations as well as utilizing a kinetic model. Water lowers the barrier to dehydration by at most 15.78 kcal mol-1 and the barrier to decarboxylation by up to 15.90 kcal mol-1. Our research also examines for the first time the effect of one water molecule on the interconversion barrier and we find that the barrier from trans to cis is not catalyzed by water due to the strong FA and water interactions. Our results highlight some instances where different binary complexes result in different decomposition pathways and even a case where one binary complex can form the same decomposition products via two distinct mechanisms. Our results provide a reliable benchmark of the FAH2O system as well as provide insight into future studies of similar atmospheric systems.

Carboxylic acids in aqueous solutions: Hydrogen bonds, hydrophobic effects, concentration fluctuations, ionization, and catalysis

In the frequency range between 100 kHz and 2 GHz, ultrasonic absorption spectra have been measured for a series of carboxylic acids from formic to enanthic acid, including constitutional isomers. Also investigated have been the spectra for mixtures with water of short-chain formic, acetic, propionic, butyric, and isobutyric acid, in each case covering the complete composition range. The neat carboxylic acids feature two Debye-type relaxation terms with relaxation times between 5.6 and 260 ns as well as 0.14 and 1.4 ns, respectively, at room temperature. Depending on the composition, mixtures with water reveal an additional Debye relaxation term in the intermediate frequency range (acetic acid) or a term subject to a relaxation time distribution (propionic, butyric, and isobutyric acid). The relaxations of the neat acids are assigned to the equilibrium between monomers and single-hydrogen-bonded linear dimers and between linear and twofold-hydrogen-bonded cyclic dimers. The latter equilibrium is considerably catalyzed by hydronium and carboxylate ions. Several mixtures with water indicate one of the up to three Debye relaxations to reflect the protolysis of the organic acid. The term with underlying relaxation time distribution is due to noncritical fluctuations in the local concentrations. The Debye relaxations are evaluated to yield the parameters of the relevant elementary chemical reactions, such as the rate and equilibrium constants and the isentropic reaction volumes. A comparison of the correlation length of concentration fluctuations with data for other aqueous systems confirms the idea that the hydrophobic part of the organic constituent promotes the formation of a micro-heterogeneous liquid structure, whereas the hydrophilic moiety is of minor importance in this respect. The high-frequency limiting absorption suggests the equilibrium between conformers of linear dimers to contribute to the spectra well above the frequency range of measurements.

VUV photochemistry of the H2O⋯CO complex in noble-gas matrices: formation of the OH⋯CO complex and the HOCO radical

VUV photolysis of the H₂O⋯CO complexes leads to the formation of the OH⋯CO radical–molecule complexes and trans-HOCO radicals.

Hydrogen bonding in the hydroxysulfinyl radical-formic acid-water system: A theoretical study

Quantum chemical methods have been employed to evaluate the possible configurations of the 1:1 and 1:2 HOSO-formic acid complexes and 1:1:1 HOSO-formic acid-water complexes. The first type of complex involves two H bonds, while the other two types comprise three H bonds in a ring. The complexes are relatively stable, with CBS-QB3 computed binding energies of 14.3 kcal mol(-1) , 23.4 kcal mol(-1) , and 21.1 kcal mol(-1) for the lowest-energy structures of the 1:1, 1:2, and 1:1:1 complexes, respectively. Complex formations induce a large spectral red-shift and an enhancement of the IR intensity for the H-bonded OH stretching modes relative to those in the parent monomers. TDDFT calculations of the low-lying electronic excited states demonstrate that the complexes are photochemically quite stable in the troposphere. Small spectral shifts in comparison to the free HOSO radical suggest that the radical and the complexes would not be easily distinguishable using standard UV/vis absorption spectroscopy. © 2016 Wiley Periodicals, Inc.

Acetic acid-water complex: The first observation of structures containing the higher-energy acetic acid conformer

Non-covalent interaction of acetic acid (AA) and water is studied experimentally by IR spectroscopy in a nitrogen matrix and theoretically at the MP2 and coupled-cluster with single and double and perturbative triple excitations [CCSD(T)]/6-311++G(2d,2p) levels of theory. This work is focused on the first preparation and characterization of complexes of higher-energy (cis) conformer of AA with water. The calculations show three 1:1 structures for the trans-AA⋯H2O complexes and three 1:1 structures for the cis-AA⋯H2O complexes. Two trans-AA⋯H2O and two cis-AA⋯H2O complexes are found and structurally assigned in the experiments. The two cis-AA⋯ ⋅ H2O complexes are obtained by annealing of a matrix containing water and cis-AA molecules prepared by selective vibrational excitation of the ground-state trans form. The less stable trans-AA⋯H2O complex is obtained by vibrational excitation of the less stable cis-AA⋯H2O complex. In addition, the 1:2 complexes of trans-AA and cis-AA with water molecules are studied computationally and the most stable forms of the 1:2 complexes are experimentally identified.

Deprotonation of formic acid in collisions with a liquid water surface studied by molecular dynamics and metadynamics simulations

Formic acid has a lower barrier to deprotonation at the air–water interface than in bulk liquid water.

Calorimetric and spectroscopic studies on solvation energetics for H2storage in the CO2/HCOOH system

The role of solvent interactions in H₂-storage/delivery in the carbon dioxide–formic acid couple.

Matrix-isolation studies of noncovalent interactions: more sophisticated approaches

Noncovalent interactions are crucial for many physical, chemical, and biological phenomena. Matrix isolation is a powerful method to study noncovalent interactions, including hydrogen-bonded species, and it has been extensively used in this field. However, there are difficult situations, such as in the case of species that are impossible to prepare in the gas phase. In this article, we describe some advanced approaches allowing studies of complexes that are problematic for the traditional methods. Photolysis of a suitable precursor in a matrix can lead to a large concentration of 1:1 complexes, which are otherwise very difficult to prepare (e.g., the H2O···O complex). Photolysis of species combined with annealing can lead to complexes of molecules with mobile atoms (e.g., the same H2O···O complex). Simultaneous photolysis of two species combined with annealing can produce complexes of radicals via reactions of the photogenerated complexes with mobile atoms (e.g., the H2O···HCO complex). Interaction of noble-gas (Ng) hydrides with other species is another topic (e.g., the N2···HArF complex) and very large blue shifts of the H-Ng stretching modes are normally observed for these systems. Complexes and dimers of the higher-energy conformer of formic acid have been prepared by using selective vibrational excitation of the ground-state conformer. The higher-energy conformer of formic acid can be efficiently stabilized in the complexes with strong hydrogen bonding. We also consider some problematic cases when the changes in the vibrational frequencies of the 1:1 complexes are very small (e.g., the phenol···Xe complex) and when the complex formation is prevented by strong solvation in the matrix (e.g., species in solid xenon).

Radiation-induced transformations of matrix-isolated formic acid: evidence for the HCOOH → HOCO + H channel

HOCO radicals are produced under X-ray irradiation of isolated formic acid at low temperatures, and possible mechanisms and astrochemical implications are discussed.

A theoretical study on the red- and blue-shift hydrogen bonds of cis-trans formic acid dimer in excited states

The excited states of cis-trans formic acid dimer and its monomers have been investigated by time-dependent density functional theory (TDDFT) method. The formation of intermolecular hydrogen bonds O1-H1...O2=C2 and C2-H2...O4=C1 induces bond length lengthening of the groups related to the hydrogen bond, while that of the C2-H2 group is shortened. It is demonstrated that the red-shift hydrogen bond O1-H1...O2=C2 and blue-shift hydrogen bond C2-H2...O4=C1 are both weakened when excited to the S1 state. Moreover, it is found that the groups related to the formation of red-shift hydrogen bond O1-H1...O2=C2 are both strengthened in the S1 state, while the groups related to the blue-shift hydrogen bond C2-H2...O4=C1 are both weakened. This will provide information for the photochemistry and photophysical study of red- and blue-shift hydrogen bond.

Microsolvation of the formic acid dimer — (HCOOH)2(H2O)n clusters with n = 1, . . ., 5

Density functional theory and spin-component-scaled Møller–Plesset perturbation theory calculations are used to examine the microsolvation of the formic acid dimer. The lowest energy structures with n water molecules consist of a n-water cluster, not necessarily of lowest energy, with two formic acid molecules attached to its surface by hydrogen bonds. The total number of hydrogen bonds does not correlate directly with relative stability.