Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interatomic mixing”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A new interatomic potential for mixed Mg-Al-Ga-In spinels

While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. Here, in this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO 4 structure while correctly predicting that neither MgAlInO 4 nor MgGaInO 4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimization of Thermal Conductance at Interfaces Using Machine Learning Algorithms

We report optimization of thermal transport across the interface of two different materials is critical to micro-/nanoscale electronic, photonic, and phononic devices. Although several examples of compositional intermixing at the interfaces having a positive effect on interfacial thermal conductance (ITC) have been reported, an optimum arrangement has not yet been determined because of the large number of potential atomic configurations and the significant computational cost of evaluation. On the other hand, computation-driven materials design efforts are rising in popularity and importance. Yet, the scalability and transferability of machine learning models remain as challenges in creating a complete pipeline for the simulation and analysis of large molecular systems. In this work we present a scalable Bayesian optimization framework, which leverages dynamic spawning of jobs through the Message Passing Interface (MPI) to run multiple parallel molecular dynamics simulations within a parent MPI job to optimize heat transfer at the silicon and aluminum (Si/Al) interface. We found a maximum of 50% increase in the ITC when introducing a two-layer intermixed region that consists of a higher percentage of Si. Because of the random nature of the intermixing, the magnitude of increase in the ITC varies. We observed that both homogeneity/heterogeneity of the intermixing and the intrinsic stochastic nature of molecular dynamics simulations account for the variance in ITC.

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

Energy-efficient, Large-scale Molecular Dynamics Simulations via Hardware- and Algorithm-level Optimization

This work aims to develop a framework for energy-efficient computing that will enable molecular dynamics (MD) simulations of large-scale phenomena with atomic precision and simultaneously remove computational bottlenecks limiting the speed of MD simulations. We seek to implement such an approach through the development of surrogate models for the interatomic force calculation combined with the use of mixed numerical precision formats. For a model system of neutral atoms (only pairwise interactions), significant force calculation efficiency improvements were achieved, without detrimental effects on atomic structures or average energies, using single precision, by developing a surrogate model (deep neural network), and by quantizing this surrogate model. For a model system of charged atoms, the reciprocal-space calculation of electrostatic interactions was identified as the main bottleneck, and the development of a surrogate model should be pursued to achieve an estimated one-order-of-magnitude additional speedup.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Suppressed phase segregation for triple-junction perovskite solar cells

The tunable band gaps and facile fabrication of perovskites make them attractive for multi-junction photovoltaics. However, light-induced phase segregation limits their efficiency and stability: this occurs in wide band gap (> 1.65 eV) I/Br mixed perovskite absorbers, and becomes even more acute in the top cells of triple-junction solar photovoltaics that requires a fully 2.0 eV band gap absorber. We report herein that lattice distortion in I/Br mixed perovskites is correlated with the suppression of phase segregation, generating an increased ion migration energy barrier arising from the decreased average interatomic distance between A-site cation and iodide. Using a ~2.0 eV Rb/Cs mixed-cation inorganic perovskite with large lattice distortion in the top subcell, we fabricated all-perovskite triple-junction solar cells and achieved an efficiency of 24.3% (23.3% certified quasi-steady-state efficiency) with an open-circuit voltage of 3.21 V. This is, to our knowledge, the first reported certified efficiency for perovskite-based triple-junction solar cells. The triple-junction devices retain 80% of their initial efficiency following 420 hours of operation at the maximum power point.

14 SOLAR ENERGY↗
DOE OSTI2022

A scalable solid-state nanoporous network with atomic-level interaction design for carbon dioxide capture

Carbon capture and sequestration reduces carbon dioxide emissions and is critical in accomplishing carbon neutrality targets. Here, we demonstrate new sustainable, solid-state, polyamine-appended, cyanuric acid–stabilized melamine nanoporous networks (MNNs) via dynamic combinatorial chemistry (DCC) at the kilogram scale toward effective and high-capacity carbon dioxide capture. Polyamine-appended MNNs reaction mechanisms with carbon dioxide were elucidated with double-level DCC where two-dimensional heteronuclear chemical shift correlation nuclear magnetic resonance spectroscopy was performed to demonstrate the interatomic interactions. We distinguished ammonium carbamate pairs and a mix of ammonium carbamate and carbamic acid during carbon dioxide chemisorption. The coordination of polyamine and cyanuric acid modification endows MNNs with high adsorption capacity (1.82 millimoles per gram at 1 bar), fast adsorption time (less than 1 minute), low price, and extraordinary stability to cycling by flue gas. This work creates a general industrialization method toward carbon dioxide capture via DCC atomic-level design strategies.

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

