Computational Study on the Kinetics and Mechanisms for the Unimolecular Decomposition of Formic and Oxalic Acids

Liberation of carbon monoxide from formic acid mediated by molybdenum oxyanions

Multistage mass spectrometry experiments, isotope labelling and DFT calculations were used to explore whether selective decarbonylation of formic acid could be mediated by molybdate anions [(MoO3)x(OH)]- (x = 1 and 2) via a formal catalytic cycle involving two steps. In step 1, both molybdate anions undergo gas-phase ion-molecule reactions (IMR) with formic acid to produce the coordinated formates [(MoO3)x(O2CH)]- and H2O. In step 2, both coordinated formates [(MoO3)x(O2CH)]- undergo decarbonylation under collision-induced dissociation (CID) conditions to reform the molybdate anions [(MoO3)x(OH)]- (x = 1 and 2), thus closing a formal catalytic cycle. In the case of [MoO3(O2CH)]- an additional decarboxylation channel also occurs to yield [MoO3(H)]-, which is unreactive towards formic acid. The reaction between [Mo18O3(18OH)]- and formic acid gives rise to [Mo18O3(O2CH)]- highlighting that ligand substitution occurs without 18O/16O exchange between the coordinated 18OH ligand and HC16O2H. The reaction between [(MoO3)x(OD)]- (x = 1 and 2) and DCO2H initially produces [(MoO3)x(OH)]- (x = 1 and 2), indicating that D/H exchange occurs. DFT calculations were carried out to investigate the reaction mechanisms and energetics associated with both steps of the formal catalytic cycle and to better understand the competition between decarbonylation and decarboxylation, which is crucial in developing a selective catalyst. The CO and CO2 loss channels from the monomolybdate anion [MoO3(O2CH)]- have similar barrier heights which is in agreement with experimental results where both fragmentation channels are observed. In contrast, the dimolybdate anion is more selective, since the decarbonylation pathway of [(MoO3)2(O2CH)]- is both kinetically and thermodynamically favoured, which agrees with experimental observations where the CO loss channel is solely observed.

Direct conversion of lignin to functionalized diaryl ethers via oxidative cross-coupling

Efficient valorization of lignin, a sustainable source of functionalized aromatic products, would reduce dependence on fossil-derived feedstocks. Oxidative depolymerization is frequently applied to lignin to generate phenolic monomers. However, due to the instability of phenolic intermediates, repolymerization and dearylation reactions lead to low selectivity and product yields. Here, a highly efficient strategy to extract the aromatic monomers from lignin affording functionalized diaryl ethers using oxidative cross-coupling reactions is described, which overcomes the limitations of oxidative methods and affords high-value specialty chemicals. Reaction of phenylboronic acids with lignin converts the reactive phenolic intermediates into stable diaryl ether products in near-theoretical maximum yields (92% for beech lignin and 95% for poplar lignin based on the content of β-O-4 linkages). This strategy suppresses side reactions typically encountered in oxidative depolymerization reactions of lignin and provides a new approach for the direct transformation of lignin into valuable functionalized diaryl ethers, including key intermediates in pharmaceutical and natural product synthesis.

Near thermal, selective liberation of hydrogen from formic acid catalysed by copper hydride ate complexes

A near thermal two-step catalytic cycle for the selective release of hydrogen from formic acid by mononuclear cuprate anions was revealed using multistage mass spectrometry experiments, deuterium labelling and DFT calculations. In gas-phase ion-molecule reactions, mononuclear copper hydride anions [(L)Cu(H)]^- (where L = H^-, O₂CH^-, BH₄^- and CN^-) were found to react with formic acid (HCO₂H) to yield [(L)Cu(O₂CH)]^- and H₂. The copper formate anions [(L)Cu(O₂CH)]^- can decarboxylate via collision-induced dissociation (CID) to reform the copper hydride [(L)Cu(H)]^-, thereby closing the two-step catalytic cycle. Analogous labelling experiments with d₁-formic acid (DCO₂H) reveal that the decarboxylation process also occurs spontaneously. A kinetic study was carried out to provide further insights into the species involved in this reaction. Energetics from density functional theory (DFT) calculations show that the key decarboxylation step can occur without CID, thus in support of experimental observations.

Probing the pyrolysis of ethyl formate in the dilute gas phase by synchrotron radiation and theory

