SEARCH · Search NASA

Results for “cluster expansion”

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 OSTI2024

Towards accurate prediction of configurational disorder properties in materials using graph neural networks

Abstract The prediction of configurational disorder properties, such as configurational entropy and order-disorder phase transition temperature, of compound materials relies on efficient and accurate evaluations of configurational energies. Previous cluster expansion methods are not applicable to configurationally-complex material systems, including those with atomic distortions and long-range orders. In this work, we propose to leverage the versatile expressive capabilities of graph neural networks (GNNs) for efficient evaluations of configurational energies and present a workflow combining attention-based GNNs and Monte Carlo simulations to calculate the disorder properties. Using the dataset of face-centered tetragonal gold copper without and with local atomic distortions as an example, we demonstrate that the proposed data-driven framework enables the prediction of phase transition temperatures close to experimental values. We also elucidate that the variance of the energy deviations among configurations controls the prediction accuracy of disorder properties and can be used as the target loss function when training and selecting the GNN models. The work serves as a fundamental step toward a data-driven paradigm for the accelerated design of configurationally-complex functional material systems.

Chemistry↗

DOE OSTI2026

Optimal invariant sets for atomistic machine learning

The representation of atomic configurations for machine learning models has led to numerous sets of descriptors. However, many descriptor sets are incomplete and/or functionally dependent. Incomplete sets cannot faithfully represent atomic environments. Yet complete constructions often suffer from a high degree of functional dependence, where some descriptors are functions of others. These redundant descriptors do not improve discrimination between atomic environments. We employ pattern recognition techniques to remove dependent descriptors to produce the smallest possible set that satisfies completeness. We apply this in two ways: First, we refine an existing description, the atomic cluster expansion. Second, we augment an incomplete construction, yielding a new message-passing neural network architecture that can recognize up to 5-body patterns. This architecture shows strong accuracy on state-of-the-art benchmarks while retaining low computational cost. Our results demonstrate the utility of this strategy to optimize descriptor sets across a range of descriptors and application datasets.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗

DOE OSTI2020

Prediction of structure and cation ordering in an ordered normal-inverse double spinel

Abstract Spinels represent an important class of technologically relevant materials, used in diverse applications ranging from dielectrics, sensors and energy materials. While solid solutions combining two “single spinels” have been explored in a number of past studies, no ordered “double” spinels have been reported. Based on our first principles computations, here we predict the existence of such a double spinel compound MgAlGaO 4 , formed by an equimolar mixing of MgAl 2 O 4 normal and MgGa 2 O 4 inverse spinels. After studying the details of its atomic and electronic structure, we use a cluster expansion based effective Hamiltonian approach with Monte Carlo simulations to study the thermodynamic behavior and cation distribution as a function of temperature. Our simulations provide strong evidence for short-ranged cation order in the double spinel structure, even at significantly elevated temperatures. Finally, an attempt was made to synthesize the predicted double spinel compound. Energy Dispersive X-ray Spectrometry and X-ray diffraction Rietveld refinements were performed to characterize the single-phase chemical composition and local configurational environments, which showed a favorable agreement with the theoretical predictions. These findings suggest that a much larger number of compounds can potentially be realized within this chemical space, opening new avenues for the design of spinel-structured materials with tailored functionality.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study

