Publications

Abstract not provided.

Cost-effective and reliable hydrogen compression remains a challenging barrier in the widespread adoption of hydrogen as an energy carrier. The prevailing technology of mechanical compression suffers from several drawbacks, some of which can be addressed by nonmechanical compression strategies (e.g., electrochemical or metal hydride-based thermal compression). Thermally driven metal hydride compression strategies typically rely on multistage metal hydride-based compressors; however, discovering or optimizing low-stability metal hydrides that can pressurize hydrogen upward of 1000 bar is difficult, both with respect to computational predictions and experimental validation. Here, we (1) demonstrate that simple machine learning-derived design rules can inform the rational design of alloying strategies yielding low-stability hydrides, (2) validate their experimental pressure–composition–temperature (PCT) isotherms up to 875 bar, and (3) utilize a dynamic system-level model of a metal hydride compressor design to evaluate their performance under realistic operating conditions. Importantly, this analysis yields predicted operational efficiencies of both 2-stage (90–875 bar) and 3-stage (20–875 bar) metal hydride compressors to enable further evaluation of this technology and its techno-economic outlook.

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., O(102)GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the O(100)GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va)3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. In this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

A series of high entropy AB2-type Ti2-xZrxMnCrFeNi alloys (x = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2) were synthesized to investigate their potential for hydrogen storage and chemical compression. The influence of the Ti/Zr ratio was explored in terms of structural, microstructural and thermodynamic properties. The storage capacity together with the reaction enthalpy and entropy changes of the synthesized high entropy alloys were compared to predictions from Machine Learning (ML) to investigate changes in these properties across the explored composition space. The results revealed that a decreasing Zr content consistently lowered the hydride formation enthalpy and increased the plateau pressure from 8 to >90 bar H2 at 25 °C, in good agreement with ML predictions. Selected compositions (x = 1.0 and 1.2) demonstrated reversible hydrogen storage capability over 150 cycles, with capacities of 1.34–1.40 wt % H2 and remarkable reaction kinetics (<4 min) at ambient temperature. These experimental and computational findings highlight the potential of this Laves-HEA system as tuneable, stable, and cost-effective materials suitable for long-term operations in stationary hydrogen storage and compression applications.

A trilateral agreement has been finalized involving research institutions in Korea, Japan, and the U.S. The project partners are Sandia, LLNL, KIST, KAIST, and AIST. The project title is “Structure-Property Relationships in Metal Alloys for Hydrogen Storage and Processing.” Funding for the U.S. portion of the effort is through NNSA; the PI is Vitalie Stavila. The overall objective of this project is to identify detailed structure-property relationships governing hydrogen separation, purification, storage, and compression in compositionally complex metal alloys.

Machine learning (ML) models in the materials sciences that are validated by overly simplistic cross-validation (CV) protocols can yield biased performance estimates for downstream modeling or materials screening tasks. This can be particularly counterproductive for applications where the time and cost of failed validation efforts (experimental synthesis, characterization, and testing) are consequential. We propose a set of standardized and increasingly difficult splitting protocols for chemically and structurally motivated CV that can be followed to validate any ML model for materials discovery. Among several benefits, this enables systematic insights into model generalizability, improvability, and uncertainty, provides benchmarks for fair comparison between competing models with access to differing quantities of data, and systematically reduces possible data leakage through increasingly strict splitting protocols. Performing thorough CV investigations across increasingly strict chemical/structural splitting criteria, local vs. global property prediction tasks, small vs. large datasets, and structure vs. compositional model architectures, some common threads are observed; however, several marked differences exist across these exemplars, indicating the need for comprehensive analysis to fully understand each model's generalization accuracy and potential for materials discovery. For this we provide a general-purpose, featurization-agnostic toolkit, MatFold, to automate reproducible construction of these CV splits and encourage further community use in model benchmarking.

There is growing interest in material candidates with properties that can be engineered beyond traditional design limits. Compositionally complex oxides (CCO), often called high entropy oxides, are excellent candidates, wherein a lattice site shares more than four cations, forming single-phase solid solutions with unique properties. However, the nature of compositional complexity in dictating properties remains unclear, with characteristics that are difficult to calculate from first principles. Here, compositional complexity is demonstrated as a tunable parameter in a spin-transition oxide semiconductor La1− x(Nd, Sm, Gd, Y)x/4CoO3, by varying the population x of rare earth cations over 0.00≤ x≤ 0.80. Across the series, increasing complexity is revealed to systematically improve crystallinity, increase the amount of electron versus hole carriers, and tune the spin transition temperature and on-off ratio. At high a population (x = 0.8), Seebeck measurements indicate a crossover from hole-majority to electron-majority conduction without the introduction of conventional electron donors, and tunable complexity is proposed as new method to dope semiconductors. First principles calculations combined with angle resolved photoemission reveal an unconventional doping mechanism of lattice distortions leading to asymmetric hole localization over electrons. Thus, tunable complexity is demonstrated as a facile knob to improve crystallinity, tune electronic transitions, and to dope semiconductors beyond traditional means.

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, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb)75Cr25-xFex 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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.

