Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic dynamics”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Facet-Dependent Cold Welding of Au Nanorods Revealed by Liquid Cell Transmission Electron Microscopy

Cold welding of metals at the nanoscale has been demonstrated to play a significant role in bottom-up manufacturing and self-healing processes of nanostructures and nanodevices. However, the welding mechanism at the nanoscale is not well understood. In this study, a comprehensive demonstration of the cold welding process of gold nanorods with different modes is presented through in situ liquid cell transmission electron microscopy. The experimental results and molecular dynamics simulations reveal that the nanorods are welded through the facet-dependent atomic surface diffusion and rearrangement along {100} facets. The density functional theory calculations indicate that the preferred coalescence of two {100} surfaces is thermodynamically favorable. Unlike the prevalent “oriented attachment” in the nanoparticle coalescence, the misalignment of nanorod orientations and local stresses can induce grain boundaries and stacking faults in the welded interface.

42 ENGINEERING↗
DOE OSTI2024

Self-Inhibition Phenomena in Cu 3 Pt Oxidation by CO 2

Here, this study investigates the oxidation behavior of Cu 3 Pt(100) in CO 2 using a combination of ambient-pressure X-ray photoelectron spectroscopy, mass spectroscopy, and density functional theory modeling. Our in situ measurements reveal the simultaneous oxidation and reduction of Cu 2 O due to the opposing effects of atomic oxygen and CO generated from dissociative CO 2 adsorption, leading to a dynamic equilibrium state of simultaneously occurring redox reactions. Complementary atomistic calculations elucidate the inhibitory effects of subsurface Pt enrichment and the counteracting roles of CO 2 and CO in surface oxidation and reduction. These results provide mechanistic insights into the dissociative pathway of CO 2 molecules and dynamic evolution of surface composition and reactivity of Cu-based alloy catalysts in CO 2 -rich environments, with broader implications for tuning gas–surface reactions by manipulating gas reactants or solid surface composition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Activation of Lattice and Adatom Oxygen by Highly Stable Ceria-Supported Cu Single Atoms

Requiring catalysts to be both active yet stable over long periods of time under variable reaction conditions including high and low temperatures is a daunting challenge due to the almost mutual exclusivity of these constraints. Using CO oxidation as a probe reaction, in this work we demonstrate that thermally stable single atom copper catalysts prepared by high temperature synthesis (atom trapping) on ceria can achieve this feat by allowing modulation of the Cu charge state through facile charge transfer between active site and support. This provides the catalysts with an ability to activate either lattice or adatom oxygen atoms, accessing additional reaction channels depending on the reaction conditions. Such adaptability allows dynamic response of catalysts enabling them to remain active under variable reaction conditions. The inherent stability of the catalyst arises from the enhanced strength of the Cu-O interactions from high temperature synthesis that exist even when Cu oxidation state varies, which retards sintering and deactivation. As we show here, one can circumvent the dilemma of designing catalysts that are simultaneously active and stable by matching the redox properties of the active site and support and establishing an environmental adaptability around the active sites.

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

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗
DOE OSTI2025

Atomically resolved edges and defects in lead halide perovskites

Although edges and defects constitute only a small fraction of crystalline materials, they exert an outsized impact on a material′s properties. Organic–inorganic halide perovskites are promising next-generation semiconductor materials with superior cost effectiveness and interesting optoelectronic properties. However, clear images of their edges have remained challenging to obtain owing to their extreme sensitivity. Using truly high-speed ultralow-dose four-dimensional scanning transmission electron microscopy with dose fractionation, we perform ptychography at, to our knowledge, the lowest-dose atomic resolution to date, revealing not only the detailed atomic structure of the edges of a halide perovskite but also their structural dynamics. Here, a majority methylammonium (MA) and iodine (I) edge termination is observed in methylammonium lead iodide (MAPbI 3 ), and the damage rate of its edges and internal defects is found to depend on the concentration and type of vacancies present, with a preponderance of I vacancies in particular correlating with higher rates of damage.

4D-STEM↗
DOE OSTI2024

Cu site differentiation in tetracopper(I) sulfide clusters enables biomimetic N 2 O reduction

