Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Leveraging Exchange Kinetics for the Synthesis of Atomically Precise Porous Catalysts

Metal-organic frameworks (MOFs) have attracted significant attention as porous catalyst platforms due to the synthetic modularity of these materials and the diversity of lattice-confined catalytic active sites that are readily embedded within periodic crystalline frameworks. MOFs offer platforms to heterogenize molecular catalysts, stabilize novel coordination motifs, and leverage confinement effects in catalysis. Crystallinity allows diffraction-based methods to be employed in the characterization of these catalysts. Access to crystalline MOFs typically requires reversible construction of the metal-ligand (M-L) bonds that connect the SBUs, which provides a mechanism to anneal defects during crystallization. While the required M-L bond reversibility is often promoted by synthesis at elevated temperature, access to crystalline materials based on either transition metals with characteristic slow exchange kinetics or highly basic donor ligands remains a synthetic challenge. Here in this paper, we highlight synthetic strategies that leverage M-L exchange kinetics to access MOFs based on kinetically inert ions and extensions of these strategies to the assembly of atomically precise multimetallic materials.

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

Experimental Investigation of the Thermal Decomposition Pathways and Kinetics of TATB by Isotopic Substitution

Real-time measurements of the product gases arising from the thermal decomposition of triamino-trinitro benzene (TATB), its deuterated analogue, and plastically bonded TATB (LX-17) are presented in this study. Gas-phase decomposition products are identified by IR absorption spectroscopy. The frequency shifts in rovibrational spectra due to isotopic substitution and the change in rate of formation of decomposition products due to the kinetic-isotope-effect (KIE) help elucidate the decomposition pathways. The formation of H 2 O precedes other molecules (e.g., HCN, HNCO) during decomposition. After the concentrations of HCN and HNCO molecules reach a peak, their amounts gradually decrease. The concentrations of the other decomposition products (e.g., NH 3 and CO 2 ) rapidly rise after an induction period, which is attributed to the presence of autocatalytic reactions. The trends of chemical evolution are similar for all the samples, but their kinetic behaviors are different. This indicates the rates of consistent pathways are changed during thermal decomposition. Additionally, the kinetics of deuterated TATB decomposition is slower than that of unsubstituted TATB due to the KIE (kH/kD~1.41). The rate of LX-17 decomposition is slightly lower than unsubstituted TATB (k TATB /k LX-17 ~1.15). The KIE is more pronounced during the early stage of decomposition, which is attributed to the first steps of TATB decomposition involving water formation (i.e., H vs D transfer). The KIE slows down the formation of all gases, including those lacking hydrogen (e.g., CO 2 ). These results suggest the TATB thermal decomposition mechanism might involve a series of pathways rather than a set of independent and parallel reactions.

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

Thermal Decomposition Kinetics of 4,6‐Diamino‐5,7‐dinitro‐benzo‐furazan

This experimental study investigated the thermal decomposition kinetics of 4,6-diamino-5,7-dinitro-benzo-furazan (referred to as F1 hereafter)—an important decomposition product of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB—a prototypical insensitive high explosive). Simultaneous differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), and mass spectrometry (MS) measurements were employed to determine the decomposition kinetics of F1 and to track the evolution of product gases. The DSC profiles were measured at 10 different heating rates between 0.025°C/min and 10°C/min. The measured exotherms were influenced by F1 melting at heating rates above 0.25°C/min, and corresponding changes in decomposition enthalpy and TGA mass-loss-rate profiles indicated a transition from solid-to-gas decomposition to an increasing contribution from liquid-to-gas decomposition. Analysis of low-heating-rate DSC data between 0.025°C/min and 0.17°C/min with the extended Prout–Tompkins model yielded an activation energy of 305 kJ/mol for solid-to-gas F1 decomposition, higher than previous values inferred from TATB decomposition models involving F1. This study provides the first direct experimental determination of the energy barrier for F1 decomposition. MS measurements showed that the major gaseous products matched species previously reported for TATB decomposition (e.g., CO 2 , HCN, C 2 N 2 , etc.), with water identified as the dominant product. Furthermore, these results provide important experimental constraints for improving chemical kinetics models of TATB decomposition and for predicting the reactivity, stability, and safety of TATB-based high explosives under long-term aging conditions and abnormal thermal environments.

