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 955 records · Page 53

Monte Carlo simulation of energy deposition by low-energy electrons in molecular hydrogen

A set of detailed atomic cross sections has been used to obtain the spatial deposition of energy by 1-20-eV electrons in molecular hydrogen by a Monte Carlo simulation of the actual trajectories. The energy deposition curve (energy per distance traversed) is quite peaked in the forward direction about the entry point for electrons with energies above the threshold of the electronic states, but the peak decreases and broadens noticeably as the electron energy decreases below 10 eV (threshold for the lowest excitable electronic state of H2). The curve also assumes a very symmetrical shape for energies below 10 eV, indicating the increasing importance of elastic collisions in determining the shape of the curve, although not the mode of energy deposition.

Heaps, M. G.↗

Determination of atomic oxygen fluence using spectrophotometric analysis of infrared transparent witness coupons for long duration exposure tests

Atomic oxygen degradation is one of several major threats to the durability of spaceborne systems in low Earth orbit. Ground-based simulations are conducted to learn how to minimize the adverse effects of atomic oxygen exposure. Assessing the fluence of atomic oxygen in test chambers such as a plasma asher over long periods of time is necessary for accurate determination of atomic oxygen exposure. Currently, an atomic oxygen susceptible organic material such as Kapton is placed next to samples as a witness coupon and its mass loss is monitored and used to determine the effective atomic oxygen fluence. However, degradation of the Kapton witness coupons occurs so rapidly in plasma ashers that for any long term test many witness coupons must be used sequentially in order to keep track of the fluence. This necessitates opening vacuum to substitute fresh coupons. A passive dosimetry technique was sought to monitor atomic oxygen exposure over longer periods without the need to open the plasma asher to the atmosphere. This paper investigates the use of spectrophotometric analysis of durable IR transparent witness coupons to measure atomic oxygen exposure for longer duration testing. The method considered would be conductive to making in situ measurements of atomic oxygen fluence.

Podojil, Gregg M.↗

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atomic-Scale Investigation of Oxidation at the Black Phosphorus Surface

Black phosphorus (BP) exhibits extraordinary electronic properties that are desirable for a wide variety of electronic and optoelectronic applications. However, applications of BP are hindered by its rapid degradation in ambient conditions. Despite significant advances that have been made in understanding the degradation mechanism, no consensus has yet been reached on how BP oxidation occurs at the atomic scale as experimental studies have been mostly restricted to averaged effects of degradation over a micron- to millimeter-sized region. Here, BP oxidation is investigated using scanning tunneling microscopy/spectroscopy (STM/S). Introducing O 2 gas to the BP surface in ultrahigh vacuum at a pressure of 10 –5 mbar for 1 min creates two new types of defects on the surface. We identify these defects as dangling atomic oxygen and phosphorus multivacancies using density functional theory simulations. In addition to the structural changes to the surface, the electronic structure is also drastically altered by the introduction of oxygen. The 300 meV band gap of BP is lifted due to dosing. This change in the electronic structure is reversible through STM tip manipulation. As a result, these are the first experimental results showing the atomic-scale oxidation of BP, an important step toward understanding the degradation process

2D materials↗

Quantum simulation of an extended Dicke model with a magnetic solid

Abstract The Dicke model describes the cooperative interaction of an ensemble of two-level atoms with a single-mode photonic field and exhibits a quantum phase transition as a function of light–matter coupling strength. Extending this model by incorporating short-range atom–atom interactions makes the problem intractable but is expected to produce new physical phenomena and phases. Here, we simulate such an extended Dicke model using a crystal of ErFeO 3 , where the role of atoms (photons) is played by Er 3+ spins (Fe 3+ magnons). Through terahertz spectroscopy and magnetocaloric effect measurements as a function of temperature and magnetic field, we demonstrated the existence of a novel atomically ordered phase in addition to the superradiant and normal phases that are expected from the standard Dicke model. Further, we elucidated the nature of the phase boundaries in the temperature–magnetic-field phase diagram, identifying both first-order and second-order phase transitions. These results lay the foundation for studying multiatomic quantum optics models using well-characterized many-body solid-state systems.

Marquez Peraca, Nicolas↗

Probing High-Pressure Structural Evolution in Polyurea with In Situ Energy-Dispersive X-ray Diffraction and Molecular Dynamics Simulations