The thermal decomposition of the atmospheric constituent ethyl formate was studied by coupling flash pyrolysis with imaging photoelectron photoion coincidence (iPEPICO) spectroscopy using synchrotron vacuum ultraviolet (VUV) radiation at the Swiss Light Source (SLS). iPEPICO allows photoion mass-selected threshold photoelectron spectra (ms-TPES) to be obtained for pyrolysis products. By threshold photoionization and ion imaging, parent ions of neutral pyrolysis products and dissociative photoionization products could be distinguished, and multiple spectral carriers could be identified in several ms-TPES. The TPES and mass-selected TPES for ethyl formate are reported for the first time and appear to correspond to ionization of the lowest energy conformer having a cis (eclipsed) configuration of the O=C(H)-O-C(H₂ )-CH₃ and trans (staggered) configuration of the O=C(H)-O-C(H₂ )-CH₃ dihedral angles. We observed the following ethyl formate pyrolysis products: CH₃ CH₂ OH, CH₃ CHO, C₂ H₆ , C₂ H₄ , HC(O)OH, CH₂ O, CO₂ , and CO, with HC(O)OH and C₂ H₄ pyrolyzing further, forming CO + H₂ O and C₂ H₂  + H₂ . The reaction paths and energetics leading to these products, together with the products of two homolytic bond cleavage reactions, CH₃ CH₂ O + CHO and CH₃ CH₂  + HC(O)O, were studied computationally at the M06-2X-GD3/aug-cc-pVTZ and SVECV-f12 levels of theory, complemented by further theoretical methods for comparison. The calculated reaction pathways were used to derive Arrhenius parameters for each reaction. The reaction rate constants and branching ratios are discussed in terms of the residence time and newly suggest carbon monoxide as a competitive primary fragmentation product at high temperatures.

Enabling simultaneous redox transformation of toxic chromium(VI) and arsenic(III) in aqueous media-A review

Simultaneous conversion of most harmful As(III) and Cr(VI) to their less toxic counterparts is environmentally desirable and cost-effective. It has been confirmed that simultaneous oxidation of As(III) to As(V) and reduction of Cr(VI) to Cr(III) can occur via free radical or mediated electron transfer processes. While Cr(VI) is reduced by reacting with H^•, e_aq^-, photoelectron directly or undergoing ligand exchange with H₂O₂ and SO₃^2-, As(III) is oxidized by HO^•, SO₄^•-, O₂^•-, and holes (h⁺) in free radical process. The ability to concentrate Cr and As species on heterogeneous interface and conductivity determining the co-conversion efficiency in mediated electron transfer process. Acidity has positive effect on these co-conversion, while mediated electron transfer process is not much affected by dissolved oxygen (O₂). Organic compounds (e.g., oxalate, citrate and phenol) commonly favor Cr(VI) reduction and inhibit As(III) oxidation. To better understand the trends in the existing data and to identify the knowledge gaps, this review elaborates the complicated mechanisms for co-conversion of As(III) and Cr(VI) by various methods. Some challenges and prospects in this active field are also briefly discussed.

The influence of the environment in chemical reactivity: the HCOOH formation from the H2O + CO reaction

The reaction between carbon monoxide and water was studied occurring in an aerosol medium rich in methanol. This environment is plausible for the primitive and prebiotic Earth atmosphere. The chemical environment is expressed in terms of dielectric constant (ε) and the chemical system was modeled employing the polarizable continuum model (PCM). The main results were acquired from calculations employing the M06-2X density functional for the electronic structure calculations and the canonical variational theory with small curvature tunneling for the chemical kinetic calculations. The rise of ε affects both the thermochemistry and the kinetics of the reaction, increasing the barrier height and decreasing the rate constant for the reaction occurring at room temperature. For example, the rate constant at 300 K is 5-10× 10^- 53 cm³ ⋅molecule^- 1 ⋅s^- 1 for low dielectric constant (ε < 3) and around 2-4× 10^- 53 cm³ ⋅molecule^- 1 ⋅s^- 1 for ε between 7 and 40. Our results indicate that the ε variation allows a fine tuning to the rate of the reaction.

Direct Dynamics Simulations of the 3CH2 + 3O2 Reaction at High Temperature

Direct dynamics simulations with the M06/6-311++G(d,p) level of theory were performed to study the ³CH₂ + ³O₂ reaction at 1000 K temperature on the ground state singlet surface. The reaction is complex with formation of many different product channels in highly exothermic reactions. CO, CO₂, H₂O, OH, H₂, O, H, and HCO are the products formed from the reaction. The total simulation rate constant for the reaction at 1000 K is (1.2 ± 0.3) × 10^-12 cm³ molecule^-1 s^-1, while the simulation rate constant at 300 K is (0.96 ± 0.28) × 10^-12 cm³ molecule^-1 s^-1. The simulated product yields show that CO is the dominant product and the CO:CO₂ ratio is 5.3:1, in good comparison with the experimental ratio of 4.3:1 at 1000 K. On comparing the product yields for the 300 and 1000 K simulations, we observed that, except for CO and H₂O, the yields of the other products at 1000 K are lower at 300 K, showing a negative temperature dependence.

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.

