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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamic formation of preferentially lattice oriented, self trapped hydrogen clusters

Abstract A series of MD and DFT simulations were performed to investigate hydrogen self-clustering and retention in tungsten. Using a newly develop machine learned interatomic potential, spontaneous formation of hydrogen platelets was observed after implanting low-energy hydrogen into tungsten at high fluxes and temperatures. The platelets formed along low miller index orientations and neighboring tetrahedral and octahedral sites and could grow to over 50 atoms in size. High temperatures above 600 K and high hydrogen concentrations were needed to observe significant platelet formation. A critical platelet size of six hydrogen atoms was needed for long term stability. Platelets smaller than this were found to be thermally unstable within a few nanoseconds. To verify these observations, characteristic platelets from the MD simulations were simulated using large-scale DFT. DFT corroborated the MD results in that large platelets were also found to be dynamically stable for five or more hydrogen atoms. The LDOS from the DFT simulated platelets indicated that hydrogen atoms, particularly at the periphery of the platelet, were found to be at least as stable as hydrogen atoms in bulk tungsten. In addition, electrons were found to be localized around hydrogen atoms in the platelet itself and that hydrogen atoms up to 4.2 Å away within the platelet were found to share charge suggesting that the hydrogen atoms are interacting across longer distances than previously suggested. These results reveal a self-clustering mechanisms for hydrogen within tungsten in the absence of radiation induced or microstructural defects that could be a precursor to blistering and potentially explain the experimentally observed high hydrogen retention particularly in the near surface region.

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

Dynamic Tuning of Moiré Excitons in a WSe 2 /WS 2 Heterostructure via Mechanical Deformation

Moiré superlattices in van der Waals (vdW) heterostructures form by stacking atomically thin layers on top of one another with a twist angle or lattice mismatch. The resulting moiré potential leads to a strong modification of the band structure, which can give rise to exotic quantum phenomena ranging from correlated insulators and superconductors to moiré excitons and Wigner crystals. Here, in this study, we demonstrate the dynamic tuning of moiré potential in a WSe 2 /WS 2 heterostructure at cryogenic temperature. We utilize the optical fiber tip of a cryogenic scanning near-field optical microscope (SNOM) to locally deform the heterostructure and measure its near-field optical response simultaneously. The deformation of the heterostructure increases the moiré potential, which leads to a red shift of the moiré exciton resonances. We observe the interlayer exciton resonance shifts up to 20 meV, while the intralayer exciton resonances shift up to 17 meV.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Revealing atomic-to-nanoscale oxidation mechanisms of metallic materials

Oxidation and corrosion are the leading causes of degradation and failure of metallic materials. Future alloy development requires the incorporation of corrosion resistance into alloy design and processing, and this begins with attaining a fundamental insight into dynamic reaction mechanisms and kinetics between a metal and aggressive environments. With recent advances in environmental transmission electron microscopy, there have been increased efforts to apply this approach to atomically probe the microscopic mechanisms that govern the oxidation and corrosion behavior of metallic materials. Consequently, fundamental insights have been obtained in understanding the underlying processes of oxidation and passivation, including active sites, surface restructuring, oxide/metal interface dynamics, and microstructure and phase evolution. Additionally, we discuss on-going and future developments that are expected to significantly advance the field of high-temperature oxidation and corrosion.

25 ENERGY STORAGE↗
DOE OSTI2022

The effect of Cr alloying on defect migration at Ni grain boundaries

Mass transport along grain boundaries in alloys depends not only on the atomic structure of the boundary, but also its chemical make-up. In this work, we use molecular dynamics to examine the effect of Cr alloying on interstitial and vacancy-mediated transport at a variety of grain boundaries in Ni. We find that, in general, Cr tends to reduce the rate of mass transport, an effect which is greatest for interstitials at pure tilt boundaries. However, there are special scenarios in which it can greatly enhance atomic mobility. Cr tends to migrate faster than Ni, though again this depends on the structure of the grain boundary. Further, grain boundary mobility, which is sometimes pronounced for pure Ni grain boundaries, is eliminated on the time scales of our simulations when Cr is present. We conclude that the enhanced transport and grain boundary mobility often seen in this system in experimental studies is the result of non-equilibrium effects and is not intrinsic to the alloyed grain boundary. These results provide new insight into the role of grain boundary alloying on transport that can help in the interpretation of experimental results and the development of predictive models of materials evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗
DOE OSTI2026

Physics-informed machine learning exploration of Na storage mechanisms in disordered carbon

