Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamics”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Tightening the thermodynamic uncertainty relations with null-entropy events: What we learn when nothing happens

Fluctuation theorems establish that thermodynamic processes at the microscale can occasionally result in negative entropy production. At the microscale, another distinct possibility becomes more likely: processes where no entropy is produced overall. In this work, we explore the constraints imposed by such null-entropy events on the fluctuations of thermodynamic currents. By incorporating the probability of null-entropy events, we obtain tighter bounds on finite-time thermodynamic uncertainty relations derived from fluctuation theorems. We validate this framework using an example of a qudit SWAP engine.

FOS: Physical sciences↗
DOE OSTI2022

Development Of Thermodynamic and Kinetic Simulation Tools and Testing Procedures for Enhanced Durability of Concrete Containing Industrial By-Products

This project developed screening tools that enable evaluation of alternative cementitious binders that create concretes to significantly reduce energy and emissions while remaining cost competitive on both initial and long-term costs. The team began with the viewpoint that acceptance of new cementitious binder products has a substantially greater chance of successful implementation when capital investment is not excessive and the end product has customers that have experience and resources to use this. As such, a binder system that is based on portland cement with blended using industrial by-products (alternative cementitious materials) has the potential for dramatic and meaningful impact. The team has focused on developing implementable solutions in specifications and current practice. This however requires three main factors: 1) ability to screen byproducts and alternative materials for success, 2) ability to ‘treat’ materials chemically to enhance kinetics, and 3) ability to provide predictions of performance of both binders and concrete from first principles. The project developed/refined a state of the art and scientifically based screening test for SCM called the pozzolanic reactivity test. The team developed kinetic models to simulate these materials as well as experimental approaches to alter selected reactions. Simulation tools were developed that enable the performance of concrete to be predicted based on the chemistry and reactivity of the cement and alternative SCM. Specifically, this project: • enhanced the kinetic reactivity models for use in multi-scale computational programs that use thermodynamics to predict reaction products. 2 • developed scaling models to extend thermodynamic modeling to link these models with pore structure. This enables strength, transport property, and coupled transport prediction. • developed tools to predict performance of cementitious materials using the pozzolanic reactivity test and chemical composition. The predicted properties are consistent with AASHTO R101 and the CEB-FIP model code and can be measured using associated test procedures. • demonstrated mechanical and fracture based modeling tools that thermodynamic predictions and inputs to predict concrete service life. These results have been used to demonstrate the value of enabling specifications to include ASTM C 595 cement as well ASTM C150 cement. In addition, these products are being used to expedite the evaluation of alternative SCM to aide in determining which materials have potential value and what ‘compositions’ of blended cements merit further investment.

42 ENGINEERING↗
DOE OSTI2022

Development of a Thermodynamic Database for Corrosion in Chloride MSRs

Control of the salt redox condition is the salient method for corrosion inhibition in molten salts. For chloride salt-fueled molten salt reactors of primary concern is the U(IV)/U(III) ratio, influenced directly by the Cl2 activity in the salt. Although thermodynamic assessments of MSR-relevant salts abound, none have included consideration of UCl4 which is essential to correctly describe the salt redox condition. This work reports the assessment of the NaCl-KCl-MgCl2 base salt with UCl3-UCl4 and corrosion chlorides involving Cr, Fe, and Ni as minor components and the resultant thermodynamic models. The effort utilized available thermodynamic measurements and phase equilibria complemented by differential scanning calorimetry determinations of phase equilibria. Example calculations explore the temperature and composition space to understand its influence on corrosion. The foregoing models are incorporated into an update of the Molten Salt Thermal Properties Database – Thermochemical which is publicly available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Defect Thermodynamics and Transport Properties of Perovskite and Fluorite Materials for Solid-Oxide and Proton Conducting Oxide Cells Evaluated Based on Density Functional Theory Modeling