Kinetic and mechanism study of UV/pre-magnetized-Fe0/oxalate for removing sulfamethazine

In this study, UV irradiated photochemical reactions of oxalate (Ox) with premagnetized-Fe⁰ (pre- Fe⁰) as the catalyst was used to degrade sulfamethazine (SMT). Magnetic field promoted the release of iron ion from Fe⁰ thus enhanced SMT and Ox removal in UV/pre- Fe⁰/Ox process. X-ray photoelectron spectroscopy demonstrated that the presence of UV and Ox promoted the transformation of Fe³⁺ to Fe²⁺ on Fe⁰, which enhanced the surface bound •OH (•OH_surf) generation. Ox inhibited the formation of iron (hydro)xides and enhanced the hydroxylation of Fe⁰ surface. •OH_surf was mainly responsible for SMT removal (44%), while UV direct photolysis and •OH in the solution both caused around 28% SMT removal. The process with Ox exhibited much higher efficiency in SMT degradation than that added with H₃PO₄, citric acid and ethylenediaminetetraacetic acid, which greatly expanded the chelate-modified Fenton processes and their treatment efficiency.

Role of Ligand in the Selective Production of Hydrogen from Formic Acid Catalysed by the Mononuclear Cationic Zinc Complexes [(L)Zn(H)]+ (L=tpy, phen, and bpy)

A series of zinc-based catalysts was evaluated for their efficiency in decomposing formic acid into molecular hydrogen and carbon dioxide in the gas phase using quadrupole ion trap mass spectrometry experiments. The effectiveness of the catalysts in the series [(L)Zn(H)]⁺ , where L=2,2':6',2''-terpyridine (tpy), 1,10-phenanthroline (phen) or 2,2'-bipyrydine (bpy), was found to depend on the ligand used, which turned out to be fundamental in tuning the catalytic properties of the zinc complex. Specifically, [(tpy)Zn(H)]⁺ displayed the fastest reaction with formic acid proceeding by dehydrogenation to produce the zinc formate complex [(tpy)Zn(O₂ CH)]⁺ and H₂ . The catalysts [(L)Zn(H)]⁺ are reformed by decarboxylating the zinc formate complexes [(L)Zn(O₂ CH)]⁺ by collision-induced dissociation, which is the only reaction channel for each of the ligands used. The decarboxylation reaction was found to be reversible, since the zinc hydride complexes [(L)Zn(H)]⁺ react with carbon dioxide yielding the zinc formate complex. This reaction was again substantially faster for L=tpy than L=phen or bpy. The energetics and mechanisms of these processes were modelled using several levels of density functional theory (DFT) calculations. Experimental results are fully supported by the computational predictions.

Electrophilicity of oxalic acid monomer is enhanced in the dimer by intermolecular proton transfer

We have analyzed the effect of excess electron attachment on the network of hydrogen bonds in the oxalic acid dimer (OA)₂. The most stable anionic structures may be viewed as complexes of a neutral hydrogenated moiety HOA˙ coordinated to an anionic deprotonated moiety (OA-H)^-. HOA˙ acts as a double proton donor and (OA-H)^- as a double proton acceptor. Thus the excess electron attachment drives intermolecular proton transfer. We have identified several cyclic hydrogen bonded structures of (OA)₂^-. Their stability has been analyzed in terms of the stability of the involved conformers, the energetic penalty for deformation of these conformers to the geometry of the dimer, and the two-body interaction energy between the deformed HOA˙ and (OA-H)^-. There are at least seven isomers of (OA)₂^- with stabilization energies in the range of 1.26-1.39 eV. These energies are dominated by attractive two-body interaction energies. The anions are vertically bound electronically by 3.0-3.4 eV and adiabatically bound by at least 1.6 eV. The computational predictions are consistent with the anion photoelectron spectrum of (OA)₂^-. The spectrum consists of a broad feature, with an onset of 2.5 eV and spanning to 4.3 eV. The electron vertical detachment energy (VDE) is assigned to be 3.3 eV.

Combined Ab Initio, Kinetic Modeling, and Shock Tube Study of the Thermal Decomposition of Ethyl Formate

