Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temperature effect”

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 163 records · Page 9

Lignin Extraction and Condensation as a Function of Temperature, Residence Time, and Solvent System in Flow-through Reactors

Solvolytic extraction of lignin from biomass is a critical step in lignin-first biorefining, including the reductive catalytic fractionation (RCF) process. Key to optimal RCF processing is the ability to rapidly extract lignin from biomass at high delignification extents and transfer the lignin molecules to a catalyst surface in a time frame that minimizes lignin condensation reactions. Here, we use a flow-through reactor to study the effects of temperature (175–250 °C), residence time (9 to 36 min), and solvent composition (methanol and methanol–water) on lignin extraction and condensation. We evaluated three metrics at each condition: total delignification, delignification rate, and extent of condensation, the latter measured by a decrease in monomer yield for batch hydrogenolysis reactions of solvolysis liquor compared to batch RCF reactions. We observe that delignification is predominantly determined by temperature, while residence time dictates the lignin condensation extent. Moreover, the extent of both extraction and condensation increased in the methanol–water solvent system compared to that in the methanol system. Lignin extracted in methanol is stable up to 18-min residence times at or below 225 °C, while a majority of the lignin extracted in methanol–water is condensed with a 9-min residence time at 200 °C. These results can inform reactor designs and solvent selection for lignin-first biorefining processes that aim to physically separate the biomass and catalyst.

biorefining↗

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗

An Acid-Free, Temperature-Based Cation Contamination Removal Strategy for PEM Water Electrolysis

It is widely understood that the durability and reliability of polymer electrolyte membrane (PEM) water electrolyzers are heavily dependent on feedwater purity, with cation contaminants that originate from incomplete water purification and balance of plant materials significantly harming electrolyzer performance. However, contamination remains a challenge and a common cause of failure at the stack level, indicating the need for strategies to recover the performance of contaminated cells. In this study, we investigate the effects of temperature on the uptake, electrochemical impacts, and removal of contaminant calcium and iron cations. Lower operating temperatures increase the sensitivity of the cell performance to contaminant cations, while also decreasing cation uptake and promoting contaminant removal. Computational charge transfer modelling shows that lower temperature increases the concentration of contaminant at the cathode and facilitates their removal from the cell. By testing single cells under scenarios designed to mimic stack temperature dynamics, we investigate low-temperature operation as an approach to stack-relevant contaminant recovery. Together, these results demonstrate that the low-temperature recovery approach is a promising approach for acid-free contamination recovery for PEM water electrolysis to promote stack reliability and durability.

08 HYDROGEN↗

Temperature-induced hexagonal–orthorhombic phase transition in lutetium ferrite nanoparticles

The x-ray diffraction, Raman, and infrared spectroscopies and magnetic measurements were used to explore the correlated changes of the structure, lattice dynamics, and magnetic properties of the LuFeO3 nanoparticles, which appear in dependence on their sintering temperature. We revealed a gradual substitution of the hexagonal phase by the orthorhombic phase in the nanoparticles, with sintering temperature increasing from 700 to 1100 °C. The origin and stability of the hexagonal phase in the LuFeO3 nanoparticles are of the special interest, because the nanoparticles in the phase can be a room-temperature multiferroic with a weak ferromagnetic and pronounced structural and ferroelectric long-range ordering. The antiferromagnetic and nonpolar orthorhombic phase is more stable in the bulk LuFeO3. To define the ranges of the hexagonal phase stability, we determine the bulk and interface energy densities of different phases from the comparison of the Gibbs model with experimental results. Using effective parameters of the Gibbs model, we predict the influence of size effects and temperature on the structural and polar properties of the LuFeO3 nanoparticles. Analysis of the obtained results shows that the combination of the x-ray diffraction, Raman and infrared spectroscopies, magnetic measurements, and theoretical modeling of structural and polar properties allows us to establish the interplay between the phase composition, lattice dynamics, and multiferroic properties of the LuFeO3 nanoparticles prepared under different conditions.

Materials Science↗

Cold air quench control of local crystallization environment in fully air-processed carbon-based perovskite solar cells

