Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic equilibrium”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adsorption Thermodynamics for Process Simulation

Adsorption has rapidly evolved in recent decades and is an established separation technology extensively practiced in gas separation industries and others. However, rigorous thermodynamic modeling of multicomponent adsorption equilibrium remains elusive, and industrial practitioners rely heavily on expensive and time-consuming trial-and-error pilot studies to develop adsorption units. Here, this article highlights the need for rigorous adsorption thermodynamic models and the limitations and deficiencies of existing models such as the extended Langmuir isotherm, dual-process Langmuir isotherm, and adsorbed solution theory. It further presents a series of recent advances in the generalization of the classical Langmuir isotherm of single-component adsorption by deriving an activity coefficient model to account for the adsorbed phase adsorbate–adsorbent interactions, substituting adsorbed phase adsorbate and vacant site concentrations with activities, and extending to multicomponent competitive adsorption equilibrium, both monolayer and multilayer. Requiring a minimum set of physically meaningful model parameters, the generalized Langmuir isotherm for monolayer adsorption and the generalized Brunauer–Emmett–Teller isotherm for multilayer adsorption address various thermodynamic modeling challenges including adsorbent surface heterogeneity, isosteric enthalpies of adsorption, BET surface areas, adsorbed phase nonideality, adsorption azeotrope formation, and multilayer adsorption. Also discussed is the importance of quality adsorption data that cover sufficient temperature, pressure, and composition ranges for reliable determination of the model parameters to support adsorption process simulation, design, and optimization.

09 BIOMASS FUELS↗
DOE OSTI2025

Selection of solvents for integrated CO 2 absorption and electrochemical reduction systems

Abstract Solvent‐based electrochemical CO 2 reduction (CO 2 R) enables the production of chemicals or fuels using CO 2 from a preceding absorption process. Employing previously tested CO 2 capture solvents does not ensure their suitability for either CO 2 R or integrated CO 2 absorption‐reduction. We propose solvent selection criteria that include the CO 2 solubility, kinetic constant, ionic conductivity, concentration of the bicarbonate, carbamate, and solvent cation in the CO 2 ‐loaded solution, and sustainability indicators. They are implemented for solvent selection (a) from novel, aqueous mixtures of N ‐methylcyclohexylamine (MCA) with piperazine (PZ), 2‐amino‐2‐methyl‐1‐propanol (AMP), potassium hydroxide (KOH), and potassium chloride (KCl) and (b) from aqueous monoethanolamine (MEA), AMP, KOH, MCA, and PZ solutions. Versions of a modified Kent‐Eisenberg model for strong bases, carbamate, and non‐carbamate‐forming amine solutions are developed and parameterized through experimental equilibrium measurements. CO 2 R experimental results are presented for solutions of KOH and MCA + KOH, as these indicate desired trade‐offs for CO 2 absorption and reduction.

Amines↗
DOE OSTI2024

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-exchange Resin at Initial pH ~3 – 7 and Temperature 25 – 55 °C

This work presents systematic aqueous-phase adsorption equilibrium, kinetics, and column breakthrough measurements with three key biointermediates that are common compounds in many bioprocesses. Adsorption equilibrium experiments were carried out with acetic acid, butyric acid, and lactic acid on a commercial ion-exchange resin, Amberlite IRN-78, at wide ranges of acid concentration (8-500 mmol/L), initial pH (similar to 3-7), and temperature (25-55 degrees C), simulating the effluent characteristics from different fermenter operations. The kinetics and column breakthrough experiments were conducted at an initial pH of 6 and a concentration of 200 mmol/L. The equilibrium study shows a higher loading at the initial pH < pK(a) and a lower loading at the initial pH > pK(a). Overall removal varies between 16 and 99% depending on the initial pH, temperature, and organic acid concentration and type. The study further indicates monolayer adsorption at the equilibrium pH > 10 and multilayer adsorption at the equilibrium pH < 6. The thermodynamic modeling of adsorption isotherm data was carried out using Langmuir and Freundlich isotherms. IRN-78 presents fast adsorption kinetics as the maximum loading was attained in <= 10 min and nearly the same breakthrough time for all three organic acids involved in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Quantum Thermodynamics of Nonequilibrium Processes in Lattice Gauge Theories

A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2024

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Electrothermal Phenomena in Ferroelectrics

Physical properties of lead-zirconate-titanate (PZT) ceramics change according to the initial electric poling process and electrical boundary conditions. This paper reports the electrothermal, piezothermal, and piezoelectric coupling phenomena in ferroelectrics from thermodynamics viewpoints, in particular, thermal property differences between unpoled and poled PZT’s in the poling direction for open circuit and short circuit conditions. We propose a new terminology, “secondary electrothermal” coupling factor kλ, which is analogous to the electromechanical coupling factor k, relating the elastic compliances under short- and open-circuit conditions, in order to explain the fact that the short-circuit condition exhibited the larger thermal diffusivity than the open-circuit condition. On the other hand, the unpoled specimen exhibits the lowest thermal diffusivity. This tutorial paper was authored for providing comprehensive knowledge on equilibrium and time-dependent thermodynamics in ferroelectrics.