The potential energy surfaces (PESs) and reaction rate constants of the unimolecular decomposition of ethyl formate (EF) were investigated using high-precision theoretical methods at the CCSD(T)/CBS(T-Q)//M06-2X/6-311++G(d,p) level of theory. The calculated PESs of EF dissociation and molecular decomposition reactions indicate that the intramolecular H-shift to produce formic acid and ethylene is the dominant decomposition pathway. A detailed chemical kinetic mechanism for EF pyrolysis was constructed by incorporating the important reactions of EF and its radicals into an existing mechanism previously developed for small methyl esters. The updated mechanism was first used to reproduce CO, CO₂, and H₂O concentration time histories during EF pyrolysis in the shock tube reported by Ren et al. [ Ren , W. , Mitchell Spearrin , R. , Davidson , D. F. , and Hanson , R. K. J. Phys. Chem. A 2014 , 118 , 1785 - 1798 ]. The rate of production and sensitivity analyses show that the competing dehydration and decarboxylation channels of the intermediate formic acid control the final product yields of EF pyrolysis. The EF mechanism was further validated against the shock tube data of OH, CO, CO₂, and H₂O time histories measured during EF oxidation (equivalence ratio Φ = 1.0) at 1331-1615 K and 1.52-1.74 atm. This revised EF mechanism captured all of the species' time histories over the entire temperature range. Such modeling capability was due to the more accurate rate constants of EF reactions determined by high-precision theoretical calculations and a high-fidelity C₀-C₂ basis mechanism.

Evaluation of the photooxidation efficiency of As(III) applying the UVC/oxalate technique

In this study, the photooxidation capacity of UVC/Oxalate (Ox) was evaluated using As(III) as a typical pollutant. The results show that the direct oxidation amount of As(III) induced by UVC in water was negligible, but the presence of Ox remarkably accelerated the oxidation rate of As(III). Under UVC irradiation, 50 μM As(III) can be completely oxidized to As(V) in the case of Ox concentration above 300 μM within 60 min. As(III) oxidation was found greatly related with the photodecomposition of Ox. Much more Ox can be mineralized in more acidic solution. At the same time, the photooxidation of As(III) was significantly favored at decreased initial pH from 8.0 to 3.0. In this reaction system, the role of oxygen was indispensable for Ox photodecomposition and As(III) photooxidation, which can be ascribed to its special roles as a precursor of reactive superoxide and an electron acceptor. In oxygen-present atmosphere, the in situ production of H₂O₂ was detected during the photolysis of Ox and its photolysis product, i.e., OH primarily contributed to the oxidation of As(III). However, the photodecomposition of Ox and photooxidation of As(III) were significantly inhibited in the anaerobic environment. In general, the homogeneous photolysis of Ox in many commonly practiced UVC oxidation processes can be also proposed as a supplementary method of generating highly oxiditive species in aerobic condition.

Ligand-induced substrate steering and reshaping of [Ag2(H)](+) scaffold for selective CO2 extrusion from formic acid

Metalloenzymes preorganize the reaction environment to steer substrate(s) along the required reaction coordinate. Here, we show that phosphine ligands selectively facilitate protonation of binuclear silver hydride cations, [LAg₂(H)]⁺ by optimizing the geometry of the active site. This is a key step in the selective, catalysed extrusion of carbon dioxide from formic acid, HO₂CH, with important applications (for example, hydrogen storage). Gas-phase ion-molecule reactions, collision-induced dissociation (CID), infrared and ultraviolet action spectroscopy and computational chemistry link structure to reactivity and mechanism. [Ag₂(H)]⁺ and [Ph₃PAg₂(H)]⁺ react with formic acid yielding Lewis adducts, while [(Ph₃P)₂Ag₂(H)]⁺ is unreactive. Using bis(diphenylphosphino)methane (dppm) reshapes the geometry of the binuclear Ag₂(H)⁺ scaffold, triggering reactivity towards formic acid, to produce [dppmAg₂(O₂CH)]⁺ and H₂. Decarboxylation of [dppmAg₂(O₂CH)]⁺ via CID regenerates [dppmAg₂(H)]⁺. These gas-phase insights inspired variable temperature NMR studies that show CO₂ and H₂ production at 70 °C from solutions containing dppm, AgBF₄, NaO₂CH and HO₂CH.

Designing catalysts and understanding the influence of ligands for particular transformations remains a highly challenging task. Here, the authors show that bisphosphine ligands can alter the geometry of the active site in silver catalysts, driving protonation and ultimately extrusion of carbon dioxide from formic acid.

An unusual co-crystal [(μ2-dcpm)Ag2(μ2-O2CH)(η2-NO3)]2·[(μ2-dcpm)2Ag4(μ2-NO3)4] and its connection to the selective decarboxylation of formic acid in the gas phase

Binuclear silver cluster with μ₂-formate and μ₂-dcpm catalyses the decomposition of HCO₂H. This structural motif is present in the crystal.