Density functional theory based defect thermodynamic modeling was performed to determine the effect of humidity and H2/O2 gas pressure on various defect chemistry and transport properties of perovskite and fluorite oxides for solid-oxide and proton-conducting-oxide cell applications, with inclusion of the electronic-conducting oxides (as electrodes) and insulating oxides (as electrolytes). Automatic defect generation workflow and first-principles charged defect analysis were implemented on NETL Joule supercomputer for modeling defect equilibria and transport properties of insulating oxides as electrolytes in SOCs and proton-conducting ceramic cells. A GNU Octave defect model subroutines were developed to facilitate defect modeling of electronic conducting oxides in a wide range of operating conditions guided by modeling and experiments. The model includes the hydride defect formation reaction under reducing conditions and allows to incorporate nonstoichiometry effects on the defect thermodynamic parameters. The developed model serves as a platform to facilitate fundamental understanding of the defect thermodynamics in SOC oxide materials and can be used as a novel tool in computational materials screening for SOC and other energy applications.

Lee, Yueh-Lin↗
DOE OSTI2022

Thermodynamic Model and Initial Experimental Investigation of Air Dehumidification through Electrically Charged Vapor Capturing Electrostatic Droplets

The removal of water vapor from the air to reduce relative humidity is a well known indoor environmental comfort requirement. Common dehumidification approaches require a substantial amount of energy and usually involve the cooling of atmospheric humid air below its dew point or the use of absorbent/adsorbent materials to extract water vapor out of the air. More recently, researchers investigated the effect of electrostatic forces for enhancing water vapor condensation with the goal to reduce the energy consumption associated with dehumidification. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. In this thesis, a broad theoretical investigation to study the electrostatically enhanced condensation processes was carried out. These processes consist of the use of highly charged particles, preferably highly charged water droplets to attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. The electrical charge promotes the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation on curved surfaces and, consequently, the equilibrium between evaporation and condensation is shifted towards condensation. This investigation resulted in the development of a thermodynamics model, which was able to predict an effective size range of the charged droplets for optimal dehumidification, under ideal conditions. The range resulted in about 2 to 4 m in diameter, while the electrical charge was kept to the maximum limit predicted by the Rayleigh model. A sensitivity analysis on the variation of the size and the charge of the charged droplets was also performed. In terms of dehumidification rates, when six electrospray heads, i.e. the maximum number of heads tested in the experimental campaign, were considered, the model predicted very limited rates, way lower than the target of this research, i.e. a dehumidification rate in terms of relative humidity of 5 % with an air flow rate of 5 cfm. While ways to increase this rate existed, their implementation proved difficult for this initial investigation. According to the requirements predicted by the model, the use of electrosprays appeared the most suitable solution for the production of small but highly charged droplets. The electrospray features and operational modes were studied in detail and an electrospray assembly was designed to be tested in this initial experimental investigation. For this initial investigation, an experimental system was designed and built to validate the thermodynamics model and with the goal to achieve a dehumidification rate of 5 % in terms of relative humidity with an air flow rate of 5 cfm. The experimental system consisted of a wind tunnel test apparatus, dew point sensors, thermocouples for temperature measurements, pressure sensors and other devices. The experimental system was equipped with a computerized data acquisition and storage system. The system was suitable to evaluate the air water content differential before and after the test section where the electrospray heads were installed. The air water content differential was evaluated in terms of the ∆ω/ω_1 ratio, to eliminate the dependence on the dry bulb temperature along the test apparatus. This work presented initial experimental data for electrostatically enhanced dehumidification processes. Working parameters, such as air and water flow rates, high voltage potentials and polarity, deionized water types with different electrical conductivity, were varied to find the best combination for an improved dehumidification. With the set up used and the conditions considered, it was never possible to achieve a 5 % dehumidification rate with 5 cfm. In general, the dehumidification rate was always limited and lower than 1 % for all the air flow rates considered and mainly within the uncertainty of the dew point sensors. These results confirmed the already limited predictions of the thermodynamics model, for ideal conditions and the interesting and promising results obtained in some experimental investigations available in the open literature were not achievable in this thesis. The more probable reasons were the higher air flow rates considered and the complete absence of the use of cooling power, which was mainly implemented to facilitate the vapor condensation.