4,6-Diamino-5,7-dinitro-benzo-furazan↗
DOE OSTI2021

Lessons Learned in Employing Data Analytics to Predict Oxidation Kinetics and Spallation Behavior of High-Temperature NiCr-Based Alloys

Machine learning (ML) can offer many advantages in predicting material properties over traditional materials development methods based solely on limited experimental investigations or physical-based simulations with the capability to reduce development cost, risk, and time. However, so far, limited efforts have been made to predict alloy oxidation kinetics and spallation behavior via ML due to the lack of consistently measured and sufficient experimental data and the inherent complexity in oxidation behavior of multicomponent high-temperature alloys. A previous study reported the ability of ML to predict oxidation kinetics of NiCr-based alloys as a function of alloy composition and operating conditions. Here, the performance of a ML model in predicting rate constants and spallation probability was evaluated in light of the roles of the data distribution of the experimental dataset (data analytics), the alloy composition, the exposure environment and the chosen oxidation approach to extracting kinetic values from the measured mass changes (but using either a simple parabolic law or a statistical cyclic oxidation model). Potential strategies to improve the predictions and enhance the extrapolative capability of the previously trained model will be discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Microstructural Coarsening Kinetics and Mechanical Property Changes in Long-Term Aged Sn–Pb–Sb Solder Joints

Tin-lead-antimony (50Sn–47Pb–3Sb wt.%) soldered assemblies were mechanically tested approximately 30 years after initial production and found to have solder joints of reduced strength. The microstructure of this solder alloy exhibits a ternary eutectic structure with Sn-rich, Pb-rich, and SnSb phases. Accelerated aging was performed to evaluate solder microstructural coarsening and associated strength of laboratory solder joints to correlate these properties to the “naturally aged” solder joints. Isothermal aging was conducted at room temperature, 55, 70, 100, and 135 °C and aging times that ranged from 0.1 to 365 days. The coarsening kinetics of the Pb-rich phase were determined through optical microscopy and image analysis methods established in previous studies on binary Sn–Pb solder. A kinetic equation was developed with time exponent n of 0.43 and activation energy of 24000 J/mol, suggesting grain boundary diffusion or other fast diffusion pathways controlling the microstructural evolution. Compression testing and Vickers microhardness showed significant strength loss within the first 20–30 days after soldering; then, the microstructure and mechanical properties changed more slowly over long periods of time. Further, by combining accelerated aging data and the microstructure-based kinetics, strength predictions were made that match well with the properties of the actual soldered assemblies naturally aged for 30 years. However, aging at the highest temperature of 135 °C produced anomalous behavior suggesting that extraneous aging mechanisms are active. Therefore, data obtained at this temperature or higher should not be used. Overall, the combined microstructural and mechanical property methods used in this study confirmed that the observed reduction in strength of ~ 30-year-old solder joints can be accounted for by the microstructural coarsening that takes place during long-term solid-state aging.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O 2 and H 2 O 2 ), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R 2 > 0.95). Kinetic observations also enabled a direct comparison between O 2 and H 2 O 2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

3-hydroxypropionic acid↗
DOE OSTI2022

Addition of transient kinetics capabilities to an infrared reflection absorption spectroscopy system through synchronized gas pulsing and data acquisition

Surface science methodologies for understanding thermodynamic aspects of surface processes are at an advanced level. However, instrumentation and approaches for extracting kinetic parameters from elementary steps are far less accessible. In this work, we present an approach combining the use of a fast gas pulsing valve synchronized with data acquisition to enable surface transient kinetics studies using infrared reflection absorption spectroscopy. This methodology applies to the study of reversible processes and borrows concepts and ideas from molecular beam scattering and temporal analysis of products. Here, a temporal resolution of ~67 ms is achieved, and this is illustrated through the study of CO adsorption and desorption on a Pd(111) crystal in the presence and absence of background O 2 . The same approach can be extended to other surface spectroscopies, such as X-ray photoelectron spectroscopy, to obtain spectra with high temporal resolution and signal-to-noise ratio and enable future multimodal surface transient kinetic studies aiming at elucidating reaction mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reduction kinetics of hercynite redox materials for solar thermochemical water splitting