Uchino, Kenji (ORCID:0000000316517912)↗
DOE OSTI2023

Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI

The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Harnessing dimethyl ether with ultra-low-grade heat for scaling-resistant brine concentration and fractional crystallization

Solvent-driven separations may enable scalable concentration of hypersaline brines, supporting a circular resource economy from the extraction of lithium and rare earth elements from spent battery and magnet leachates. This work analyses a novel solvent-driven water extraction (SDWE) system employing dimethyl ether (DME) and ultra-low-grade heat for brine concentration and fractional crystallization. SDWE exploits DME’s unique properties: (1) a low dielectric constant that promotes water solubility over charged solutes by a factor of 10 3 , and (2) a high volatility that facilitate efficient DME reconcentration with ultra-low-grade heat. The techno-economic viability of SDWE is assessed with a computational framework that encompasses a liquid–liquid separator and a solvent concentrator. We integrate the extended universal quasichemical model with the virial equation of state to predict the compositions of the complex three-phase DME-water mixture at vapor–liquid and liquid–liquid equilibrium. Subsequently, we optimize the thermodynamic and economic performance of SDWE, by controlling the interstage flash pressure, heat source temperature, and the number of concentrating stages. DME-based SDWE concentrates an input saline feed to 5.5 M and regenerates over 99 % of the DME using ultra-low-grade heat below 50 °C, with a DME/water selectivity ratio of 125. Here our calculations reveal that optimal performance is achieved at interstage flash pressures of 0.4 – 0.5 bar for heat source temperatures between 323–373 K, with improved exergetic efficiencies at lower temperatures. At a heat source temperature of 323 K and an interstage pressure of 0.489 bar, DME-driven SDWE achieves an optimal thermodynamic efficiency of 20.5 % and a projected specific cost of US$ 1.93 m -3 . These specific costs suggest that SDWE is competitive with commercialized thermal distillation technologies, while mitigating the traditional risks associated with scaling in heat and mass exchangers with hypersaline brines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Quantum sensing of paramagnetic analytes by nanodiamonds in levitated microdroplets and aqueous solutions

Nanodiamonds (ND) hosting negatively charged nitrogen-vacancy (NV-) color centers have received attention for applications in magnetic field, electric field, chemical, and bio-sensing. The versatility of these probes is their excellent room-temperature optical and spin properties, along with their small size, functionalized surfaces and resistance to bleaching, making them ideal as nanoscopic sensors in picoliter volumes (e.g. single cells, but also microcompartments and aerosols). For quantitative ND-NV- sensing of paramagnetic analytes in such contexts, however, there remains an incomplete understanding of how factors related to the aqueous phase environment control detection efficiency. To address this, optically detected magnetic resonance (ODMR) is measured in bulk macroscale solutions and single levitated microdroplets as a function of Gd+3 concentration (340 nM to 1.5 mM), nanodiamond size, pH, competitor ions, and ligands. The ODMR response to [Gd+3] is found to be nonlinear, and pH, ND and sample volume dependent; indicating the detection of Gd+3 requires efficient adsorption of the analyte to the diamond surface. Langmuir adsorption isotherms embedded in a quantitative photophysical model links the ODMR response to adsorption thermodynamics of Gd+3. The equilibrium constant for Gd+3 adsorption to a carboxylated ND surface is determined to be (1 ± 0.5) x 105 M-1 corresponding to a free energy of adsorption of (-28 ± 1) kJ mol-1. These results provide general insight into how complex aqueous and microscale environments impact nanodiamond based quantum sensing modalities, and portend their application as quantitative chemical sensors in microenvironments.

Brown, Emily K↗
DOE OSTI2025

Specificity and tunability of efflux pumps: A new role for the proton gradient?

Efflux pumps that transport antibacterial drugs out of bacterial cells have broad specificity, commonly leading to broad spectrum resistance and limiting treatment strategies for infections. It remains unclear how efflux pumps can maintain this broad spectrum specificity to diverse drug molecules while limiting the efflux of other cytoplasmic content. We have investigated the origins of this broad specificity using theoretical models informed by the experimentally determined structural and kinetic properties of efflux pumps. We developed a set of mathematical models describing operation of efflux pumps as a discrete cyclic stochastic process across a network of states characterizing pump conformations and the presence/absence of bound ligands and protons. These include a minimal three-state model that lends itself to clear analytic calculations as well as a five-state model that relaxes some of the simpler model’s most strict assumptions. We found that the pump specificity is determined not solely by the drug affinity to the pump–as is commonly assumed–but it is also directly affected by the periplasmic pH and the transmembrane potential. Therefore, changes to the proton concentration gradient and voltage drop across the membrane can influence how effective the pump is at extruding a particular drug molecule. Furthermore, we found that while both the proton concentration gradient across the membrane and the transmembrane potential contribute to the thermodynamic force driving the pump, their effects on the efflux enter not strictly in a combined proton motive force. Rather, they have two distinguishable effects on the overall throughput. These results highlight the unexpected effects of thermodynamic driving forces out of equilibrium and illustrate how efflux pump structure and function are conducive to the emergence of multidrug resistance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Operational Limits of Redox Metal Oxides Performing Thermochemical Water Splitting