Carbon-based perovskite solar cells (C-PSCs) present a low-cost route to efficient photovoltaic technology. Gas quenching is an essential process commonly used for solvent extraction in scalable module fabrication using solution processes, but film quality remains highly sensitive to the localized environmental processing temperature, especially during intermediate phase transitions before perovskite film high temperature annealing. In this study, we introduce a cold air quench strategy to precisely control the local crystallization temperature that simultaneously promotes uniform crystallization and facilitates strain relaxation in fully air-processed C-PSCs. It is found that implementing dry air at 10 °C as the quenching gas yields a champion power conversion efficiency (PCE) of 20.52 %, with preferential (110) orientation and reduced, homogenized in-plane strain. This newly developed technique of a cost-effective low temperature air quenching process also enhances carrier lifetime, reduces interface recombination, and improves charge extraction. Furthermore, this work presents cold air quenching as a scalable, economic method to control perovskite crystallization and strain in C-PSC fabrication, advancing the industrial viability of high-performance perovskite modules.

14 SOLAR ENERGY↗

Degradation and Accelerated Recovery of Surface Passivation in n+ Poly-Si/SiOx Passivating Contacts for TOPCon Solar Cells

We studied the surface degradation and recovery of fired poly-Si/SiOx passivating contacts during subsequent dark and illuminated annealing. We report on an industrially-viable path for accelerated recovery of surface passivation. The degradation is influenced by the type of doping in the poly-Si. Phosphorus doped n+ poly-Si/SiOx contacts show degradation followed by recovery, undoped poly-Si/SiOx contacts only show recovery during annealing. Boron doped p+ poly-Si/SiOx contacts show neither degradation nor improvement with annealing. Both degradation and recovery are thermally-activated processes and are completely reversible and cyclic in nature. The activation energy of degradation and recovery in dark for n+ poly-Si/SiOx contacts are 1.27 and 1.33 eV respectively. Dark annealing at elevated temperatures is effective for complete recovery but takes a long time (~30 min at 350 degrees C) due to higher activation energies. Annealing under 7.5 Suns of illumination lowers the activation energy for degradation and recovery to 0.88 and 0.90 eV, respectively. Using this data, we have developed an industrially viable post-firing treatment for accelerated recovery of TOPCon cells by annealing them at elevated temperatures and under intense illumination for a few minutes.

14 SOLAR ENERGY↗

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

From viruses to protists: temperature response of the neglected components of microbial controls on peatland nutrient cycling

The response of plant-associated microbial communities to rising temperatures likely plays a key role in global Carbon (C) and Nitrogen (N) biogeochemical cycles. Mosses of the genus Sphagnum, in particular, sequester 25% of all terrestrial C as recalcitrant peat. However, their ability to sequester C is mediated by interactions with microbes that fix N and pass it to the moss to grow in otherwise oligotrophic systems. Two important biotic controls on this processes are likely to respond to rising temperatures –predation by protists and infections by viruses– but whether and how this might occur is not well understood. This proposal addressed these questions using a combination of genomics and mathematical modeling with both field data from whole-ecosystem warming experiments, and laboratory-based microcosm experiments. We have discovered 1) mechanistic links between cell traits and thermal performance in protists, which 2) enabled predictions from populations to ecosystems, 3) we have identified eco-phenotypic feedbacks that determine protist–bacterial predator–prey dynamics, 4) have documented concomitant compositional and trait shifts in protist communities with warming in the field, 5) revealed consistent interactive effects of temperature and nutrients on bacterial and protist communities; 6) characterized the Sphagnum virome, and 7) made modeling predictions on its implications for thermal responses of microbial communities in peatlands globally. This award also substantially strengthened U.S. scientific capacity and workforce development.

Gibert, Jean Philippe [Duke University]↗

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a high-temperature Inconel 625 heat exchanger by model design and binder jetting additive manufacturing

A nickel-based Inconel 625 superalloy heat exchanger for high-temperature applications was developed via binder jetting additive manufacturing. The material properties were characterized first on printed and sintered parts. Two sintering temperatures were used to investigate the effects of the temperature on densification and microstructure. For channel geometry design, the heat transfer capabilities of the heat exchanger were optimized for the cross-section geometry of fluid flow channels in a counterflow configuration, and a stress analysis was conducted to investigate the effects of channel geometry. After headers were incorporated into the heat exchanger as a one-piece component, a prototype was printed. A complete depowdering was achieved via compressed air blowing, and sintering was performed with the developed profile. Ultimately, a one-piece heat exchanger with a nearly full density was obtained. Heat transfer tests were performed on this unit, and the results were compared with those from the simulations. In this study, systematic processes were developed for an additively manufactured Inconel-based heat exchanger for high-pressure, high-temperature heat transfer applications.

36 MATERIALS SCIENCE↗

Spontaneously formed phonon frequency combs in van der Waals solid CrGeTe 3 and CrSiTe 3

