SEARCH · Search NASA

Results for “constituency”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A New Method for Polishing Laminated Dissimilar Metals

Here, this paper presents a novel method for preparing metallographic samples of laminated dissimilar metals, specifically a layered structure of zirconium (Zr), vanadium (V), and tungsten (W), which pose a challenge due to their differing crystal structures and hardness levels. The work was motivated by the need to develop a cladding for tungsten targets in spallation neutron sources, where a zirconium face is applied over an interlayer of either vanadium or titanium to prevent the formation of brittle intermetallic compounds such as ZrW 2 . Traditional polishing techniques prove inadequate because each constituent metal requires a different polishing approach based on its unique material properties. The methodology outlined includes a sequential process involving conductive epoxy mounting, successive grinding and polishing with specific abrasives, and a final chemical-mechanical polish using a colloidal silica-based solution. Optical microscopy confirmed that this multi-step process successfully produced scratch-free surfaces and minimized surface relief across the different metal interfaces. Although developed for refractory metals, this approach shows promise for other layered structures comprised of transition metals, offering a meticulous yet effective technique for high-quality metallographic sample preparation.

Newberry, Daniel [Oak Ridge National Laboratory (O↗

DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗

DOE OSTI2026

A modified scheil approach for nucleation-dependent solidification pathways

As-solidified microstructures of near-eutectic alloys often contain multiple primary phases that are not expected from equilibrium phase diagrams. Such microstructures are caused by cooling-rate-dependent solidification pathways, a factor not captured by the Scheil–Gulliver model or variations thereof. Here, we present a model and algorithm that incorporate the critical nucleation undercooling for each solid phase into the Scheil–Gulliver model. We hypothesize that the non-equilibrium microstructure formation is primarily governed by a nucleation-competition mechanism. This mechanism accounts for both stable/metastable phase selection and primary-phase formation within eutectic regions driven by asymmetric nucleation barriers. The model is validated against a hypereutectic Al-Fe alloy, where it successfully reproduces the observed microstructural constituents, revealing the key dependencies of solidification microstructure on nucleation kinetics. Applicability to multicomponent systems is demonstrated through a hypereutectic Al–Fe–Si ternary alloy, where the model successfully predicts divorced eutectic microstructures and the associated oscillatory solidification pathways along univariant lines. As a result, the proposed framework establishes a nucleation-dependent computational approach for interpreting and predicting solidification microstructures.

Alloy design↗

DOE OSTI2026

Pulsed laser synthesis of mesoporous metal chalcogenide thin films

Mesoporous films of the metal chalcogenide B-FeSe were grown on MgO substrates by KrF pulsed laser deposition (PLD) in an argon background. At 100 mTorr, gated intensified charge-coupled device imaging and ion probe measurements showed that the plasma plume responsible for crystal growth initially comprised three components, with distinct expansion velocities. Plume interactions with the substrate heater and ablation target gave rise to complex dynamics, including collisions between the charged leading edge—rebounding between the substrate and the target—and slower-moving species in the plume interior. Film growth was dominated by species with kinetic energies ≤0.5 eV/atom. X-ray reflectivity revealed that films grown in this environment—with a substrate temperature of 350 ° C, a laser fluence of 1.0 J cm −2 , and a 7.5 mm 2 spot area—formed a porous framework with 15% porosity. Atomic force microscopy showed surface features that suggest pore sizes below 100 nm. X-ray diffraction indicated that the porous films were epitaxial with respect to the substrate and likely grew by oriented-attachment of gas-phase molecular clusters or very small nanoparticles, in contrast to the conventional epitaxy of vacuum films from atomic constituents. The in-plane orientation of the mesoporous films was B-FeSe [100]||[110] MgO, attributed to the soft landing of pre-formed crystallites on the MgO substrates, where protruding Se rows of B-FeSe aligned with corrugations of the MgO surface. In conclusion, this work implies that growth of candidate electrocatalyst materials by PLD in inert gas background may allow mesoporous frameworks with a single crystallographic orientation that expose specific crystal facets for electrochemical reactions and active site engineering.

cluster epitaxy↗

DOE OSTI2025

A worldwide aerosol phenomenology: Elemental and organic carbon in PM 2.5 and PM 10

Elemental carbon (EC), organic carbon (OC), and particulate matter (PM) concentrations in the inhalable (PM 10 ) and fine (PM 2.5 ) size fractions are measured worldwide, albeit with different analytical methods. These measurements from many researchers were collected and analyzed for Africa, America, Asia, and Europe for 2012–2019. EC/PM, OC/PM, and OC/EC ratios were examined based on region, site type, and season to infer potential sources and impacts. These analyses demonstrate that carbonaceous materials are important PM constituents throughout the world. Mean EC/PM ratios were lowest in PM 10 in Sahelian Africa and Europe (∼0.01), highest (>0.07) in PM 2.5 at urban sites in North America, South America, and Japan. Mean OC/PM ratios were lowest in PM 10 in the Sahel (∼0.06) and in PM 2.5 in China and Thailand (0.10), and highest in central and eastern Europe (∼0.3) and North America (∼0.4). OC/EC ratios were elevated in western and northern Europe, and at regional background sites in North America. EC/PM increased with PM 10 in Thailand, while OC/PM increased with higher PM mass in Thailand, India, and North America, highlighting the specific contribution of carbonaceous aerosols to PM pollution in these regions. At European and North American background sites, OC/EC ratios increased with PM mass. Higher OC/EC ratios in dry periods indicate influence of wildfires, prescribed burns, and secondary aerosol formation. Elevated wintertime EC/PM ratios coincide with residential heating in temperate climate zones.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗

DOE OSTI2025

Physical model and experimental validation of a high temperature proton exchange membrane electrochemical hydrogen pump cell for efficient single-stage extraction of low concentration hydrogen gas

There is interest in valorization of existing natural gas infrastructure to facilitate the co-transportation of hydrogen via blending of hydrogen gas initially at limited concentrations of 1–20 vol% H2 and to subsequently extract hydrogen at fuel cell quality standards (SAE J2719/ISO14687-2). High temperature proton exchange membrane electrochemical hydrogen pump (HT-PEM EHP) based on phosphoric acid doped polybenzimidazole (PA-PBI) exhibits good performance at elevated temperatures (>120 °C), which provides desirable tolerance to non-methane natural gas constituents that are problematic for lower temperature based EHP. To better understand the suitability of the HT-PEM EHP for such gas separation processes, a two-dimensional model of EHP based on PA-PBI was developed. The model is validated for several relevant operating conditions and across cells with differing amounts of phosphoric acid content in the electrodes. Operando micro x-ray computed tomography (CT) imaging of an HT-PEM EHP was used to further validate physical parameters and assumptions of the model. The impacts of pressure, relative humidity of the anode feed, and concentration of feed gases on separation performance are investigated. This study shows that a specific energy of separation of 5.1 kWh/kg H2 at a hydrogen recovery factor (HRF) of 50 vol% can be achieved in a single stage with the EHP, producing fuel cell quality hydrogen purity of 99.99 vol% H2 from a 2 vol% H2/CH4 feed blend, while pressurizing the product H2 at a pressure ratio of 1.3 relative to feed pressure.

Stansberry, John M↗

DOE OSTI2025

Component Identification of Solid Biomass Fuels Using Reflected Light Microscopy: Interlaboratory Study 2

As nations transition toward sustainable energy systems, biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

09 BIOMASS FUELS↗

DOE OSTI2025

Protecting honey bees through microbiome engineering

Honey bees are indispensable insects. Their pollination services support modern agriculture and natural ecosystems. Managed honey bee colonies face increasing threats to their survival, ranging from environmental stressors that include agrochemicals to infestations of arthropod pests and infections with microbial pathogens. Like humans, honey bees have a native gut microbiome that supports their health. However, the bee gut microbiome has a simpler composition than the gut microbiome of mammals, and its main constituent bacterial species can be easily cultured outside of the host. Furthermore, this experimental tractability and the need for new methods for protecting hive health have made honey bees a testbed for synthetic microbiomes augmented with probiotic bacteria and engineered DNA. Here, we discuss the natural benefits of bee gut bacteria, recent progress in genetically modifying these bacteria, and how symbiont-mediated RNA interference and other microbiome engineering approaches can boost bee immunity and suppress bee pathogens and parasites. Finally, we discuss how emerging methods for microbiome engineering and biocontainment could be applied to honey bees and used to address challenges in translating these proof-of-principle achievements into safe and effective technologies for field applications at scale.

Biological and medical sciences↗

DOE OSTI2025

Screening green solvents for multilayer plastic film recycling processes

Multilayer (ML) plastic films are essential packaging materials that help protect products from diverse external factors; however, only 5% of all ML films are recycled in the United States. Solvent-based technologies are a promising alternative for recycling ML films because they enable recovery of constituent polymer resins. For example, the Solvent Targeted Recovery and Precipitation (STRAPTM) process sequentially dissolves and separates polymer components using a series of targeted solvent washes. A crucial design aspect of this process is the impact of selected solvents on human health and on the environment. Here, this work introduces a computational framework that integrates molecular modeling, process modeling, techno-economic analysis (TEA), and life-cycle analysis (LCA) to quickly screen green solvents for solvent-based ML recycling processes. Initial screening for solvents based on selectivity is performed by estimating temperature-dependent solubilities using molecular-scale models. Subsequent screening uses basic estimates of energy use and octanol-water partition coefficients (logP) as key measures of health, safety, and environmental hazards. Detailed process modeling, TEA, and LCA are used on a reduced set of promising solvents identified in early screening steps to more accurately determine how solvent selection and associated operating conditions impact overall economics and environmental impacts. The framework is used for the identification of green solvents (from a database of 1,000 solvents) that separate an industrial ML film composed of polyethylene (PE), ethylene vinyl alcohol (EVOH), and polyethylene terephthalate (PET). Our analysis shows the effectiveness of the framework and reveals fundamental trade-offs between solvent greenness, solubility, and economics. Our work emphasizes the importance of taking a holistic systems view during solvent design and aims to inform the development of new processes for ML film recycling and the identification of new ML films that are easier to recycle.

economics↗

DOE OSTI2025

A fast computational framework for the design of solvent-based plastic recycling processes

Multicomponent plastics cannot be processed using mechanical recycling technologies, hindering efforts to deal with plastic waste. Multicomponent plastics include multilayer plastic films, which are widely used for food and healthcare packaging. Multilayer films combine several layers (potentially dozens) of different polymers to protect products from external factors (e.g., oxygen, water, temperature, shock, and light). Solvent-based separation processes have emerged as a promising alternative to recycle these complex materials. For instance, the Solvent-Targeted Recovery and Precipitation (STRAP TM ) process uses sequential solvent washes to selectively dissolve and separate constituent polymers from multicomponent plastic waste, including films. STRAP TM process design (separation sequence, type of solvents, and operating conditions) changes significantly depending on the design of the multilayer plastic film (e.g., number, types, and proportions of polymers). The ability to quickly quantify the economic and environmental benefits of diverse STRAP TM process designs is essential to accelerate the development of sustainable recycling processes and more recyclable multilayer film products. In this work, we present a fast computational framework that integrates molecular-scale models, process modeling, and techno-economic and life cycle analysis to quickly evaluate STRAP TM designs. The computational framework is general and can be used to study the processing of complex multilayer plastic waste streams that contain many layers. Furthermore, we highlight the different uses of the framework via targeted case studies.

Computational framework↗

DOE OSTI2025

Machine learning-assisted identification of potential sources of bias in measurements of prompt-fission neutron spectra

Unrecognized sources of uncertainty (USU) can bias the reported mean and/or covariance of experimental nuclear data. These biases, in turn, can propagate through evaluated nuclear data to application simulations or may poorly inform nuclear theory that is fitted to the experimental data. Such unknown sources of bias must be tied to the inherent physical constituents of the measurements such as the characteristics of a detector response or a background reduction technique. Here, in this article, a sparse Bayesian learning model is used to support experts in their efforts to identify and characterize USU in experimental prompt fission neutron spectra (PFNS) for spontaneous fissioning of 252 Cf by linking observed biases to features of the measurement system. Three different bias components were found. The first acts as a verification case for the algorithm as it identifies a bias coming from a well-known source related to the use of 6 Li in the neutron detection system. The second two cases demonstrate how this method can benefit the evaluation of experimental nuclear data by identifying, quantifying, and relating unknown biases to potential causes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Influence of strain-rate on the response of elastomeric architected materials

Architected materials have shown substantial promise in impact mitigation and protective applications, and there has accordingly been great interest in better characterizing their response at elevated strain rates due to impact. There remains ambiguity regarding the contribution of inertial and material responses to strain rate sensitivity, and, in particular, when these effects begin to gain dominance in the impact response of an architected material. The response of soft polymer architected materials as a function of strain rate, in particular, has been little investigated. We characterize the experimental impact response of four soft polymer architected lattice geometries across varying strain rates in the intermediate strain rate regime (∼10 3 s −1 ) using split-Hopkinson pressure bar loading and high speed video characterization of the resulting deformation fields. In conclusion, our results highlight the interplay of influence between constituent material, lattice geometry, length scale, and strain rate in determining the onset of significant inertia effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Hybrid power plants: An effective way of decreasing loss-of-load expectation

Diversifying variable renewable resources by combining wind, solar photovoltaic, and battery assets in a hybrid power plant can increase renewable energy usage efficiency and improve system flexibility, particularly in distributed energy systems. However, the resilience impact of these systems, particularly outage mitigation, can be difficult to quantify due to uncertainty in resource, energy demand, and outage occurrence. Here, this study outlines a framework to quantify the incremental benefit of hybrid power plant assets for reducing loss-of-load expectation during random outage events. Hybrid power plant performance during outages (considering varying duration and severity) is simulated using a Monte Carlo methodology to reflect uncertainty associated with renewable resource, load demand, and outage timing. Results demonstrate the additional incremental value from increasingly hybrid designs, in which relative capacities of wind, solar photovoltaic, and storage assets contribute to lower loss-of-load expectation than the constituent technologies would alone. The value of added wind or solar capacity increases as the plant composition approaches an equal split. The value of added battery capacity depends on the outage duration and severity, but the first 50 MWh of added storage capacity is the most valuable for reducing the loss-of-load expectation for all plant designs.

14 SOLAR ENERGY↗

DOE OSTI2026

A collective trigger for widespread planetesimal formation revealed by accretion ages

The formation of planetesimals was an integral part of the cascading series of processes that built the terrestrial planets. To illuminate planetesimal formation, here we develop a refined thermal evolution model to calculate the formation ages of meteorite parent planetesimals. This model includes chemical reactions and phase changes during heating, as well as natural variations in the proportions of the constituent phases of these planetesimals. We find that the parent bodies of non-carbonaceous (NC) and carbonaceous (CC) iron meteorites start forming at very similar times ( ∼ 0.95 Myr after calcium-aluminium-rich inclusion [CAI] formation) and occupy overlapping time windows. NC and CC chondrite parent bodies formed later during non-overlapping periods. We combine these ages with proportions of isotopic end-members we recover from mixing models to construct records of motion throughout the protoplanetary disk. These records argue that NC and CC material traversed the barrier in the disk after ∼ 0.95 Myr after CAI formation. The onset of this motion coincided with planetesimal formation, indicating that the phenomenon that drove motion also triggered planetesimal formation. We argue that this feature also served as the semi-permeable barrier in the disk. Although its identity is uncertain, the effects this phenomenon had on the timing of planetesimal formation and motion through the disk can now serve as constraints on models of disk evolution. Models that reproduce these effects would elucidate the nature and implications of this phenomenon, which is key to unlocking a holistic model of terrestrial planet building.

Astronomy and AstroPhysics↗

DOE OSTI2024

Microwave-Assisted dry reforming of toluene as a model tar compound using low-cost iron catalyst for syngas clean-up

Gasification of waste feedstock such as biomass, waste plastics suffers from high tar yields during hydrogen-rich syngas production. The presence of tars result in lower quality syngas, lower syngas yields, reactor blockage, reactor down time, and costly maintenance. Therefore, removal of tar from syngas during gasification is essential. Catalytic reforming of tars via in-situ syngas cleanup is an effective way of mitigating tars. This work explores the possibility of microwave-assisted catalytic dry reforming of tars for syngas production. Further, due to its chemical complexity, toluene, which is one of the main tar constituents, could be used as a model tar compound. Toluene dry reforming was studied using the Fe/Al 2 O 3 catalyst under CO 2 under the temperature range of 400–700 °C. The toluene reforming reaction was conducted using microwave and conventional thermal reactors. Under microwave irradiation, CO 2 and toluene conversions are boosted to 80% at 500 °C. Hydrogen and carbon monoxide yields were approximately five times and ten times higher in the microwave reactor at 500 °C, respectively, compared to the productions obtained in the conventional fixed-bed reactor at 700 °C. Filamentous carbon was also produced as a valuable side product to improve the economy of this process and such value-added carbon was only observed in the microwave reactor. Three reaction pathways were observed during microwave reaction: the toluene decomposition produces an initial hydrogen and carbon deposit on the catalyst; the formation of methane and benzene suggests toluene hydrodemethylation as a secondary reaction; and toluene hydrogenolysis forms light alkanes such as methane, and through reforming reaction under CO 2 to syngas.

10 SYNTHETIC FUELS↗

DOE OSTI2025

Mineralogical and isotopic constraints on early, high-temperature events and reservoirs recorded in the interior of a Type B Ca-Al-rich inclusion from the reduced CV3 chondrite Vigarano

A coordinated mineralogical and oxygen and Al-Mg isotopic study of a Type B Ca-Al-rich inclusion (CAI) from the reduced CV3 chondrite Vigarano was carried out using electron microprobe analysis, transmission electron microscopy, and secondary ion mass spectrometry. This CAI, a once molten igneous object, shows heterogeneous oxygen isotopic compositions among its constituent mineral phases. Spinel is uniformly 16 O-rich with Δ 17 O ≤ −23 ‰. The Δ 17 O values of Al,Ti-diopside in the CAI core range from −11 ‰ to −18 ‰, which are positively correlated with its TiO 2 contents, whereas diopside in the Wark–Lovering (WL) rim is 16 O-rich with Δ 17 O = −23 ‰. Melilite in the CAI mantle shows Δ 17 O values ≤ −10 ‰ only when Åk <15 but becomes more 16 O-poor (Δ 17 O = −3 ‰ to −9 ‰) when Åk >15 . These correlated chemical and oxygen isotopic variations were likely established during crystallization of melilite and Al,Ti-diopside from a partial melt, while spinel preserves the original 16 O-rich composition of the CAI precursor. We infer that a partial melt was isotopically evolving from 16 O-rich to 16 O-poor during melilite crystallization, then back to 16 O-rich during Al,Ti-diopside crystallization via exchange with different gas reservoirs of 16 O-poor and 16 O-rich compositions during heating event(s). Our Al-Mg isotopic measurements of the CAI core and mantle define a single isochron with an inferred initial 26 Al/ 27 Al ratio of (4.93 ± 0.18) × 10 −5 , indistinguishable from that of the WL rim. This indicates that multiple high-temperature events and oxygen isotope exchange with isotopically distinct gas reservoirs occurred rapidly during the CAI formation.

Al-Mg isotopic systematic↗

DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