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 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Unravelling the origin of reaction-driven aggregation and fragmentation of atomically dispersed Pt catalyst on ceria support

Metal–support interaction plays a crucial role in governing the stability and activity of atomically dispersed platinum catalysts on ceria support. The migration and aggregation of platinum atoms during the catalytic reaction leads to the redistribution of active sites. In this study, by utilizing a multimodal characterization scheme, we observed the aggregation of platinum atoms at high temperatures under reverse water gas shift reaction conditions and the subsequent fragmentation of platinum clusters, forming “single atoms” upon cooling. Theoretical simulations of both effects uncovered the roles of carbon monoxide binding on perimeter Pt sites in the clusters and hydrogen coverage in the aggregation and fragmentation mechanisms. Finally, this study highlights the complex effects of adsorbate and supports interactions with metal sites in Pt/ceria catalysts that govern their structural transformations under in situ conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Experimental and theoretical examination of shock-compressed copper through the fcc to bcc to melt phase transitions

Recent studies show a face-centered cubic (fcc) to body-centered cubic (bcc) transformation along the shock Hugoniot for several metals (i.e., Cu, Au, and Ag). Here, we combine laser-shock compression of Cu foils on nanosecond timescales with in situ x-ray diffraction (XRD) to examine the microstructural changes with stress. We study the fcc phase and the phase transition from fcc to bcc (pressures greater than 180 GPa). Textural analysis of the azimuthal intensities from the XRD images is consistent with transformation into the bcc phase through the Pitsch-distortion mechanism. We use embedded atom model molecular dynamics simulations to determine the stability of the bcc phase in pressure–temperature space. Finally, our results indicate that the bcc phase is stabilized only at high temperatures and remains stable at pressures greater than 500 GPa.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Structural evolution in Au- and Pd-based metallic glass forming liquids and the case for improved molecular dynamics force fields