Optical phonon engineering through nonlinear effects has been utilized in ultrafast control of material properties. However, nonlinear optical phonons typically exhibit rapid decay due to strong mode-mode couplings, limiting their effectiveness in temperature or frequency sensitive applications. Here we report the observation of long-lived nonlinear optical phonons through the spontaneous formation of phonon frequency combs in the van der Waals material CrXTe 3 (X=Ge, Si) using high-resolution Raman scattering. Unlike conventional optical phonons, the highest A g mode in CrGeTe 3 splits into equidistant, sharp peaks forming a frequency comb that persists for hundreds of oscillations and survives up to 200K. These modes correspond to localized oscillations of Ge 2 Te 6 clusters, isolated from Cr hexagons, behaving as independent quantum oscillators. Introducing a cubic nonlinear term to the harmonic oscillator model, we simulate the phonon time evolution and successfully replicate the observed comb structure. Similar frequency comb behavior is observed in CrSiTe 3 , demonstrating the generalizability of this phenomenon. Our findings demonstrate that Raman scattering effectively probes high-frequency nonlinear phonon modes, offering insight into the generation of long-lived, tunable phonon frequency combs with potential applications in ultrafast material control and phonon-based technologies.

36 MATERIALS SCIENCE↗

Measurements of the electron energy distribution function in partially magnetized low temperature plasmas

While Langmuir probes (LPs) are relatively simple and inexpensive plasma diagnostics for the electron density, temperature, and the electron energy distribution function (EEDF), the interpretation of the measured current–voltage (I–V) characteristic is complicated considerably by the presence of a magnetic field. In regimes where the electron gyroradius is comparable to the probe radius, the electron flux to the probe surface is retarded by reduced mobility across field lines and can no longer be described by a thermal model. Predicting the current collected by the probe in these regimes requires accurate estimates of the plasma diffusion coefficients, which are usually difficult to obtain. In this work, we measure electron energy distribution functions in E×B plasmas with magnetized electrons and non-magnetized ions in argon and krypton gases, using both a LP and laser Thomson scattering (LTS) at various magnetic fields. Using the LTS measurements to provide a robust benchmark for comparison, we compare existing theories describing the flux to the probe under magnetized conditions. We find that even when the electron gyroradius associated with the effective electron temperature is small compared to the probe radius, the EEDF computed using classical probe theory is still robust at energies higher than the energy at which the gyroradius becomes larger than the probe size. For plasmas that are not strongly non-Maxwellian, we formulate a method to extract robust density and temperature measurements using physics-informed fitting techniques to analyze EEDFs computed using classical theory.

