Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Cu Evolution over Bimetallic Cu‐Y/Beta Zeolite Under H 2 and Ethanol Atmospheres: Unveiling the Role of Diatomic Metal–Metal Interactions

Understanding the dynamic evolution of Cu species under varying environmental conditions is critical for addressing challenges related to the activity and the stability of copper‐based catalysts in thermo‐, photo‐, and electrocatalysis. However, metal–metal interactions between dual single atoms and their effects on Cu evolution after exposure to different environmental molecules remain underexplored. Herein, we synthesized bimetallic Cu‐Y/Beta catalysts with dual single‐atom Cu and Y sites and monometallic Cu‐Beta catalysts with isolated Cu sites in dealuminated Beta zeolites. By varying Cu and Y compositions, diatomic interactions were studied under H 2 and ethanol atmospheres. With 6 wt% Y loading, approximately 0.4 wt% of Cu species in Cu‐Y/Beta remained partially oxidized as Cu(I) after reduction in pure H 2 at 350 °C, in contrast to the full transition to metallic Cu observed in Cu‐Beta. Combining X‐ray absorption spectroscopy with kinetic studies revealed that metallic Cu became the predominant species after reduction with H 2 as Cu loading increased from 0.4 to 1.7 wt%, quadrupling the initial ethanol dehydrogenation rate and demonstrating the dominant role of Cu(0) sites. In conclusion, scanning transmission electron microscopy and density functional theory simulations indicated spatial proximity between dual single‐atom Cu and Y sites and elucidated Cu speciation controlled by diatomic interactions.

Copper↗
DOE OSTI2019

Interaction of Methane with Calcium Ferrite in the Chemical Looping Partial Oxidation Application: Experimental and DFT Study

In this study, experimental studies using mass spectroscopy, thermogravimetric analysis, in situ X-ray diffraction (XRD), and theoretical density functional theory calculations were carried out to investigate the mechanism of methane partial oxidation over a CaFe 2 O 4 oxygen carrier. The XRD analysis revealed the formation of the calcium ferrite phases during methane reduction occurred in the sequence CaFe 2 O 4 →CaFe 3 O 5 + Ca 2 Fe 2 O 5 →FeO + Fe + CaO→Fe + CaO. Each of these phases was modeled using periodic density functional theory to study the activity of these phases for the reactivity with CH 4 , CO, and H 2 along with the adsorbed surface species to help elucidate the effect that phase changes during reduction have on the chemical looping products CO and H 2 . Mass spectroscopic analysis revealed the formation of H 2 O and CO 2 occurred prior to the formation of CO and H 2 , and these products may be due to the removal of lattice oxygen bound in the discrete primary phase, which contributes to a small portion (~16%) of the total transferable oxygen in the solid. The density functional theory (DFT) results suggest that the gaseous CO and H 2 are not thermodynamically favorable to react with the Ca 2 Fe 2 O 5 where the release of CO and H 2 led to an increase in selectivity for these products. The DFT results suggest that some of the adsorbed atomic hydrogen combined to form adsorbed H 2 and then gaseous H 2 on the CaFe 2 O 4 surface, while the carbon deposition and some adsorbed atomic hydrogen react to form the CaFe 3 O 5 and Ca 2 Fe 2 O 5 phases simultaneously. DFT simulations suggest that both the gaseous H 2 and adsorbed atomic H are less thermodynamically favorable to drive the deep reduction of Ca 2 Fe 2 O 5 to CaO and Fe0. The deep reduction of the calcium ferrite therefore favors the generation of syngas products (CO and H 2 ) through mechanisms of methane dehydrogenation and selective oxidation of carbon to CO. The significant difference in mechanisms for partial oxidation of methane to syngas with the CaFe 2 O 4 carrier can be correlated to the differences in the oxide phases present during the reduction process, concentration of oxygen species, and bonding in the lattice structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Materials selection for long life in low earth orbit - A critical evaluation of atomic oxygen testing with thermal atom systems

The use of thermal atom test methods as a materials selection and screening technique for low-earth orbit (LEO) spacecraft is critically evaluated. The chemistry and physics of thermal atom environments are compared with the LEO environment. The relative reactivities of a number of materials determined in thermal atom environments are compared with those observed in LEO and in high-quality LEO simulations. Reaction efficiencies (cu cm/atom) measured in a new type of thermal atom apparatus are one-thousandth to one ten-thousandth those observed in LEO, and many materials showing nearly identical reactivities in LEO show relative reactivities differing by as much as a factor of eight in thermal atom systems. A simple phenomenological kinetic model for the reaction of oxygen atoms with organic materials can be used to explain the differences in reactivity in different environments. Certain speciic thermal atom test environments can be used as reliable materials screening tools.

