Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2024

Fundamental Thermodynamic, Kinetic, and Mechanical Properties of Lithium and Its Alloys

Lithium alloying reactions are beneficial in promoting uniform plating and stripping of lithium metal in all-solid-state batteries. First-principles calculations are performed to predict thermodynamic, kinetic, and mechanical properties of lithium and several important Li–M alloys (M = Mg, Ag, Zn, Al, Ga, In, Sn, Sb, and Bi). While the Li–Mg binary system forms a solid solution, most other lithium–metal alloys prefer stoichiometric intermetallic compounds with common local motifs that enable fast Li diffusion. Lithium and Li-rich alloys exhibit an unusually flat energy landscape along paths that connect BCC to close-packed structures like FCC and HCP, with important implications for mechanical properties. Very low migration barriers for Li diffusion that rival those of superion conductors are predicted, both in pure Li and in Li–M intermetallics. However, vacancy concentration, which is crucial for substitutional diffusion, is predicted to be low in metallic Li and most Li–M intermetallics. Compounds such as B32 LiAl and LiGa as well as D0 3 Li 3 Sb and Li 3 Bi exhibit structural vacancies at higher ends of their voltage windows, which together with low migration barriers leads to exceptionally high Li mobilities. In the Li–Mg solid solution, the addition of Mg is found to decrease the vacancy tracer diffusion coefficient by an order of magnitude.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Potential bioavailability of representative pyrogenic organic matter compounds in comparison to natural dissolved organic matter pools

Abstract. Pyrogenic organic matter (PyOM) from wildfires impacts river corridors globally and is widely regarded as resistant to biological degradation. Though recent work suggests PyOM may be more bioavailable than historically perceived, estimating bioavailability across its chemical spectrum remains elusive. To address this knowledge gap, we assessed potential bioavailability of representative PyOM compounds relative to ubiquitous dissolved organic matter (DOM) with a substrate-explicit model. The range of potential bioavailability of PyOM was greater than natural DOM; however, the predicted thermodynamics, metabolic rates, and carbon use efficiencies (CUEs) overlapped significantly between all OM pools. Compound type (e.g., natural versus PyOM) had approximately 6-fold less impact on predicted respiration rates than simulated carbon and oxygen limitations. Within PyOM, the metabolism of specific chemistries differed strongly between unlimited and oxygen-limited conditions – degradations of anhydrosugars, phenols, and polycyclic aromatic hydrocarbons (PAHs) were more favorable under oxygen limitation than other molecules. Notably, amino sugar-like, protein-like, and lignin-like PyOM had lower carbon use efficiencies relative to natural DOM of the same classes, indicating potential impacts in process-based model representations. Overall, our work illustrates how similar PyOM bioavailability may be to that of natural DOM in the river corridor, furthering our understanding of how PyOM may influence riverine biogeochemical cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Fundamental Effects of Al and Ta on Microstructure and Phase Transformations in the Al–Cr–Mo–Ta–Ti Refractory Complex Concentrated Alloy System

The effect of aluminum and tantalum concentrations on a refractory metal complex concentrated alloy is reported, particularly with respect to their effect on microstructure and phase composition of the alloy in cast and annealed form. Alloys with an equiatomic composition, (AlCrMoTaTi), an aluminum-lean composition (Al 0.75 CrMoTaTi), and a tantalum-lean composition (AlCrMoTa 0.75 Ti) are produced via arc melting. The alloys exhibit multiphase structures, confirmed by X-ray diffraction, microstructural characterization, and thermal analysis. The minor off-equiatomic adjustments of aluminum and tantalum in this alloy system did not drastically alter the prevalence of the Cr–Ta-based Laves phase. Correlations between thermodynamic predictions and observed phase transformations via thermal analysis are improved upon refinement of calculations removing impractical intermediate phases. As a result, experimental findings provide information for the refinement of thermodynamic modeling and deliver additional insight into the optimization of alloy compositions within this five-component system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermodynamic modeling of aqueous lithium salt solutions with association electrolyte nonrandom two-liquid activity coefficient model