Copper clusters feature prominently in both metalloenzymes and synthetic nanoclusters that mediate catalytic redox transformations of gaseous small molecules. Such reactions are critical to biological energy conversion and are expected to be crucial parts of renewable energy economies. However, the precise roles of individual metal atoms within clusters are difficult to elucidate, particularly for cluster systems that are dynamic under operating conditions. Here, we present a metal site-specific analysis of synthetic Cu 4 (μ 4 -S) clusters that mimic the Cu Z active site of the nitrous oxide reductase enzyme. Leveraging the ability to obtain structural snapshots of both inactive and active forms of the synthetic model system, we analyzed both states using resonant X-ray diffraction anomalous fine structure (DAFS), a technique that enables X-ray absorption profiles of individual metal sites within a cluster to be extracted independently. Using DAFS, we found that a change in cluster geometry between the inactive and active states is correlated to Cu site differentiation that is presumably required for efficient activation of N 2 O gas. More precisely, we hypothesize that the Cu δ+ ∙∙∙Cu δ- pairs produced upon site differentiation are poised for N 2 O activation, as supported by computational modeling. These results provide an unprecedented level of detail on the roles of individual metal sites within the synthetic cluster system and how those roles interplay with cluster geometry to impact the reactivity function. We expect this fundamental knowledge to inform understanding of metal clusters in settings ranging from (bio)molecular to nanocluster to extended solid systems involved in energy conversion.

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

Nonlinear X-ray Optics

Novel X-ray free-electron lasers (XFELs) are capable of producing radiation with unprecedented properties. Specifically, their ultrashort coherent X-ray pulses can generate extreme peak intensities where ordinary rules of light-matter interaction may no longer apply, and nonlinear processes start to become important. XFELs have made it possible to observe some of the most important coherent non-sequential nonlinearities, including X-ray-optical sum frequency generation (SFG), X-ray second harmonic generation (XSHG), two-photon absorption (TPA) and nonlinear two-photon Compton scattering (2PCS). In particular, the 2PCS experiment has led to unexpected results, namely the observation of a substantial anomalous red shift in the energy of photons generated by nonlinear two-photon X-ray Compton scattering in beryllium. The energy shift is in addition to the predicted nonlinear Compton shift and was neither expected from extrapolations from linear X-ray interactions nor from nonlinear effects at optical wavelengths. Our results suggest a novel nonlinear scattering mechanism that can only be observed at high fields at X-ray wavelengths, where the photon energy is a significant fraction of the electron rest mass (and thus a significant momentum transfer to the electron occurs during the scattering) and where the interaction is not dominated by the quiver motion of the electron in the field (ponderomotive potential) as is usually the case at longer (optical) wavelengths. The goal of this research project was the investigation of fundamental nonlinear X-ray-matter interactions. As these processes are typically weak, part of this research is the exploration of more efficient nonlinear X-ray processes by studying novel mechanisms including solid-state effects, the assistance of optical radiation and novel phase-matching techniques. These experiments require extremely large X-ray field strengths, such as produced by XFELs. Nonlinear X-ray effects can lead to important applications in many research fields. For example, it could lead to instantaneous plasma diagnostics for materials in extreme conditions, as a method for combining atomic-scale structural sensitivity with chemical specificity, or as a probe of electron dynamics in solids on the attosecond (1 as = 10 -18 s) time- and Ångstrom (10 -10 m) length scales. In addition, a detailed understanding of nonlinear X-ray matter interactions is crucial for understanding and extending the limits of structural determination at high intensity such as in serial femto-crystallography. The results could have a profound impact on advanced XFEL light sources such as the LCLS-II.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Electron Correlation and High-Temperature Superconductivity

Strong electron correlation plays a central role in the high-temperature superconductivity (HTSC) of cuprates. However, to date, research has focused only on its role in spin dynamics and related effects, even though it is becoming increasingly clear that spin alone may not be sufficient to create HTSC. Here, we discuss a possible role of electron correlation in the Bose–Einstein condensation (BEC) of Cooper pairs. Recently, we succeeded in observing dynamic electron correlation via inelastic X-ray scattering through results presented in real space. We discovered that electron correlations are strongly modified in the plasmon, proving that electron dynamics significantly affect electron correlation. Earlier, we found that in 4 He, the atom–atom distance in the BE condensate is 10% longer than that in the non-condensate. These results suggest the possibility that the reduction in electron-repulsion energy upon BEC is driving T c to high values. Thus, electron correlation itself could be the origin of the HTSC phenomenon.

Bose–Einstein condensation↗
DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Breaking the mold: Overcoming the time constraints of molecular dynamics on general-purpose hardware

The evolution of molecular dynamics (MD) simulations has been intimately linked to that of computing hardware. For decades following the creation of MD, simulations have improved with computing power along the three principal dimensions of accuracy, atom count (spatial scale), and duration (temporal scale). Since the mid-2000s, computer platforms have, however, failed to provide strong scaling for MD, as scale-out central processing unit (CPU) and graphics processing unit (GPU) platforms that provide substantial increases to spatial scale do not lead to proportional increases in temporal scale. Important scientific problems therefore remained inaccessible to direct simulation, prompting the development of increasingly sophisticated algorithms that present significant complexity, accuracy, and efficiency challenges. While bespoke MD-only hardware solutions have provided a path to longer timescales for specific physical systems, their impact on the broader community has been mitigated by their limited adaptability to new methods and potentials. In this work, we show that a novel computing architecture, the Cerebras wafer scale engine, completely alters the scaling path by delivering unprecedentedly high simulation rates up to 1.144 M steps/s for 200 000 atoms whose interactions are described by an embedded atom method potential. This enables direct simulations of the evolution of materials using general-purpose programmable hardware over millisecond timescales, dramatically increasing the space of direct MD simulations that can be carried out. In this paper, we provide an overview of advances in MD over the last 60 years and present our recent result in the context of historical MD performance trends.

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

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2021