Computational Analysis of the Energetic Stability of High-Entropy Structures of a Prototypical Lanthanide-Based Metal–Organic Framework

High-entropy materials are characterized by their complex compositions, typically comprising five or more elements in near-equiatomic proportions. Applying this concept to metal ions in metal−organic frameworks (MOFs) has paved the way for exploring a new class of high-entropy MOFs. While the compositional strategy of high-entropy materials leverages configurational entropy to aid thermodynamic stability, it also poses significant analytical challenges due to the vast compositional landscape and diverse phases that these materials can adopt. We present a computational study of several complexities associated with selecting potential high-entropy versions of a prototype lanthanidebased MOF. We compute the energetics of metal mixing of these heterometallic MOFs using density functional theory (DFT) and machine learning interatomic potential (MLIP) methods. The use of MLIP methods allows a systematic exploration of the convex hull of thermodynamically stable MOF structures containing up to 5 distinct metals.

Chemical structure↗
DOE OSTI2025

Machine-Learning-Guided Insights into Solid-Electrolyte Interphase Conductivity: Are Amorphous Lithium Fluorophosphates the Key?

Despite decades of study, the identity of the dominant Li + -conducting phase within the inorganic SEI of Li-ion batteries remains unresolved. While the mosaic model describes LiF/Li 2 O/Li 2 CO 3 nanocrystallites within a disordered matrix, these crystalline phases inherently offer limited ionic conductivity. Growing evidence suggests that interfaces, grain boundaries, and amorphous phases may instead host the primary fast-ion pathways. Using diffusion-based generative structure prediction and machine-learning interatomic potentials (MLIPs), we investigate lithium difluorophosphate (LiPO 2 F 2 ), a key mixed-anion decomposition product of phosphorus- and fluorine-containing electrolytes. We identify a stable crystalline polymorph and demonstrate that the amorphous counterpart is conductive, with projected room-temperature σ ≈ 0.18 mS cm –1 and E a ≈ 0.40 eV. Here, this enhancement stems from structural disorder flattening the Li site-energy landscape and a low formation energy for Li-interstitial defects, which supplies additional mobile carriers. We propose amorphous mixed-anion Li-P-O-F phases as a promising conducting medium in the SEI, offering a specific target for engineering improved battery interfaces.