Polyurea, an elastomer with a phase-segregated microstructure, has been proven as an effective coating in defense applications. To gain a more complete understanding of the high-pressure atomic-level morphology of these phases and to validate molecular dynamics (MD) simulations, multi-angle energy-dispersive X-ray diffraction experiments were performed in situ up to pressures of ~6 GPa at room temperature. Structure factors were obtained and compared to MD simulations with an average error of less than 5% between major peak positions. The first sharp diffraction peak shifted from 4.56 Å to lower d-spacing with pressure, indicating compression between hard segments. This was further supported by the behavior of a peak at ~3.86 Å from the pair distribution function (PDF), suspected to represent π-stacking and separation between soft segments. Compression within the hard segments themselves is minimal as low-r peaks in the PDF are not greatly affected by pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aerodynamic-wave break-up of liquid sheets in swirling airflows and combustor modules

Experimental mean drop diameter data were obtained for the atomization of liquid sheets injected axially downstream in high velocity swirling and nonswirling airflow. Conventional simplex pressure atomizing fuel nozzles and splash type fuel injectors were studied under simulated combustor inlet airflow conditions. A general empirical expression relating recirprocal mean drop diameter to airstream mass velocity was obtained and is presented. The finest degree of atomization, i.e., the highest value of the coefficient C, was obtained with swirl can combustor modules (C = 15) as compared with pressure atomizing nozzles (C = 12).

Ingebo, R.↗

Aerodynamic-wave break-up of liquid sheets in swirling airflows and combustor modules

Experimental mean drop diameter data were obtained for the atomization of liquid sheets injected axially downstream in high velocity swirling and nonswirling airflow. Conventional simplex pressure atomizing fuel nozzles and splash type fuel injectors were studied under simulated combustor inlet airflow conditions. A general empirical expression relating reciprocal mean drop diameter to airstream mass velocity was obtained and is presented. The finest degree of atomization, i.e., the highest value of the coefficient C, was obtained with swirl can combustor modules (C = 15) as compared with pressure atomizing nozzles (C = 12). Previously announced in STAR as N83-23545

Ingebo, R.↗

Optical Time Transfer for Bistatic SAR Spacecraft

A spacecraft-to-spacecraft optical time-transfer simulation has been developed as a tool for informing NASA’s Surface Deformation and Change (SDC) mission architecture. The SDC mission will combine radar images from multiple spacecraft to improve understanding of the Earth’s sea-level and landscape changes. Spacecraft must be precisely synchronized in order to create sharp radar images. Simulation of multiple spacecraft time-synchronizing via laser communication can inform technology choices of a mission by providing a picosecond-precision level estimate of clock error. This timing and ranging simulation has been combined with a radar system performance analysis pipeline. The simulated timing errors are used to predict performance of bistatic SAR systems in the presence of oscillator noise and time synchronization in accuracy. This analysis includes both analytic approximation equations from existing literature, and a numerical radar simulation to extract key system performance parameters like phase error and signal-to-noise ratio (SNR)degradation. Precision time-transfer techniques facilitate the accurate synchronization of clocks between any combination of terminals. Most time-transfer technology for comparing two clocks at different terminals use radio frequencies (RF) to measure the time delay between the sending and receiving of signals. Laser technology offers the capability to transmit high data rates with systems that are of smaller size and lower power than comparable RF systems. The clocks on independent spacecraft will have some phase and frequency errors between them that result in clock drift. The two clock models that are included in this bi-directional MATLAB simulation are a cesium-based Chip-Scale Atomic Clock (CSAC) and a rubidium-based Miniature Atomic Clock (MAC). The CSAC has flown as hardware for small satellite missions such as the University of Florida’s CHOMPTT mission. A study of example orbits, including that of NASA NASA-ISRO Synthetic Aperture Radar Mission (NISAR) mission, and lasing rates demonstrate the impact of flight configuration parameters on the synchronization error between two spacecraft. The MATLAB timing simulation uses a Runge-Kutta 4th-order method to propagate spacecraft orbits and computes the light-travel time estimate between them. The simulation outputs the estimated range and estimated clock error based on a user-defined spacecraft cluster configuration. The radar simulation and analytic approximations are applied to evaluate a potential future NASA bistatic SAR constellation architecture. In the proposed architecture, satellites follow each other in the same orbit at 800 km altitude, with a 210 km baseline. We also baseline the CSAC as an ultra stable oscillator, and use NASA’s NISAR for baseline radar system parameters to compute a clock-system introduced phase error of 5.6 degrees without synchronization by frequent time transfer. We build on this base case with a sensitivity analysis of radar performance over a proposed range of constellation and radar system parameters. With this analysis pipeline, we comment on which radar parameters should or should not be changed to minimize synchronization requirements. This analysis technique could be extended or modified to evaluate the timing requirements of other geometries for other future multistatic SAR missions, or other interferometric satellite missions.