Morcelli, Stefano↗
DOE OSTI2019

A Thermodynamic Model Toward the Comprehension of Ferrous Burden Softening and Melting Using FactSage Macro-Processing

The softening and melting properties of iron-bearing materials play a decisive role in the formation of the cohesive zone, which greatly affects blast furnace gas flow distribution and heat-transfer efficiency. To improve the understanding regarding the evolution of condensed phases during reduction, softening, and melting, a thermodynamic model has been developed using FactSage™ thermodynamic software and macro-processing. The model was constructed using a series of equilibrium stages and splitters to determine flow directions of streams and to consider kinetic inhibitions. For all iron-bearing materials studied, the methodology proposed for modeling was capable of obtaining reduction degrees in very good agreement to its experimental data. The evolution of solid phases was qualitatively comparable to the available literature, with Fayalite, Kirschsteinite, Melilite, and FeO as the main solid phases in equilibrium before slag formation. The comparison between the profiles of the calculated slag mass fraction and the experimental pressure drop showed a close relation between these properties. Moreover, the level of heterogeneity of each raw material may play a significant role in the interpretation of its results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chemical degradation of fly ash blended concrete with the seasonal variation of rainwater in a radioactive waste repository: A thermodynamic modeling approach

Highlights: • The thermodynamic model for the cementitious system was successfully developed. • The chemical degradation of concrete was greatly affected by rainwater types. • Inorganic carbon species in rainwater reduced the chemical resistance of concrete. • The durability of concrete to rainwater decreased with increasing fly ash content. This study presents a long-term degradation behavior of fly ash blended concrete in a vault type low- and intermediate-level radioactive waste (LILW) repository by thermodynamic equilibrium calculations using PHREEQC combined with CEMDATA18 database. Since rainwater plays a predominant role as leachate for the concrete waste container in the repository, the effect of chemical properties of seasonal rainwater on the chemical degradation of concrete was analyzed. In addition, the impact of the blending with fly ash on the chemical degradation was evaluated through a gradual replacement of ordinary Portland cement (OPC) by fly ash. Regardless of the replacement level, autumn rainwater showed the greatest influence on the concrete degradation owing to the presence of inorganic carbon species, which induce the carbonation. The resistance of concrete to chemical degradation by rainwater was persistently reduced due to the decreasing formation of portlandite and calcium silicate hydrate (C-S-H) with increasing replacement level.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Temperature screening and cross-field impurity accumulation from a thermodynamic perspective

In a variety of different systems, high-Z ion species show a marked tendency to accumulate in regions of high plasma density. It has previously been suggested that the apparent universality of this behavior could be explained thermodynamically, in terms of the maximum-entropy state attainable when the system must obey an ambipolarity condition. Furthermore, the previous treatment did not allow for the possibility of temperature gradients. Here, tools from non-equilibrium (Onsager) thermodynamics are used to show that ambipolarity continues to play a key role in producing this behavior in the presence of temperature gradients, and to recover well-known temperature screening effects that appear in these cases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

PET-FBA: A lightweight enzyme allocation and thermodynamics-constrained flux analysis approach to explore Escherichia coli metabolic adaptation to intracellular acidification