Emergent Ergodicity at the Transition between Many-Body Localized Phases

Strongly disordered systems in the many-body localized (MBL) phase can exhibit ground state order in highly excited eigenstates. The interplay between localization, symmetry, and topology has led to the characterization of a broad landscape of MBL phases ranging from spin glasses and time crystals to symmetry protected topological phases. Understanding the nature of phase transitions between these different forms of eigenstate order remains an essential open question. In this work, we conjecture that no direct transition between distinct MBL orders can occur in one dimension; rather, an ergodic phase always intervenes. Motivated by recent advances in Rydberg-atom-based quantum simulation, we propose an experimental protocol where the intervening ergodic phase can be diagnosed via the dynamics of local observables.

1-dimensional spin chains↗
DOE OSTI2022

Elucidating the Atomic Structures of the Gel Layer Formed during Aluminoborosilicate Glass Dissolution: An Integrated Experimental and Simulation Study

The altered glasses produced during aqueous dissolution of silicate and borosilicate glasses are among the most complex structures to understand at the atomic level due to their amorphous nature, random porosity and various levels of hydration. In this study, we gained insights of the complex atomic structures of altered aluminoborosilicate glasses by combining a range of experimental and computational approaches. The altered glasses were prepared by dissolution of three glasses with varying level of alumina in acid for 7 days. A comprehensive set of experimental [elemental analysis, high-energy X-ray diffraction, 29 Si and 27 Al solid-state nuclear magnetic resonance (NMR), and O 1s X-ray photoelectron spectroscopy (XPS)] and modeling (molecular dynamics (MD) simulations using non-reactive and reactive force fields) were used to study the atomic structures of these altered glasses. Elemental analysis showed that most of the B in the pristine glasses was leached into the solution and was not contained in the altered glass. The 29 Si and 27 Al solid-state NMR spectra revealed that the altered glasses have more polymerized silicate network as compared to those in the pristine glasses due to the reformation of linkages among Si and Al oxygen polyhedral in the altered glasses. The bridging and non-bridging (or hydroxyl O) atoms in the altered glasses were also quantified from their O 1s XPS spectra. Atomic structure models of the altered glasses were constructed by using MD simulations using the reactive force field based on the compositional information obtained from experiments. Various pore structures were generated using the charge-scaling (CS) method by using different initial densities and CS temperatures; the best CS parameters of each alternated glass were then determined by comparing with the experimental structure factors obtained from high energy X-ray diffraction. Pore and the atomic structures and vibrational properties around these pore surfaces were analyzed. Furthermore, these results from this comprehensive study thus provides a realistic insight into the pore morphology, atomic structure, and vibrational properties of altered glasses.

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

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Resonant interaction between phonons and PbTe/PbSe (001) misfit dislocation networks

Here this work aims at a quantitative and mechanistic understanding of the dynamic process of the phonon-dislocation interaction in PbTe/PbSe (001) heterostructures using the Concurrent Atomistic-Continuum (CAC) method as the simulation tool. The misfit dislocation network and the atomic-scale dislocation core structure obtained in the simulations are found to agree reasonably well with the experimental observations of the PbTe/PbSe (001) interface. Through visualizing the dynamic interaction between phonons and dislocations, as well as quantifying the dislocation vibration amplitude, the phonon energy transmission, and the thermal resistance of the misfit interfaces, this work has illustrated and quantified two mechanisms for phonon-dislocation interaction: (1) phonon scattering by the strain field of dislocations, and (2) phonon scattering by dislocations that vibrate via the local modes of a dislocation network; the latter, leads to resonant phonon-dislocation interaction, which is manifested as local maxima of out-of-phase vibration of the atoms on the two sides of the slip plane, leading to local minima of the energy transmission in the heterostructure that contains one interface. The local vibrational modes are found to be excited only by shear stress induced by transverse phonons. Among various resonant modes, the one with the lowest frequency has the strongest effect. This work has also demonstrated the collective motion of dislocations under ultrafast phonon pulses. In addition, the dynamic properties of the misfit dislocation network localized within one interface are found to be significantly altered by the presence of misfit dislocations at other interfaces, thus further confirming the cooperative dynamic nature of the motion of dislocations and phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emergent Glassy Behavior in a Kagome Rydberg Atom Array

Here, we present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids—as well as the crossover towards a trivial paramagnetic phase—are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

74 ATOMIC AND MOLECULAR PHYSICS↗