Surface Deformation and Change↗

Learning local equivariant representations for large-scale atomistic dynamics

Abstract A simultaneously accurate and computationally efficient parametrization of the potential energy surface of molecules and materials is a long-standing goal in the natural sciences. While atom-centered message passing neural networks (MPNNs) have shown remarkable accuracy, their information propagation has limited the accessible length-scales. Local methods, conversely, scale to large simulations but have suffered from inferior accuracy. This work introduces Allegro, a strictly local equivariant deep neural network interatomic potential architecture that simultaneously exhibits excellent accuracy and scalability. Allegro represents a many-body potential using iterated tensor products of learned equivariant representations without atom-centered message passing. Allegro obtains improvements over state-of-the-art methods on QM9 and revMD17. A single tensor product layer outperforms existing deep MPNNs and transformers on QM9. Furthermore, Allegro displays remarkable generalization to out-of-distribution data. Molecular simulations using Allegro recover structural and kinetic properties of an amorphous electrolyte in excellent agreement with ab-initio simulations. Finally, we demonstrate parallelization with a simulation of 100 million atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

Strongly Anharmonic Phonons and Their Role in Superionic Diffusion and Ultralow Thermal Conductivity of Cu 7 PSe 6

We report the quest for advanced superionic materials requires understanding their complex atomic dynamics, but detailed studies of the interplay between lattice vibrations and ionic diffusion remain scarce. Here inelastic and quasielastic neutron scattering measurements in the superionic argyrodite Cu 7 PSe 6 are reported, combined with molecular dynamics (MD) based on ab initio and machine-learned potentials, providing critical insights into the atomistic mechanisms underlying fast ion conduction. The results reveal how long-range Cu diffusion is limited by intercluster hopping, controlled by selective anharmonic phonons of the crystalline framework. Further, the Green–Kubo simulations reproduce the ultralow lattice thermal conductivity and identify contributions from mobile ions, phonons, and their cross-correlations. The mode resolved analysis shows that the thermal conductivity is dominated by low-energy acoustic phonon modes of the overall crystal framework. The analysis of mode-resolved spectral functions further show that vibrational modes with significant Cu contributions are strongly damped, corresponding to the breakdown of associated phonon quasiparticles. These results highlight the importance of strongly anharmonic effects in superionic systems, in which the traditional quasiharmonic phonon picture is insufficient, and pave the way toward combining machine-learning accelerated simulations with neutron scattering experiments to rationalize the complex atomic dynamics underlying ionic and thermal transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗

Exploring model complexity in machine learned potentials for simulated properties

Abstract Machine learning (ML) enables the development of interatomic potentials with the accuracy of first principles methods while retaining the speed and parallel efficiency of empirical potentials. While ML potentials traditionally use atom-centered descriptors as inputs, different models such as linear regression and neural networks map descriptors to atomic energies and forces. This begs the question: what is the improvement in accuracy due to model complexity irrespective of descriptors? We curate three datasets to investigate this question in terms of ab initio energy and force errors: (1) solid and liquid silicon, (2) gallium nitride, and (3) the superionic conductor Li $$_{10}$$ 10 Ge(PS $$_{6}$$ 6 ) $$_{2}$$ 2 (LGPS). We further investigate how these errors affect simulated properties and verify if the improvement in fitting errors corresponds to measurable improvement in property prediction. By assessing different models, we observe correlations between fitting quantity (e.g. atomic force) error and simulated property error with respect to ab initio values. Graphical abstract

Rohskopf, A. (ORCID:0000000227128296)↗

Atomic Oxygen Effects on Spacecraft Materials