Escherichia coli employs diverse strategies to adapt to acidic environments that disrupt enzyme activity and the thermodynamic feasibility of essential reactions. To understand the impact of pH stress on cell metabolism, we present the PET-FBA (pH-, Enzyme protein allocation-, and Thermodynamics-constrained Flux Balance Analysis) framework. PET-FBA extends genome-scale modeling by integrating enzyme protein costs and reaction Gibbs free energy changes. Additionally, by incorporating pH-dependent enzyme kinetics in response to intracellular acidification, this framework enables the simulation of E. coli's metabolic adjustments across varying external pH levels. The model's accuracy is validated by comparing in silico growth simulations with experimental measurements under both anaerobic and aerobic conditions, as well as in silico gene knockouts of essential genes. By explicitly incorporating pH effects, our model accurately replicates the metabolic shift towards lactate production as the primary fermentation product at low pH in anaerobic conditions. This shift is only predicted when enzyme kinetics are dynamically adjusted as a function of pH. Further analysis revealed that this shift can be attributed to the reduced protein efficiency of the acetyl-CoA branch compared to lactate dehydrogenase under acidic stress, which then becomes crucial for maintaining NAD regeneration and cell growth at low pH. Furthermore, we identified strategies for enhancing cell growth under acidic anaerobic conditions by improving the enzyme activity of lactate dehydrogenase and pyruvate formate lyase, which increases NAD production efficiency and reduces enzyme protein allocation costs. Designed as a lightweight yet versatile framework, PET-FBA enables efficient genome-scale metabolic analysis. Using E. coli as a model system, our framework provides a systematic approach to understanding metabolic responses to environmental stress, pinpointing key metabolic bottlenecks, and identifying potential targets for strain optimization.

42 ENGINEERING↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Probing the Effect of Electrode Thermodynamics on Reaction Heterogeneity in Thick Battery Electrodes

Thick electrodes present a viable strategy for enhancing energy density and reducing manufacturing costs of lithium-ion batteries. However, reaction heterogeneity during cycling compromises their rate capability and cycle life. While this nonuniformity is commonly attributed to sluggish charge transport, it is demonstrated here that the thermodynamic properties of the electrode material play an equally critical role. Through combined X-ray fluorescence microscopy and absorption near-edge structure spectroscopy, reaction distributions in LiFePO 4 (LFP) and LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC) thick electrodes with matched porosity and tortuosity are compared. LFP electrodes develop pronounced depth-oriented state-of-charge (SOC) gradients that worsen with increasing discharge rates, whereas NMC maintains much more uniform SOC distributions under such conditions. This difference originates from their distinct SOC dependence of equilibrium potentials and is quantifiable through a dimensionless “reaction uniformity” number. Intriguingly, LFP thick electrodes also exhibit lateral SOC variations that strengthen during slow discharge. In conclusion, the enhanced reaction uniformity in NMC correlates with better active material utilization and slower capacity fade than LFP, highlighting electrode thermodynamics as a key design consideration for thick electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Roles of Solvent in the Catalytic Hydrogen Release from Liquid Organic Hydrogen Carriers: Chemical, Thermodynamical and Technological Aspects

A Liquid Organic Hydrogen Carrier (LOHC) enables the storage and transport of hydrogen at ambient pressures and temperatures in a safe and convenient form using current infrastructure. However, it is challenging to directly compare reactivity and selectivity for hydrogen release, especially when comparing the catalytic efficiencies of neat LOHCs to highly diluted LOHCs in different solvents, reaction conditions, and catalysts. This work evaluates the role of solvents in catalysis and quantifies the energy efficiency of the overall process. The presence of solvent dilutes the volumetric density of available hydrogen, but may be necessary to achieve optimal catalysts stability, reactivity, and product selectivity. With respect to the reaction conditions as determined by thermodynamics, solvents with higher vapor pressures than that of the carrier can cause the erroneous impression of a more favorable reaction equilibrium. Concerning energy efficiency, solvents can result in increased energy demand for hydrogen release as the inert solvent must be heated to reaction temperatures required for release of H 2 from the LOHC. Further, this work recommends that investigations of catalyst reactivity should be carried out at different ratios of solvent to LOHC to understand how the reactivity changes and what the implications are for maximizing energy density and catalyst stability and reactivity. Investigations should also consider how these implications will affect the technical needs of applications intended for the LOHC system. Based on the results of this study, it is advised to focus research activities on LOHC systems with a gravimetric solvent content below about 50% as the thermodynamic disadvantages become very pronounced beyond this threshold.

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

Ab Initio Modeling on The Thermodynamic and Temperature-Dependent Elastic Properties of Subsystems of The FCC FeNiCoCr Medium Entropy Alloys (MEAs)

