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.

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗

DOE OSTI2023

Topology driven and soft phonon mode enabled Na-ion diffusion in quaternary chalcogenides, Na 3 ZnGaX 4 (X = S, and Se)

The compounds Na 3 ZnGaX 4 (X = S, Se) are potential solid electrolyte materials in sodium-based batteries, which have certain advantages over oxide materials and have shown significant ionic conductivity at ambient temperature. In this paper, we bring out atomic-level features of the diffusion process in these new materials using the microscopic techniques of inelastic neutron scattering (INS), quasielastic neutron scattering (QENS), and ab initio molecular dynamics (AIMD) simulations. The insights obtained from these techniques are unique and not available from other macroscopic experiments. Neutron scattering experiments have been performed at temperatures from 100 to 700 K. The simulations have been carried out up to 900 K. We have calculated the phonon spectra and the space–time correlation functions and found good agreement with the results of the neutron scattering experiments. The simulations enable detailed analysis of the atomic-site dependent dynamical information. We observe low-energy phonon modes of ~6 meV involving the vibrations of certain Na atoms in the lattice. This reveals that the Na at 32g Wyckoff sites (Na2) has sufficiently shallow potential among the two available crystallographic sites. This shallow potential facilitates diffusion. Furthermore, the specific structural topology of the network of interconnected zig-zag chains of the Na2 atomic sites provides the low-barrier energy pathways for diffusion. A small fraction of vacancy defects appears essential for diffusion. We further observe that the Na2 atoms undergo jump-like diffusion to the vacant next or the 2nd next neighbour sites at ~4 Å. While the QENS experiments reveal the jump-like diffusion and its time scale, detailed analysis of the AIMD simulations shows that the jumps appear mostly along zig-zag chains of the Na2 sites in the tetragonal ab-plane, as well as between the chains along the c-axis. In conclusion, the net diffusion is essentially 3-dimensional, with little anisotropy despite the anisotropy of the tetragonal crystal structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Reactor wall effects in Si–Cl 2 –Ar atomic layer etching

This work complements our previous manuscript [J. Vac. Sci. Technol. A41, 062602 (2023)] where predictions from molecular dynamics (MD) simulations of silicon–chlorine–argon (Si–Cl 2 –Ar) atomic layer etching (ALE) are compared to experiments. When etch product distributions for atomic chlorine (Cl) and silicon chlorides were initially compared to optical emission spectroscopy (OES) signals, it appeared that there was a discrepancy between the MD predictions and experimental results at higher ion fluences. Experiments showed a relatively long period of nearly constant Cl-containing etch products released from the ion-bombarded surface (referred to as the “plateau”) but this effect was not observed in MD simulations. In this report, we demonstrate that the “plateau” observed in the OES signals is most likely due to the desorption of Cl-containing etch products from the walls of the reactor and subsequent adsorption on the Si substrate. Experiments varying the gas residence time in the chamber while keeping incoming gas concentrations and pressure constant support this interpretation. We also conducted experiments with an additional Ar-only flow in the chamber to reduce the concentration of Cl-containing species on the chamber walls. For both sets of flow modification experiments, we observe results consistent with the hypothesis that Cl-containing species desorbing from chamber walls are a significant cause of the observed discrepancy between MD predictions and experimental observations. If the measured OES signals are corrected for this “additional” source of Cl-containing species at the surface, the MD predictions and measured OES signals are in excellent agreement. This further supports the predictive capability of MD simulations to accurately capture the relevant physical and chemical processes in plasma-assisted ALE processes. We provide an order of magnitude estimate of the required density of Cl-containing species that would account for the additional etch products observed. Finally, we discuss the implications of this effect on ALE in plasma nanofabrication.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Dataset of simulated vibrational density of states and X-ray diffraction profiles of mechanically deformed and disordered atomic structures in Gold, Iron, Magnesium, and Silicon

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

NASA NTRSDate not supplied

Nitric Oxide Concentrations from Solar Maximum to Solar Minimum Years Derived from SABER