Sodium-ion batteries are a cost-effective, sustainable alternative to lithium-ion systems for large-scale energy storage. However, optimizing sodium storage in carbon-based anodes with microstructural complexity and atomic disorder remains a major challenge. The intrinsic inhomogeneity of these materials produces diverse local environments, making it difficult for conventional methods to predict and control ion dynamics. Hard carbon (HC) anodes, composed of ranges of ordered-to-disordered graphitic and amorphous nanodomains, offer tunable ion storage and rate capacity, yet rationale design remains a challenge due to poorly understood correlation between local atomic feature and ion transport mechanism. Here, to address this challenge, we introduce a data-driven framework that integrates validated machine-learned interatomic potentials, large-scale molecular dynamics simulations, and machine learning to elucidate sodium transport mechanisms as a function of carbon and sodium loading densities. By computing per-ion structural descriptors and applying unsupervised learning, we identify distinct diffusion modes governed by microscopic features. Supervised analysis and correlation mapping then establish quantitative links between these transport regimes and processing variables such as bulk carbon density and sodium content. This physics-informed approach establishes quantitative structure–transport relationships and offers actionable design principles for engineering high-performance HC anodes.

Data-driven framework↗
DOE OSTI2025

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

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

Atomic fluctuations in electronic materials revealed by dephasing

The microscopic origin and timescale of the fluctuations of the energies of electronic states has a significant impact on the properties of interest of electronic materials, with implication in fields ranging from photovoltaic devices to quantum information processing. Spectroscopic investigations of coherent dynamics provide a direct measurement of electronic fluctuations. Modern multidimensional spectroscopy techniques allow the mapping of coherent processes along multiple time or frequency axes and thus allow unprecedented discrimination between different sources of electronic dephasing. Exploiting modern abilities in coherence mapping in both amplitude and phase, we unravel dissipative processes of electronic coherences in the model system of CdSe quantum dots (QDs). The method allows the assignment of the nature of the observed coherence as vibrational or electronic. The expected coherence maps are obtained for the coherent longitudinal optical (LO) phonon, which serves as an internal standard and confirms the sensitivity of the technique. Fast dephasing is observed between the first two exciton states, despite their shared electron state and common environment. This result is contrary to predictions of the standard effective mass model for these materials, in which the exciton levels are strongly correlated through a common size dependence. In contrast, the experiment is in agreement with ab initio molecular dynamics of a single QD. Electronic dephasing in these materials is thus dominated by the realistic electronic structure arising from fluctuations at the atomic level rather than static size distribution. The analysis of electronic dephasing thereby uniquely enables the study of electronic fluctuations in complex materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulation, Structural, Thermal and Mechanical Properties of the FeTiTaVW High Entropy Alloy

Developing new materials to be applied in extreme environments is an opportunity and a challenge for the future. High entropy alloys are new materials that seem promising approaches to work in nuclear fusion reactors. In this work, FeTaTiVW high entropy alloys were developed and characterized with Molecular Dynamic and Hybrid Molecular Dynamic Monte Carlo simulations. The simulation results show that phase separation originates a lower potential energy per atom and a high level of segregation compared to those of a uniform solid solution. Moreover, the experimental diffractogram of the milled powder shows the formation of a body-centred cubic-type structure and the presence of TiO2. In addition, the microstructure of the consolidated material evidenced three phases: W-rich, Ti-rich, and a phase with all the elements. This phase separation observed in the microstructure agrees with the Hybrid Molecular Dynamic Monte Carlo simulation. Moreover, the consolidated material’s thermal conductivity and specific heat are almost constant from 25 °C to 1000 °C, and linear expansion increases with increasing temperature. On the other hand, specific heat and thermal expansion values are in between CuCrZr and W values (materials chosen for the reactor walls). The FeTaTiVW high entropy alloy evidences a ductile behaviour at 1000 °C. Therefore, the promising thermal properties of this system can be attributed to the multiple phases and systems with different compositions of the same elements, which is exciting for future developments.

Martins, Ricardo (ORCID:0000000177685704)↗
DOE OSTI2024

XFEL SASE pulses can enhance time-dependent observables

Abstract X-ray free electron lasers (XFELs) have emerged as powerful sources of short and intense x-ray pulses. We propose a simple and robust procedure which takes advantage of the inherent stochasticity of self-amplified stimulated emission (SASE) pulses to enhance the time-resolution and signal strength of the recorded data. Notably, the proposed method is able to enhance the average signal without knowledge of the signal strength of individual shots. Simple metrics for the probe pulses are introduced, such as an effective pulse duration applicable to SASE pulses characterised in the time domain using e.g. an X-band transverse cavity. The approach is evaluated using simulated and real pulse data in the context of ultrafast electron dynamics in a molecule. Utilising H 2 as a model system, we demonstrate the efficacy of the method theoretically, successfully enhancing the predicted nonresonant ultrafast x-ray scattering signal associated with electron dynamics. The method presented is broadly applicable and offers a general strategy for enhancing time-dependent observables at XFELs.

data analysis↗
DOE OSTI2024

Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry

Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550–600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.

Polymer Science↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning approach to pattern recognition in nuclear dynamics from the ab initio symmetry-adapted no-core shell model