The stability, phonon spectra, thermodynamic, and temperature-dependent elasticity of subsystems of the FCC FeNiCoCr MEAs are systematically studied by the ab initio approach. Especially, the quasi-harmonic approximation (QHA) and the innovative Zentropy theory were utilized to predict the thermodynamic properties and elastic properties of FeNi, NiCo, FeNiCo, and FeNiCoCr MEAs with the consideration of magnetic transition. With the ensemble of the partition function based on the multiplicity of each magnetic configuration, the current work successfully reproduced the Curie temperature and the Schottky anomaly of heat capacity of these four MEAs purely based on the ab initio input. Meanwhile, the elastic properties of these alloys at finite temperatures are also successfully predicted with the consideration of magnetic transition. The overall results agree well with the available experimental data and CALPHAD prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
DOE OSTI2025

Distinctive features of structural evolution and thermodynamic response in wide-bandgap semiconductors driven by intense electronic excitation

Radiation-tolerant material selection requires balancing lattice rigidity, defect dynamics, and electronic stability, as shown by covalent SiC outperforming ionic Ga 2 O 3 and GaN under extreme environments. Responding to intense electronic excitation, irradiation-driven phase segregation (β → δ/κ in Ga 2 O 3 ) and core–shell track (disordered structure in GaN), accompanied by elemental redistribution, contrastingly, exceptional radiation tolerance manifested by comparatively minimal lattice distortion (0.17 % strain variation) was demonstrated in SiC. These differential responses are primarily attributed to two fundamental mechanisms: (i) thermodynamic driving forces governing defect migration and phase separation, and (ii) the synergistic effects of robust covalent bonding composition coupled with efficient defect recombination processes. Here, the stronger electron–phonon (e-ph) coupling in Ga 2 O 3 (4.34 × 1018 W m −3 K −1 ) and GaN (3.55 × 10 18 W m −3 K −1 ) enhances lattice energy deposition, triggering thermal spikes (ΔT ≫ T m ) and structural transition behaviors, whereas weaker e-ph coupling in SiC (3.69 × 10 18 W m −3 K −1 ), relatively high thermodynamic parameters and efficient energy dissipation suppress thermal spikes to maintaining lattice integrity. The photoresponse degradation driven by enhanced radiative recombination is dominant in N-doped SiC, while V-doped systems achieve defect-mediated photoconduction optimization characterized by abrupt current transitions, matching fluorescence yield evolutions, and directly connecting defect engineering to optoelectronic performance.

Intense electronic excitation↗
DOE OSTI2022

Thermodynamic modeling of the Pd–Zn system with uncertainty quantification and its implication to tailor catalysts

Pd–Zn intermetallic catalysts show encouraging combinations of activity and selectivity on well-defined active site ensembles. Thermodynamic description of the Pd–Zn system, delineating phase boundaries and enumerating site occupancies within intermediate alloy phases, is essential to determining the ensembles of Pd–Zn atoms as a function of composition and temperature. Combining the present extensive first-principles calculations based on density functional theory (DFT) and available experimental data, the Pd–Zn system was remodeled using the CALculation of PHAse Diagrams (CALPHAD) approach. High throughput modeling tools with uncertainty quantification, i.e., ESPEI and PyCalphad, were incorporated in the phase analysis. Here, the site occupancies across the γ phase composition region were given special attention. A four-sublattice model was used for the γ phase owing to its four Wyckoff positions, i.e., the outer tetrahedral (OT) site 8c, the inner tetrahedral (IT) site 8c, the octahedral (OH) site 12e, and the cuboctahedral (CO) site 24g. The site fractions of Pd and Zn calculated from the present thermodynamic model show the occupancy preference of Pd in the OT and OH sublattices in agreement with experimental observations. The force constants obtained from DFT-based phonon calculations further supports the tendency of Pd occupying the OH sublattice compared with the IT and CO sublattice. The catalytic ensembles changing from Pd monomers (Pd 1 ) to trimers (Pd 3 ) on the surface of γ phase are attributed to the increase of Pd occupancy in the OH sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2026

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↗