High-entropy alloys (HEAs) represent an interesting alloying strategy that can yield exceptional performance properties needed across a variety of technology applications, including hydrogen storage. Examples include ultrahigh volumetric capacity materials (BCC alloys → FCC dihydrides) with improved thermodynamics relative to conventional high-capacity metal hydrides (like MgH2), but still further destabilization is needed to reduce operating temperature and increase system-level capacity. In this work, we demonstrate efficient hydride destabilization strategies by synthesizing two new Al0.05(TiVNb)0.95-xMox (x = 0.05, 0.10) compositions. We specifically evaluate the effect of molybdenum (Mo) addition on the phase structure, microstructure, hydrogen absorption, and desorption properties. Both alloys crystallize in a bcc structure with decreasing lattice parameters as the Mo content increases. The alloys can rapidly absorb hydrogen at 25 °C with capacities of 1.78 H/M (2.79 wt %) and 1.79 H/M (2.75 wt %) with increasing Mo content. Pressure-composition isotherms suggest a two-step reaction for hydrogen absorption to a final fcc dihydride phase. The experiments demonstrate that increasing Mo content results in a significant hydride destabilization, which is consistent with predictions from a gradient boosting tree data-driven model for metal hydride thermodynamics. Furthermore, improved desorption properties with increasing Mo content and reversibility were observed by in situ synchrotron X-ray diffraction, in situ neutron diffraction, and thermal desorption spectroscopy.

The objective of this project was to evaluate material- and chemical-based solutions for hydrogen storage in rail applications as an alternative to high-pressure hydrogen gas and liquid hydrogen. Three use cases were assessed: yard switchers, long-haul locomotives, and tenders. Four storage options were considered: metal hydrides, nanoporous sorbents, liquid organic hydrogen carriers, and ammonia, using 700 bar compressed hydrogen as a benchmark. The results suggest that metal hydrides, currently the most mature of these options, have the highest potential. Storage in tenders is the most likely use case to be successful, with long-haul locomotives the least likely due to the required storage capacities and weight and volume constraints. Overall, the results are relevant for high-impact regions, such as the South Coast Air Quality Management District, for which an economical vehicular hydrogen storage system with minimal impact on cargo capacity could accelerate adoption of fuel cell electric locomotives. The results obtained here will contribute to the development of technical storage targets for rail applications that can guide future research. Moreover, the knowledge generated by this project will assist in development of material-based storage for stationary applications such as microgrids and backup power for data centers.

Hydrogen diffusion in metals and alloys plays an important role in the discovery of new materials for fuel cell and energy storage technology. While analytic models use hand-selected features that have clear physical ties to hydrogen diffusion, they often lack accuracy when making quantitative predictions. Machine learning models are capable of making accurate predictions, but their inner workings are obscured, rendering it unclear which physical features are truly important. To develop interpretable machine learning models to predict the activation energies of hydrogen diffusion in metals and random binary alloys, we create a database for physical and chemical properties of the species and use it to fit six machine learning models. Our models achieve root-mean-squared errors between 98-119 meV on the testing data and accurately predict that elemental Ru has a large activation energy, while elemental Cr and Fe have small activation energies. By analyzing the feature importances of these fitted models, we identify relevant physical properties for predicting hydrogen diffusivity. While metrics for measuring the individual feature importances for machine learning models exist, correlations between the features lead to disagreement between models and limit the conclusions that can be drawn. Instead grouped feature importance, formed by combining the features via their correlations, agree across the six models and reveal that the two groups containing the packing factor and electronic specific heat are particularly significant for predicting hydrogen diffusion in metals and random binary alloys. This framework allows us to interpret machine learning models and enables rapid screening of new materials with the desired rates of hydrogen diffusion.

The hydrogen sorption properties of single-phase bcc (TiVNb)100-xCrx alloys (x = 0-35) are reported. All alloys absorb hydrogen quickly at 25 °C, forming fcc hydrides with storage capacity depending on the Cr content. A thermodynamic destabilization of the fcc hydride is observed with increasing Cr concentration, which agrees well with previous compositional machine learning models for metal hydride thermodynamics. The steric effect or repulsive interactions between Cr-H might be responsible for this behavior. The cycling performances of the TiVNbCr alloy show an initial decrease in capacity, which cannot be explained by a structural change. Pair distribution function analysis of the total X-ray scattering on the first and last cycled hydrides demonstrated an average random fcc structure without lattice distortion at short-range order. If the as-cast alloy contains a very low density of defects, the first hydrogen absorption introduces dislocations and vacancies that cumulate into small vacancy clusters, as revealed by positron annihilation spectroscopy. Finally, the main reason for the capacity drop seems to be due to dislocations formed during cycling, while the presence of vacancy clusters might be related to the lattice relaxation. Having identified the major contribution to the capacity loss, compositional modifications to the TiVNbCr system can now be explored that minimize defect formation and maximize material cycling performance.

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. It will help to tackle future materials discovery problems in clean energy and beyond.

Abstract not provided.

The ability to rapidly screen material performance in the vast space of high entropy alloys is of critical importance to efficiently identify optimal hydride candidates for various use cases. Given the prohibitive complexity of first principles simulations and large-scale sampling required to rigorously predict hydrogen equilibrium in these systems, we turn to compositional machine learning models as the most feasible approach to screen on the order of tens of thousands of candidate equimolar high entropy alloys (HEAs). Critically, we show that machine learning models can predict hydride thermodynamics and capacities with reasonable accuracy (e.g. a mean absolute error in desorption enthalpy prediction of ∼5 kJ molH2−1) and that explainability analyses capture the competing trade-offs that arise from feature interdependence. We can therefore elucidate the multi-dimensional Pareto optimal set of materials, i.e., where two or more competing objective properties can't be simultaneously improved by another material. This provides rapid and efficient down-selection of the highest priority candidates for more time-consuming density functional theory investigations and experimental validation. Various targets were selected from the predicted Pareto front (with saturation capacities approaching two hydrogen per metal and desorption enthalpy less than 60 kJ molH2−1) and were experimentally synthesized, characterized, and tested amongst an international collaboration group to validate the proposed novel hydrides. Additional top-predicted candidates are suggested to the community for future synthesis efforts, and we conclude with an outlook on improving the current approach for the next generation of computational HEA hydride discovery efforts.