The Role of Hydrogen Bonding in the Decomposition of H2CO3 in Water: Mechanistic Insights from Ab Initio Metadynamics Studies of Aqueous Clusters

13 C NMR as an analytical tool for the detection of carbonic acid and pKa determination

NMR spectroscopy has become a standard technique in studies both on carbon capture and storage. 13 C NMR allows the detection of two peaks for carbonated aqueous samples: one for CO2(aq) and another one for the species H2 CO3 , HCO3 - , and CO3 2- -herein collectively named Hx CO3 x-2 . The chemical shift of this second peak depends on the molar fraction of the three species in equilibrium and has been used to assess the equilibrium between HCO3 - and CO3 2- . The detection of H2 CO3 at low pH solutions is hindered, because of the concurrent liberation of CO2 when the medium is acidified. Herein, a valved NMR tube facilitates the detection of the Hx CO3 x-2 peak across a wide pH range, even at pH 1.8 where the dominant species is H2 CO3 . The method employed the formation of frozen layers of NaH13 CO3 and acid solutions within the tube, which are mixed as the tube reaches room temperature. At this point, the tube is already securely sealed, preventing any loss of CO2 to the atmosphere. A spectrophotometry approach allowed the measurement of the actual pH inside the pressurized NMR tube. The chemical shift for H2 CO3 was determined as 160.33 ± 0.03 ppm, which is in good agreement with value obtained by DFT calculations combined with Car-Parrinello molecular dynamics. The H2 CO3 pKa value determined by the present method was 3.41 ± 0.03, for 15% D2 O aqueous medium and 0.8 mol/L ionic strength. The proposed method can be extended to studies about analogs such as alkyl carbonic and carbamic acids.

Hydrogen bonds-triggered differential extraction efficiencies for bifenthrin by three polymeric ionic liquids with varying anions based on FT-IR spectroscopy

Herein, we fabricated three imidazolium-based polymeric ionic liquids (PILs) with different anions (P[VEIM]BF4, P[VEIM]PF6 and P[VEIM]Br), and analyzed their differential extraction efficiencies for bifenthrin through H-bonding induced effects. Three PILs all presented an irregular block structure with rough surface and lower specific-surface area (SSA, 11.2-18.7 m2 g-1) than carbon-based nanomaterials. They formed hydrogen bonds with free-water molecules in the lattice of PILs, including C2,4,5-H⋯O-H, Br⋯H-O-H⋯Br, O-H⋯Br, C2,4,5-H⋯F-P, P-F⋯H-O-H⋯F-P, C2,4,5-H⋯F-B and B-F⋯H-O-H⋯F-B. After extraction, the O-H stretching-vibration peak was prominently intensified, whereas the C-H bond varied slightly concomitant with reduced B-F and P-F vibration. Theoretically, the C-H vibration should become more intense in the C4,5-H⋯H2O and C2-H⋯H2O bonds after extraction in contrast to before extraction. These contrary spectral changes demonstrated that the hydrogen bonds between cations in the PILs and free-water molecules were broken after extraction, yielding the H-bonding occurrence between bifenthrin and H-O-H in the lattice. As a time indicator for the free-water binding and releasing process, the highest slope for the plot of I t /I 0 against time implied that the shortest time was required for P[VEIM]PF6 to reach an adsorption equilibrium. Overall, the strong hydrophobicity, small SSA and electrostatic-repulsion force for P[VEIM]PF6 are all not conducive to its efficient adsorption. Beyond our anticipation, P[VEIM]PF6 provided the highest extraction recovery for bifenthrin up to 92.4% among three PILs. Therefore, these data lead us to posit that the above high efficiency results from the strongest H-bonding effect between P[VEIM]PF6 and bifenthrin. These findings promote our deep understanding of PILs-triggered differential efficiency through a H-bonding induced effect.

Large Presence of Carbonic Acid in CO2-Rich Aqueous Fluids under Earth's Mantle Conditions

The chemistry of carbon in aqueous fluids at extreme pressure and temperature conditions is of great importance to Earth's deep carbon cycle, which substantially affects the carbon budget at Earth's surface and global climate change. At ambient conditions, the concentration of carbonic acid in water is negligible; therefore, aqueous carbonic acid was simply ignored in previous geochemical models. However, by applying extensive ab initio molecular dynamics simulations at pressure and temperature conditions similar to those in Earth's upper mantle, we found that carbonic acid can be the most abundant carbon species in aqueous CO₂ solutions at ∼10 GPa and 1000 K. The mole percent of carbonic acid in total dissolved carbon species increases with increasing pressure along an isotherm, while its mole percent decreases with increasing temperature along an isobar. In CO₂-rich solutions, we found significant proton transfer between carbonic acid molecules and bicarbonate ions, which may enhance the conductivity of the solutions. The effects of pH buffering by carbonic acid may play an important role in water-rock interactions in Earth's interior. Our findings suggest that carbonic acid is an important carbon carrier in the deep carbon cycle.

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide Salt on the Deprotonation of Formic Acid in Aqueous Solution

The deprotonation of formic acid is investigated using metadynamics in tandem with Born-Oppenheimer molecular dynamics simulations. We compare our findings for formic acid in pure water with previous studies before examining formic acid in aqueous solutions of lithium bromide. We carefully consider different definitions for the collective variable(s) used to drive the metadynamics, emphasizing that the variables used must include all of the possible reactive atoms in the system, in this case carboxylate oxygens and water hydrogens. This ensures that all the various possible proton exchange events can be accommodated and the collective variable(s) can distinguish the protonated and deprotonated states, even over rather long ab initio simulation runs (ca. 200-300 ps). Our findings show that the formic acid deprotonation barrier and the free energy of the deprotonated state are higher in concentrated lithium bromide, in agreement with the available experimental data for acids in salt solution. We show that the presence of Br^- in proximity to the formic acid hydroxyl group effectively inhibits deprotonation. Our study extends previous work on acid deprotonation in pure water and at air-water interfaces to more complex multicomponent systems of importance in atmospheric and marine chemistry.