Palladium-atom catalyzed formic acid decomposition and the switch of reaction mechanism with temperature

We carefully calculated the mechanism of one-atom model and its poisoned species, PdCO, as formic acid decomposition catalysts.

Reaction dynamics of O(¹D) + HCOOD/DCOOH investigated with time-resolved Fourier-transform infrared emission spectroscopy

We investigated the reaction dynamics of O((1)D) towards hydrogen atoms of two types in HCOOH. The reaction was initiated on irradiation of a flowing mixture of O3 and HCOOD or DCOOH at 248 nm. The relative vibration-rotational populations of OH and OD (1 ≦ v ≦ 4, J ≤ 15) states were determined from time-resolved IR emission recorded with a step-scan Fourier-transform spectrometer. In the reaction of O((1)D) + HCOOD, the rotational distribution of product OH is nearly Boltzmann, whereas that of OD is bimodal. The product ratio [OH]/[OD] is 0.16 ± 0.05. In the reaction of O((1)D) + DCOOH, the rotational distribution of product OH is bimodal, but the observed OD lines are too weak to provide reliable intensities. The three observed OH/OD channels agree with three major channels of production predicted with quantum-chemical calculations. In the case of O((1)D) + HCOOD, two intermediates HOC(O)OD and HC(O)OOD are produced in the initial C-H and O-D insertion, respectively. The former undergoes further decomposition of the newly formed OH or the original OD, whereas the latter produces OD via direct decomposition. Decomposition of HOC(O)OD produced OH and OD with similar vibrational excitation, indicating efficient intramolecular vibrational relaxation, IVR. Decomposition of HC(O)OOD produced OD with greater rotational excitation. The predicted [OH]/[OD] ratio is 0.20 for O((1)D) + HCOOD and 4.08 for O((1)D) + DCOOH; the former agrees satisfactorily with experiments. We also observed the v3 emission from the product CO2. This emission band is deconvoluted into two components corresponding to internal energies E = 317 and 96 kJ mol(-1) of CO2, predicted to be produced via direct dehydration of HOC(O)OH and secondary decomposition of HC(O)O that was produced via decomposition of HC(O)OOH, respectively.

Facile C(sp(2))-C(sp(2)) bond cleavage in oxalic acid-derived radicals

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. Many questions remain about the catalytic mechanism of OxDC although it has been proposed that the reaction proceeds via substrate-based radical intermediates. Using coupled cluster theory combined with implicit solvation models we have examined the effects of radical formation on the structure and reactivity of oxalic acid-derived radicals in aqueous solution. Our results show that the calculated solution-phase free-energy barrier for C-C bond cleavage to form CO2 is decreased from 34.2 kcal/mol for oxalic acid to only 9.3 kcal/mol and a maximum of 3.5 kcal/mol for the cationic and neutral oxalic acid-derived radicals, respectively. These studies also show that the C-C σ bonding orbital of the radical cation contains only a single electron, giving rise to an elongated C-C bond distance of 1.7 Å; a similar lengthening of the C-C bond is not observed for the neutral radical. This study provides new chemical insights into the structure and stability of plausible intermediates in the catalytic mechanism of OxDC, and suggests that removal of an electron to form a radical (with or without the concomitant loss of a proton) may be a general strategy for cleaving the unreactive C-C bonds between adjacent sp(2)-hybridized carbon atoms.

Experimental and modeling study of the thermal decomposition of C3-C5 ethyl esters behind reflected shock waves

The thermal decomposition of three ethyl esters, ethyl formate (C3H6O2), ethyl acetate (C4H8O2), and ethyl propanoate (C5H10O2), was studied behind reflected shock waves using laser absorption to measure concentration time-histories of H2O, CO2, and CO. Experimental conditions covered temperatures of 1301-1636 K, pressures of 1.48-1.72 atm, and reactant concentrations of 2000 ppm in argon. Recently developed mid-infrared laser diagnostics for H2O (2.5 μm), CO2 (4.3 μm), and CO (4.6 μm) provide orders-of-magnitude greater detectivity compared to previous near-infrared absorption sensors. The experimental results have highlighted significant differences among these three ethyl esters: negligible CO2 production during ethyl formate pyrolysis, quite slow CO formation rate during ethyl acetate pyrolysis, and nearly equal formation rate of H2O, CO2, and CO during ethyl propanoate pyrolysis. Detailed kinetic modeling was performed to understand the destruction pathways of these three ethyl esters with different alkyl chain lengths. Rate of production and sensitivity analyses were also carried out to interpret the experimental results and to identify the key reactions affecting experimental results.