Nitric oxide (NO) is of particular importance in the upper atmosphere nominally the region between 100 km and 250 km, as it plays the role in cooling and regulating the energy budget of the system. The NO cooling rates have been routinely derived for the past two decades from the measurements made by Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite (TIMED) satellite since 2002. However, a comprehensive long-term dataset of global thermospheric NO concentration as a function of local time is currently unavailable. Given that solar activity significantly influences NO abundance on timescales ranging from days to decades, analyzing the behavior of NO concentration during solar storm and solar cycle can provide valuable insights into thermospheric energy balance, improving storm-time and climatological knowledge. In this study, we used newly calculated NO cooling rates and simulations of neutral kinetic temperature and atomic oxygen number density from the MSIS 2.1 model to derive 20-year NO concentration data. The global annual mean NO concentration shows a strong correlation with the solar cycle, with high concentrations and upper boundaries up to 250 km during solar maximum years and lower concentrations and upper boundaries up to ~230 km during solar minimum years. We also compared NO concentrations obtained from SABER measurements and MSIS 2.1 simulations during solar maximum (2014) and solar minimum (2009) years and found good agreement between the two during both day and night.

Ningchao Wang↗

NASA NTRS1992

Atomic oxygen effects on thin film space coatings studied by spectroscopic ellipsometry, atomic force microscopy, and laser light scattering

The University of Nebraska is currently evaluating Low Earth Orbit (LEO) simulation techniques as well as a variety of thin film protective coatings to withstand atomic oxygen (AO) degradation. Both oxygen plasma ashers and an electron cyclotron resonance (ECR) source are being used for LEO simulation. Thin film coatings are characterized by optical techniques including Variable Angle Spectroscopic Ellipsometry, Optical spectrophotometry, and laser light scatterometry. Atomic Force Microscopy (AFM) is also used to characterize surface morphology. Results on diamondlike carbon (DLC) films show that DLC degrades with simulated AO exposure at a rate comparable to Kapton polyimide. Since DLC is not as susceptible to environmental factors such as moisture absorption, it could potentially provide more accurate measurements of AO fluence on short space flights.

Synowicki, R. A.↗

DOE OSTI2020

Deep learning for NLTE spectral opacities

Computer simulations of high energy density science experiments are computationally challenging, consisting of multiple physics calculations including radiation transport, hydrodynamics, atomic physics, nuclear reactions, laser–plasma interactions, and more. To simulate inertial confinement fusion (ICF) experiments at high fidelity, each of these physics calculations should be as detailed as possible. However, this quickly becomes too computationally expensive even for modern supercomputers, and thus many simplifying assumptions are made to reduce the required computational time. Much of the research has focused on acceleration techniques for the various packages in multiphysics codes. In this work, we explore a novel method for accelerating physics packages via machine learning. The non-local thermodynamic equilibrium (NLTE) package is one of the most expensive calculations in the simulations of indirect drive inertial confinement fusion, taking several tens of percent of the total wall clock time. We explore the use of machine learning to accelerate this package, by essentially replacing the physics calculation with a deep neural network that has been trained to emulate the physics code. Overall, we demonstrate the feasibility of this approach on a simple problem and perform a side-by-side comparison of the physics calculation and the neural network inline in an ICF Hohlraum simulation. We show that the neural network achieves a 10× speed up in NLTE computational time while achieving good agreement with the physics code for several quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Imaging multiphoton ionization dynamics of CH 3 I at a high repetition rate XUV free-electron laser

XUV multiphoton ionization of molecules is commonly used in free-electron laser experiments to study charge transfer dynamics. However, molecular dissociation and electron dynamics, such as multiple photon absorption, Auger decay, and charge transfer, often happen on competing time scales, and the contributions of individual processes can be difficult to unravel. We experimentally investigate the Coulomb explosion dynamics of methyl iodide upon core–hole ionization of the shallow inner-shell of iodine (4d) and classically simulate the fragmentation by phenomenologically introducing ionization dynamics and charge transfer. Under our experimental conditions with medium fluence and relatively long XUV pulses (~75 fs), we find that fast Auger decay prior to charge transfer significantly contributes to the charging mechanism, leading to a yield enhancement of higher carbon charge states upon molecular dissociation. Furthermore, we argue for the existence of another charging mechanism for the weak fragmentation channels leading to triply charged carbon atoms. This study shows that classical simulations can be a useful tool to guide the quantum mechanical description of the femtosecond dynamics upon multiphoton absorption in molecular systems.

