The fate of carbon dioxide in water-rich fluids under extreme conditions

Platinum Carbonates in Aqueous Fluids under Extreme Conditions

Platinum is among the rarest elements on the planet, and the understanding of its formation and transport through aqueous fluids in the Earth, at high pressure and temperature, may help in the identification of new deposits. While complexation of platinum with sulfides, chlorides, and hydroxyl has been the topic of numerous investigations, the interaction of Pt and carbonates in aqueous fluids under pressure remains largely unexplored. Here, we present extensive first-principles molecular simulations of Pt (bi)carbonates at conditions (1 GPa, 1000 K and 11 GPa, 1000 K) relevant to the Earth crust and upper mantle and we predict how the metal speciation varies as a function of pressure and how it depends on its oxidation state. Furthermore, we compute Raman spectra and identify vibrational signatures that may be used to characterize the varied species in solutions. Our simulations provide valuable inputs to the Deep Earth Water model.

Unveiling hidden reaction kinetics of carbon dioxide in supercritical aqueous solutions

Dissolution of CO2 in water followed by the subsequent hydrolysis reactions is of great importance to the global carbon cycle, and carbon capture and storage. Despite numerous previous studies, the reactions are still not fully understood at the atomistic scale. Here, we combined ab initio molecular dynamics (AIMD) simulations with Markov state models to elucidate the reaction mechanisms and kinetics of CO2 in supercritical water both in the bulk and nanoconfined states. The integration of unsupervised learning with first-principles data allows us to identify complex reaction coordinates and pathways automatically instead of a priori human speculation. Interestingly, our unbiased modeling found an unknown pathway of dissolving CO2(aq) under graphene nanoconfinement, involving the pyrocarbonate anion [C2O[Formula: see text](aq)] as an intermediate state. The pyrocarbonate anion was previously hypothesized to have a fleeting existence in water; however, our study reveals that it is a crucial reaction intermediate and stable carbon species in the nanoconfined solutions. We even observed the formation of pyrocarbonic acid [H2C2O5(aq)], which was unknown in water, in our AIMD simulations. The unexpected appearance of pyrocarbonates is related to the superionic behavior of the confined solutions. We also found that carbonation reactions involve collective proton transfer along transient water wires, which exhibits concerted behavior in the bulk solution but proceeds stepwise under nanoconfinement. The first-principles Markov state models show substantial promise for elucidating complex reaction kinetics in aqueous solutions. Our study highlights the importance of large oxocarbons in aqueous carbon reactions, with great implications for the deep carbon cycle and the sequestration of CO2.

Synthesis and Stability of Biomolecules in C-H-O-N Fluids under Earth's Upper Mantle Conditions

How life started on Earth is an unsolved mystery. There are various hypotheses for the location ranging from outer space to the seafloor, subseafloor, or potentially deeper. Here, we applied extensive ab initio molecular dynamics simulations to study chemical reactions between NH3, H2O, H2, and CO at pressures (P) and temperatures (T) approximating the conditions of Earth's upper mantle (i.e., 10-13 GPa, 1000-1400 K). Contrary to the previous assumptions that large organic molecules might readily disintegrate in aqueous solutions at extreme P-T conditions, we found that many organic compounds formed without any catalysts and persisted in C-H-O-N fluids under these extreme conditions, including glycine, ribose, urea, and uracil-like molecules. Particularly, our free-energy calculations showed that the C-N bond is thermodynamically stable at 10 GPa and 1400 K. Moreover, while the pyranose (six-membered ring) form of ribose is more stable than the furanose (five-membered ring) form at ambient conditions, we found that the formation of the five-membered-ring form of ribose is thermodynamically more favored at extreme conditions, which is consistent with the exclusive incorporation of β-d-ribofuranose in RNA. We have uncovered a previously unexplored pathway through which the crucial biomolecules could be abiotically synthesized from geofluids in the deep interior of Earth and other planets, and these formed biomolecules could potentially contribute to the early stage of the emergence of life.

Nanoconfinement facilitates reactions of carbon dioxide in supercritical water

