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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Charge transfer rate constants for the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation

Charge transfer rate constants were calculated for the carotenoid-porphyrin-C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor–acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi’s golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Furthermore, within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC 60 dissolved in tetrahydrofuran.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Quantum membrane phases in synthetic lattices of cold molecules or Rydberg atoms

Here, we calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semisynthetic three-dimensional configuration—one synthetic dimension plus a two-dimensional real-space optical lattice or periodic microtrap array—using the stochastic Green's function quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite-size scalings, we show that there is a second-order transition to a low-temperature phase in which two-dimensional “sheets” form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction V, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large V as well as a quantum phase transition—a critical value V c below which the transition temperature vanishes.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

Spatiotemporal dynamics of fast electron heating in solid-density matter via XFEL

Abstract High-intensity, short-pulse lasers are crucial for generating energetic electrons that produce high-energy-density (HED) states in matter, offering potential applications in igniting dense fusion fuels for fast ignition laser fusion. High-density targets heated by these electrons exhibit spatially non-uniform and highly transient conditions, which have been challenging to characterize due to limitations in diagnostics that provide simultaneous high spatial and temporal resolution. Here, we employ an X-ray Free Electron Laser (XFEL) to achieve spatiotemporally resolved measurements at sub-micron and femtosecond scales on a solid-density copper foil heated by laser-driven fast electrons. Our X-ray transmission imaging reveals the formation of a solid-density hot plasma localized to the laser spot size, surrounded by Fermi degenerate, warm dense matter within a picosecond, and the energy relaxation occurring within the hot plasma over tens of picoseconds. These results validate 2D particle-in-cell simulations incorporating atomic processes and provide insights into the energy transfer mechanisms beyond current simulation capabilities. This work significantly advances our understanding of rapid fast electron heating and energy relaxation in solid-density matter, serving as a key stepping stone towards efficient high-density plasma heating and furthering the fields of HED science and inertial fusion energy research using intense, short-pulse lasers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Assembly and Analysis of Cell-Scale Membrane Envelopes

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms-a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗

DOE OSTI2023

Multiscale simulations of novel lithium electrolytes for improved processability and performance of solid-state batteries

In FY23, we performed atomic-scale simulations and developed analysis tools to understand the dynamic behaviors of Li ions and cation/anion components and predict mechanical properties in a new class of Li-ion solid electrolytes (SEs) developed by TRINA. Specifically, we revealed the transport characteristics of Li ions in anion sublattice, rotational dynamics of anion and cation components, and the effect of cation and anion on Li transport as well as interactions between cation and anions. The results imply that the transport characteristics in these soft materials may be affected by processing conditions significantly, which will offer the necessary insights to optimize these future SEs and accelerate the deployment of practical and easily processable Solid-state batteries (SSBs).

25 ENERGY STORAGE↗

NASA NTRS1994

Computer Modeling Of Atomization

Improved mathematical models based on fundamental principles of conservation of mass, energy, and momentum developed for use in computer simulation of atomization of jets of liquid fuel in rocket engines. Models also used to study atomization in terrestrial applications; prove especially useful in designing improved industrial sprays - humidifier water sprays, chemical process sprays, and sprays of molten metal. Because present improved mathematical models based on first principles, they are minimally dependent on empirical correlations and better able to represent hot-flow conditions that prevail in rocket engines and are too severe to be accessible for detailed experimentation.

Giridharan, M.↗

NASA NTRS1993

Flight- and ground-test correlation study of BMDO SDS materials: Phase 1 report

Chung, Shirley Y.↗

NASA NTRS1995

BMDO materials testing in the EOIM-3 experiment