The high charge density of lithium ion and the resulting strong association phenomena make thermodynamic modeling of aqueous lithium electrolyte solutions extremely challenging. In this study, the association electrolyte nonrandom two-liquid activity coefficient model of Lin et al. (AIChE J. 2022, 68(2), e17422) is utilized to correlate and predict thermodynamic properties and solubility behavior of aqueous single electrolyte solutions of LiCl, LiBr, LiI, and LiNO 3 , and their mixed electrolyte solutions. Capturing self-association of water, cross-association of ion and water for hydration, and cross-association of cation and anion for ion-pairing, the association model accurately represents the literature experimental data up to saturation concentrations and at the temperature ranging from 263 K to 523 K. Here, this study further investigated the effect of anions of the lithium salts, and re-confirmed that the order of solution non-ideality as LiI > LiBr > LiCl > LiNO 3 because the anions with stronger association strengths are more likely to form ion pairs and thus lower the mean ionic activity coefficients.

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

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Temperature-Dependent Thermal and Mechanical Properties of a Wire Arc Additively Manufactured Low Transformation Temperature Steel

Recent research has studied the use of low transformation temperature (LTT) martensite steel as feedstock for wire arc additive manufacturing (WAAM) and low tensile residual stresses or compressive residual stresses were detected in the printed walls. These residual stress states help to improve printed product properties such as fatigue strength and corrosion resistance. However, the thermal and mechanical properties of WAAM printed LTT martensite steel walls are largely unknown. In this work, a printed LTT martensite steel was characterized for its thermal, metallurgical, and mechanical behavior at room and elevated temperatures. The temperature-dependent specific heat capacity, thermal expansion, atomic lattice spacing, and tensile properties were measured during both heating and cooling and related to observed microstructural features and computational thermodynamics predictions. These results revealed a large hysteresis in the martensitic transformation, with a martensite start temperature of 240 °C and austenite start temperature of 680 °C. Additional thermal cycles and specimen orientation did not affect the printed specimen austenite and martensite transformations. However, it was observed that the printed metal may exhibit tempering embrittlement at about 350 °C but further studies are needed to confirm that. Further, these results suggest that a temperature control of 250 °C to 350 °C during WAAM is needed to maximize the stress reduction potential of the LTT250 martensite steel. Opportunities for future implementation of LTT martensite steels and optimization of additive manufacturing process conditions are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Destabilizing high-capacity high entropy hydrides via earth abundant substitutions: From predictions to experimental validation

The vast chemical space of high entropy alloys (HEAs) makes trial-and-error experimental approaches for materials discovery intractable and often necessitates data-driven and/or first principles computational insights to successfully target materials with desired properties. In the context of materials discovery for hydrogen storage applications, a theoretical prediction-experimental validation approach can vastly accelerate the search for substitution strategies to destabilize high-capacity hydrides based on benchmark HEAs, e.g. TiVNbCr alloys. Here, in this study, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb) 75 Cr 25-x Fe x system. The as-prepared alloys crystallize in a single-phase bcc lattice for limited Fe content x < 7, while larger Fe content favors the formation of a secondary C14 Laves phase intermetallic. Short range order for alloys with x < 7 can be well described by a random distribution of atoms within the bcc lattice without lattice distortion. Hydrogen absorption experiments performed on selected alloys validate the predicted thermodynamic destabilization of the corresponding fcc hydrides and demonstrate promising lifecycle performance through reversible absorption/desorption. This demonstrates the potential of computationally expedited hydride discovery and points to further opportunities for optimizing bcc alloy ↔ fcc hydrides for practical hydrogen storage applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of elevated temperature exposure on the residual stress state and microstructure of PVD Cr coatings on SiC investigated via in situ X-ray diffraction and transmission electron microscopy