Characterization of a trans-trans Carbonic Acid-Fluoride Complex by Infrared Action Spectroscopy in Helium Nanodroplets

The high Lewis basicity and small ionic radius of fluoride promote the formation of strong ionic hydrogen bonds in the complexation of fluoride with protic molecules. Herein, we report that carbonic acid, a thermodynamically disfavored species that is challenging to investigate experimentally, forms a complex with fluoride in the gas phase. Intriguingly, this complex is highly stable and is observed in abundance upon nanoelectrospray ionization of an aqueous sodium fluoride solution in the presence of gas-phase carbon dioxide. We characterize the structure and properties of the carbonic acid–fluoride complex, F^–(H₂CO₃), and its deuterated isotopologue, F^–(D₂CO₃), by helium nanodroplet infrared action spectroscopy in the photon energy range of 390–2800 cm^–1. The complex adopts a C_2v symmetry structure with the carbonic acid in a planar trans–trans conformation and both OH groups forming ionic hydrogen bonds with the fluoride. Substantial vibrational anharmonic effects are observed in the infrared spectra, most notably a strong blue shift of the symmetric hydrogen stretching fundamental relative to predictions from the harmonic approximation or vibrational second-order perturbation theory. Ab initio thermostated ring-polymer molecular dynamics simulations indicate that this blue shift originates from strong coupling between the hydrogen stretching and bending vibrations, resulting in an effective weakening of the OH···F^– ionic hydrogen bonds.

A DFT/PCM-based methodology for predicting solvolytic reactivities of organic carbonates

The possibility of employing quantum mechanical computations to predict solvolysis rates of benzhydryl aryl/alkyl (Ar/R) carbonates and to determine nucleofugalities of various Ar/R carbonate leaving groups in terms of N_f values is examined. Since unassisted S_N1 transition states of neutral substrates cannot be optimized by using implicit solvation models, a model reaction that includes anchimerically assisted heterolysis of 2-oxyethyl Ar/R carbonates is utilized to determine the relative reactivities of both benzhydryl Ar/R carbonate substrates and Ar/R carbonate leaving groups. Very good linear correlations have been obtained between activation free energies of the model reaction, calculated by using the M06-2X method in conjunction with the IEFPCM solvation model, and activation free energies in the literature for solvolysis of the corresponding benzhydryl Ar/R carbonates in a given solvent. The slopes of close to unity demonstrate that calculated and measured relative reactivities of the Ar/R carbonate leaving groups are practically the same. Very good agreement between experiment and theory has enabled extending the nucleofugality (N_f) scale established by Mayr and co-workers (Acc. Chem. Res., 2010, 43, 1537-1549) with numerous new Ar/R carbonate leaving groups.

Reaction Mechanism for Direct Proton Transfer from Carbonic Acid to a Strong Base in Aqueous Solution I: Acid and Base Coordinate and Charge Dynamics

Protonation by carbonic acid H2CO3 of the strong base methylamine CH3NH2 in a neutral contact pair in aqueous solution is followed via Car-Parrinello molecular dynamics simulations. Proton transfer (PT) occurs to form an aqueous solvent-stabilized contact ion pair within 100 fs, a fast time scale associated with the compression of the acid-base hydrogen-bond (H-bond), a key reaction coordinate. This rapid barrierless PT is consistent with the carbonic acid-protonated base pKa difference that considerably favors the PT, and supports the view of intact carbonic acid as potentially important proton donor in assorted biological and environmental contexts. The charge redistribution within the H-bonded complex during PT supports a Mulliken picture of charge transfer from the nitrogen base to carbonic acid without altering the transferring hydrogen's charge from approximately midway between that of a hydrogen atom and that of a proton.

Ammonia as an efficient catalyst for decomposition of carbonic acid: a quantum chemical investigation

Electronic structure calculations using M06-2X, MP2 and CCSD(T) methods have been employed to show ammonia as an efficient catalyst for decomposition of carbonic acid.

Solvation of CO2 in water: effect of RuBP on CO2 concentration in bundle sheath of C4 plants

An understanding of the temperature-dependence of solubility of carbon dioxide (CO2) in water is important for many industrial processes. Voluminous work has been done by both quantum chemical methods and molecular dynamics (MD) simulations on the interaction between CO2 and water, but a quantitative evaluation of solubility remains elusive. In this work, we have approached the problem by considering quantum chemically calculated total energies and thermal energies, and incorporating the effects of mixing, hydrogen bonding, and phonon modes. An overall equation relating the calculated free energy and entropy of mixing with the gas-solution equilibrium constant has been derived. This equation has been iteratively solved to obtain the solubility as functions of temperature and dielectric constant. The calculated solubility versus temperature plot excellently matches the observed plot. Solubility has been shown to increase with dielectric constant, for example, by addition of electrolytes. We have also found that at the experimentally reported concentration of enzyme RuBP in bundle sheath cells of chloroplast in C4 green plants, the concentration of CO2 can effectively increase by as much as a factor of 7.1-38.5. This stands in agreement with the observed effective rise in concentration by as much as 10 times.

Molecular dynamics simulations predict an accelerated dissociation of H2CO3 at the air–water interface

Ab initio molecular dynamics simulations of carbonic acid (H₂CO₃) at the air–water interface yield a lower dissociation barrier than in bulk water.