Devin, E. G. [Princeton Plasma Physics Laboratory ↗

Study of phase decoherence in GeSn (8%) through measurements of the weak antilocalization effect

Alloying germanium with tin offers a means to modulate germanium's electronic structure, enabling a greater degree of control over quantum properties such as the retention of the phase or spin of the electron wave. However, the extent to which the presence of high dopant concentrations in GeSn alters these quantum behaviors is poorly understood. Here, we investigate the role of dopant concentrations on phase coherence through measurements of the weak antilocalization (WAL) effect at temperatures between 30 mK and 10 K in p-GeSn (8%) thin films, which were doped to a series of carrier densities on the order of 10 12 cm –2 . Phase coherence and spin–orbit lengths were extracted from the magnetoconductivities using the 2D Hikami–Larkin–Nagaoka model. Phase coherence lengths peaked at 577, 593, and 737 nm for the low-, mid-, and high-density samples, while upper limits on the spin–orbit lengths of less than 25 nm were relatively independent of carrier density and temperature. The phase coherence lengths increased as the temperature decreased but changed only minimally with carrier density, contrary to common models of temperature-dependent inelastic scattering. Saturation of the phase coherence lengths occurred below 600 mK. Based on these findings, intrinsically generated inelastic scattering mechanisms such as two-level systems or impurity band scattering likely contribute to phase decoherence in these alloys. Our results provide insight into the inelastic scattering mechanisms of GeSn, while suggesting a need for further investigation into phase decoherence mechanisms in doped group-IV alloys.

42 ENGINEERING↗

A Simple Route for Open Fluidic Devices with Particle Walls

Open fluidics, allowing liquid in a flow channel to interact with the external environment, is a revolutionary concept. However, fabricating a highly stable open fluidic device of arbitrary complexity, while maintaining reconfigurability, is still a challenge. This is achieved by the use of a patterned substrate and liquids that are covered with functional, readily available hydrophobic particles, providing great flexibility in the construction and use of open fluidic structures. Decorated with a coating of modified carbon nanotubes (CNTs) to encapsulate the fluids, the study capitalizes on the photothermal characteristics of CNTs to fabricate a device to probe the effects of temperature on tumor chemotherapy. The strategy substantially increases the availability and potential use of open fluidic devices.

Liu, Heng↗

Dynamic Diketoenamine Crosslinking Unlocks Vinyl Polymer Vitrimers From β ‐Triketone Chemistry

Covalent adaptable networks (CANs) offer a compelling strategy to unite the mechanical robustness of thermosets with the reprocessability of thermoplastics, yet achieving simultaneous durability, processability, and true recyclability remains challenging. We introduce diketoenamine (DKE) vitrimers derived from β-triketone methacrylate monomers and demonstrate how rational monomer design dictates network processability, viscoelasticity, and recyclability. By systematically varying the spacer length between the β-triketone (TK) moiety and the polymer backbone, we identify a key structure–property relationship that dictates vitrimer behavior. Networks bearing TK pendants minimally displaced from the backbone suppress creep but exhibit limited stress relaxation, whereas extended spacers yield lower glass transition temperatures, higher effective crosslink densities, and efficient stress dissipation, enabling optical transparency and reprocessability. Extending this platform to ultra-high molecular-weight prepolymers introduces physical entanglements as secondary crosslinks, further enhancing dimensional stability without compromising processability. Both mechanical and chemical recycling validate the closed-loop circularity of these materials. Furthermore, these results establish TK methacrylates as a versatile platform for designing high-performance vitrimers that integrate durability, reprocessability, and true closed-loop recyclability.

closed-loop recyclability↗

Screening and Development of Sacrificial Cathode Additives for Lithium‐Ion Batteries

Abstract This work presents a computational screening approach to identify Li‐rich transition‐metal oxide sacrificial cathode additives and provides experimental validation of antifluorite‐structured Li 6 MnO 4 as a potential candidate. Initial attempts to synthesize this compound result in low purity (≤40% by weight) owing to close thermodynamic competition with Li 2 O and MnO at low temperature. However, it is shown that a much higher purity of 85% by weight can be achieved by combining Li excess with rapid cooling from high temperature, which effectively stabilizes the Li 6 MnO 4 phase. The synthesized product delivers a high irreversible Li release capacity that exceeds 700 mAh g −1 by utilizing combined Mn oxidation (Mn 2+/3+ and Mn 3+/4+ ) and O oxidation. These results demonstrate that Li 6 MnO 4 may therefore be useful as a potential sacrificial cathode additive in Li‐ion batteries and motivate further investigation of other structurally‐related compounds. While attempts were made to synthesize two additional compounds among computationally screened candidates, it was not successful to experimentally realize the two candidates. The difficulty of experimental realization of the newly predicted materials remains a challenge and it is suggested that more efforts need to be devoted to developing computational techniques to precisely predict synthesizability and propose potential synthetic routes of the predicted materials.

Kim, Haegyeom↗

Ab initio thermodynamics of Ni and Co incorporation in Mg hydroxide, carbonate, and hydroxycarbonate minerals

Ni and Co are critical elements needed for modern technologies, and a better understanding of the ability of Mg-based minerals to incorporate these elements would benefit strategy development for Ni and Co recovery from mafic and ultramafic deposits. Here, in this work, we performed density functional theory (DFT) calculations of Ni and Co incorporation in six potential products of the carbonation of mafic and ultramafic silicates: brucite (Mg(OH) 2 ), magnesite (MgCO 3 ), nesquehonite (MgCO 3 ⸱3H 2 O), lansfordite (MgCO 3 ⸱5H 2 O), artinite (Mg 2 CO 3 (OH) 2 ⸱3H 2 O), and hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ⸱4H 2 O). The DFT results were used in an ab initio thermodynamics framework to explore the pH 2 O–pCO 2 conditions at which the Mg-based minerals were predicted to be thermodynamically stable and to quantify the Gibbs free energy of Ni and Co substitution at Mg sites. Among the six Mg-based minerals, brucite and magnesite were predicted to have the lowest Ni and Co substitution free energy. An analysis of the effect of temperature indicated that, at low temperature (<100 K), brucite more readily accommodated Ni and Co, while, at higher temperature (>335 K), magnesite more favorably incorporated Ni and Co. Between 100 K and 335 K, Ni was predicted to preferentially substitute for Mg in brucite and Co for Mg in magnesite, thus leading to a driving force for separating Ni and Co in conditions where brucite and magnesite both form. Insights gained in this work could therefore help select experimental conditions that either promote or inhibit incorporation of these critical elements into Mg-based mineral phases.

Critical elements↗