Optics↗

DOE OSTI2022

Bulk and boundary quantum phase transitions in a square Rydberg atom array

Motivated by recent experimental realizations of exotic phases of matter on programmable quantum simulators, we carry out a comprehensive theoretical study of quantum phase transitions in a Rydberg atom array on a square lattice, with both open and periodic boundary conditions. In the bulk, we identify several first-order and continuous phase transitions by performing large-scale quantum Monte Carlo simulations and develop an analytical understanding of the nature of these transitions using the framework of Landau-Ginzburg-Wilson theory. Remarkably, we find that under open boundary conditions, the boundary itself undergoes a second-order quantum phase transition, independent of the bulk. These results explain recent experimental observations and provide important insights into both the adiabatic state preparation of novel quantum phases and quantum optimization using Rydberg atom array platforms.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Structure of disordered TiO 2 phases from ab initio based deep neural network simulations

Amorphous TiO 2 (a-TiO 2 ) is widely used in many fields, ranging from photo-electrochemistry to bio-engineering, hence detailed knowledge of its atomic structure is of scientific and technological interest. Here we use an ab initio-based deep neural network potential (DP) to simulate large scale atomic models of crystalline and disordered TiO 2 with molecular dynamics. Our DP reproduces the structural properties of all (11) TiO 2 crystalline phases, predicts the densities and structure factors of molten and amorphous TiO 2 with only a few percent deviation from experiments, and describes the pressure dependence of the amorphous structure in agreement with recent observations. Furthermore, it can be extended to model additional structures and compositions and can be thus of great value in the study of TiO 2 -based (nano-)materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Understanding the phase transformation mechanisms that affect the dynamic response of Fe-based microstructures at the atomic scales

Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer↗

NASA NTRSDate not supplied

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Gabriel Plummer↗

DOE OSTI2020

On the Stark Effect of the O I 777-nm Triplet in Plasma and Laser Fields

The O I 777-nm triplet transition is often used for plasma density diagnostics. It is also employed in nonlinear optics setups for producing quasi-comb structures when pumped by a near-resonant laser field. Here, we apply computer simulations to situations of the radiating atom subjected to the plasma microfields, laser fields, and both perturbations together. Our results, in particular, resolve a controversy related to the spectral line anomalously broadened in some laser-produced plasmas. The importance of using time-dependent density matrix is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Effects of elevated-temperature deposition on the atomic structure of amorphous Ta 2 O 5 films

Brownian thermal noise as a result of mechanical loss in optical coatings will become the dominant source of noise at the most sensitive frequencies of ground-based gravitational-wave detectors. Experiments found, however, that a candidate material, amorphous Ta 2 O 5 , is unable to form an ultrastable glass and, consequently, to yield a film with significantly reduced mechanical loss through elevated-temperature deposition alone. X-ray scattering PDF measurements are carried out on films deposited and subsequently annealed at various temperatures. Inverse atomic modeling is used to analyze the short and medium range features in the atomic structure of these films. Furthermore, in silico deposition simulations of Ta 2 O 5 are carried out at various substrate temperatures and an atomic level analysis of the growth at high temperatures is presented. It is observed that upon elevated-temperature deposition, short range features remain identical, whereas medium range order increases. After annealing, however, both the short and medium range orders of films deposited at different substrate temperatures are nearly identical. A discussion on the surface diffusion and glass transition temperatures indicates that future pursuits of ultrastable low-mechanical-loss films through elevated temperature deposition should focus on materials with a high surface mobility, and/or lower glass transition temperatures in the range of achievable deposition temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1987

Image Processing

A new spinoff product was derived from Geospectra Corporation's expertise in processing LANDSAT data in a software package. Called ATOM (for Automatic Topographic Mapping), it's capable of digitally extracting elevation information from stereo photos taken by spaceborne cameras. ATOM offers a new dimension of realism in applications involving terrain simulations, producing extremely precise maps of an area's elevations at a lower cost than traditional methods. ATOM has a number of applications involving defense training simulations and offers utility in architecture, urban planning, forestry, petroleum and mineral exploration.

Source record↗