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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Computational chemistry

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Arnold, J. O.↗

DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗

NASA NTRS1988

Oxidation and protection of fiberglass-epoxy composite masts for photovoltaic arrays in the low Earth orbital environment

Fiberglass-epoxy composites are considered for use as structural members for the mast of the space station solar array panel. The low Earth orbital environment in which space station is to operate is composed mainly of atomic oxygen, which has been shown to cause erosion of many organic materials and some metals. Ground based testing in a plasma asher was performed to determine the extent of degradation of fiberglass-epoxy composites when exposed to a simulated atomic oxygen environment. During exposure, the epoxy at the surface of the composite was oxidized, exposing individual glass fibers which could easily be removed. Several methods of protecting the composite were evaluated in an atomic oxygen environment and with thermal cycling and flexing. The protection techniques evaluated to date include an aluminum braid covering, an indium-tin eutectic and a silicone based paint. The open aluminum braid offered little protection while the CV-1144 coating offered some initial protection against atomic oxygen, but appears to develop cracks which accelerate degradation when flexed. Coatings such as the In-Sn eutectic may provide adequate protection by containing the glass fibers even though mass loss still occurs.

Rutledge, Sharon K.↗

DOE OSTI2025

Study of self-assembly behavior and ionic conductivity of conjugated liquid crystals with T-shaped facial-polyphilic structure

The unique self-assembly of liquid crystals (LCs), combined with their potential application as organic semiconductors, has become a focus of recent research. Here, a joint experimental and computational study of the self-assembly and ionic conduction was carried out on a series of T-shaped conjugated LCs consisting of three incompatible components. By extending the EOn side-chain length, several experimental evaluations confirmed a decrease of the order-disorder transition temperature, while coarse-grained simulations revealed a structural evolution from a smectic phase to a columnar phase. Ionic conductivity of these molecules was achieved by adding Li salt, leading to a maximum conductivity of 1.1 × 10 −3 siemens per centimeter observed at 120°C. All-atom simulations were performed to examine the Li-ion solvation environment and to evaluate the intrachain and interchain Li-ion hopping mechanisms. The molecule with a long EOn side chain was found to generate a densely distributed network of Li-ion solvation sites, which can facilitate effective interchain hopping to promote ion transport.