Low Earth orbital (LEO) atomic oxygen cannot only erode the external surfaces of polymers on spacecraft, but can cause degradation of surfaces internal to components on the spacecraft where openings to the space environment exist. Although atomic oxygen attack on internal or interior surfaces may not have direct exposure to the LEO atomic oxygen flux, scattered impingement can have can have serious degradation effects where sensitive interior surfaces are present. The effects of atomic oxygen erosion of polymers interior to an aperture on a spacecraft is simulated using Monte Carlo computational techniques. A 2-dimensional model is used to provide quantitative indications of the attenuation of atomic oxygen flux as a function of distance into a parallel walled cavity. The degree of erosion relative is compared between the various interior locations and the external surface of an LEO spacecraft.

Banks, Bruce A.↗

Scattered Atomic Oxygen Effects on Spacecraft Materials

Low Earth orbital (LEO) atomic oxygen cannot only erode the external surfaces of polymers on spacecraft, but can cause degradation of surfaces internal to components on the spacecraft where openings to the space environment exist. Although atomic oxygen attack on internal or interior surfaces may not have direct exposure to the LEO atomic oxygen flux scattered impingement can have serious degradation effects where sensitive interior surfaces are present. The effects of atomic oxygen erosion of polymer interior to an aperture on a spacecraft is simulated using Monte Carlo computational techniques. A 2-dimensional model is used to provide quantitative indications of the attenuation of atomic oxygen flux as a function of distance into a parallel walled cavity. The degree of erosion re1ative is compared between the various interior locations and the external surface of a LEO spacecraft.

Banks, Bruce A.↗

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sputtering Erosion in the Ion Thruster

During the first phase of this research, the sputtering yields of molybdenum by low energy (100 eV and higher) xenon ions were measured by using the methods of secondary neutral mass spectrometry (SNMS) and Rutherford backscattering spectrometry (RBS). However, the measured sputtering yields were found to be far too low to explain the sputtering erosions observed in the long-duration tests of ion thrusters. The only difference between the sputtering yield measurement experiments and the ion thruster tests was that the later are conducted at high ion fluences. Hence, a study was initiated to investigate if any linkage exists between high ion fluence and an enhanced sputtering yield. The objective of this research is to gain an understanding of the causes of the discrepancies between the sputtering rates of molybdenum grids in an ion thruster and those measured from our experiments. We are developing a molecular dynamics simulation technique for studying low-energy xenon ion interactions with molybdenum. It is difficult to determine collision sequences analytically for primary ions below the 200 eV energy range where the ion energy is too low to be able to employ a random cascade model with confidence and it is too high to have to consider only single collision at or near the surface. At these low energies, the range of primary ions is about 1 to 2 nm from the surface and it takes less than 4 collisions on the average to get an ion to degrade to such an energy that it can no longer migrate. The fine details of atomic motion during the sputtering process are revealed through computer simulation schemes. By using an appropriate interatomic potential, the positions and velocities of the incident ion together with a sufficient number of target atoms are determined in small time steps. Hence, it allows one to study the evolution of damages in the target and its effect on the sputtering yield. We are at the preliminary stages of setting up the simulation program.

Ray, Pradosh K.↗

Protective coatings for composite tubes in space applications

Protective coatings for graphite/epoxy (Gr/Ep) tubular structures for a Manned Space Station truss structure were evaluated. The success of the composite tube truss structure depends on its stability to long-term exposure to the Low Earth Orbit (LEO) environment with particular emphasis placed on atomic oxygen. Concepts for protectively coating Gr/Ep tubes include use of inorganic coated metal foils and electroplating. These coatings were applied to Gr/Ep tubes and then subjected to simulated LEO environmnet to evaluate survivability of coatings and coated tubes. Evaluation included: atomic oxygen resistance, changes in optical properties and adhesion, abrasion resistancem surface preparation required, coating uniformity, and formation of microcracks in the Gr/Ep tubes caused by thermal cycling. Program results demonstrated that both phosphoric and chromic acid anodized Al foil provided excellent adhesion to Gr/Ep tubes and exhibited stable optical properties when subjected to simulated LEO environment. The SiO2/Al coatings speuttered onto Al foils also resulted in an excellent protective coating. Electroplated Ni exhibited unaccepatble adhesion loss to Gr/Ep tubes during atomic oxygen exposure.

Dursch, Harry W.↗