The NASA Evaluation of Oxygen Interactions with Materials-3 (EOIM-3) experiment served as a testbed for a variety of materials that are candidates for Ballistic Missile Defense Organization (BMDO) space assets. The materials evaluated on this flight experiment were provided by BMDO contractors and technology laboratories. A parallel ground-based exposure evaluation was conducted using the Fast Atom Sample Tester (FAST) atomic-oxygen simulation facility at Physical Sciences, Inc. The EOIM-3 flight materials were exposed to an atomic oxygen fluence of approximately 2.3 x 10(exp 20) atoms/sq cm. The ground-based exposure fluence of 2.0 - 2.5 x 10(exp 20) atoms/sq cm permits direct comparison with that of the flight-exposed specimens. The results from the flight test conducted aboard STS-46 and the correlative ground-based exposure are summarized here. A more detailed correlation study is presented in the JPL Publication 93-31 entitled 'Flight-and Ground-Test Correlation Study of BMDO SDS Materials: Phase 1 Report'. In general, the majority of the materials survived the AO environment with their performance tolerances maintained for the duration of the exposure. Optical materials, baffles, and coatings performed extremely well as did most of the thermal coatings and tribological materials. A few of the candidate radiator, threat shielding, and structural materials showed significant degradation. Many of the coatings designed to protect against AO erosion of sensitive materials performed this function well.

Chung, Shirley Y.↗

DOE OSTI2025

Efficiently Measuring 𝑑-Wave Pairing and Beyond in Quantum Gas Microscopes

Understanding the mechanism of high-temperature superconductivity is among the most important problems in physics, one for which quantum simulation can provide new insights. However, it remains challenging to characterize superconductivity in existing cold-atom quantum simulation platforms. Here, in this work, we introduce a protocol for measuring a broad class of observables in fermionic quantum gas microscopes, including long-range superconducting pairing correlations (after a repulsive-to-attractive mapping). The protocol only requires global controls followed by site-resolved particle number measurements—capabilities that have been already demonstrated in multiple experiments—and is designed by analyzing the Hilbert space of fermions on two sites. The protocol is sample-efficient and we further optimize our pulses for robustness to experimental imperfections such as lattice inhomogeneity. Our Letter introduces a general tool for manipulating quantum states on optical lattices, enhancing their ability to tackle problems such as high-temperature superconductivity.

Hubbard model↗

DOE OSTI2024

Atomic Structure and Dynamics of Unusual and Wide-Gap Phase-Change Chalcogenides: A GeTe 2 Case

Brain-inspired computing, reconfigurable optical metamaterials, photonic tensor cores, and many other advanced applications require next-generation phase-change materials (PCMs) with better energy efficiency and a wider thermal and spectral range for reliable operations. Germanium ditelluride (GeTe 2 ), with higher thermal stability and a larger bandgap compared to current benchmark PCMs, appears promising for THz metasurfaces and the controlled crystallization of atomically thin 2D materials. Using high-energy X-Ray diffraction supported by first-principles simulation, the atomic structure in semiconducting pulsed laser deposition films and metallic high-temperature liquids is investigated. The results suggest that the structural and chemical metastability of GeTe 2 , leading to disproportionation into GeTe and Te, is related to high internal pressure during a semiconductor–metal transition, presumably occurring in the supercooled melt. Similar phenomena are expected for canonical GeS 2 and GeSe 2 under high temperatures and pressures.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

E-beam-enhanced solid-state mechanical amorphization of α-quartz: Reduced deformation barrier via localized excess electrons as network modifiers

Under hydrostatic pressure, α-quartz (α-SiO 2 ) undergoes solid-state mechanical amorphization wherein the interpenetration of [SiO 4 ] 2- tetrahedra occurs and the material loses crystallinity. This phase transformation requires a high hydrostatic pressure of 14GPa because the repulsive forces resulting from the ionic nature of the Si–O bonds prevent the severe distortion of the atomic configuration. Herein, we experimentally and computationally demonstrate that e-beam irradiation changes the nature of the interatomic bonds in α-quartz and enhances the solid-state mechanical amorphization at nanoscale. Specifically, during in situ uniaxial compression, a larger permanent deformation occurs in α-quartz submicron pillars compressed during e-beam irradiation than in those without e-beam irradiation. Microstructural analysis reveals that the large permanent deformation under e-beam irradiation originates from the enhanced mechanical amorphization of α-quartz and the subsequent viscoplastic deformation of the amorphized region. Further, atomic-scale simulations suggest that the delocalized excess electrons introduced by e-beam irradiation move to highly distorted atomic configurations and alleviate the repulsive force, thus reducing the barrier to the solid-state mechanical amorphization. In conclusion, these findings deepen our understanding of electron–matter interactions and can be extended to new glass forming and processing technologies at nano- and microscale.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Atomic oxygen stimulated outgassing

