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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Ground radiation tests and flight atomic oxygen tests of ITO protective coatings for Galileo Spacecraft

Radiation simulation tests (protons and electrons) were performed along with atomic oxygen flight tests aboard the Shuttle to space qualify the surface protective coatings. The results, which contributed to the selection of indium-tin-oxide (ITO) coated polyester as the material for the thermal blankets of the Galileo Spacecraft, are given here. Two candidate materials, polyester and Fluorglas, were radiation-tested to determine changes at simulated Jovian radiation levels. The polyester exhibited a smaller weight loss (2.8) than the Fluorglas (8.8 percent). Other changes of polyester are given. During low-earth orbit, prior to transit to Jupiter, the thermal blankets would be exposed to atomic oxygen. Samples of uncoated and ITO-coated polyesters were flown on the Shuttle. Qualitative results are given which indicated that the ITO coating protected the underlying polyester.

Bouquet, Frank L.↗
DOE OSTI2026

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The Importance of Surface Oxygen for Lithiation and Morphology Evolution during Calcination of High-Nickel NMC Cathodes

Nanoscale morphology has a direct impact on the performance of materials for electrochemical energy storage. Despite this importance, little is known about the evolution of primary particle morphology nor its effect on chemical pathways during synthesis. In this work, operando characterization is combined with atomic-scale and continuum simulations to clarify the relationship between morphology of cathode primary particles and their lithiation during calcination of LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC-811). This combined approach reveals a key role for surface oxygen adsorption in facilitating the lithiation reaction by promoting metal diffusion and oxidation, and simultaneously providing surface sites for lithium insertion. Furthermore, oxygen surface termination is shown to increase the activation energy for sintering, leading to smaller primary particle sizes at intermediate temperatures. Smaller particles provide both shorter diffusion lengths for lithium incorporation and increased surface site density for lithium insertion. These insights provide a foundation for more tailored syntheses of cathode materials with optimized performance characteristics.

25 ENERGY STORAGE↗
DOE OSTI2023

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure–function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Finally, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Novel effects of grain size and ion implantation on grain boundary segregation in ion irradiated austenitic steel

Nuclear reactor environments are extreme in nature due to the combination of exposure to corrosive coolant, mechanical stresses, and neutron irradiation damage. This combination of environments can lead to irradiation induced/assisted precipitation, segregation, and embrittlement coupled with enhanced creep and stress corrosion cracking. While nanostructured materials are of interest in a variety of applications due to their unique properties including high mechanical strength, for nuclear applications nanostructured materials have the added benefit of irradiation resistance due to the high volume density of sinks for irradiation-induced defects. In this study, we show unique radiation-induced segregation (RIS) behavior in nanocrystalline 304 stainless steel (SS304). Typical RIS in austenitic steels includes Ni enrichment and Cr depletion within the grain boundaries. In nanostructured SS304 this RIS is reduced. Furthermore, interestingly, some grain boundaries in the nanocrystalline SS304 in this study were found to be enriched in Cr after Fe 2+ ion irradiation at 500 °C up to 50 displacements per atom (dpa). To understand this unique observation, lattice-based atomic kinetic Monte Carlo simulations were performed to investigate the effects of grain size, thermal segregation, as well as injected interstitials (introduced as Fe self-ions). The results indicate that the Cr enrichment at grain boundaries may be caused by combined effects of grain size, self-ion injection and preferential Cr diffusion via the interstitial mechanism. In conclusion, the findings also suggest a possible difference in RIS between ion and neutron irradiations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fluctuation–dissipation relation for a quantum Brownian oscillator in a parametrically squeezed thermal field

In this paper we study the nonequilibrium evolution of a quantum Brownian oscillator, modeling the internal degree of freedom of a harmonic atom or an Unruh–DeWitt detector, coupled to a nonequilibrium and nonstationary quantum field bath and inquire whether a fluctuation–dissipation relation (FDR) can exist after/if it approaches equilibration. This is a nontrivial issue because a squeezed field bath cannot reach equilibration and yet, as this work shows, the system oscillator indeed can, which is a necessary condition for FDRs. We discuss three different settings: (A) The bath field essentially remains in a squeezed thermal state throughout, whose squeeze parameter is a mode- and time-independent constant. This situation is often encountered in quantum optics and quantum thermodynamics. (B) The bath field is initially in a thermal state, but is subjected to a parametric process leading to mode- and time-dependent squeezing. This scenario is encountered in cosmology and dynamical Casimir effects. The squeezing in the bath in both types of processes will affect the oscillator’s nonequilibrium evolution. We show that at late times it approaches equilibration and this stationarity condition warrants the existence of a FDR. The trait of squeezing is marked by the oscillator’s effective equilibrium temperature, and the proportionality factor in the FDR is only related to the stationary component of the noise kernel of the bath field. Setting (C) is more subtle: A finite system–bath coupling strength can set the oscillator in a squeezed state even though the bath field is stationary and does not engage in any parametric process. The squeezing of the system in this case is in general time-dependent but becomes constant when the internal dynamics is fully relaxed. We begin with comments on the broad range of physical processes involving squeezed thermal baths and end with some remarks on the significance of FDRs in capturing the essence of quantum backreaction in nonequilibrium and stochastic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic X-ray scattering (IXS) measurements — a technique for probing vibrations in soft and biological materials — are now possible with sub-meV resolution, permitting direct observation of low energy, phonon-like modes in lipid membranes. Here, in this study, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (L o ) and liquid-disordered (L d ) phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm -1 in the L o system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represnt the fast, collective dynamics in the L o phase. It’s optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provide molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on molecular dynamics simulations and later by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the L o substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular-Level Behavior of Imidazolium-Based Ionic Liquid Mixtures