Zhong, Peichen [University of California, Berkeley↗
NASA NTRS1985

Statistical mechanics of light elements at high pressure. VIII - Thomas-Fermi-Dirac theory for binary mixtures of H with He, C, and O

We present three-dimensional Thomas-Fermi-Dirac calculations of lattice mixing energies of hydrogen with carbon and oxygen atoms, respectively. The results are used to derive effective interatomic potentials for use in liquid-state mixture calculations. We then use the potentials to derive analytic expressions for binary mixture-free energies and to map out the phase diagrams of mixtures of hydrogen with, respectively, helium, carbon, and oxygen, over a pressure range of about 5 to about 10 to the 3rd Mbar. Within this pressure range, all three of the latter elements are found to have unlimited solubility in metallic hydrogen over a temperature range which lies above their pure-element melting temperatures, and which includes likely interior temperatures in the Jovian planets.

Hubbard, W. B.↗
NASA NTRS2011

Lattice Thermal Conductivity of Ultra High Temperature Ceramics ZrB2 and HfB2 from Atomistic Simulations

Atomistic Green-Kubo simulations are performed to evaluate the lattice thermal conductivity for single crystals of the ultra high temperature ceramics ZrB2 and HfB2 for a range of temperatures. Recently developed interatomic potentials are used for these simulations. Heat current correlation functions show rapid oscillations which can be identified with mixed metal-Boron optical phonon modes. Agreement with available experimental data is good.

Lawson, John W.↗
DOE OSTI2022

Flexible machine-learning interatomic potential for simulating structural disordering behavior of Li 7 La 3 Zr 2 O 12 solid electrolytes

Batteries based on solid-state electrolytes, including Li 7 La 3 Zr 2 O 12 (LLZO), promise improved safety and increased energy density; however, atomic disorder at grain boundaries and phase boundaries can severely deteriorate their performance. Machine-learning (ML) interatomic potentials offer a uniquely compelling solution for simulating chemical processes, rare events, and phase transitions associated with these complex interfaces by mixing high scalability with quantum-level accuracy, provided that they can be trained to properly address atomic disorder. To this end, we report the construction and validation of an ML potential that is specifically designed to simulate crystalline, disordered, and amorphous LLZO systems across a wide range of conditions. The ML model is based on a neural network algorithm and is trained using ab initio data. Performance tests prove that the developed ML potential can predict accurate structural and vibrational characteristics, elastic properties, and Li diffusivity of LLZO comparable to ab initio simulations. As a demonstration of its applicability to larger systems, we show that the potential can correctly capture grain boundary effects on diffusivity, as well as the thermal transition behavior of LLZO. Here these examples show that the ML potential enables simulations of transitions between well-defined and disordered structures with quantum-level accuracy at speeds thousands of times faster than ab initio methods.

25 ENERGY STORAGE↗
DOE OSTI2026

User Manual - HydraGNN v5.0: Distributed Implementation of Multi-Tasking Graph Neural Networks

This document serves as the user manual for HydraGNN v5.0, a scalable graph neural network (GNN) architecture for simultaneous prediction of multiple target properties using multi-task learning (MTL). This version of HydraGNN has been developed primarily to support the development, training, and deployment of predictive graph-based deep learning (DL) models for atomistic materials modeling. HydraGNN is templated over 13 message-passing policies, including invariant models (GIN, PNA, PNAPlus, GAT, MFC, CGCNN, SAGE, SchNet, DimeNet) and equivariant models (EGNN, PNAEq, PAINN, MACE), and supports distributed training via distributed data parallelism (DDP), DeepSpeed, and Fully Sharded Data Parallelism (FSDP) on leadership-class supercomputers. Although HydraGNN can be applied to problems beyond atomistic materials modeling, its current use is confined to homogeneous graphs. Additional capabilities include machine-learned interatomic potentials with energy-conserving forces, General, Powerful, and Scalable Graph Transformer (GraphGPS) global attention, periodic boundary conditions, hyperparameter optimization, mixed-precision training, and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗
DOE OSTI2023

What dictates soft clay-like lithium superionic conductor formation from rigid salts mixture

Abstract Soft clay-like Li-superionic conductors, integral to realizing all-solid-state batteries, have been recently synthesized by mixing rigid-salts. Here, through computational and experimental analysis, we clarify how a soft clay-like material can be created from a mixture of rigid-salts. Using molecular dynamics simulations with a deep learning-based interatomic potential energy model, we uncover the microscopic features responsible for soft clay-formation from ionic solid mixtures. We find that salt mixtures capable of forming molecular solid units on anion exchange, along with the slow kinetics of such reactions, are key to soft-clay formation. Molecular solid units serve as sites for shear transformation zones, and their inherent softness enables plasticity at low stress. Extended X-ray absorption fine structure spectroscopy confirms the formation of molecular solid units. A general strategy for creating soft clay-like materials from ionic solid mixtures is formulated.

25 ENERGY STORAGE↗
DOE OSTI2024

Probing degradation at solid-state battery interfaces using machine-learning interatomic potential

Solid-state batteries featuring fast ion-conducting solid electrolytes are promising next-generation energy storage technologies, yet challenges remain for practical deployment due to electro-chemo-mechanical instabilities at solid-solid interfaces. These interfaces, which include homogeneous/internal interfaces such as grain boundaries (GBs) and heterogeneous/external interfaces between solid-electrolyte and electrode materials, can impede Li-ion transport, deteriorate performance, and eventually lead to cell failure. Here, in this study, we leverage large-scale molecular simulations, enabled by validated machine-learning interatomic potentials, to directly probe the onset of interfacial degradation at the garnet Li 7 La 3 Zr 2 O 12 (LLZO) solid-electrolyte/LiCoO 2 (LCO) cathode interface. By surveying different interfacial geometries and compositions, it is found that Li-deficient interfaces can lead to severe interfacial disordering with cation mixing and Co interdiffusion from LCO into LLZO. By contrast, Li-sufficient interfaces are less disordered, although elemental segregation with local ordering is observed. As a consequence of Co interdiffusion, Co-rich regions are formed at the GBs of LLZO due to cation segregation and trapping effects. This behavior is independent of the GB tilting axis, degree of disorder at the GBs, and Co concentration, which implies Co clustering at GBs is a general phenomenon in polycrystalline LLZO and can dictate its overall transport and mechanical properties. Our findings elucidate the underlying fundamental mechanisms that give rise to experimentally observed physicochemical properties and provide guidelines for interface design that can mitigate interfacial degradation and improve cycling performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Plutonium oxide melt structure and covalency

Advances in nuclear power reactors include the use of mixed oxide fuel, containing uranium and plutonium oxides. The high-temperature behaviour and structure of PuO 2-x above 1,800 K remain largely unexplored, and these conditions must be considered for reactor design and planning for the mitigation of severe accidents. Here, in this study, we measure the atomic structure of PuO 2-x through the melting transition up to 3,000 +/- 50 K using X-ray scattering of aerodynamically levitated and laser-beam-heated samples, with O/Pu ranging from 1.57 to 1.76. Liquid structural models consistent with the X-ray data are developed using machine-learned interatomic potentials and density functional theory. Molten PuO 1.76 contains some degree of covalent Pu-O bonding, signalled by the degeneracy of Pu 5f and O 2p orbitals. The liquid is isomorphous with molten CeO 1.75 , demonstrating the latter as a non-radioactive, non-toxic, structural surrogate when differences in the oxidation potentials of Pu and Ce are accounted for. These characterizations provide essential constraints for modelling pertinent to reactor safety design. The molten structure of plutonium oxide-a component of mixed oxide nuclear fuels-is measured, showing some degree of covalent bonding. Its atomic structure is similar to that of cerium oxide, which could be a non-radioactive structural surrogate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Design of Sodium Chalcohalide Solid Electrolytes with Mixed Anions for All‐Solid‐State Sodium‐Ion Batteries

Solid-state sodium-ion batteries (SSNIBs) have emerged as a promising alternative to lithium-ion systems for grid-scale energy storage, owing to sodium's abundance and the improved safety of solid-state designs. Among various solid-state electrolytes (SSEs), halide-based Na + SSEs offer high electrochemical stability but are limited by low ionic conductivity and poor thermal stability. Herein, a novel class of sodium hafnium chalcohalide SSEs is reported with a dual-anion (S 2− /Cl − ) framework, with a high ionic conductivity of 4.5 × 10 −4 S cm −1 . The incorporation of sulfur enhances Na⁺ mobility by reducing the migration barrier through increased anion polarizability and expanded diffusion pathways. Additionally, S 2− contributes to stronger interatomic bonding, leading to higher cohesive energy density, improved thermal stability, and mechanical robustness. These SSEs exhibit minimal sulfur oxidation and excellent chemical/electrochemical interface stability with different cathode materials, such as O3-layered NaNi 1/3 Fe 1/3 Mn 1/3 O 2 , P2/O3 layered Na 0.85 Mn 0.5 Ni 0.4 Fe 0.1 O 2 , and Na 3 V 2 (PO 4 ) 3 cathodes. As a result, SSNIBs with P2/O3 layered Na 0.85 Mn 0.5 Ni 0.4 Fe 0.1 O 2 employing the sodium hafnium chalcohalide SSEs demonstrate outstanding cycling performance, achieving a capacity retention of 88.5% after 200 cycles at 0.1 C. This study establishes a new design strategy for high-performance SSEs, demonstrating that mixed-anion frameworks offer a viable route to overcome the intrinsic limitations of single-anion electrolytes in next-generation SSNIBs.

DFT calculation↗