The passive Long Duration Exposure Facility (LDEF) Experiment A0034, Atomic Oxygen Simulated Outgassing, consisted of two identical one-sixth tray modules, exposing selected thermal control coatings to atomic oxygen and the combined space environment on the leading edge and, for reference, to the relative wake environment on the trailing edge. Optical mirrors were included adjacent to the thermal coatings for deposition of outgassing products. Ultraviolet grade windows and metal covers were provided for additional assessment of the effects of the various environmental factors. Preliminary results indicate that orbital atomic oxygen is both a degrading and a optically restorative factor in the thermo-optical properties of selected thermal coatings. There is evidence of more severe optical degradation on collector mirrors adjacent to coatings that were exposed to the RAM-impinging atomic oxygen. This evidence of atomic oxygen stimulated outgassing is discussed in relation to alternative factors that could affect degradation. The general effects of the space environment on the experiment hardware as well as the specimens are discussed.

Linton, Roger C.↗

NASA NTRS1989

Study of indium tin oxide films exposed to atomic axygen

A qualitative simulation of the effects of atomic oxygen has been conducted on indium tin oxide (ITO) films prepared by dc sputtering onto room-temperature substrates, by exposing them to an RF-excited oxygen plasma and characterizing the resulting changes in optical, electrical, and structural properties as functions of exposure time with ellipsometry, spectrophotometry, resistivity, and X-ray measurements. While the films thus exposed exhibit reduced resistivity and optical transmission; both of these effects, as well as partial crystallization of the films, may be due to sample heating by the plasma. Film resistivity is found to stabilize after a period of exposure.

Snyder, Paul G.↗

NASA NTRS1992

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) will assist in understanding the mechanisms involved, and will lead to improved reliability in predicting in-space durability of materials based on ground laboratory testing. A computational simulation of atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of assumed mechanistic behavior of atomic oxygen and results of both ground laboratory and LDEF data, a predictive Monte Carlo model was developed which simulates the oxidation processes that occur on polymers with applied protective coatings that have defects. The use of high atomic oxygen fluence-directed ram LDEF results has enabled mechanistic implications to be made by adjusting Monte Carlo modeling assumptions to match observed results based on scanning electron microscopy. Modeling assumptions, implications, and predictions are presented, along with comparison of observed ground laboratory and LDEF results.

Banks, Bruce A.↗

DOE OSTI2022

Multiscale water drop contact angles at selected silica surfaces

In this paper, multiscale advancing contact angles for glycerol/water drops at silica surfaces are reported for millidrops, submicron-drops, and nanodrops. Selected silica surfaces were muscovite, silicon, and talc. The contact angles for millidrops (1–2 mm) were determined by the traditional sessile drop technique. For submicron-drops (0.1–1.0 μm), a hollow tip Atomic Force Microscope (AFM) procedure was used. The contact angles for nanodrops (~7 nm) were examined from Molecular Dynamics (MD) simulation. The results were compared to evaluate the effect of drop size on the contact angle. In the case of the hydrophobic talc surface, the 75° advancing contact angle did not vary significantly with drop size. For the hydrophilic muscovite surface, the water drop wet the surface and an advancing contact angle of about 10° was found for the millidrops and submicron-drops. However, for the MD simulated nanodrops, attachment and spreading of the ~7 nm drop created a 2D film of molecular dimensions, the contact angle of which was difficult to define and varied from 0° to 17°. Perhaps of equal interest from the MD simulation results was that the spreading of the glycerol/water nanodrop at the muscovite surface resulted in crystallographic directional transport of water molecules to the extremities of the 2D film. Such separation and segregation left the center of the film with an increased concentration of glycerol. Based on these results, the line tension, which has been found in other investigations to account for contact angle decrease with a decrease in drop size, does not seem to be a significant factor in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