The development of high-entropy alloys (HEAs) focuses on exploring compositional regions in multi-component systems with all alloy elements in equal or near-equal atomic concentrations. Initially it was based on the main idea that high mixing configurational entropy contributions to the alloy free energy could promote the formation of a single solid solution phase. By using the ab-initio based Cluster Expansion (CE) Hamiltonian model constructed for the quinary bcc Cr–Ta–Ti–V–W system in combination with Monte Carlo (MC) simulations, we show that the phase stability and chemical short-range order (SRO) of the equiatomic quinary and five sub-quaternary systems, as well as their derivative alloys, can dramatically change the order–disorder transition temperatures (ODTT) as a function of alloy compositions. In particular, it has been found, that the equiatomic quaternary Ta–Ti–V–W and Cr–Ta–Ti–W alloys had the lowest order–disorder transition temperature (500 K) among all the analysed equiatomic compositions. In all investigated alloy systems, the strongest chemical ordering has been observed between Cr and V, which led to the conclusion that decreasing the concentration of either Cr or V might be beneficial in terms of decreasing the ODTT. It also predicts that increasing concentration of Ti significantly decreases the ODTT. Our analysis of chemical SRO as a function of alloy composition allows to understand the microstructure evolution of HEAs as a function of temperature in excellent agreement with available experimental observations. Importantly, our free energy of mixing and SRO calculations predict that the origin of precipitates formed by Cr- and V-rich in the sub-quaternary Cr–Ta–V–W system is driven by the thermodynamics. The modelling results are in an excellent agreement with experimental observation of Cr and V segregation in the W 0.38 Ta 0.36 Cr 0.15 V 0.11 alloy which in turns shows an exceptional radiation resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Ab initio study of short-range ordering in vanadium-based disordered rocksalt structures

Disordered rocksalt Li-excess (DRX) compounds are attractive new cathode materials for Li-ion batteries as they contain resource-abundant metals and do not require the use of cobalt or nickel. Understanding the delithiation process and cation short-range ordering (SRO) in DRX compounds is essential to improve these promising cathode materials. Herein, we use first-principles calculations along with the cluster-expansion approach to model the disorder in DRX Li 2-x VO 3 , 0 < x < 1. We discuss the SRO of Li in tetrahedral and octahedral sites, and the order in which Li delithiates and V oxidizes with respect to local environments. We reveal that the number of nearest-neighbor V dictates the order of delithiation from octahedral sites and that V is oxidized in a manner that minimizes the electrostatic interactions among V. In conclusion, our results provide valuable insight for tailoring the performance of V-based DRX cathode materials in general by controlling the SRO features that reduce energy density.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗

DOE OSTI2025

Structural origin of disorder-induced ion conduction in NaFePO 4 cathode materials

Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.

Christensen, Rasmus↗

DOE OSTI2020

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Finite-temperature kinetic ferromagnetism in the square-lattice Hubbard model

While the exact phase diagram of the Fermi-Hubbard model remains poorly understood despite decades of progress, nearly 60 years ago Nagaoka proved that a single dopant in an otherwise halffilled Hubbard system can bring about ferromagnetism through kinetic means. The phenomenon was recently observed with ultracold atoms in triangular optical lattices. Here, we explore the kinetic ferromagnetism within the square lattice Hubbard model and its strong-coupling counterpart, the t − J model, at finite temperatures in the thermodynamic limit via numerical linked-cluster expansions. We find evidence of ferromagnetic Nagaoka polarons at dopings up to ∼ 30% away from half filling for a variety of interaction strengths and at temperatures as low as 0.2 of the hopping energy. We map out the boundaries of this phase through analyzing various correlation functions.

