SEARCH · Engineering Papers

Results for “Atomic-scale 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗

DOE OSTI2025

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Condensed matter imaging↗

DOE OSTI2025

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗

DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Atomically Revealing Bulk Point Defect Dynamics in Hydrogen‐Driven γ‐Fe 2 O 3 → Fe 3 O 4 → FeO Transformation

Understanding how point defects in the bulk govern redox transformations is essential for advancing hydrogen-based metal production and designing high-performance oxide materials. This study reveals the atomic-scale mechanisms driving hydrogen-induced reduction of γ-Fe 2 O 3 to Fe 3 O 4 , focusing on how bulk vacancy dynamics dictate structural evolution and reaction kinetics. A key finding is the pronounced contrast in defect behavior between the two oxides: in γ-Fe 2 O 3 , intrinsic Fe vacancies promote oxygen vacancy clustering, destabilizing the local lattice and driving nanopore formation. In contrast, Fe 3 O 4 exhibits a higher oxygen vacancy formation energy and lacks intrinsic Fe vacancies, suppressing vacancy aggregation and maintaining a dense, pore-free structure. This divergence governs distinct reduction pathways—γ-Fe 2 O 3 undergoes an interface-reaction-limited transformation confined to the γ-Fe 2 O 3 /Fe 3 O 4 boundary, while Fe 3 O 4 supports a uniform increase in oxygen vacancy concentration, enabling bulk-phase reduction to lower-oxide FeO. Integrated in situ electron microscopy and density functional theory modeling uncover a vacancy-mediated mechanism, where synergistic cation-anion vacancy dynamics steer microstructure evolution and phase progression. These insights highlight the critical role of vacancy dynamics in controlling oxide reactivity and offer a pathway toward vacancy engineering to enhance reduction kinetics in hydrogen metallurgy and to tailor porosity, reactivity, and structural resilience in oxide-based catalysts and energy materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Acidities of MgO surface sites: implications for the formation mechanism of Mg(OH) 2

The hydroxylation of periclase (MgO) to brucite (Mg(OH) 2 ) is thought to be an important intermediate step when using MgO to capture CO 2 from the atmosphere. However, the mechanism of hydroxylation of MgO to form Mg(OH) 2 is poorly understood. In this work, we used atomic-scale density functional tight binding simulations coupled with the metadynamics rare event method to analyze the surface chemistry of MgO and the acid dissociation equilibrium constants (pK a ) of its surface sites. The method and parameters were validated by calculating the pK a for hydroxylation of the first shell water bound to aqueous Mg 2+ ion. The pK a value derived using a probabilistic method was 12.3, which is in fair agreement with the accepted value of 11.4, with the difference between them equal to a ∼5 kJ mol −1 error in the calculations. We then extended these pK a calculations to probe the hydroxylation reactions of the surface sites of the MgO(100)–water interface, arriving at pK a s of 5.4 to deprotonate terminal water molecules bound to the surface magnesium sites (η-OH 2 or 〉MgOH 2 ), and 13.9 to deprotonate hydroxylated bridging oxygen sites (μ 5 -oxo or 〉O). Hydroxide (OH − ) adsorption on the surface was also probed and found to be less thermodynamically favorable than deprotonation of the terminal water molecule. The plausibility of the computed pK a s was verified using an activity-based speciation model and compared to pH measurements of water equilibrated with MgO nanoparticles and single crystals. The model predicted a solution pH of 7.1 when surface sites buffered and the pH of 12.0 when MgO dissolution dominated. These are close to the experimental initial solution pHs of 7–7.5 and the long term pHs of ∼10.5. The similarity suggests that the calculated pK a values from the DFTB+/metadynamics simulations are plausible and that these methods can be a useful tool to probe reaction mechanisms involving covalent bonds.