The reactions of CO2 in water under extreme pressure-temperature conditions are of great importance to the carbon storage and transport below Earth's surface, which substantially affect the carbon budget in the atmosphere. Previous studies focus on the CO2(aq) solutions in the bulk phase, but underground aqueous solutions are often confined to the nanoscale, and nanoconfinement and solid-liquid interfaces may substantially affect chemical speciation and reaction mechanisms, which are poorly known on the molecular scale. Here, we apply extensive ab initio molecular dynamics simulations to study aqueous carbon solutions nanoconfined by graphene and stishovite (SiO2) at 10 GPa and 1000 ~ 1400 K. We find that CO2(aq) reacts more in nanoconfinement than in bulk. The stishovite-water interface makes the solutions more acidic, which shifts the chemical equilibria, and the interface chemistry also significantly affects the reaction mechanisms. Our findings suggest that CO2(aq) in deep Earth is more active than previously thought, and confining CO2 and water in nanopores may enhance the efficiency of mineral carbonation.

Photoelectron spectra of water and simple aqueous solutions at extreme conditions

Determining the electronic structure of aqueous solutions at extreme conditions is an important step towards understanding chemical bonding and reactions in water under pressure (P) and at high temperature (T)....

Carbon Dioxide under Earth Mantle Conditions: From a Molecular Liquid through a Reactive Fluid to Polymeric Regimes

In both its gaseous and condensed forms, carbon dioxide has an ever-increasing impact on Earth's chemistry and human life and activities. However, many aspects of its high-pressure phase diagram remain unclear. In this work, we present a complete structural characterization of carbon dioxide fluids under geological conditions using extensive ab initio molecular dynamics simulations throughout a wide pressure and temperature range, corresponding to Earth's lower mantle. We identify and describe four different disordered regimes, including two polymeric forms and two molecular ones, all within the geothermal conditions of the lower mantle. At pressures below 40 GPa, we find that the molecular liquid becomes very reactive above 2000 K: the C-O double bond routinely breaks, resulting in small and transient chains composed of CO₂ units and frequently leading to an exchange of oxygen atoms between molecules. At higher pressures, in addition to the polymeric fluid previously reported at 3000 K, we find a polymeric system with glass-like behavior at lower temperatures, suggesting a complex interplay between kinetics and stability.

Water-Gas Shift Reaction Produces Formate at Extreme Pressures and Temperatures in Deep Earth Fluids

The water-gas shift reaction is one of the most important reactions in industrial hydrogen production and plays a key role in Fischer-Tropsch-type synthesis, which is widely believed to generate hydrocarbons in the deep carbon cycle but is little known at extreme pressure-temperature conditions found in the Earth's upper mantle. Here, we performed extensive ab initio molecular dynamics simulations and free energy calculations to study the water-gas shift reaction. We found the direct formation of formic acid from CO and supercritical water at 10-13 GPa and 1400 K without any catalyst. Contrary to the common assumption that formic acid or formate is an intermediate product, we found that HCOOH is thermodynamically more stable than the products of the water-gas shift reaction above 3 GPa and at 1000-1400 K. Our study suggests that the water-gas shift reaction may not happen in the Earth's upper mantle, and formic acid or formate may be an important carbon carrier in reducing environments, participating in many geochemical processes in deep Earth.

Solvation of simple ions in water at extreme conditions

The interaction of ions and water at high pressure and temperature plays a critical role in Earth and planetary science yet remains poorly understood. Aqueous fluids affect geochemical properties ranging from water phase stability to mineral solubility and reactivity. Here, we report first-principles molecular dynamics simulations of mono-valent ions (Li⁺, K⁺, Cl^-) as well as NaCl in liquid water at temperatures and pressures relevant to the Earth's upper mantle (11 GPa, 1000 K) and concentrations in the dilute limit (0.44-0.88 m), in the regime of ocean salinity. We find that, at extreme conditions, the average structural and vibrational properties of water are weakly affected by the presence of ions, beyond the first solvation shell, similar to what was observed at ambient conditions. We also find that the ionic conductivity of the liquid increases in the presence of ions by less than an order of magnitude and that the dielectric constant is moderately reduced by at most ∼10% at these conditions. Our findings may aid in the parameterization of deep earth water models developed to describe water-rock reactions.

Dissociation of salts in water under pressure

The investigation of salts in water at extreme conditions is crucial to understanding the properties of aqueous fluids in the Earth. We report first principles (FP) and classical molecular dynamics simulations of NaCl in the dilute limit, at temperatures and pressures relevant to the Earth’s upper mantle. Similar to ambient conditions, we observe two metastable states of the salt: the contact (CIP) and the solvent-shared ion-pair (SIP), which are entropically and enthalpically favored, respectively. We find that the free energy barrier between the CIP and SIP minima increases at extreme conditions, and that the stability of the CIP is enhanced in FP simulations, consistent with the decrease of the dielectric constant of water. The minimum free energy path between the CIP and SIP becomes smoother at high pressure, and the relative stability of the two configurations is affected by water self-dissociation, which can only be described properly by FP simulations.

