Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic Materials Modeling”

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.

1,810 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The impact of chemistry on anion migration in bixbyite-structured lanthanide oxides

This manuscript describes atomistic calculations of oxygen vacancy and interstitial migration in bixbyite structured lanthanide oxides. We examine two types of compounds, one in which only one type of lanthanide cation is present and a second class in which two lanthanides are present in a 3:1 ratio as dictated by the symmetry of the bixbyite crystal structure. Using temperature accelerated dynamics and the nudged elastic band method, we quantify the role of chemistry on the energy barriers for the most important pathways for both vacancy and interstitial migration. We then analyze the impact of these variations on the overall diffusivity of each defect, quantifying the contribution of each pathway using the theory of kinosons. We find that vacancy mobility can vary by as much as three orders of magnitude through changes in chemistry at 500 K. Changes in interstitial mobility are more modest but can still vary by an order of magnitude. This points to the ability to tune the mass transport characteristics of these compounds through appropriate choices in chemistry. We have also included supplementary information containing the atomic structures of the relevant pathways.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗
DOE OSTI2024

Tailoring Ion Transport in Li 3‐3y Ho 1+y Cl 6‐x Br x via Transition‐Metal Free Structural Planes and Charge Carrier Distribution

Abstract Localized atomistic disorder in halide‐based solid electrolytes (SEs) can be leveraged to boost Li + mobility. In this study, Li + transport in structurally modified Li 3 HoCl 6 , via Br − introduction and Li + deficiency, is explored. The optimized Li 3‐3 y Ho 1+ y Cl 6‐ x Br x achieves an ionic conductivity of 3.8 mS cm −1 at 25 °C, the highest reported for holmium halide materials. 6,7 Li nuclear magnetic resonance and relaxometry investigations unveil enhanced ion dynamics with bromination, attaining a Li + motional rate neighboring 116 MHz. X‐ray diffraction analyses reveal mixed‐anion‐induced phase transitions with disproportionate octahedral expansions and distortions, creating Ho‐free planes with favorable energetics for Li + migration. Bond valence site energy analysis highlights preferred Li + transport pathways, particularly in structural planes devoid of Ho 3+ blocking effects. Molecular dynamics simulations corroborate enhanced Li + diffusion with Br − introduction into Li 3 HoCl 6 . Li‐Ho electrostatic repulsions in the (001) plane presumably drive Li + diffusion into the Ho‐free (002) layer, enabling rapid intraplanar Li + motion and exchange between the 2d and 4h sites. Li 3‐3 y Ho 1+ y Cl 6‐ x Br x also demonstrates good battery cycling stability. These findings offer valuable insights into the intricate correlations between structure and ion transport and will help guide the design of high‐performance fast ion conductors for all‐solid‐state batteries.

Chemistry↗
DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

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

Direct Observation of Elusive (DTBM‐SEGPHOS)CuH Monomer Enables Mechanistic Insights Into Hydrocupration, Aggregation, and Dynamics of Alkene Functionalization Catalysis

The bulky diphosphine DTBM-SEGPHOS is widely employed in CuH-catalyzed transformations as it provides remarkably active catalyst systems. The transient (DTBM-SEGPHOS)CuH monomer (LCuH) is the often-invoked active species. However, its instability has prevented spectroscopic characterization and mechanistic elucidation, hindering mechanistic understanding. We report low-temperature NMR spectroscopic characterization of LCuH, enabling quantitative kinetic analysis of the stoichiometric hydrocupration and catalytic hydroboration of cyclopentene, as well as the structural identification of two CuH clusters. LCuH inserts cyclopentene at −43°C, reaffirming its high reactivity toward olefins. LCuH deactivates to form L 2 Cu 3 H 3 and L 2 Cu 4 H 4 clusters, in which LCuH dimerization initiates aggregation. Kinetic analysis of reactions of unactivated alkenes indicates that competing on-cycle alkene hydrocupration and LCuH dimerization impact performance, as catalyst deactivation and turnover occur on comparable timescales. Structure–activity analysis using atomistic simulations shows that the steric profile of DTBM-SEGPHOS increases the CuH dimerization barrier by ∼7.7 kcal mol−1 compared to that of SEGPHOS, rationalizing the unique ability of DTBM-SEGPHOS to stabilize a reactive monomer for hydrocupration of broader alkene substrates. These findings illustrate the fundamental design principle that steric control of aggregation governs CuH catalyst performance, explaining both the exceptional activity of (DTBM-SEGPHOS)CuH and the limitations imposed by competing deactivation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Extending mARC II Arc-Jet Test Duration via Design and Implementation of Model System Cooling Sleeve

The mARC II is a 30 kW arc-jet facility at NASA Ames Research Center used to generate high-enthalpy flows for low-cost thermal protection system (TPS) technology development. Sustained operation of downstream instrumentation and material samples is constrained by thermal loading transmitted through the arc-jet test environment, limiting achievable run times and experimental throughput. This work presents the design, integration, and validation of a cooling sleeve implemented on the sweep arm drive motor feedthrough to mitigate thermal accumulation during testing. The addition of the cooling sleeve is a simple, robust upgrade that translates directly into enhanced facility capability by supporting longer run durations, reduced turnaround time, and higher throughput.

numerical simulations↗
NASA NTRSDate not supplied

Extending mARC II Arc-Jet Test Duration via Design and Implementation of Model System Cooling Sleeve

The mARC II is a 30 kW arc-jet facility at NASA Ames Research Center used to generate high-enthalpy flows for low-cost thermal protection system (TPS) technology development. Sustained operation of downstream instrumentation and material samples is constrained by thermal loading transmitted through the arc-jet test environment, limiting achievable run times and experimental throughput. This work presents the design, integration, and validation of a cooling sleeve implemented on the sweep arm drive motor feedthrough to mitigate thermal accumulation during testing. The addition of the cooling sleeve is a simple, robust upgrade that translates directly into enhanced facility capability by supporting longer run durations, reduced turnaround time, and higher throughput.

sting arm↗