Understanding the molecular-level organization of anions and cations in a mixture of different ionic liquids (ILs) is challenging, but it can be very important for designing solvents and IL-based membrane materials. In this work, we perform all-atom molecular dynamics (MD) simulations to understand the intermolecular interactions and structuring of both neat and mixed IL compositions containing 1-benzyl-3-methylimidazolium [ZMIM + ], 1-butyl-3-methylimidazolium [BMIM + ], bistriflimide [Tf 2 N – ], and tetrafluoroborate [BF 4 – ]. We observe very strong underlying molecular-level migrations in these mixtures, and this tends to mainly involve the anion species. As compared to the systems with a single anion, the mixed anions lead to large increases/decreases in the radial distribution function peak heights, whereas the cation mixing is relatively benign. The [BF 4 – ] anions have a strong tendency to displace the [Tf 2 N – ] anions in the mixtures, and this is seen with respect to both [BMIM + ] and [ZMIM + ].

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Ab initio molecular dynamics simulation of Nd3+ incorporation in calcite

Trivalent rare-earth elements (REEs) readily incorporate in Ca-bearing minerals such as calcite (CaCO3), a phenomenon with important implications for the use of REEs as biogeochemical tracers as well as for environmental remediation. Despite substantial efforts, determining the incorporation modes of trivalent REEs in calcite has remained challenging due to their presence in multiple coordination environments. This is particularly problematic for the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. In this work, ab initio molecular dynamics (AIMD) simulations of Nd3+ incorporated in calcite were performed to resolve its incorporation modes. Nd3+ and Ca2+ have very similar ionic radii but the additional positive charge requires a charge compensation scheme (CCS). Five CCSs were considered in AIMD simulations, and the resulting atomic trajectories were then used to compute EXAFS spectra for direct fitting to a published experimental spectrum. The results of the fits indicated that Nd3+ incorporates at the Ca2+ site in calcite in a mixture of six- and seven-fold coordination environments and through three main incorporation modes: associated with a protonated nearest-neighbor Ca2+ vacancy, incorporated as NdOH2+, and charge compensated by a remote species (on the spatial scale probed by EXAFS). The incorporation of Nd3+ in calcite via multiple modes was consistent with previous work, but the AIMD-EXAFS approach allowed for resolving the nature, proportion, and structure of the different coordination environments.

Kerisit, Sebastien N.↗
DOE OSTI2025

Corrosion resistant hot melt adhesive to bind metals

Hot melt adhesives (HMAs) play an important role in many industries, and their demand is expected to grow. HMAs don't require any solvents, and their application results in the formation of strong bonds with the substrate upon cooling within seconds. These properties differentiate them from liquid glues and make them preferable for practical application. Currently, commercial HMAs are used in bonding lightweight materials such as paper, polymers, and cartons, and have limited usage in areas necessitating the bonding of heavier objects like metals. Here, in this study, we report a design and testing of versatile platform comprising an ion-coordinating polymer and ionic fillers for performance optimization and understanding of structure–property relationships, enabling the rational design of HMAs with improved adhesion to metal surfaces. All-atom Molecular dynamics (MD) simulations and various characterization methods are used to elucidate the adhesion mechanism in model composite system containing polyethylene oxide mixed with chemically diverse salt particles. The maximum adhesion strength is found in composites with Al(OH) 3 and FeCl 3 fillers. Interestingly, the presence of Al(OH) 3 also provides a multifunctional anticorrosion property as measured electrochemically using the Tafel method. The discovered path to formulations with improved adhesion to metal surfaces constitutes an important step toward advancing HMAs for use in the structural and semi-structural metal work domain.

Adhesive and cohesive forces↗
DOE OSTI2020

Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al

Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An improved semi-global intrinsic kinetics model for high temperature carbon oxidation

Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO 2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. Here, in this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N 2 and CO 2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗
DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