Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stack Sampling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Weak Lensing Mass Calibration of the ACT DR5 Galaxy Clusters with the DES Year 3 Weak Lensing Data

We use weak gravitational lensing measurements from Year 3 Dark Energy Survey data to calibrate the masses of 443 galaxy clusters selected via the Sunyaev-Zel'dovich effect from Atacama Cosmology Telescope Data Release 5 maps of the cosmic microwave background. We incorporate redshift and SZ measurements for individual clusters into a hierarchical model for the stacked lensing signals and perform Bayesian analyses to constrain the hydrostatic mass bias of the clusters. Our treatment of systematic uncertainties includes a prescription for measuring and accounting for the weak lensing boost factor, consideration of a miscentering effect, as well as marginalization over uncertainties in the source galaxy photometric redshift distributions and shear calibration. The resultant constraints on the normalization of the mass-observable relation have a precision of approximately 7%, with the mean WL halo mass of M $_{500c}$ = 5.4 × 10$^{14}$ M $_{⊙}$. We measure the bias between the true cluster mass and the mass estimated from the SZ signal based on an X-ray-calibrated scaling relation assuming hydrostatic equilibrium, to be 1 - b = 0.74$^{+0.06}$ $_{-0.05}$ over the full sample. When splitting the clusters into high (z = 0.43-0.70) and low (z = 0.15-0.43) redshift bins, we measure 1 - b = 0.58$^{+0.06}$ $_{-0.05}$ and 0.81$^{+0.08}$ $_{-0.06}$, respectively. When introducing additional freedom in redshift and mass to the hydrostatic bias model, we find that 1 - b decreases with redshift (with the power law of -1.8$^{+0.5}$ $_{-0.6}$, 99.95% confidence), consistent with findings from other recent studies, while we do not find any significant trend in mass. We also demonstrate that our result is robust against various systematics such as a scale cut, priors on baryonic and miscentering parameters, and degree of scatter in mass-observable relation. The weak-lensing mass calibration presented in this study will be a useful tool for using the ACT clusters as probes of astrophysics, and as a step towards using their abundance as a cosmological probe.