Liu, Ziwei [Cornell Univ., Ithaca, NY (United Stat↗

NASA NTRS1993

Monte Carlo modeling of atomic oxygen attack of polymers with protective coatings on LDEF

Characterization of the behavior of atomic oxygen interaction with materials on the Long Duration Exposure Facility (LDEF) assists in understanding of the mechanisms involved. Thus the reliability of predicting in-space durability of materials based on ground laboratory testing should be improved. A computational model which simulates atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of an assumed mechanistic behavior of atomic oxygen interaction based on in-space atomic oxygen erosion of unprotected polymers and ground laboratory atomic oxygen interaction with protected polymers, prediction of atomic oxygen interaction with protected polymers on LDEF was accomplished. However, the results of these predictions are not consistent with the observed LDEF results at defect sites in protected polymers. Improved agreement between observed LDEF results and predicted Monte Carlo modeling can be achieved by modifying of the atomic oxygen interactive assumptions used in the model. LDEF atomic oxygen undercutting results, modeling assumptions, and implications are presented.

Banks, Bruce A.↗

NASA NTRS1991

Preparation of atomic oxygen resistant polymeric materials

Polyphenyl quinoxalines (PPQs) are an important family of high performance polymers that offer good chemical and thermal stability coupled with excellent mechanical properties. These aromatic heterocyclic polymers are potentially useful as films, coatings, adhesives, and composite materials that demand stability in harsh environments. Our approach was to prepare PPQs with pendent siloxane groups using the appropriate chemistry and then evaluate these polymers before and after exposure to simulated atomic oxygen. Either monomer, the bis(o-diamine)s or the bis(alpha-diketone)s can be synthesized with a hydroxy group to which the siloxane chain will be attached. Several novel materials were prepared.

Tortorelli, Victor J.↗

DOE OSTI2022

Ingrained: An Automated Framework for Fusing Atomic-Scale Image Simulations into Experiments

To fully leverage the power of image simulation to corroborate and explain patterns and structures in atomic resolution microscopy, an initial correspondence between the simulation and experimental image must be established at the outset of further high accuracy simulations or calculations. Furthermore, if simulation is to be used in context of highly automated processes or high-throughput optimization, the process of finding this correspondence itself must be automated. In this work, "ingrained," an open-source automation framework which solves for this correspondence and fuses atomic resolution image simulations into the experimental images to which they correspond, is introduced. Here, the overall "ingrained" workflow, focusing on its application to interface structure approximations, and the development of an experimentally rationalized forward model for scanning tunneling microscopy simulation are described.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Exploring Cr and molten salt interfacial interactions for molten salt applications

Molten salts play an important role in various energy-related applications such as high-temperature heat transfer fluids and reaction media. However, the extreme molten salt environment causes the degradation of materials, raising safety and sustainability challenges. A fundamental understanding of material–molten salt interfacial evolution is needed. This work studies the transformation of metallic Cr in molten 50/50 mol% KCl–MgCl 2 via multi-modal in situ synchrotron X-ray nano-tomography, diffraction and spectroscopy combined with density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. Notably, in addition to the dissolution of Cr in the molten salt to form porous structures, a δ-A15 Cr phase was found to gradually form as a result of the metal–salt interaction. This phase change of Cr is associated with a change in the coordination environment of Cr at the interface. DFT and AIMD simulations provide a basis for understanding the enhanced stability of δ-A15 Cr vs. bcc Cr, by revealing their competitive phase thermodynamics at elevated temperatures and probing the interfacial behavior of the molten salt at relevant facets. This study provides critical insights into the morphological and chemical evolution of metal–molten salt interfaces. Finally, the combination of multimodal synchrotron analysis and atomic simulation also offers an opportunity to explore a broader range of systems critical to energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Photoluminescence, photophysics, and photochemistry of the V B - defect in hexagonal boron nitride

Extensive photochemical and spectroscopic properties of the V B - defect in hexagonal boron nitride are calculated, concluding that the observed photoemission associated with recently observed optically detected magnetic resonance is most likely of (1) 3 E'' → (1) 3 A' 2 origin. Rapid intersystem crossing from the defect's triplet to singlet manifolds explains the observed short excited-state lifetime and very low quantum yield. New experimental results reveal smaller intrinsic spectral bandwidths than previously recognized, interpreted in terms of spectral narrowing and zero-phonon-line shifting induced by the Jahn-Teller effect. Different types of computational methods are applied to map out the complex triplet and singlet defect manifolds, including the doubly ionized formulation of the equation-of-motion coupled-cluster theory that is designed to deal with the open-shell nature of defect states, and mixed quantum-mechanics/molecular-mechanics schemes enabling 5763-atom simulations. Finally, two other energetically feasible spectral assignments from among the singlet and triplet manifolds are considered, but ruled out based on inappropriate photochemical properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Atomistic Explanation of the Dramatically Improved Oxygen Reduction Reaction of Jagged Platinum Nanowires, 50 Times Better than Pt

Pt is the best catalyst for the oxygen reduction reactions (ORRs), but it is far too slow. Huang and co-workers showed that dealloying 5 nm Ni 7 Pt 3 nanowires (NW) led to 2 nm pure Pt jagged NW (J-PtNW) with ORRs 50 times faster than Pt/C. They suggested that the undercoordinated surface Pt atoms, mechanical strain, and high electrochemically active surface area (ECSA) are the main contributors. We report here multiscale atomic simulations that further explain this remarkably accelerated ORR activity from an atomistic perspective. We used the ReaxFF reactive force field to convert the 5 nm Ni 7 Pt 3 NW to the jagged 2 nm NW. We applied quantum mechanics to find that 14.4% of the surface sites are barrierless for O ads + H 2 O ads → 2OH ads , the rate-determining step (RDS). The reason is that the concave nature of many surface sites pushes the OH bond of the H 2 O ads close to the O ads , leading to a dramatically reduced barrier. We used this observation to predict the performance improvement of the J-PtNW relative to Pt (111). Assuming every surface site reacts independently with this predicted rate leads to a 212-fold enhancement at 298.15 K, compared to 50 times experimentally. The atomic structures of the active sites provide insights for designing high-performance electrocatalysts for ORR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Observations of xenon spectra on the EAST x-ray crystal spectrometer for high-temperature plasma diagnostics

The Xe 44+ 2.7203 Å line, which has been proposed as one of the diagnostic lines for the x-ray imaging crystal spectrometer on ITER, is observed on the EAST tokamak together with its several satellite lines. The observations are made under high electron temperature (T e ) conditions (core T e > 5 keV). Most of the observed xenon lines are identified by comparing the experiment results with the atomic simulation results. The first ion temperature measurements made by the xenon spectra on EAST are also reported in this article. These xenon spectra observations contribute to the justification for using xenon as the diagnostic impurity in x-ray crystal spectrometers in future reactor-scale high-temperature plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Modelling Thomson Scattering in a Hydrogen Plasma at Stellar Interior Conditions Using the Hypernetted‐Chain Approach

Under the extreme conditions found in small stars, where electron degeneracy and Coulomb coupling are significant, accurate modeling of Thomson scattering is crucial for determining opacity, a primary quantity for stellar energy transport. We use hypernetted-chain calculations, incorporating quantum pseudopotentials and electron-exchange effects to obtain the electron–electron static structure factor to calculate the Thomson scattering transport cross-section for conditions prevailing in the interior of small stars. These results are compared to those from average-atom simulations and analytical calculations. Our findings support laboratory astrophysics experiments aimed at benchmarking opacity models for stellar interiors, particularly for red dwarf stars, and help to bridge theoretical models with observations.

ICF↗

DOE OSTI2021

Normal and anomalous self-healing mechanism of crystalline calcium silicate hydrates

The origin of different stability of crystalline calcium silicate hydrates was investigated. The tobermorite crystal has been used as an analog of cement hydrate that is being mostly manufactured material on earth. Normal tobermorite is thermally unstable and transforms to amorphous at low pressure. Meanwhile, anomalous tobermorite with high Al content does not significantly transform under high pressure or high temperature. Conducted X-ray absorption spectroscopy explains the weak stability of normal tobermorite which was originally hypothesized by the role of zeolitic Ca ions in the cavities of silicate chains. Atomic simulations reproduced the experimentally observed trend of pressure behavior once the ideal structures were modified to account for the Al content as well as the chain defects. The simulations also suggested that the stability of tobermorite under stress could be rationalized as a self-healing mechanism in which the structural instabilities were accommodated by a global sliding of the CaO layer.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Integrating Machine Learning Potential and X-ray Absorption Spectroscopy for Predicting the Chemical Speciation of Disordered Carbon Nitrides

Precise determination of atomic structural information in functional materials holds transformative potential and broad implications for emerging technologies. Spectroscopic techniques, such as X-ray absorption near-edge structure (XANES), have been widely used for material characterization; however, extracting chemical information from experimental probes remains a significant challenge, particularly for disordered materials. We present an integrated approach that combines atomic simulations, data-driven techniques, and experimental measurements to investigate chemical speciation of amorphous carbon nitride systems as a case study. Here, we discuss the development of machine learning potentials that can efficiently explore the vast configuration space of amorphous carbon nitrides. By employing statistical methods, this structural database enables the elucidation of the most representative local structures and how they evolve with chemical compositions and density. Density functional theory simulations are used to establish a correlation between the local structure and spectroscopic signatures, which then serve as the basis for interpreting and extracting chemical content from experimental data. Although our framework is specifically demonstrated for XANES and carbon nitrides, the approach described herein is readily adaptable as applied to other experimental characterization probes and materials classes.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

In this work, we have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Å 2 , a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Å2). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S CD ) profiles calculated from the simulations were in relatively good agreement with S CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Å 2 resulted in good agreement between the simulation and experimental scattering form factors, but not with S CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