Solar thermochemical hydrogen production is an attractive technology that stores intermittent solar energy in the form of chemical bonds. Efficient operation requires the identification of a redox-active metal oxide (MO x ) material that can achieve high conversion of water to hydrogen at minimal energy input. Water splitting occurs by consecutive reduction and reoxidation reactions of MO x . MO x is reduced to MO x-δ and, in the second step, is reoxidized by water recovering the initial MO x and generate H 2 . The material must reduce at temperatures achievable in concentrated solar receiver/reactors, while maintaining a thermodynamic driving force to split water. At equilibrium, extent of reduction depends on temperature and oxygen partial pressure, and in this analysis, a set of thermodynamic properties, namely, enthalpy and entropy of oxygen vacancy formation, is sufficient to represent MO x . Herein, a method to easily classify materials based on these thermodynamic properties under any condition of oxygen partial pressure and temperature is presented. This method is based on fundamental thermodynamic principles and is applicable for any redox material with known thermodynamic properties. Despite the simplicity of the method, it is believed that this analysis will support future research in targeting thermodynamic properties of redox-active metal oxides.

water splitting↗
DOE OSTI2024

Catalytic Reduction of Esters over Zirconia-Supported Metal Catalysts

Esters are often produced as unwanted byproducts during the catalytic upgrading of ethanol to diesel fuel precursors through Guerbet coupling. Removal of esters from the product stream is important to prevent the loss of downstream catalyst activity from ester-derived carboxylic acids. In this work, we studied ester hydrogenolysis to the parent alcohols as a viable route for enhanced diesel fuel production. Specifically, we investigated the reduction of hexyl acetate in butanol over ZrO 2 -supported Ni, Co, Cu, Rh, Pd, and Pt catalysts, where Cu/ZrO 2 was the most selective catalyst for the hydrogenolysis of hexyl acetate into hexanol and ethanol. Thermodynamic analysis reveals that a 90% alcohol yield can be obtained at 200 °C, 30 bar, and a relatively high H 2 :hexyl acetate molar ratio of 480:1. Experimentally, an alcohol yield of 88% yield was obtained with a 10 wt % Cu/ZrO 2 catalyst at these conditions with a residence time of 5.4 h kg cat kmol gas –1 . Catalytic tests on the support revealed that ZrO 2 catalyzes the transesterification reaction between hexyl acetate and butanol. However, only the Cu sites can catalyze the hydrogenolysis of the esters into the final alcohols. We developed a kinetic model for our experimental results, which shows that the transesterification and hydrogenolysis reactions run at two different timescales, the former being 10 times faster than the latter. Data regression has been used to develop a model to predict the mole fraction distribution of ester hydrogenolysis products over a wide range of contact times. Cu/ZrO 2 loses half its catalytic activity after 80 h of time on stream. Modeling of deactivation data reveals that the ZrO 2 support conserves a residual activity due to external active sites, while active sites over the Cu surface deactivate at different rates. Furthermore, the catalytic conversion of esters into their parent alcohols is relevant to the production of surrogate liquid fuels since alcohols can be bimolecularly dehydrated to produce a blend of ethers with diesel fuel-like properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

LLEPE

LLEPE (liquid-liquid equilibrium parameter estimator) is a tool for estimating thermodynamic parameters like standard enthalpies, standard molar volumes and Pitzer coefficients, by combining Gibbs minimization with regression to multiphase equilibrium isotherm data. It also includes visualization tools to access data quality and prediction model performance.

ILOEJE, CHUKWUNWIKE↗
DOE OSTI2021

Mass action model of solution activity via speciation by solvation and ion pairing equilibria

Abstract Solutes and their concentrations influence many natural and anthropogenic solution processes. Electrolyte and solution models are used to quantify and predict such behavior. Here we present a mechanistic solution model based on mass action equilibria. Solvation and ion pairing are used to model speciated solute and solvent concentrations such that they correlate to a solution’s vapor pressure (solvent activity) according to Raoult’s law from dilute conditions to saturation. This model introduces a hydration equilibrium constant ( K ha ) that is used with either an ion dissociation constant ( K id ) or a hydration modifier ( m ) with an experimentally determined ion dissociation constant, as adjustable parameters to fit vapor–liquid equilibrium data. The modeled solvation equilibria are accompanied by molecular dynamics (MD) studies that support a decline in the observed degree of solvation with increased concentration. MD calculations indicate this finding is a combination of a solvent that solvates multiple solutes, and changes in a solute’s solvation sphere, with the dominant factor changing with concentration. This speciation-based solution model is lateral to established electrostatics-based electrolyte theories. With its basis in mass action, the model can directly relate experimental data to the modeled solute and solvent speciated concentrations and structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2024

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