Solar thermochemical water splitting (STWS) using a hercynite redox cycle is a promising technology for producing renewable H 2 . Here, the reduction kinetics of hercynite (FeAl 2 O 4 ) is evaluated using thermogravimetric and XRD analyses. The results indicate that as-prepared hercynite materials undergo reduction via two different reaction mechanisms. The reaction first proceeds by a nucleation and growth mechanism (AE1), followed by a third-order kinetic model (F3). XRD analyses show the occurrence of superstoichiometric oxygen in the spinel structure of FeAl 2 O 4+δ in the second reaction mechanism, which indicates the formation of cationic vacancies. TGA and XRD analyses confirm that hercynite maintains its spinel structure when the materials are thermally reduced and oxidized with steam. Differential scanning calorimetry (DSC) performed on FeAl 2 O 4 particles shows two heats of reaction in agreement with the two observed mechanisms. This is the first kinetics study to report on the thermal reduction of hercynite.

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

A Transient-Response methodology based on experiments and modeling for Cu-Redox Half-Cycle kinetic analysis on a Cu-SSZ-13 SCR catalyst

In this report the Standard SCR reaction catalyzed by Cu-SSZ-13 is a redox process consisting of a reduction half cycle (RHC) and an oxidation half cycle (OHC) that cycle the active Cu sites between the Cu(II) and Cu(I) states. In the current work, a transient-response methodology consisting of experimental transient response Cu-redox (TRCR) measurements and kinetic modeling was developed for detailed study of individual SCR Cu-redox half cycles. The TRCR protocol allows quantification of the reducible Cu density, Cu(II)-Cu(I) partitioning, relative native RHC and OHC rates, and limiting half cycle during SCR. The half-cycle kinetics are studied over a wide (170–450 °C) temperature, and in differential segments along the catalyst length using spatially resolved capillary inlet mass-spectrometry (SpaciMS). The protocol alone provides two independent measures indicating that OHC increases faster than RHC with temperature, and that at all but the lowest temperature SCR is RHC limited. Introducing the transient kinetic models allows half-cycle reaction pathways, orders and activation energies to be studied and determined; e.g., RHC involves surface and gas-phase NH 3 routes, and OHC with O 2 involves dimer formation. Application of the full methodology provides further confirmation and quantification of the protocol conclusions, and specifically the first measurements of RHC and OHC activation energies; specifically for parallel RHC surface and gas-phase routes, and experimentally obtained OHC activation energy vs. first-principle calculations in the literature. Beyond insights related to the degreened catalyst studied here, the methodology provides a broadly available approach for quantifying how aging processes selectively impact RHC and OHC.

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

Insights into co-pyrolysis of polyethylene terephthalate and polyamide 6 mixture through experiments, kinetic modeling and machine learning

The non-isothermal pyrolysis of polyethylene terephthalate (PET), polyamide 6 (PA6), and their mixtures was studied in a thermogravimetric analyzer at different heating rates. Temperature of maximum decomposition (T max ) decreased by 25–45 °C and 35–55 °C for the PET:PA6 mixtures (3:1, 1:1, 1:3) compared to PET and PA6, respectively. The kinetic analysis was initially carried out using isoconversional method. However, the dependency of activation energy on conversion was observed for the co-pyrolysis of PET and PA6 that suggested the occurrence of multi-step reactions in the mixtures. Distributed activation energy model (DAEM) was used in this study to describe the multistep reactions occurring during pyrolysis of PET:PA6 mixtures. Here, in this work, a four-parallel reaction DAEM was developed to describe the pyrolysis kinetics of PET:PA6 mixtures. The apparent mean activation energies (E o ) for PET, PA6, and mixtures varied in the range of 244–255, 140–215, and 138–255 kJ mol –1 , respectively. The mass loss profiles of PET and PA6 mixtures were also modeled using artificial neural network (ANN). Out of 155 ANN models, the best prediction was made by ANN511 with R 2 greater than 0.997 for both test and unseen data. The interaction effects observed through TGA experiments and subsequent kinetic analysis were further assessed in terms of product composition using analytical pyrolysis coupled with gas chromatograph/mass spectrometer (Py-GC/MS). Co-pyrolysis of PET and PA6 resulted in the formation of new aromatic compounds with nitrogen-containing functional groups, which were not detected when PET or PA6 were pyrolyzed individually.

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

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2021