Shin, T. [Carnegie Mellon U.] (ORCID:0000000263895↗
DOE OSTI2026

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗
DOE OSTI2026

H I Depletion Begins Well Beyond the Virial Radius: A FAST Stacking Study of 36 Galaxy Clusters to 5 × R 200

Abstract We present a stacking study of the neutral atomic hydrogen (H i ) content in and around 36 local galaxy clusters at z < 0.07, using a combination of the FAST All Sky H i survey (FASHI) and the extensive spectroscopic catalog mainly from the Dark Energy Spectroscopic Instrument (DESI). We employ spectral stacking techniques to probe the average H i mass and HI-to-stellar mass ratio ( M HI / M * ) for member galaxies down to stellar masses of M * ∼ 10 9 M ⊙ , spanning a projected cluster-centric distance of up to 5 R 200 . Our analysis reveals a pronounced environmental effect; both M HI and M HI / M * decrease steadily toward the cluster center, dropping by ∼0.5 dex on average from the outskirts to the core. Crucially, we find that M HI / M * of galaxies remain lower than the field galaxies even at the 5 R 200 . This provides direct, statistical evidence for substantial gas stripping and preprocessing in the cluster outskirts, likely occurring in infalling groups and large-scale filaments. By further splitting the sample by g − r color, we show that the H i deficiency persists at fixed galaxy color; even the bluest cluster members exhibit ∼0.5 dex lower M HI / M * than field galaxies of similar color, reflecting environmental effects on the cold gas reservoir prior to full optical transformation. The total H i mass within clusters and their outskirts agrees broadly with predictions from cosmological simulations. Our results underscore the critical role of the extended cluster environment in quenching galaxies by depleting their cold gas reservoirs well before they enter the dense cluster core.

Cheng, Cheng [Chinese Academy of Sciences South Am↗
DOE OSTI2025

Strength, deformation, and the fcc–hcp phase transition in condensed Kr and Xe to the 100 GPa pressure range

The rare gas solids exhibit systematic differences in crystal structure, phase transition conditions, bond strength, and other physical properties. The physical properties of heavy rare gas solids krypton and xenon are modified by the martensitic phase transition from face-centered cubic to hexagonal close packed structure over a broad pressure range. Crystal structure, strength, and plastic deformation of krypton and xenon have been investigated at 300 K using compression in the diamond-anvil cell with synchrotron angle-dispersive x-ray diffraction and complementary ruby fluorescence spectroscopy for Xe. Stacking faults indicative of the fcc–hcp phase transition are observed at pressures at and above 1.23 ± 0.05 and 1.9 ± 0.6 GPa in Kr and Xe, respectively. The transition remains incomplete in both solids to pressures greater than 100 GPa. Strength determined from stress measurements in Pt and ruby standards at pressures up to 111 GPa and complemented by observations of strain and texture measurements obtained by x-ray diffraction in the radial geometry to 100 GPa indicates similar or higher strength than Ar at all conditions, with significant stiffening at 15–20 GPa. Radial diffraction data reveal the persistence of broad highly textured fcc diffraction lines to 101 GPa in Xe, suggesting that the axial measurements may underestimate the metastable persistence of the fcc phase due to biased sampling of hcp crystallites resulting from preferred crystallite orientation. Kr and Xe are compared with He, Ne, and Ar for a systematic understanding of physical properties and phase equilibria of rare gas solids.

Compressive stress↗
DOE OSTI2025

Sustainable Shape Memory Elastomers with Reduced Melt Viscosity and Enhanced Stiffness

Melt reactive processing of lignin with nitrile rubber is a promising approach to synthesizing shape memory materials. The strong intramolecular interactions in lignin macromolecular structures, caused by π–π stacking in aromatic rings and hydrogen bonding, often result in large phase separation or low miscibility with rubbers. In this study, we investigated the chemical and molecular characteristics, as well as the stiffness and complex viscosity, of modified kraft lignin melt-reacted with an acrylonitrile/butadiene copolymer containing 41% acrylonitrile (NBR41). To enhance the macromolecular compatibility of kraft lignin with NBR41, kraft lignin was cross-linked with poly(propylene glycol) diglycidyl ether (PPDE) and trimethylolpropane triglycidyl ether (TTE), both rich in epoxy reactive groups capable of forming chemical bonds with hydroxyl and carboxyl groups. Here, our findings demonstrate that the modification of kraft lignin with PPDE and TTE resulted in significantly increased stiffness of the composites. The elastic modulus of NBR41-Kraft lignin-PPDE and NBR41-Kraft lignin-TTE increased by 82 and 162%, respectively. Both the yield strength and Young’s modulus of these two samples showed dramatic improvements. Specifically, the yield strength and Young’s modulus of NBR41-Kraft lignin-TTE increased nearly 4 and 3-fold, respectively, compared to the control sample. Interestingly, despite significant improvements in mechanical properties, the viscosity of NBR41-Kraft lignin-PPDE was substantially lower than that of the control sample. At 210 °C and an angular frequency of 1 rad/s, the complex viscosity of NBR41-Kraft lignin was approximately 100.25 ± 4.77 kPa·s, while that of NBR41-Kraft lignin-PPDE was significantly lower at 56 ± 0.93 kPa·s. These findings were validated through Fourier transform infrared spectroscopy, scanning electron microscopy, dynamic mechanical analysis, thermal characterization, rheological tests, and quasi-elastic neutron scattering techniques.

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

Mechanistic insights into rubidium ion adsorption at the quartz (101) surface from quantum chemical metadynamics

Ion sorption extent and mechanism depend in part on the mineral surface termination, which can be highly complex. Variations in surface functional groups, particularly with defect density and surface roughness, influence mineral reactivity towards solutes. In this work, we investigate the adsorption of a rubidium cation (Rb + ) at pristine and defect quartz (101) surface sites using well-tempered metadynamics, based on simulations with the quantum chemical density-functional tight-binding (DFTB) method. We compare the relative energetics of Rb + adsorption across selected sites for each surface, with nanosecond-level sampling, highlighting similarities between vicinal and geminal silanol sites. We find that the positive Rb + partial atomic charge can increase by as much as 0.5e as it approaches the surface, with implications for modulation of ion adsorption strength and extent at the quartz (101) surface with surface vacancies, silanol coverage, and charge.

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

Design and Commissioning of a Deuterium-Tritium Gas Delivery System for Muon Catalyzed Fusion in a Diamond Anvil Cell

We report the design, commissioning, and operation of deuterium-deuterium (DD) and deuterium-tritium (DT) gas delivery systems developed to load a diamond anvil cell (DAC) beam target for muon-catalyzed fusion (muCF). The DAC approach enables DT fuel to be compressed to GPa pressures at more than twice the liquid density and heated from cryogenic temperatures through 500 K, opening access to a substantially expanded parameter range for muCF kinetics and yield measurements. In this approach, DT is cryo-condensed to a liquid in a minichamber and then compressed in the DAC using a helium-driven pneumatic membrane, achieving high pressures in a millimeter-scale DT sample volume. A DD gas delivery system was designed and used to validate the experimental apparatus, measure the gas quantities needed for filling, develop operational experience, and collect kinetics and yield data with DD targets. The DT gas delivery system adds tritium-specific capabilities for inventory minimization, secondary containment, and activity monitoring. The DT system integrates depleted uranium storage beds and a liquid helium cryogenic condenser used for pressure building and cryopumping. High-purity delivery is provided by a rapid-response palladium permeator. The system is housed in a helium-atmosphere glovebox held at negative pressure with continuous cleanup. We present the process and instrumentation design, a failure modes and effects analysis (FMEA), and data from the experiment's in situ Raman spectrometer, which provides direct confirmation of target loading and composition through the optically clear diamond anvils. The 2024 and 2025 DT campaigns achieved repeatable target fills and operation with no measurable tritium releases to the stack, demonstrating safe, high-purity DT loading at novel density-temperature conditions for muCF studies.

Koukina, Elena [Acceleron Fusion]↗
DOE OSTI2025

Highly Tunable Valley Polarization of Potential-Trapped Moiré Excitons in WSe 2 /WS 2 Heterojunctions

Moiré superlattices created by stacking atomic layers of transition metal dichalcogenide semiconductors have emerged as a class of fascinating artificial photonic and electronic materials. An appealing attribute of these structures is the inheritance of the valley degree of freedom from the constituent monolayers. Recent studies show evidence that the valley polarization of the moiré excitons is highly tunable. In heterojunctions of WSe 2 /WS 2 , marked improvement in valley polarization is observed by increasing optical excitation power, a behavior that is quite distinct from the monolayers, and lacks a clear understanding so far. Here, in this Letter, we show that this highly tunable valley property arises from filling of the moiré superlattice, which provides an intriguing mechanism for engineering these quantum opto-valleytronic platforms. Our data further demonstrate that the long-range electron-hole exchange interaction, despite being significantly weakened in the junctions, is the dominant source of moiré exciton intervalley scattering at low population. Using magnetic field tuning, we quantitatively determine the exchange interaction strength to be 0.03 and 0.24 meV for 0° and 60° twisted samples, respectively, in our experiments, about 1 order of magnitude weaker than that in the monolayers.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A new data-driven map predicts substantial undocumented peatland areas in Amazonia

Tropical peatlands are among the most carbon-dense terrestrial ecosystems yet recorded. Collectively, they comprise a large but highly uncertain reservoir of the global carbon cycle, with wide-ranging estimates of their global area (441 025–1700 000 km 2 ) and below-ground carbon storage (105–288 Pg C). Substantial gaps remain in our understanding of peatland distribution in some key regions, including most of tropical South America. Here we compile 2413 ground reference points in and around Amazonian peatlands and use them alongside a stack of remote sensing products in a random forest model to generate the first field-data-driven model of peatland distribution across the Amazon basin. Our model predicts a total Amazonian peatland extent of 251 015 km 2 (95th percentile confidence interval: 128 671–373 359), greater than that of the Congo basin, but around 30% smaller than a recent model-derived estimate of peatland area across Amazonia. The model performs relatively well against point observations but spatial gaps in the ground reference dataset mean that model uncertainty remains high, particularly in parts of Brazil and Bolivia. For example, we predict significant peatland areas in northern Peru with relatively high confidence, while peatland areas in the Rio Negro basin and adjacent south-western Orinoco basin which have previously been predicted to hold Campinarana or white sand forests, are predicted with greater uncertainty. Similarly, we predict large areas of peatlands in Bolivia, surprisingly given the strong climatic seasonality found over most of the country. Very little field data exists with which to quantitatively assess the accuracy of our map in these regions. Data gaps such as these should be a high priority for new field sampling. This new map can facilitate future research into the vulnerability of peatlands to climate change and anthropogenic impacts, which is likely to vary spatially across the Amazon basin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗
DOE OSTI2024

Velocity reconstruction in the era of DESI and Rubin/LSST. II. Realistic samples on the light cone

Reconstructing the galaxy peculiar velocity field from the distribution of large-scale structure plays an important role in cosmology. On one hand, it gives us an insight into structure formation and gravity; on the other, it allows us to selectively extract the kinetic Sunyaev-Zel’dovich (kSZ) effect from cosmic microwave background maps. In this work, we employ high-accuracy synthetic galaxy catalogs on the light cone to investigate how well we can recover the velocity field when utilizing the three-dimensional spatial distribution of the galaxies in a modern large-scale structure experiment such as the Dark Energy Spectroscopic Instrument (DESI) and the Rubin Observatory Legacy Survey of Space and Time. In particular, we adopt the standard technique used in baryon acoustic oscillation analysis for reconstructing the Zel’dovich displacements of galaxies through the continuity equation, which yields a first-order approximation to their large-scale velocities. We investigate variations in the number density, bias, mask, area, redshift noise, and survey depth, as well as modifications to the settings of the standard reconstruction algorithm. Since our main goal is to provide guidance for planned kSZ analysis between DESI and the Atacama Cosmology Telescope, we apply velocity reconstruction to a faithful representation of DESI spectroscopic and photometric targets. We report the cross-correlation coefficient between the reconstructed and the true velocities along the line of sight. For the DESI Y1 spectroscopic survey, we expect the correlation coefficient to be r ≈ 0.64, while for a photometric survey with δ z /(1+z) = 0.02, as is approximately the case for the Legacy Survey used in the target selection of DESI galaxies, r shrinks by half to r ≈ 0.31. Here, we hope the results in this paper can be used to inform future kSZ stacking studies and other velocity reconstruction analyses planned with the next generation of cosmology experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Precision Kinematic Sunyaev--Zel'dovich Measurements Across Halo Mass and Redshift with DESI DR2 and ACT DR6: Part II. Bright Galaxy Survey and Emission-Line Galaxies

We present the first high-significance spectroscopic stacked kinetic Sunyaev-Zel'dovich (kSZ) measurements of circumgalactic gas profiles for both Bright Galaxy Survey (BGS) and Emission Line Galaxy (ELG) tracers, combining DESI Data Release 2 with ACT Data Release 6. Using reconstructed line-of-sight velocities from the DESI galaxies and high-resolution ACT temperature maps, we detect the kSZ signal at high significance, reaching signal-to-noise ratios of up to $\sim$9 for BGS and $\sim$7.5 for ELGs in optimal stellar-mass selections. Together with the LRG measurements presented in Paper I, these constitute the most significant kSZ detections from any spectroscopic survey to date. We perform the analysis in both real and harmonic space, obtaining consistent results. By splitting both tracers into stellar-mass bins, we study the scaling of the kSZ amplitude with galaxy properties. Combining the kSZ measurements with ACT Data Release 6 (DR6) CMB lensing maps enables a joint calibration of the galaxy-halo connection and the gas fractions of host halos. For the BGS galaxies, we observe low gas fractions around the virial radius relative to standard expectations, likely attributable to active galactic nuclei (AGN) activity. We find some evidence for higher-mass halos retaining a larger fraction of their baryons, consistent with more efficient feedback in lower-mass systems. For the ELG sample, dominated by blue, star-forming galaxies, we provide the first detection of the gas distribution in ELG host halos. The ELGs appear to exhibit relatively high gas fractions, which points to the possibility of weaker feedback (due to e.g. low AGN and supernova feedback activity) at their mass scale. Finally, we present generalized Navarro-Frenk-White (GNFW) fits to the harmonic-space measurements, providing a compact parametrization of gas profiles for forward modeling in large-scale structure analyses.

Hadzhiyska, B. [Cambridge U.; Cambridge U., Inst. ↗
DOE OSTI2025

Galaxy cluster profiles: a Gaussian mixture model approach to halo miscentering

Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center defined as the minimum in the potential well. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys. Our Python package for the Gaussian mixture model is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model.

Bayesian reasoning↗
DOE OSTI2025

Atomic-Scale Behavior of Radiation-Resistant ZnO under High-Energy Electron Bombardment

Understanding the atomic structure and defect characteristics of ZnO thin films is crucial for optimizing their electronic properties and performance in advanced applications. Here, we investigate the atomic structure and defect characteristics of atomic layer deposition (ALD)-grown ZnO thin films by using aberration-corrected scanning transmission electron microscopy (STEM). Atomic-resolution imaging identifies prevalent stacking faults, dipole disorder, and various grain boundary types, which are believed to influence the electronic properties of ZnO. Additionally, real-time electron beam exposure experiments demonstrate structural transformations, including crystal growth and surface rearrangements. These findings provide insights into the growth mechanisms of ALD ZnO under high-energy electron irradiation conditions, an important finding for the use of polycrystalline ZnO wide bandgap semiconductors in space-like conditions. In conclusion, our results underscore the capability of STEM in directly visualizing and quantifying atomic-scale defects and beam-induced transformations in radiation-resistant ZnO.

Defects↗
DOE OSTI2025

Inhibition of Reaction Layer Formation on MgO(100) by Doping with Trace Amounts of Iron

Despite extensive research on MgO’s reactivity in the presence of CO 2 under various conditions, little is known about whether impurities incorporated into the solid, such as iron, enhance or impede hydroxylation and carbonation reactions. The purity of the MgO required for the successful implementation of MgO looping as a direct air capture technology affects the deployment costs. With this motivation, we tested how incorporated iron impacts MgO (100) reactivity and passivation layer formation under ambient conditions by using atomic force microscopy, electron microscopy, and synchrotron-based X-ray scattering. Based on electron microprobe analysis, our MgO samples were 0.5 wt % iron, and Mössbauer spectroscopy results indicated that 70% of the iron is present as Fe(II). We find that even these low levels of iron dopants impeded both the hydroxylation at various relative humidities (10%, 33%, 75%, and >95%) and carbonation in CO 2 (33%, 75%, and >95%) on the (100) surface. Crystalline reaction products were formed. Reaction layers on the sample were easily removed by exposing the sample to deionized water for 2 min. Overall, our findings demonstrate that the presence of iron dopants slows the reaction rate of MgO, indicating that MgO without incorporated iron is preferable for mineral looping applications.

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

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

58 GEOSCIENCES↗
DOE OSTI2025

Reaction-Induced Fracturing of Porous Carbonate Rocks during Volume-Increasing Replacement by Witherite

Volume-increasing replacement reactions can lead either to fracturing of the parent phase or to the formation of cohesive layers that passivate further reaction, but the factors that drive one outcome or the other are not understood. In this experimental study, we investigated the volume-increasing replacement of carbonate rocks (consisting of CaCO 3 and CaMg(CO 3 ) 2 ) of different porosity by witherite (BaCO 3 ). Samples were characterized using scanning electron microscopy, Raman spectroscopy, small-angle neutron/X-ray scattering, X-ray tomography, and scanning transmission electron microscopy. We observed the formation of a witherite reaction layer and witherite formation within pores and along grain boundaries. Despite this being a volume-increasing replacement reaction, newly formed witherite was porous, potentially allowing further replacement. Filled fractures were observed in the low-porosity carbonates, whereas witherite formed within pores in high-porosity carbonates. We conclude that fracturing of the parent phase versus passivation is contingent on the initial microstructure of the parent with an optimal degree of porosity required for fracturing.

Calcite↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