Pure Cr coatings are being considered as corrosion-resistant layers for SiC-SiC f composite fuel cladding in light water reactors. The residual stress state and interdiffusion behavior of Cr coatings deposited on SiC substrates using bipolar high-power impulse magnetron sputtering (B-HiPIMS) and direct current impulse magnetron sputtering (DCMS) were studied using in situ X-ray diffraction experiments at 300 °C, 700 °C, and 1200 °C. After elevated temperature exposure, evolution of the coating microstructure and coating/substrate interdiffusion was characterized by transmission electron microscopy. The compressive residual stress of the dense B-HiPIMS coating at room temperature was gradually relieved at 300 °C, while the stress transitioned to a tensile state at 700 °C. Grain growth of the Cr coating and a Cr–SiC interdiffusion layer consisting of newly formed Cr–C and Cr–Si phases were observed in the B-HiPIMS coating after exposure at 700 °C. The phases in the interdiffusion layer were identified to be Cr 7 C 3 , Cr 3 Si, and Cr 5 Si 3 C x experimentally and by thermodynamic predictions. For the DCMS coating, however, changes in residual stress and microstructure at the elevated temperatures were limited by the inherent porosity and inhomogeneous structure of the as-deposited coating. At 1200 °C, both coating types interacted rapidly with the underlying SiC substrate in the early stages of the experiment. The results from this study suggest that the characteristics of the Cr coatings do not significantly change at the reactor operating temperature of 300 °C, while changes in microstructure and stress state are expected in off-normal conditions of 700 °C and 1200 °C. Finally the results from this study suggest the B-HiPIMS coating is best-suited as a protective coating for SiC fuel cladding.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Partitioning of tramp elements Cu and Si in a Ni-based superalloy and their effect on creep properties

An alloy's processing history, including melting, remelting or the choice of stock material, affects its purity and eventually involves tramp element pickup or retention. Herein, variants of a novel Ni-based superalloy were manufactured with different levels of purity. The so-called low purity alloys contained 0.138 wt. % Cu and 0.019 wt. % Si while the Cu and Si levels were below x-ray fluorescence (XRF) detection limits of 0.003 and 0.010 wt. %, respectively, in the high purity ingots. Atom-probe tomography (APT) was carried out and revealed Si partitioning at the following interfaces: grain boundaries, MC carbide/γ, M 3 B 2 boride/$γ$ and M 3 B 2 boride/$γ$'. Copper was found to primarily segregate to the $γ$' precipitates. An average of 2.4 × decrease in creep life and 4.3 × decrease in creep ductility was measured in the low purity alloys, which was attributed to the embrittlement caused by Si segregation to grain boundaries. Furthermore, the positive effect of B on the creep properties was mitigated by the presence of Si. Thermodynamic predictions for the matrix and γ' precipitate compositions represented the trends observed experimentally although the extent of preferential partitioning lacks accuracy. Monte Carlo simulations were performed to describe the partitioning of Cu and Si atoms to either $γ$ or $γ$' phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Air separation and N 2 purification with Ba 0.15 Sr 0.85 FeO 3-$δ$ via a two-step thermochemical process

Thermochemical air separation to produce high-purity N 2 was demonstrated in a vertical tube reactor via a two-step reduction–oxidation cycle with an A-site substituted perovskite Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585). BSF1585 particles were synthesized and characterized in terms of their chemical, morphological, and thermophysical properties. A thermodynamic cycle model and sensitivity analysis using computational heat and mass transfer models of the reactor were used to select the system operating parameters for a concentrating solar thermal-driven process. Thermal reduction up to 800 °C in air and temperature-swing air separation from 800 °C to minimum temperatures between 400 and 600 °C were performed in the reactor containing a 35 g packed bed of BSF1585. The reactor was characterized for dispersion, and air separation was characterized via mass spectrometry. Gas measurements indicated that the reactor produced N 2 with O 2 impurity concentrations as low as 0.02 % for > 30 min of operation. Additionally, a parametric study of air flow rates suggested that differences in observed and thermodynamically predicted O 2 impurities were due to imperfect gas transport in the bed. Temperature swing reduction/oxidation cycling experiments between 800 and 400 °C in air were conducted with no statistically significant degradation in N 2 purity over 50 cycles.

14 SOLAR ENERGY↗
DOE OSTI2025

Electrochemical Oxidation in Garnet-Type Solid Electrolyte by Formation of Point Defects