Koontz, S. L.↗
NASA NTRS1989

Arc jet diagnostics tests

Two objectives were addressed during a 10 week 1988 NASA/ASEE summer faculty fellowship at the Johnson Space Center Atmospheric Reentry Materials Structures Evaluation Facility (ARMSEF). These objectives were the evaluation of mass spectrometry for the measurement of atomic and molecular species in an arc jet environment, and the determination of atomic recombination coefficients for reaction cured glass (RCG) coated high temperature surface insulation (HRSI) materials subjected to simulated reentry conditions. Evaluation of mass spectrometry for the measurement of atomic and molecular species provided some of the first measurements of point compositions in arc jet tunnel environments. A major objective of this project centered around the sampling residence time. A three staged vacuum sampling system pulled the molecules and atoms from the arc jet to a quadrupole ionization mass spectrometer in 400 milliseconds. Conditions investigated included a composition survey across the nozzle exit at 3 cm z-distance from the nozzle exit for 3 different currents. Also, a point composition survey was taken around a shock created by the presence of a blunt body.

Willey, Ronald J.↗
DOE OSTI2021

Challenges in solving structures from radiation-damaged tomograms of protein nanocrystals assessed by simulation

Structure-determination methods are needed to resolve the atomic details that underlie protein function. X-ray crystallography has provided most of our knowledge of protein structure, but is constrained by the need for large, well ordered crystals and the loss of phase information. The rapidly developing methods of serial femtosecond crystallography, micro-electron diffraction and single-particle reconstruction circumvent the first of these limitations by enabling data collection from nanocrystals or purified proteins. However, the first two methods also suffer from the phase problem, while many proteins fall below the molecular-weight threshold required for single-particle reconstruction. Cryo-electron tomography of protein nanocrystals has the potential to overcome these obstacles of mainstream structure-determination methods. Here, a data-processing scheme is presented that combines routines from X-ray crystallography and new algorithms that have been developed to solve structures from tomograms of nanocrystals. This pipeline handles image-processing challenges specific to tomographic sampling of periodic specimens and is validated using simulated crystals. The tolerance of this workflow to the effects of radiation damage is also assessed. The simulations indicate a trade-off between a wider tilt range to facilitate merging data from multiple tomograms and a smaller tilt increment to improve phase accuracy. Since phase errors, but not merging errors, can be overcome with additional data sets, these results recommend distributing the dose over a wide angular range rather than using a finer sampling interval to solve the protein structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Optimization of algorithmic errors in analog quantum simulations

Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. Furthermore, this computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2004

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗
DOE OSTI2020

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

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

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗
DOE OSTI2024

A Simulation Study of Quantum Clock Synchronization Using Teleportation

An important requirement in implementing distributed computing and sensing application is the synchronization of clocks at various locations. The Internet relies on the Network Time Protocol (NTP), which synchronizes clocks with accuracy in the order of milliseconds. More recently an ensemble of atomic clocks is used for navigation based on GPS. These clocks are highly accurate and provide time with very low uncertainty. Even so, many physics experiments such as distributed LIGO-based systems may require more accurate clock synchronization that is achievable using quantum entanglement. This requires the deployment of a network of quantum clocks synchronized by exploiting entangled atomic clock qubits. In this paper, we carry out a simulation study of synchronizing a network of quantum clocks interconnected by a fiber plant that supports the ESnet; the latter is used to support the classical communication needed for teleportation. We consider an existing protocol for synchro-nizing the atomic clock qubits that relies on the GHZ states. To assess the performance of the protocol we developed a discrete-event simulation of the network using IBM Qiskit framework for underlying quantum gate operations and measurements. The simulation results shed light on the resources required in terms of the number entangled qubits and the time needed to achieve the synchronization of different number of nodes in ESnet.

Kiran, Mariam↗
DOE OSTI2025

CuXASNet: Rapid and accurate prediction of copper L-edge x-ray absorption spectra using machine learning