Salts in water at extreme conditions play a fundamental role in determining the properties of the Earthʼs mantle constituents. Here the authors shed light on ion-water and ion-ion interactions for NaCl dissolved in water at conditions relevant to the Earthʼs upper mantle by molecular dynamics simulations.

Carbon dioxide, bicarbonate and carbonate ions in aqueous solutions under deep Earth conditions

We investigate the effect of pressure, temperature and acidity on the composition of water-rich carbon-bearing fluids under thermodynamic conditions that correspond to the Earth's deep crust and upper mantle. Our first-principles molecular dynamics simulations provide mechanistic insight into the hydration shell of carbon dioxide, bicarbonate and carbonate ions, and into the pathways of the acid/base reactions that convert these carbon species into one another in aqueous solutions. At temperatures of 1000 K and higher, our simulations can sample the chemical equilibrium of these acid/base reactions, thus allowing us to estimate the chemical composition of diluted carbon dioxide and (bi)carbonate ions as a function of acidity and thermodynamic conditions. We find that, especially at the highest temperature, the acidity of the solution is essential to determine the stability domain of CO₂vs. HCO₃^-vs. CO₃^2-.

A first principles method to determine speciation of carbonates in supercritical water

The determination of the speciation of ions and molecules in supercritical aqueous fluids under pressure is critical to understanding their mass transport in the Earth's interior. Unfortunately, there is no experimental technique yet available to directly characterize species dissolved in water at extreme conditions. Here we present a strategy, based on first-principles simulations, to determine ratios of Raman scattering cross-sections of aqueous species under extreme conditions, thus providing a key quantity that can be used, in conjunction with Raman measurements, to predict chemical speciation in aqueous fluids. Due to the importance of the Earth's carbon cycle, we focus on carbonate and bicarbonate ions. Our calculations up to 11 GPa and 1000 K indicate a higher concentration of bicarbonates in water than previously considered at conditions relevant to the Earth's upper mantle, with important implications for the transport of carbon in aqueous fluids in the Earth's interior.

A chemical perspective on high pressure crystal structures and properties

The general availability of third generation synchrotron sources has ushered in a new era of high pressure research. The crystal structure of materials under compression can now be determined by X-ray diffraction using powder samples and, more recently, from multi-nano single crystal diffraction. Concurrently, these experimental advancements are accompanied by a rapid increase in computational capacity and capability, enabling the application of sophisticated quantum calculations to explore a variety of material properties. One of the early surprises is the finding that simple metallic elements do not conform to the general expectation of adopting 3D close-pack structures at high pressure. Instead, many novel open structures have been identified with no known analogues at ambient pressure. The occurrence of these structural types appears to be random with no rules governing their formation. The adoption of an open structure at high pressure suggested the presence of directional bonds. Therefore, a localized atomic hybrid orbital description of the chemical bonding may be appropriate. Here, the theoretical foundation and experimental evidence supporting this approach to the elucidation of the high pressure crystal structures of group I and II elements and polyhydrides are reviewed. It is desirable and advantageous to extend and apply established chemical principles to the study of the chemistry and chemical bonding of materials at high pressure.

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.

Recent Development of CO2 Electrochemistry from Li-CO2 Batteries to Zn-CO2 Batteries