All-solid-state batteries hold greater promise for improving safety and energy density over conventional battery technology employing organic liquid electrolytes. One of the required features of a Li + conducting solid electrolyte is electrochemical stability, attained thermodynamically or kinetically, within the targeted operating voltage and temperature ranges. Therefore, understanding of the oxidative or reductive degradation mechanism is important to allow the design of stable solid electrolyte materials. This work contributes to building an understanding of the oxidative degradation mechanism in lithium solid electrolytes at cell operating conditions. Here, we have focused on resolving the oxidative decomposition mechanism of Al-doped lithium garnet Li 6.28 Al 0.24 La 3 Zr 2 O 12 (LLZO) as a state-of-the-art inorganic ceramic electrolyte. By combining experimental and computational analyses, we show that oxidation of LLZO occurs by simultaneous loss of oxygen and lithium from the structure, resulting in substoichiometric LLZO, at a moderate temperature (80 °C) and a high electrode potential (4.3 V vs Li/Li + ). Based on X-ray absorption and diffraction analyses, we find that the zirconium coordination shells in LLZO contract while the crystal structure experiences positive chemical strain upon electrochemical oxidation. The results from ex situ structural characterization of both the local structure and crystal symmetry are supported by a substoichiometric LLZO with lithium and oxygen vacancies, modeled by density functional theory (DFT) calculations. These chemical and structural changes in LLZO suppress effective lithium-ion conductivity by an order of magnitude. Formation of lithium and oxygen vacancies in LLZO upon electrochemical oxidation is different from prior thermodynamic predictions of phase decomposition of LLZO. The difference here is that the experiments were conducted at near-room temperature, which can hinder the kinetics of phase separation, and thus, the resultant LLZO solid electrolyte is still single-phase but substoichiometric in Li and O. In conclusion, these findings contribute an important degradation mechanism of the electrolyte, relevant for practical operational conditions of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Iodine Vapor Reactions with Pure Metal Wires at Temperatures of 100–139 °C in Air

In this work, I 2(g) reactions with 11 metal wires (i.e., Al, Ag, Cu, In, Mo, Nb, Ni, Pd, Pt, Sn, Ta) were evaluated at three temperatures (i.e., 100 ± 3°C, 123 ± 4°C, 139 ± 5°C) with exposure times of 24 hours at each temperature. Over these temperatures, some of the metals showed increased mass gain with higher temperatures (i.e., In, Ag, Cu), Sn showed decreased iodine capture at increased temperatures, and most metals (i.e., Mo, Nb, Ni, Pd, Pt, and Ta) showed little to no mass gain at all temperatures. Silver mordenite (AgZ) was used as a standard during these studies and showed a consistent mass gain (m% I ) of 10.6 → 12.8% over this temperature range. The values of mass of iodine captured per mass of starting metal (g g -1 ) ranged widely across the study with the highest values achieved for (in descending order) Sn-T 100°C (4.37), In-T 139°C (3.34), Sn-T 123°C (1.98), In-T 123°C (1.57), Ag-T 139°C (1.19), In-T 100°C (0.97), and Cu-T 139°C (0.71). In some cases, the metal-iodide complex was not stable at the experimental temperature and it was clear some volatility had occurred during the experiment based on discoloration in the vials. These results show that some metals can have extensive reactions with I 2(g) without showing metal-iodide preferences over metal-oxide formation based on thermodynamic predictions. It is possible that materials such as these could be implemented near nuclear facilities to getter I 2(g) in the event of a nuclear accident.

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

Itinerant Magnetism in Hydride-Synthesized CaCo 12 B 6

