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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1976

Vibration-translation energy transfer in vibrationally excited diatomic molecules

A semiclassical collision model is applied to the study of energy transfer rates between a vibrationally excited diatomic molecule and a structureless atom. The molecule is modeled as an anharmonic oscillator with a multitude of dynamically coupled vibrational states. Three main aspects in the prediction of vibrational energy transfer rates are considered. The applicability of the semiclassical model to an anharmonic oscillator is first evaluated for collinear encounters. Second, the collinear semiclassical model is applied to obtain numerical predictions of the vibrational energy transfer rate dependence on the initial vibrational state quantum number. Thermally averaged vibration-translation rate coefficients are predicted and compared with CO-He experimental values for both ground and excited initial states. The numerical model is also used as a basis for evaluating several less complete but analytic models. Third, the role of rational motion in the dynamics of vibrational energy transfer is examined. A three-dimensional semiclassical collision model is constructed with coupled rotational motion included. Energy transfer within the molecule is shown to be dominated by vibration-rotation transitions with small changes in angular momentum. The rates of vibrational energy transfer in molecules with rational frequencies that are very small in comparison to their vibrational frequency are shown to be adequately treated by the preceding collinear models.

Mckenzie, R. L.↗
NASA NTRS1985

Balloon-borne three-meter telescope for far-infrared and submillimeter astronomy

Presented are scientific objectives, engineering analysis and design, and results of technology development for a Three-Meter Balloon-Borne Far-Infrared and Submillimeter Telescope. The scientific rationale is based on two crucial instrumental capabilities: high angular resolution which approaches eight arcseconds at one hundred micron wavelength, and high resolving power spectroscopy with good sensitivity throughout the telescope's 30-micron to 1-mm wavelength range. The high angular resolution will allow us to resolve and study in detail such objects as collapsing protostellar condensations in our own galaxy, clusters of protostars in the Magellanic clouds, giant molecular clouds in nearby galaxies, and spiral arms in distant galaxies. The large aperture of the telescope will permit sensitive spectral line measurements of molecules, atoms, and ions, which can be used to probe the physical, chemical, and dynamical conditions in a wide variety of objects.

Fazio, G. G.↗
NASA NTRS1987

Collisional excitation of an asymmetric rotor, silicon dicarbide

Rotational excitation rates have been computed for the asymmetric top molecule SiC2 in collisions with low-energy He atoms. The intermolecular forces were obtained from an electron gas model, and collision dynamics were treated within the infinite-order sudden approximation. Total excitation rates, i.e., summed over final levels, are expected to be accurate to about 50 percent, and the larger state-to-state rates are likely to be within a factor of about 2 of the correct values, although some of the smaller (and less important) rates may be less accurate. These rates are also thought to reflect, within this level of accuracy, rates for excitation by collisions with H2 molecules.

Palma, Amedeo↗
NASA NTRS1994

Theoretical determination of chemical rate constants using novel time-dependent methods

The work completed within the grant period 10/1/91 through 12/31/93 falls primarily in the area of reaction dynamics using both quantum and classical mechanical methodologies. Essentially four projects have been completed and have been or are in preparation of being published. The majority of time was spent in the determination of reaction rate coefficients in the area of hydrocarbon fuel combustion reactions which are relevant to NASA's High Speed Research Program (HSRP). These reaction coefficients are important in the design of novel jet engines with low NOx emissions, which through a series of catalytic reactions contribute to the deterioration of the earth's ozone layer. A second area of research studied concerned the control of chemical reactivity using ultrashort (femtosecond) laser pulses. Recent advances in pulsed-laser technologies have opened up a vast new field to be investigated both experimentally and theoretically. The photodissociation of molecules adsorbed on surfaces using novel time-independent quantum mechanical methods was a third project. And finally, using state-of-the-art, high level ab initio electronic structure methods in conjunction with accurate quantum dynamical methods, the rovibrational energy levels of a triatomic molecule with two nonhydrogen atoms (HCN) were calculated to unprecedented levels of agreement between theory and experiment.

Dateo, Christopher E.↗
NASA NTRS2001

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. In the absence of any record of the earliest ancestors of contemporary cells, protocells, the most direct way to test our understanding of their characteristics is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs of molecules that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing and storing energy, and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following, questions about these proteins: (1) How do small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides aggregate to form membrane-spannin(y structures (e.g., channels)? (3) By what mechanisms do such aggregates perform their functions? The simulations are performed using the molecular dynamics (MD) method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated directly in real space while long-range forces are evaluated in reciprocal space, usually using a particle-mesh algorithm which is of order O(NlnN). Currently, a time step of 2 femtoseconds is typically used, thereby making studies of problems occurring on multi-nanosecond time scales (10(exp 6) - 10(exp 8) time steps) accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude. Such an extension requires both algorithmic improvements and codes scalable to a large number of parallel processors. Work in this direction is in progress. Two specific series of simulations that demonstrate how peptides self-organize and function in membranes are discussed. In one series of simulations, it was shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and, simultaneously, fold into two different helical structures, which remain in equilibrium. Once in the membrane, the peptides can readily change their orientation, especially in response to local electric fields. This structural and orientational flexibility of peptides with changing conditions may have provided a mechanism of transmitting signals between the environment and the interior of the protocell. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations demonstrate that protons move through the gate by a "shuttle" mechanism, wherein one histidine is protonated on the extracellular side and, subsequently, the proton bound on the opposite side is released.

Pohorille, Andrew↗
NASA NTRS1998

Structure of Boron Nitride Nanotubes: Tube Closing Vs. Chirality

The structure of boron nitride nanotubes is investigated using a generalized tight-binding molecular dynamics method. It is shown that dynamic relaxation results in a wavelike or "rippled" surface in which the B atoms rotate inward and the N atoms move outward, reminiscent of the surface relaxation of the III-V semiconductors. More importantly, the three different morphologies of the tube closing with flat, conical and amorphous ends, as observed in experiments, are shown to be directly related to the tube chiralities. The abundance of flat end tubes observed in experiments is, thus, shown to be an indication of the greater stability of "zig-zag" BN tubes over the "arm-chair" tubes under experimental conditions.

Srivastava, Deepak↗
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↗