An experimental and detailed kinetic modeling study of the pyrolysis and oxidation of allene and propyne over a wide range of conditions

Allene and propyne are important intermediates in the pyrolysis and oxidation of higher hydrocarbon fuels, and they are also a major source of propargyl radical formation, which can recombine into different C 6 H 6 isomers and finally produce soot. In a prior work (Panigrahy et al., “A comprehensive experimental and improved kinetic modeling study on the pyrolysis and oxidation of propyne”, Proc. Combust. Inst 38 (2021)), the pyrolysis, ignition , and laminar flame speed of propyne were investigated. To understand the kinetic features of initial fuel breakdown and oxidation of the two C 3 H 4 isomers, new measurements for allene pyrolysis and oxidation are conducted in the present paper at the same operating conditions as those studied previously for propyne. Ignition delay times of allene are measured using a high-pressure shock tube and a heated twin-opposed piston rapid compression machine in the temperature range 690– 1450 K at equivalence ratios of 0.5, 1.0 and 2.0 in ‘air’, and at pressures of 10 and 30 bar. Pyrolysis species measurements of allene and propyne are also performed using a gas chromatography integrated singlepulse shock tube in the temperature range 1000–1700 K at pressure of 2 and 5 bar. Furthermore, laminar flame speeds of allene are measured at elevated gas temperatures of 373 K at pressures of 1 and 2 bar for a wide range of equivalence ratios from 0.6 to 1.5. A newly updated kinetic mechanism developed for this study is the first model that can well reproduce all of the experimental results for both allene and propyne. It is observed that in the pyrolysis process, allene dissociates faster than propyne. Both isomers exhibit similar ignition delay times at high temperatures (>1000 K), while, at intermediate temperatures (770–1000 K) propyne is the faster to ignite, and at lower temperatures (< 770 K) allene becomes more reactive. Furthermore, laminar flame speeds for propyne are found to be slightly faster than those for allene under the conditions studied in this work

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

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

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

A wide range experimental study and further development of a kinetic model describing propane oxidation

This work presents an experimental and kinetic modeling study of propane oxidation. Ignition delay times of propane were measured in a high-pressure shock tube and in rapid compression machines in the temperature range 689 - 1700 K at equivalence ratios of 0.5, 1.0 and 2.0 in 'air', for a wide range of pressures from 20 to 90 bar. CO and H 2 O mole fraction profiles for propane oxidation were measured in a shock tube behind reflected shock waves in the temperature range 1370-1840 K at equivalence ratios of 0.5, 1.0 and 2.0 and at a pressure of approximately 1.3 atm. Moreover, propane oxidation was studied using a jet-stirred reactor coupled to a synchrotron vacuum ultraviolet photoionization mass spectrometer at low temperatures in the range 565 - 690 K and at a pressure of 1 atm. Further, this wide range of experimental datasets for propane oxidation was used to reoptimize and update our previous kinetic mechanisms, AramcoMech3.0 and NUIGMech1.1. In the current mechanism, NUIGMech1.3, the thermo-chemical parameters of all species relevant to low-temperature propane oxidation chemistry, including propyl-peroxyl, hydroperoxyl-propyl, hydroperoxyl-propyl-peroxyl, and carbonyl-hydroperoxide radicals, are updated based on newly calculated values at the CCSD(T)-F12/TZ-F12//B2PLYPD3/TZ///B2PLYP-D3/TZ level of theory. The improvements made in the thermochemical values and in the kinetic parameters for the low-temperature propane oxidation reactions in NUIGMech1.3 result in better model agreement with the new IDTs and speciation data, including carbon monoxide, formaldehyde, propene, acetaldehyde and various minor products such as ethylene, acetic acid, acrolein as well as various hydroperoxide and cyclic ether species.

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