The results of a combined experimental and computational investigation of the structural evolution of Au 81 Si 19 , Pd 82 Si 18 , and Pd 77 Cu 6 Si 17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd 77 Cu 6 Si 17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. In conclusion, a proposed order parameter distinguishing glass forming ability, 1/(S(T,q 1 )−1, shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

High-pressure and temperature neural network reactive force field for energetic materials

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive compared with electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network-based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical responses of energetic materials at extreme conditions. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and shock-loading conditions. Further, this new potential shows improved accuracy over the current state-of-the-art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock-loading conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Universal and interpretable classification of atomistic structural transitions via unsupervised graph learning

Materials processing often occurs under extreme dynamic conditions leading to a multitude of unique structural environments. These structural environments generally occur at high temperatures and/or high pressures, often under non-equilibrium conditions, which results in drastic changes in the material's structure over time. Computational techniques, such as molecular dynamics simulations, can probe the atomic regime under these extreme conditions. However, characterizing the resulting diverse atomistic structures as a material undergoes extreme changes in its structure has proved challenging due to the inherently non-linear relationship between structures as large-scale changes occur. Here, we introduce SODAS++, a universal graph neural network framework, that can accurately and intuitively quantify the atomistic structural evolution corresponding to the transition between any two arbitrary phases. We showcase SODAS++ for both solid–solid and solid–liquid transitions for systems of increasing geometric and chemical complexity, such as colloidal systems, elemental Al, rutile and amorphous TiO 2 , and the non-stoichiometric ternary alloy Ag 26 Au 5 Cu 19 . Finally, we show that SODAS++ can accurately quantify all transitions in a physically interpretable manner, showcasing the power of unsupervised graph neural network encodings for capturing the complex and non-linear pathway, a material's structure takes as it evolves.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Two-step nucleation of the Earth’s inner core

Significance Understanding the formation of the Earth’s inner core is essential to understanding the geodynamo and Earth's history. However, recent attempts to explain the initial solidification of the inner core have been unsuccessful. The supercooling necessary to form hcp iron is unrealistically large and creates the “inner core nucleation paradox.” Our work demonstrates that molten iron can crystallize via a two-step nucleation process involving the intermediate bcc phase under core conditions. This mechanism significantly reduces the required undercooling necessary to nucleate solid iron. This work also suggests that bcc and hcp iron have similar free energies at pressures near the inner core center.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Effects of interlayer spacing and applied pressure on the lanthanide transport in MoS 2 -based two-dimensional channels

Rare-earth elements (REEs) are critical to modern industry but difficult to separate due to their subtle and monotonic changes in physicochemical properties. MoS 2 -based two-dimensional (2D) materials offer novel opportunities for enhancing REE separation, exhibiting a distinct volcano-shaped transport performance distribution that peaks at Sm 3+ . However, the specific contributions of thermodynamic and kinetic factors to ion transport within 2D confinement remain unclear. In this study, we conducted a series of non-equilibrium all-atom molecular dynamics (MD) simulations to explore the effects of interlayer spacing and external pressure on the transport of lanthanide ions in Å-scale acetate functionalized 2D MoS 2 (MoS 2 -COOH) channels. We examined ion entry and permeation rates, water flux, dehydration, and binding modes. The simulation results reveal that the transport trends of lanthanide ions are jointly driven by the dehydration degree and the relative-binding strengths of ions to water and to the acetate within the 2D channels. Notably, the dehydration pattern of lanthanide ions during permeation is closely linked to kinetic factors. Overall, this study provides a detailed atomistic understanding of the mechanisms underlying lanthanide ion transport under confinement. These findings point to the significant potential for tuning confinement and chemical functionalization within Å-scale channels for more efficient REE separation.

2D channels↗

DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗

DOE OSTI2022

ezAlign

The ezAlign is aimed at clustering coarse grain simulation to find common or uncommon occurrences and convert coarse (CG) grained coordinate and topology files to atomistic formats. We use a PointNet based approach to map individual frames of simulation to points in a latent space. These points are then clustered using a variety of clustering methods. Clusters are analyzed to associate them with states in the simulation. Frames can be chosen from these clusters based on proximity to cluster centers. ezAlign takes CG coordinate and topology files and converts and outputs their corresponding atomistic formats using an alignment and relaxation procedure. ezAlign is designed to convert complex, solvated biological systems including lipid membranes with drug-like molecules using GROMACS. A GROMACS checkpoint (.cpt) file is also outputted to enable continuation simulations that retain the equilibrated atomic velocities. Independent atomistic coordinates and topologies for every molecule must already be included in ezAlign/files. A number of commonly simulated biological molecules are currently provided.

Bennett, WilliamF.↗

DOE OSTI2020

Simple and efficient algorithms for training machine learning potentials to force data

Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information than the energy alone. It may appear that training a model to this large quantity of force data would introduce significant computational costs. Actually, training to all available force data should only be a few times more expensive than training to energies alone. Here, we present a new algorithm for efficient force training, and benchmark its accuracy by training to forces from real-world datasets for organic chemistry and bulk aluminum.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Non-reactive facet specific adsorption as a route to remediation of chlorinated organic contaminants

The present work quantifies metal-contaminant interactions between palladium substrates and three salient chlorinated organic contaminants, namely trichloroethylene 1,3,5-trichlorobenzene (TCB), and 3,3′,4,4′-tetrachlorobiphenyl (PCB77). Given that Pd is one of the conventional catalytically active materials known for contaminant removal, maximizing catalytic efficiency through optimal adsorption dynamics reduces the cost of remediation of contaminants that are persistent water pollutants chronically affecting public health. Adsorption efficiency analyses from all-atom molecular dynamics (MD) simulations advance the understanding of reaction mechanisms available from density functional theory (DFT) calculations to an extractable feature scale that can fit the parametric design of supported metal catalytic systems and feed into high throughput catalyst selection. Data on residence time, site-specific adsorption, binding energies, packing geometries, orientation profiles, and the effect of adsorbate size show the anomalous behaviour of organic contaminant adsorption on the undercoordinated {110} surface as compared to the {111} and {100} surfaces. The intermolecular interaction within contaminants from molecular dynamics simulation exhibits refreshing results than ordinary single molecule density functional theory calculation. Since complete adsorption and dechlorination is an essential step for chlorinated organic contaminant remediation pathways, the presented profiles provide essential information for designing efficient remediation systems through facet-controlled palladium nanoparticles.

Guo, Hao↗

NASA NTRS1979

An analytical study of nitrogen oxides and carbon monoxide emissions in hydrocarbon combustion with added nitrogen, preliminary results

The effect of combustor operating conditions on the conversion of fuel-bound nitrogen (FBN) to nitrogen oxides NO sub x was analytically determined. The effect of FBN and of operating conditions on carbon monoxide (CO) formation was also studied. For these computations, the combustor was assumed to be a two stage, adiabatic, perfectly-stirred reactor. Propane-air was used as the combustible mixture and fuel-bound nitrogen was simulated by adding nitrogen atoms to the mixture. The oxidation of propane and formation of NO sub x and CO were modeled by a fifty-seven reaction chemical mechanism. The results for NO sub x and CO formation are given as functions of primary and secondary stage equivalence ratios and residence times.

Bittker, D. A.↗

NASA NTRS1987

Stability relationship for water droplet crystallization with the NASA Lewis icing spray

In order to produce small droplets for icing cloud simulation, high pressure air atomizing nozzles are used. For certain icing testing applications, median drop sizes as small as 5 mm are needed, which require air atomizing pressures greater than 3000 kPa. Isentropic expansion of the ambient temperature atomizing air to atmospheric pressure can result in air stream temperatures of -160 C which results in ice crystals forming in the cloud. To avoid such low temperatures, it is necessary to heat the air and water to high initial temperatures. An icing spray research program was conducted to map the temperatures below which ice crystals form. A soot slide technique was used to determine the presence of crystals in the spray.

Marek, C. John↗

NASA NTRS1988

Stability relationship for water droplet crystallization with the NASA Lewis icing spray nozzle

Marek, C. John↗

NASA NTRS1989

Laboratory investigations involving high-velocity oxygen atoms

Facilities for measuring material reactive characteristics have been under development for several years and span the atom energy range from thermal to 5 eV, the orbital collision energy. One of the high-atom energy facilities (The High Intensity/Energy Atomic Oxygen Source) capable of simulating the reactive part of LEO is described, along with results of beam characterization and preliminary material studies. The oxygen atom beam source was a continuous wave plasma produced by focusing a high-power CO2 laser through a lens system into a rare gas/molecular oxygen mixture chamber at elevated temperature. Material samples were exposed to the high velocity beam through an external feedthrough. The facility showed good stability in continued operation for more than 100 hours, producing fluences of 10 to the 21st to 10 to the 22nd atoms/sq cm. Reaction efficiencies and surface morphology have been measured for several materials at energies of 1.5 and 2.8 eV, matching with data generated from previous space flights. Activation energies for carbon and Kapton as measured in this facility were 800 cal/mole.

Leger, Lubert J.↗

NASA NTRS1991

Nonequilibrium in a low power arcjet nozzle

Emission spectroscopy measurements were made of the plasma flow inside the nozzle of a 1-kW-class arcjet thruster. The thruster propellant was a hydrogen-nitrogen mixture used to simulate fully decomposed hydrazine. Atomic electron-excitation, vibrational, and rotational temperatures were determined for the expanding plasma using relative line intensity techniques. The atomic excitation temperature decreased from 18,000 K at a location 3 mm downstream of the constrictor to 9000 K at a location 9 mm from the constrictor, while the molecular vibrational and rotational temperatures decreased from 6500 K to 2500 K and from 8000 K to 3000 K, respectively, between the same locations. The electron density, measured using hydrogen H-beta line Stark broadening, decreased from about 10 to the 15th cu cm to about 2 x 10 to the 14th cu cm during the expansion. The results show that the plasma is highly nonequilibrium throughout the nozzle, with most relaxation times equal or exceeding the particle residence time.

Myers, Roger M.↗

NASA NTRS1998

Indium Tin Oxide-Magnesium Fluoride Co-Deposited Films for Spacecraft Applications

Highly transparent coatings with a maximum sheet resistivity between 10(exp 8) and 10(exp 9) ohms/square are desired to prevent charging of solar arrays for low Earth polar orbit and geosynchronous orbit missions. Indium tin oxide (ITO) and magnesium fluoride (MgF2) were ion beam sputter co-deposited onto fused silica substrates and were evaluated for transmittance, sheet resistivity and the effects of simulated space environments including atomic oxygen (AO) and vacuum ultraviolet (VUV) radiation. Optical properties and sheet resistivity as a function of MgF2 content in the films will be presented. Films containing 8.4 wt.% MgF2 were found to be highly transparent and provided sheet resistivity in the required range. These films maintained a high transmittance upon exposure to AO and to VUV radiation, although exposure to AO in the presence of charged species and intense electromagnetic radiation caused significant degradation in film transmittance. Sheet resistivity of the as-fabricated films increased with time in ambient conditions. Vacuum beat treatment following film deposition caused a reduction in sheet resistivity. However, following vacuum heat treatment, sheet resistivity values remained stable during storage in ambient conditions.

Dever, Joycer A.↗