A new compound in the underexplored Ca–Co–B phase space has been discovered, validating high-throughput computations from the Open Quantum Materials Database, which predicted thermodynamic stability for CaCo 12 B 6 in the SrNi 12 B 6 structure type. The synthetic effects of different boron precursors and the advantages of using CaH 2 instead of Ca metal were demonstrated by the short synthesis duration and high purity of CaCo 12 B 6 , in contrast with traditional synthesis routes. Powder X-ray diffraction (PXRD) confirmed that CaCo 12 B 6 shares the SrNi 12 B 6 structure ( R $\bar{3}$m (#166), a = 9.469(4) Å, c = 7.468(2) Å, Z = 3) and is water- and air-stable. High-temperature in situ PXRD indicates that CaCo 12 B 6 is stable below 1050 K under vacuum in a silica capillary. CaCo 12 B 6 decomposes between 693 and 773 K during spark-plasma sintering. Density functional theory calculations indicate that CaCo 12 B 6 is metallic with a ferromagnetic ground state. X-ray absorption near-edge spectroscopy and Bader charge analysis indicate that Co atoms in CaCo 12 B 6 lack ionic character. Magnetometry reveals room-temperature paramagnetism with μ eff = 1.7(1)μ B per Co atom and a Weiss constant of +190(10)K. Ferromagnetic ordering occurs below 172(1)K, resulting in a saturation moment of 0.46 μ B per Co atom. Our findings demonstrate that the hydride route is a viable strategy for discovery of new ternary alkaline-earth-transition metal borides analogous to rare-earth-containing counterparts.

diffraction↗
DOE OSTI2022

Ultraviolet Photooxidation of Smectite-Bound Fe(II) and Implications for the Origin of Martian Nontronites

Clay minerals detected with orbital and in situ instruments in ancient Martian terrains constrain Mars' climate and aqueous alteration history. Early in its history, Mars experienced an atmospheric redox change and iron-bearing clay minerals may preserve the effects of that transition. Ferrous smectites, the thermodynamically predicted product of chemical weathering of basalts under anoxic conditions, may have undergone oxidation by exposure to chemical oxidants in the atmosphere or regolith, or by direct photooxidation at the surface. To assess these potential oxidation pathways, ferrous trioctahedral smectites of varying initial iron content were synthesized and subjected to oxidation by ultraviolet (UV) irradiation. Experimental UV irradiation under an anoxic atmosphere equivalent to approximately 7 years of flux on the Martian surface caused partial oxidation of smectite-bound Fe (Fe 3+ /ΣFe = 16–18%) and octahedral sheet contraction. Metal-OH vibrational bands in visible/near infrared (VNIR) reflectance spectra of oxidized smectites changed in band depth and asymmetry with higher iron content. X-ray diffraction patterns of UV irradiated samples indicate the formation of a mixed di- and trioctahedral smectite or a secondary nontronite phase, possibly on the surfaces of higher iron content smectites. Furthermore, these experiments suggest that UV irradiation is able to oxidize structurally bound iron in smectites without the presence of other chemical oxidants. Photooxidation may have influenced the mineralogy, both syndepositionally and postdepositionally, of Martian alteration assemblages formed near the surface and this process needs not be limited to one part of their formation history.

58 GEOSCIENCES↗
DOE OSTI2021

Investigations of the stability of GaAs for photoelectrochemical H 2 evolution in acidic or alkaline aqueous electrolytes

The long-term stability of p-GaAs photocathodes has been investigated for the hydrogen-evolution reaction (HER) in contact with either 1.0 M H 2 SO 4 (aq) or 1.0 M KOH(aq). Stability for the HER was evaluated using p-GaAs electrodes that were either etched or coated with active HER catalysts (Pt and CoP). Changes in surface characteristics of GaAs after exposure to electrochemical conditions were monitored by X-ray photoelectron spectroscopy (XPS), and electrode dissolution processes were evaluated by inductively coupled plasma mass spectrometry (ICP-MS). Consistent with thermodynamic predictions, after operation of the HER at pH 0 or pH 14, illuminated etched p-GaAs electrodes exhibited minimal dissolution while preserving a nearly stoichiometric surface. Electrodeposition or sputtering of Pt on the p-GaAs surface promoted the formation of excess As 0 via an interfacial reaction during the HER. The resulting non-stoichiometric As 0 -rich surface of p-GaAs/Pt electrodes caused a loss in photoactivity as well as substantial cathodic dark current. In contrast, p-GaAs electrodes coated with thin-film CoP catalysts did not display an increase in surficial As 0 after operation of the HER in acidic electrolytes. Minimization of deleterious interfacial reactions is thus critical to obtain extended stability in conjunction with high performance from p-GaAs photocathodes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