Kinetic Monte Carlo simulations of structural evolution during anneal of additively manufactured materials

Our experiments indicated that upon a post-processing anneal, an additively manufactured 316L stainless steel exhibits cubic grains rather than the conventional equiaxed grains. In this work, we have used kinetic Monte Carlo simulations to explore the origin of these cubic grains. First, we implemented a new kinetic Monte Carlo model in parallel code SPPARKS to simulate grain growth and recrystallization under a residual energy distribution. Our model incorporates physical properties and real-time, as opposed to generic properties and relative time. We further validated that our SPPARKS simulations reproduced the expected kinetic behavior of single-grain evolution. We then used the validated approach to simulate the anneal of an additively manufactured material under the same conditions used in our experiments. We found that the cubic grains can origin from a periodically varying residual energy that may be present in additively manufactured materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Biomass char gasification kinetic rates compared to data, including ash effects

This investigation provides both experimental and theoretical analyses of biomass char gasification kinetics over 1150–1350 °C in the environment of CO 2 , H 2 O and combinations of these reactants. The experimental data include continuous measurements of particle mass, external and internal temperatures, shape, and size with periodic but not continuous measurements of char porosity, pore size distribution, and pore surface area. The model describes heterogeneous char gasification for poplar wood, corn stover and switchgrass in different shapes and industrially relevant sizes (12.5 mm and smaller). Nonlinear least-squares regression produces optimized power-law model parameters for global kinetics that describe biomass char gasification with respect to both CO 2 and H 2 O separately and in combination. A single set of kinetic parameters describes reaction rates for all three fuels. Model simulations agree with measured data at all stages of char conversion even though the observed mass loss rates change significantly with conversion. In conclusion, the change in rate with conversion depends strongly on ash content, and this investigation provides a simple and theoretically based reaction rate expression that includes the effect of inert ash occupying an increasing portion of the reactive surface with increased conversion.

09 BIOMASS FUELS↗
DOE OSTI2021

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

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

Dual carbonate clumped isotopes (Δ 47 -Δ 48 ) constrains kinetic effects and timescales in peridotite-associated springs at the Cedars, Northern California

The Cedars is an area in Northern California with a chain of highly alkaline springs resulting from CO 2 -charged meteorological water interacting with a peridotite body. Serpentinization resulting from this interaction at depth leads to the sequestration of various carbonate minerals into veins accompanied by a release of Ca 2+ and OH – enriched water to the surface, creating an environment which promotes rapid precipitation of CaCO 3 at surface springs. This environment enables us to apply the recently developed Δ 47 -Δ 48 dual clumped isotope analysis to probe kinetic isotope effects (KIEs) and timescales of CO 2 transformation in a region with the potential for geological CO 2 sequestration. Here, we analyzed CaCO 3 recovered from various localities and identified significant kinetic fractionations associated with CO 2 absorption in a majority of samples, characterized by enrichment in Δ 47 values and depletion in Δ 48 values relative to equilibrium. Surface floes exhibited the largest KIEs (ΔΔ 47 : 0.163‰, ΔΔ 48 : -0.761‰). Surface floe samples begin to precipitate out of solution within the first hour of CO 2 absorption, and the dissolved inorganic carbon (DIC) pool requires a residence time of >100 h to achieve isotopic equilibria. The Δ 48 /Δ 47 slope of samples from the Cedars (-3.223 ± 0.519) is within the range of published theoretical values designed to constrain CO 2 hydrolysis-related kinetic fractionation (-1.724 to -8.330). The Δ 47 /δ 18 O slope (-0.009 ± 0.001) and Δ 47 /δ 13 C slope (-0.009 ± 0.001) are roughly consistent with literature values reported from a peridotite in Oman of -0.006 ± 0.002 and -0.005 ± 0.002, respectively. The consistency of slopes in the multi-isotope space suggests the Δ 47 -Δ 48 dual carbonate clumped isotope framework can be applied to study CO 2 -absorption processes in applied systems, including sites of interest for geological sequestration.

58 GEOSCIENCES↗