Newby, Robin C. [San José State University] (ORCID↗

DOE OSTI2024

Relativistic corrections to the correlated basis function effective nuclear Hamiltonian

We discuss the inclusion of relativistic boost corrections into the correlated basis function effective nuclear Hamiltonian, derived from a realistic model of two- and three-nucleon interactions using the formalism of correlated basis functions and the cluster expansion technique. Different procedures to take into account the effects of boost interactions are compared on the basis of the ability to reproduce the nuclear matter equation of state obtained from accurate quantum many-body calculations. Furthermore, the results of our study show that the repulsive contribution of the boost interaction significantly depends on the underlying model of the nonrelativistic potential. On the other hand, the dominant relativistic correction turns out to be the corresponding reduction of the strength of repulsive three-nucleon interactions, leading to a significant softening of the equation of state of nuclear matter at supranuclear densities.

Neutron stars & pulsars↗

DOE OSTI2024

Equation of State and Thermometry of the 2D SU ( N ) Fermi-Hubbard Model

We characterize the equation of state (EoS) of the SU ( N > 2 ) Fermi-Hubbard Model (FHM) in a two-dimensional single-layer square optical lattice. We probe the density and the site occupation probabilities as functions of interaction strength and temperature for N = 3 , 4, and 6. Our measurements are used as a benchmark for state-of-the-art numerical methods including determinantal quantum Monte Carlo and numerical linked cluster expansion. By probing the density fluctuations, we compare temperatures determined in a model-independent way by fitting measurements to numerically calculated EoS results, making this a particularly interesting new step in the exploration and characterization of the SU ( N ) FHM. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗

DOE OSTI2024

Dipolar Spin Ice Regime Proximate to an All-In-All-Out Néel Ground State in the Dipolar-Octupolar Pyrochlore Ce 2 Sn 2 O 7

The dipolar-octupolar (DO) pyrochlores, R 2 M 2 O 7 ( R = Ce , Sm , Nd ), are key players in the search for realizable novel quantum spin liquid (QSL) states as a large parameter space within the DO pyrochlore phase diagram is theorized to host QSL states of both dipolar and octupolar nature. New single crystals and powders of Ce 2 Sn 2 O 7 , synthesized by hydrothermal techniques, present an opportunity for a new characterization of the exchange parameters in Ce 2 Sn 2 O 7 using the near-neighbor X Y Z model Hamiltonian associated with DO pyrochlores. Utilizing quantum numerical linked cluster expansion fits to heat capacity and magnetic susceptibility measurements, and classical Monte Carlo calculations to the diffuse neutron diffraction of the new hydrothermally grown Ce 2 Sn 2 O 7 samples, we place Ce 2 Sn 2 O 7 ’s ground state within the ordered dipolar all-in-all-out (AIAO) Néel phase, with quantum Monte Carlo calculations showing a transition to long-range order at temperatures below those accessed experimentally. Indeed, our new neutron diffraction measurements on the hydrothermally grown Ce 2 Sn 2 O 7 powders show a broad signal at low scattering wave vectors, reminiscent of a spin ice, in striking contrast from previous powder neutron diffraction on samples grown from solid-state synthesis, which found diffuse scattering at high scattering wave vectors associated with magnetic and suggested an octupolar quantum spin ice state. We conclude that new hydrothermally grown Ce 2 Sn 2 O 7 samples host a finite-temperature proximate dipolar spin ice phase, above the expected transition to AIAO Néel order. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗

DOE OSTI2025

Efficient machine learning interatomic potentials robust for liquid and multiple solid polymorphs of NaF and KF

Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.

Crystal melting↗

DOE OSTI2025

HydraGNN v4.0

The new version of HydraGNN v4.0 provides additional core capabilities, such as: Inclusion of multi-body atomistic cluster expansion MACE, polarizable atom interaction neural network PAINN, and equivariant principal neighborhood aggregation (PNAEq) among the message passing layers supported -Inclusion of graph transformers to directly model long-range interactions between nodes that are distant in the graph topology Integration of graph transformers with message passing layers by combining the graph embedding generated by the two mechanisms, which allows for an improved expressivity of the HydraGNN architecture Improved re-implementation of multi-task learning (MTL) to allow its use for stabilized training across imbalanced, multi-source, multi-fidelity data Introduction of multi-task parallelism, a newly proposed type of model parallelism specifically for MTL architectures, which allows to dispatch different output decoding heads to different GPU devices Integration of multi-task parallelism with pre-existing distributed data parallelism to enable a 2D parallelization for distributed training Improved portability of the distributed training across Intel GPUs, which has been testes on ALCF exascale supercomputer Aurora Inclusion of 2-level fine-grained energy profilers portable across NVIDIA, AMD, and Intel GPUs to monitor the power and energy consumption associated with different functions executed by the HydraGNN code during data pre-load and training Restructuring of previous examples and inclusion of new sets of examples to illustrate the download, preprocess, and training of HydraGNN models on new large-scale open-source datasets for atomistic materials modeling (e.g., Alexandria, Transition1x, OMat24, OMol25)

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