In this work, we have developed CuXASNet, a dense neural network that predicts simulated Cu -edge x-ray absorption spectra (XAS) from atomic structures. Featurization of the Cu local environment is performed using a component of M3GNet, a graph neural network developed for predicting the potential energy surface. CuXASNet is trained on simulated spectra from FEFF9 at the multiple scattering level of theory, and can predict the and edges for Cu sites to quantitative accuracy. To validate our approach, we compare 14 experimental spectra extracted from the literature with the predictions of CuXASNet. The agreement of CuXASNet with experiments is shown by an average mean absolute error of 0.125 and an average Spearman's correlation coefficient of 0.891, which is comparable to FEFF9's values of 0.131 and 0.898 for the same metrics. As such, CuXASNet can rapidly predict a large number of -edge XAS spectra at the same accuracy as FEFF9 simulations. This can be used as a drop-in replacement for multiple scattering codes for fast screening of candidate atomic structure models of a measured system. This model establishes a general framework for Cu XAS prediction, and can be extended to more computationally expensive levels of theory and to other transition metal edges.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Atomic Oxygen Durability of Second Surface Silver Microsheet Glass Concentrators

Second surface silver microsheet glass concentrators are being developed for potential use in future solar dynamic space power systems. Traditional concentrators are aluminum honeycomb sandwich composites with either aluminum or graphite epoxy face sheets, where a reflective aluminum layer is deposited onto an organic leveling layer on the face sheet. To protect the underlying layers, a SiO2 layer is applied on top of the aluminum reflective layer. These concentrators may be vulnerable to atomic oxygen degradation due to possible atomic oxygen attack of the organic layers at defect sites in the protective and reflective coatings. A second surface microsheet glass concentrator would be inherently more atomic oxygen durable than these first surface concentrators. In addition, a second surface microsheet glass concentrator design provides a smooth optical surface and allows for silver to be used as a reflective layer, which would improve the reflectivity of the concentrator and the performance of the system. A potential threat to the performance of second surface microsheet glass concentrators is atomic oxygen attack of the underlying silver at seams and edges or at micrometeoroid and debris (MMD) impacts sites. Second surface silver microsheet glass concentrator samples were fabricated and tested for atomic oxygen durability. The samples were iteratively exposed to an atomic oxygen environment in a plasma asher. Samples were evaluated for potential degradation at fabrication seams, simulated MMD impact sites, and edges. Optical microscopy was used to evaluate atomic oxygen degradation. Reflectance was obtained for an impacted sample prior to and after atomic oxygen exposure. After an initial atomic oxygen exposure to an effective fluence of approx. 1 x 10(exp 21) atoms/cm(exp 2), oxidation of the silver at defect sites and edges was observed. Exposure to an additional approx. 1 x 10(exp 21) atoms/cm(exp 2) caused no observed increase in oxidation. Oxidation at an impact site caused negligible changes in reflectance. In all cases oxidation was found to be confined to the vicinity of the seams, impact sites, edges or defect sites. Asher to in-space atomic oxygen correlation issues will be addressed.

deGroh, Kim K.↗
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Electronic signatures of Lorentzian dynamics and charge fluctuations in lithiated graphite structures

Lithium graphite intercalation compounds (Li-GICs) are essential materials for modern day portable electronics and obtaining insights into their atomic structure and thermodynamics is of fundamental interest. Here we explore the electronic and atomic states of Li-GICs at varying degrees of Lithium loading (i.e., “staging”) by means of ab-initio molecular dynamics simulations and simulated X-ray adsorption spectroscopy (XAS). We analyze the atomic correlation functions and shows that the enhancements of the Li-ion entropy with increased staging result from Lorentzian lithium-ion dynamics and charge fluctuations, which activate low-energy phonon modes. The associated electronic signatures are modulations of the unoccupied π*/σ* orbital energy levels and unambiguous fingerprints in Carbon K-edge XAS spectra. Thus, we extend the canonical view of XAS, establishing that these “static” measurements in fact encode the signature of the thermodynamic response and relaxation dynamics of the system. This causal link between atomic structure, spectroscopy, thermodynamics, and information theory can be generally exploited to better understand stability in solid-state electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2024

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗
DOE OSTI2020

Reconstruction of effective potential from statistical analysis of dynamic trajectories

The broad incorporation of microscopic methods is yielding a wealth of information on the atomic and mesoscale dynamics of individual atoms, molecules, and particles on surfaces and in open volumes. Analysis of such data necessitates statistical frameworks to convert observed dynamic behaviors to effective properties of materials. Here, we develop a method for the stochastic reconstruction of effective local potentials solely from observed structural data collected from molecular dynamics simulations (i.e., data analogous to those obtained via atomically resolved microscopies). Using the silicon vacancy defect in graphene as a model, we apply the statistical framework presented herein to reconstruct the free energy landscape from the calculated atomic displacements. Evidence of consistency between the reconstructed local potential and the trajectory data from which it was produced is presented, along with a quantitative assessment of the uncertainty in the inferred parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