A novel machine learning approach is used to provide further insight into atomic nuclei and to detect orderly patterns amid a vast data of large-scale calculations. The method utilizes a neural network that is trained on ab initio results from the symmetry-adapted no-core shell model (SA-NCSM) for light nuclei. We show that the SA-NCSM, which expands ab initio applications up to medium-mass nuclei by using dominant symmetries of nuclear dynamics, can reach heavier nuclei when coupled with the machine learning approach. In particular, we find that a neural network trained on probability amplitudes for s- and p-shell nuclear wave functions not only predicts dominant configurations for heavier nuclei but in addition, when tested for the 20 Ne ground state, accurately reproduces the probability distribution. The non-negligible configurations predicted by the network provide an important input to the SA-NCSM for reducing ultralarge model spaces to manageable sizes that can be, in turn, utilized in SA-NCSM calculations to obtain accurate observables. The neural network is capable of describing nuclear deformation and is used to track the shape evolution along the 20-42 Mg isotopic chain, suggesting a shape coexistence that is more pronounced toward the very neutron-rich isotopes. We provide first descriptions of the structure and deformation of 24 Si and 40 Mg of interest to x-ray burst nucleosynthesis, and even of the extremely heavy nuclei such as 166,168 Er and 236 U, that build on first-principles considerations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamical Scaling as a Signature of Multiple Phase Competition in Yb 2⁢ Ti 2 ⁢O 7

Yb 2 ⁢Ti 2 ⁢O 7 is a celebrated example of a pyrochlore magnet with highly frustrated, anisotropic exchange interactions. To date, attention has largely focused on its unusual, static properties, many of which can be understood as coming from the competition between different types of magnetic order. Here, in this work, we use inelastic neutron scattering with exceptionally high energy resolution to explore the dynamical properties of Yb 2 ⁢Ti 2 ⁢O 7 . We find that spin correlations exhibit dynamical scaling, analogous to behavior found near to a quantum critical point. We show that the observed scaling collapse can be explained within a phenomenological theory of multiple-phase competition, and confirm that a scaling collapse is also seen in semiclassical simulations of a microscopic model of Yb 2 ⁢Ti 2 ⁢O 7 . These results suggest that dynamical scaling may be general to systems with competing ground states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Catching the Killer: Dynamic Disorder Design Rules for Small–Molecule Organic Semiconductors

Small-molecule organic semiconductors are promising materials for applications ranging from solar cells to medical sensors. Their nearly infinite design space means that, in theory, it is possible to tailor the material to the exact specifications of a particular application. In reality; however, design rules to improve mobility (μ) remain elusive because of the complex calculations required to understand its limiters. Transient localization theory posits that charge carriers are slowed down by the collective phonon motions, called dynamic disorder, which localize charge carriers temporarily. Traditionally, a mode-by-mode analysis is performed to try and identify “killer” modes. Herein, the flaws are demonstrated with this analysis and present a streamlined simulation workflow that simplifies simulation of dynamic disorder and enables engineering of new molecular structures based on phonon dynamics. This workflow is combined with a novel visualization technique that enables per-atom insights into how μ is limited. Furthermore, this workflow is applied and analyzed to a series of -acenes and a series of BTBT-based molecules. Then design rules are identified in each series that identify possible ways to improve the μ beyond current experimental limits using phonon engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Recent Advances in X-ray Absorption Near Edge Structure (XANES) Simulation for Catalysis: Theories and Applications

X-ray absorption spectroscopy (XAS) is a powerful spectroscopic technique for characterizing catalytic active sites. In particular, X-ray Absorption Near Edge Structure (XANES) provides insights into the oxidation state and local coordination environment of the absorbing atom, shedding light on the catalytic structure. Computational XANES simulations provide one-on-one correspondence between the molecular structure and their corresponding spectra, making them an effective approach for interpreting experimental XANES features. This work outlines various theories used to simulate XANES spectra and highlights recent advancements in computational XANES applications for homogeneous, heterogeneous, and single-atom catalysts focusing on identifying the active site structures, the site heterogeneity, and the dynamic evolution of catalysts.

Xu, Jiayi↗
DOE OSTI2026

Chemical complexity suppresses interstitial solute segregation to screw dislocation cores in the bcc NbTaTiHfZr high-entropy alloy

Recent reports suggest that interstitial solute additions to body-centered cubic (bcc) high-entropy alloys (HEAs) may enhance their mechanical properties. However, details of interactions between interstitial atoms and dislocations in these HEAs remain incompletely understood. Using first-principles calculations, we examine the energetics of C, N, and O interstitial solutes in elemental bcc metals and NbTaTiHfZr, focusing on their segregation in screw dislocation cores. We examine two types of core sites and show that the low-energy configuration and associated segregation energy depend on the transition metal group. In NbTaTiHfZr, we find that chemical complexity substantially suppresses interstitial segregation in dislocation cores, as fluctuations in bulk solute energies induce energetically favorable sites that disfavor segregation to core sites. Local chemical order further enhances this effect by promoting solute clustering away from dislocation cores. These findings reveal fundamental differences between the behavior of interstitial atoms in elemental metals and HEAs, with implications for high-temperature plasticity and dynamic strain aging.

Body-centered cubic transition metals↗