Adapa, Sai Krishna Reddy [Oak Ridge National Labor↗

DOE OSTI2026

Theory Guided Fine‐Tune of Strain Effects in Pt Ternary Alloy via Rare Earth Templating: Achieving High Performance PEMFCs Catalysts

The sluggish kinetics and insufficient durability of platinum-based catalysts remain crucial barriers limiting proton-exchange-membrane fuel cells (PEMFCs) deployment. Here, we report a theory-guided synthesis combined with rare-earth templating to realize a previously inaccessible Pt 5 Co-like phase with tailored atomic-scale strain. Guided by density functional theory (DFT) calculations, we identified that a Pt 5 Co-like sublayer can induce a unique mild compressive strain (−1.24%) to the Pt(111) shell and an optimal *OH binding energy shift (ΔE ≈ 0.11 eV). This shift positions the alloy catalyst near the apex of the oxygen reduction reaction activity volcano. This prediction guided the synthesis of ternary alloy Pt 5 (Ce)Co@Pt multilayer nanoparticles, featuring a Ce-stabilized core, a Pt 5 Co-like sublayer, and a Pt-rich shell. This catalyst demonstrates both exceptionally high activity and durability, achieving a mass activity of 2.6 A∙mg Pt −1 in rotating disk electrode testing. In fuel cell membrane electrode assembly tests, Pt 5 (Ce)Co@Pt achieves a current density of 1.9 A∙cm −2 at 0.7 V under heavy-duty vehicle conditions. Remarkably, it maintains 1.2 A∙cm −2 after 1 80 000 AST cycles, doubling the U.S. DOE 2025 target. This work demonstrates a rational design strategy that DFT-guided strain engineering integrates with rare-earth templating to advance Pt-based catalysts for fuel cell applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

The Influence of Residual Copper Catalyst on the Polymer-to-Ceramic Conversion of a Click-Derived Polycarbosilane

Preceramic polymers (PCPs) are a key class of materials for the preparation of advanced ceramic components. Following conventional polymer processing and subsequent pyrolysis, polymer-derived ceramic (PDC) fibers, monoliths, and composites offer novel engineering opportunities for aerospace applications. Despite the rapidly emerging commercial applications of PDCs, many basic scientific questions remain regarding the structural transformation of the polymer to the final ceramic and how processing variables affect PDC materials properties. PCPs prepared using a mild and scalable copper (Cu) catalyzed azide/alkyne “click” reaction have been demonstrated as a unique class of material that can bind with transition metals, ultimately producing ultrahigh temperature ceramics (UHTCs) upon pyrolysis. Unexpectedly, it was found that the removal of residual Cu catalyst from a model click-derived polycarbosilane PCP significantly reduced the ceramic yield during pyrolysis. Herein, we report the influence of Cu impurities on polycarbosilane-derived SiC ceramics and present a mechanistic understanding of the unintended role of Cu in the polymer-to-ceramic conversion process. Synchrotron-based scattering and spectroscopy methods were implemented to track the evolution of different atomic species and local structures as well as the influence of Cu toward atomic-scale structures in the final ceramics. It was determined that Cu impurities catalyze the formation of silicon carbide through the formation of an intermediate molten copper/copper silicide, ultimately increasing the ceramic yield. Overall, this contribution highlights the need to carefully assess the interplay between the chemistries needed for PCP synthesis and their potential influence on final PDCs following processing.

Ponder, Jr., James F. [Air Force Research Laborato↗

DOE OSTI2025

Nanopolysaccharide Builder: A User-Friendly Tool for Atomistic Models of Polysaccharide-Based Nanostructures

Here, we introduce Nanopolysaccharide Builder (NPB), a user-friendly software tool designed to construct polysaccharide nanostructures─mainly those based on cellulose, chitin, and chitosan─using experimental data or user-defined parameters. NPB enables the generation of cellulose and chitin allomorphs with customizable biochemical topologies and also facilitates the construction of large bundles that replicate nanostructures found in biological support systems, including plant cell walls and arthropod cuticles. The software outputs atomic Cartesian coordinates in Protein Data Bank (PDB) format and also provides atom connectivity files in PSF and PARM formats, ensuring seamless integration with major molecular dynamics (MD) engines such as NAMD, CHARMM, GROMACS, AMBER, OpenMM, and LAMMPS. Built on an interactive visualization framework, NPB features a graphical user interface (GUI) and supports both macOS and Linux operating systems. By enabling detailed atomic-scale studies of polysaccharide evolution in extracellular matrices and cell walls of algae, bacteria, fungi, and plants, NPB is poised to advance AI-guided research in sustainable chemical development and biomass utilization.

Wan, Zhangmin [Univ. of British Columbia, Vancouve↗

DOE OSTI2025

Dynamic Surface Restructuring in Cu(Au) Alloys Driven by Oxygen-Mediated Au Mobility

Alloying plays a crucial role in tuning the surface properties of metals, but the atomic-level mechanisms by which alloying elements influence surface structure dynamics under reactive conditions remain elusive. Using Cu(Au) in oxidizing environments as a model system, we reveal a dynamic oxygen-induced transformation of the topmost atomic layer into a periodically hill-and-valley morphology with reversible switching between undulated and flattened surface states. These interconversions are driven by the retreat of surface Au to the subsurface during oxygen adsorption and its resegregation to the surface upon oxygen desorption. This cyclical mobility establishes a feedback loop, allowing the surface to dynamically reconfigure in response to changes in the oxygen pressure. The results offer a broadly applicable framework for understanding atomic-scale surface restructuring in alloy systems, where differences in the chemical reactivity of alloying elements drive dynamic redistribution between surface and subsurface regions. As a result, this dynamic coupling has practical implications for designing corrosion-resistant coatings and metastable nanostructures with tunable catalytic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Atom Identification in Bilayer Moiré Materials with Gomb-Net

Moiré patterns in van der Waals bilayer materials complicate the analysis of atomic-resolution images, hindering the atomic-scale insight typically attainable with scanning transmission electron microscopy. Here, we report a method to detect the positions and identities of atoms in each of the individual layers that compose twisted bilayer heterostructures. We developed a deep learning model, Gomb-Net, which identifies the coordinates and atomic species in each layer, deconvoluting the moiré pattern. This enables layer-specific mapping of atomic positions and dopant distributions, unlike other commonly used segmentation models which struggle with moiré-induced complexity. Using this approach, we explored the Se atom substitutional site distribution in a twisted fractional Janus WS 2 -WS 2(1–x) Se 2x heterostructure and found that layer-specific implantation sites are unaffected by the moiré pattern’s local energetic or electronic modulation. In conclusion, this advancement enables atom identification within material regimes where it was not possible before, opening new insights into previously inaccessible material physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Insights into Native Single-Atom Electrocatalyst Site Structures

Single-atom electrocatalysts consisting of metal atoms embedded in a carbon matrix are promising next-generation catalysts for green hydrogen production and utilization, CO2 reduction, low-temperature CO oxidation, ammonia production, plastic decomposition, and electrochemical energy storage. The origins of activity and stability for the single-atom sites are still debatable, however, because of constrained insights into their local structure resulting from idealized models and experiments derived from a large number of individual sites. Insights into structural variations around single atomic sites are therefore critical for the continued development of these next-generation catalysts. While electron microscopy commonly provides atomic-scale information about these materials, the beam sensitivity of individual sites makes structural determination by conventional low-voltage (60 keV) techniques challenging. Here, we introduce ultralow-voltage electron ptychography, performed at 30 keV, that enables determination of the lattice structure around individual metal sites in a well-defined single-atom electrocatalyst system while essentially eliminating knock-on structural modifications. Pairing these atomic-scale, site-specific measurements with computational methods will broaden our understanding of the activity and stability of these materials, which will accelerate the development of the next generation of catalysts.

Zachman, Michael [ORNL] (ORCID:0000000319101357)↗

DOE OSTI2026

Autonomous fabrication of tailored defect structures in 2D materials using machine learning-enabled scanning transmission electron microscopy

Materials with tailored quantum properties can be engineered from atomic-scale assembly techniques, but existing methods often lack the agility and accuracy to precisely and intelligently control the manufacturing process. Here, we demonstrate a fully autonomous approach for fabricating atomic-level defects using electron beams in scanning transmission electron microscopy (STEM) that combines advanced machine learning and automated beam control. As a proof of concept, we achieved controlled fabrication of MoS-nanowire (MoS-NW) edge structures by iterative and targeted exposure of MoS 2 monolayer to a focused electron beam to selectively eject sulfur atoms, utilizing high-angle annular dark-field (HAADF) imaging for feedback-controlled monitoring of structural evolution of defects. A machine learning framework combining a random forest model and a convolutional neural network (CNN) was developed to decode the HAADF image and accurately identify atomic positions and species. This atomic-level information was then integrated into an autonomous decision-making platform, which applied predefined fabrication strategies to instruct beam control about atomic sites to be ejected. The selected sites were subsequently exposed to a localized electron beam using an FPGA-controlled scan routine with precise control over beam positioning and duration. While the MoS-NW edge structures produced exhibit promising mechanical and electronic properties, the proposed methods to build the autonomous fabrication framework is material-agnostic and can be extended to other 2D materials for the creation of diverse defect structures and heterostructures beyond Mo S2 .

Engineering↗

DOE OSTI2025

Atomic-scale imaging of frequency-dependent phonon anisotropy

Directly visualizing vibrational anisotropy in individual phonon modes is essential for understanding a wide range of intriguing optical, thermal and elastic phenomena in materials. Although conventional optical and diffraction techniques have been used to estimate vibrational anisotropies, they fall short in achieving the spatial and energy resolution necessary to provide detailed information. Here, we introduce a new form of momentum-selective electron energy-loss spectroscopy, which enables the element-resolved imaging of frequency- and symmetry-dependent vibrational anisotropies with atomic resolution. Vibrational anisotropies manifest in different norms of orthogonal atomic displacements, known as thermal ellipsoids. Using the centrosymmetric strontium titanate as a model system, we observed two distinct types of oxygen vibrations with contrasting anisotropies: oblate thermal ellipsoids below 60 meV and prolate ones above 60 meV. In non-centrosymmetric barium titanate, our approach can detect subtle distortions of the oxygen octahedra by observing the unexpected modulation of q-selective signals between apical and equatorial oxygen sites near 55 meV, which originates from reduced crystal symmetry and may also be linked to ferroelectric polarization. These observations are quantitatively supported by theoretical modelling, which demonstrates the reliability of our approach. The measured frequency-dependent vibrational anisotropies shed new light on the dielectric and thermal behaviours governed by acoustic and optical phonons. Furthermore, the ability to visualize phonon eigenvectors at specific crystallographic sites with unprecedented spatial and energy resolution opens new avenues for exploring dielectric, optical, thermal and superconducting properties.

Ferroelectrics and multiferroics↗

DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Atomic-Scale Tracking of Topological Defect Motion and Incommensurate Charge Order Melting

Charge order pervades the phase diagrams of quantum materials where it competes with superconducting and magnetic phases, hosts electronic phase transitions and topological defects, and couples to the lattice generating intricate structural distortions. Incommensurate charge order is readily stabilized in manganese oxides, where it is associated with anomalous electronic and magnetic properties, but its nanoscale structural inhomogeneity complicates precise characterization and understanding of its relationship with competing phases. Leveraging atomic-resolution variable-temperature cryogenic scanning transmission electron microscopy, we characterize the thermal evolution of charge order as it transforms from its ground state in a model manganite system. We find that mobile networks of discommensurations and dislocations generate phase inhomogeneity and induce global incommensurability in an otherwise lattice-locked modulation. Driving the order to melt at high temperatures, the discommensuration density grows, and regions of order locally decouple from the lattice periodicity. Published by the American Physical Society 2025

Schnitzer, Noah (ORCID:0000000210421146)↗

DOE OSTI2025

Identification of short-range ordering motifs in semiconductors

Chemical short-range ordering is expected to be a key factor for tuning the electronic structure of semiconductors. However, experimental evidence of short-range ordering is still lacking due to the challenge of characterizing atomic-scale ordering motifs. Here, we determined the presence of short-range order in a ternary GeSiSn semiconductor system using advanced energy-filtered four-dimensional scanning transmission electron microscopy and large-scale atomistic models generated by a machine learning neuroevolution potential of first-principles accuracy. This approach revealed preferred ordering of different atomic species with the dominant occurrence of Si–Ge–Sn triplets. Our findings not only confirmed the presence of short-range order but also directly revealed the actual atomic structure, demonstrating the potential for informed atomic order–based band engineering as a third degree of freedom beyond composition and strain tuning.

Vogl, Lilian M. [University of California Berkeley↗