Metal-CO₂ batteries with CO₂ as cathode active species give rise to opportunities to deal with energy and environmental issues simultaneously. This technology is more appealing when CO₂ is flexibly reduced to chemicals and fuels driven by surplus electricity because it represents a low-cost and controllable approach to maximized electricity utilization and value-added CO₂ utilization. Nonaqueous metal-CO₂ batteries exhibited high discharge voltage and capacity with carbon and oxalate as reduction products from CO₂ electrochemistry that lacks proton. In contrast, aqueous Zn-CO₂ batteries implemented flexible CO₂ electrochemistry for more value-added products accompanied by energy storage based on a proton-coupled electron transfer mechanism. In this Account, we have exemplified our recent results in the development of CO₂ electrochemistry from nonaqueous Li-CO₂ batteries to aqueous Zn-CO₂ batteries toward practical value-added CO₂ conversion. Aimed at the challengingly limited CO₂ electrochemistry and high cost of nonaqueous Li-CO₂ batteries, we proposed aqueous Zn-CO₂ batteries. Our previous works on nonaqueous Li-CO₂ batteries, aqueous Zn-air batteries, and aqueous CO₂ reduction electrocatalysts further shed light on battery mechanism, device construction, and electrocatalyst design. For example, bipolar membranes maintain the stability of the basic anolyte and neutral catholyte, as well as the kinetics of ion transport at the same time, forming the device base for aqueous Zn-CO₂ batteries. Moreover, in terms of the electrocatalyst catalyzing both discharge and charge reactions on the cathode, the design of multifunctional electrocatalysts is of great importance for not only CO₂ electrochemistry but also spontaneous discharge and energy efficiency of aqueous Zn-CO₂ batteries. We have explored a series of multifunctional electrocatalyst cathodes, including noble metal, transition metal, and metal-free materials, all of which facilitated CO₂ electrochemistry in aqueous Zn-CO₂ batteries with value-added carbon-based products. Meanwhile, several operating models for practical complicated situations are presented, such as rechargeable, reversible, dual-model, and solid-state batteries. Zn-CO₂ batteries with different models require different design mechanisms for electrocatalyst cathodes. Reversible aqueous Zn-CO₂ batteries with HCOOH generation were enabled by electrocatalysts capable of catalyzing the interconversion of CO₂ and HCOOH at low overpotentials, rechargeable aqueous Zn-CO₂ batteries were allowed by electrocatalysts capable of catalyzing efficient CO₂ reduction and O₂ evolution, and dual-model aqueous Zn-CO₂ batteries were realized by electrocatalysts capable of catalyzing CO₂ reduction, water oxidation, and oxygen reduction. Concluding remarks include a summary of recent CO₂ electrochemistry in metal-CO₂ batteries and a brief discussion of future challenges and opportunities for practical aqueous Zn-CO₂ batteries, such as highly reduced products and high production rate.

Ab initio spectroscopy and ionic conductivity of water under Earth mantle conditions

The phase diagram of water at extreme conditions plays a critical role in Earth and planetary science, yet remains poorly understood. Here we report a first-principles investigation of the liquid at high temperature, between 11 GPa and 20 GPa-a region where numerous controversial results have been reported over the past three decades. Our results are consistent with the recent estimates of the water melting line below 1,000 K and show that on the 1,000-K isotherm the liquid is rapidly dissociating and recombining through a bimolecular mechanism. We found that short-lived ionic species act as charge carriers, giving rise to an ionic conductivity that at 11 GPa and 20 GPa is six and seven orders of magnitude larger, respectively, than at ambient conditions. Conductivity calculations were performed entirely from first principles, with no a priori assumptions on the nature of charge carriers. Despite frequent dissociative events, we observed that hydrogen bonding persists at high pressure, up to at least 20 GPa. Our computed Raman spectra, which are in excellent agreement with experiment, show no distinctive signatures of the hydronium and hydroxide ions present in our simulations. Instead, we found that infrared spectra are sensitive probes of molecular dissociation, exhibiting a broad band below the OH stretching mode ascribable to vibrations of complex ions.

Theoretical Understanding of Mechanisms of Proton Exchange Membranes Made of 2D Crystals with Ultrahigh Selectivity

Recent reports on proton conduction across pristine graphene and hexagonal boron nitride (h-BN) provide a new avenue for the design of proton exchange membranes. The uniform pores formed by the electron clouds of two-dimensional (2D) crystals can effectively block the undesired transportation of other species thus ultrahigh selectivity can be achieved. With the aid of first-principles calculations, we investigate the proton conduction process across six kinds of intact 2D crystals, namely graphene, h-BN, β₁₂ boron sheet, χ₃ boron sheet, phosphorene, and silicene. To clarify the proton conduction mechanism, three proton penetration modes are proposed: dissociation-penetration, adsorption-penetration, and direct penetration. Based on our calculation results, for graphene and h-BN without atomic defects, they are unlikely to provide sufficient proton conductivity at room temperature when no bias potential is applied. By contrast, the β₁₂ boron sheet, χ₃ boron sheets, and silicene exhibit relatively lower proton penetration energy barriers, making them prospective candidates for future proton exchange membrane applications.