SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Effect of solar proton events in 1978 and 1979 on the odd nitrogen abundance in the middle atmosphere

Daily average solar proton flux data for 1978 and 1979 are used in a proton energy degradation scheme to derive ion pair production rates and atomic nitrogen production rates. The latter are computed in a form suitable for inclusion in an atmopheric, two-dimensional, time-dependent photochemical model. Odd nitrogen distributions are computed from the model, including atomic nitrogen production from solar protons, and are compared with baseline distributions. The comparisons show that the average effect of the solar protons in 1978 and 1979 was to cause changes in odd nitrogen only above 10 mbar and at latitudes only above about 50 deg in both hemispheres. The influence of the solar proton-produced odd nitrogen on the local abundance of odd nitrogen depends primarily on the background odd nitrogen abundance as well as the altitude and season.

Jackman, Charles H.↗

DOE OSTI2026

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗

NASA NTRS1991

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields, volume 1

The following subject areas are covered: the development of detailed nonequilibrium radiation models for molecules along with appropriate models for atoms; the inclusion of nongray radiation gasdynamic coupling in the VSL (Viscous Shock Layer) code; the development and evaluation of various electron-electronic energy models; and an examination of the effects of shock slip.

Carlson, Leland A.↗

DOE OSTI2018

In Silico Prediction of Structural Properties of a Racemic Porous Organic Cage Crystal

Porous organic cage (POC) solids are porous materials made up of individual porous molecules held together by noncovalent forces. While many examples of POCs have been made in homochiral forms (i.e., crystals in which every molecule has the same chirality), crystallization of POCs that contain mixtures of chiral molecules can also yield useful properties. For example, our previous work has demonstrated that a POC crystal [cage crystal 3 (CC3)-racemic] synthesized using a racemic mixture of diaminocyclohexane has improved sorption properties and stability than a similar POC crystal synthesized with homochiral diaminocyclohexane. The main purpose of this paper is to predict the structure of CC3-racemic with atomic detail because it is challenging to fully resolve the structure of this racemic crystal experimentally since only subtle differences exist between the racemic structure and the known crystal structure of homochiral CC3. Here, we introduce an in silico prediction method that combines electronic structure calculations and atomistic calculations to predict the structure of CC3-racemic. We first enumerate types of cage molecules that can be present in this material and establish their concentrations. A key observation from these calculations is that CC3-racemic is not made up of only CC3-R and CC3-S molecules and that an additional class of heterochiral cages is also present. By studying the packing energy of cage pairs in CC3-racemic, a lattice model representation of the racemic crystal is developed and used in MC simulations to assess the structure of extended crystals. By expanding the lattice model into atomic detail, fully detailed CC3-racemic crystal models are obtained. These models provide the most thorough description available to date of the composition and cage packing of this interesting material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Toward machine learning interatomic potentials for modeling uranium mononitride

Uranium mononitride (UN) is a promising accident-tolerant fuel because of its high fissile density and high thermal conductivity. In this study, we developed the first machine learning interatomic potentials for reliable atomic-scale modeling of UN at finite temperatures. We constructed a training set using density functional theory (DFT) calculations that was enriched through an active learning procedure, and two neural network potentials were generated. Both potentials successfully reproduce key thermophysical properties of interest, such as temperature-dependent lattice parameter, specific heat capacity, and bulk modulus. We also evaluated the energy of stoichiometric defect reactions and defect migration barriers and found close agreement with DFT predictions, demonstrating that our potentials can be used for modeling defects in UN. Additional tests provide evidence that our potentials are reliable for simulating diffusion, noble gas impurities, and radiation damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2013

Intrinsic to extrinsic phonon lifetime transition in a GaAs–AlAs superlattice

We have measured here the lifetimes of two zone-center longitudinal acoustic phonon modes, at 320 and 640 GHz, in a 14 nm GaAs/2 nm AlAs superlattice structure. By comparing measurements at 296 and 79 K we separate the intrinsic contribution to phonon lifetime determined by phonon–phonon scattering from the extrinsic contribution due to defects and interface roughness. At 296 K, the 320 GHz phonon lifetime has approximately equal contributions from intrinsic and extrinsic scattering, whilst at 640 GHz it is dominated by extrinsic effects. These measurements are compared with intrinsic and extrinsic scattering rates in the superlattice obtained from first-principles lattice dynamics calculations. The calculated room-temperature intrinsic lifetime of longitudinal phonons at 320 GHz is in agreement with the experimentally measured value of 0.9 ns. The model correctly predicts the transition from predominantly intrinsic to predominantly extrinsic scattering; however the predicted transition occurs at higher frequencies. Our analysis indicates that the 'interfacial atomic disorder' model is not entirely adequate and that the observed frequency dependence of the extrinsic scattering rate is likely to be determined by a finite correlation length of interface roughness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more stable bio-oil product. HDO involves co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Reducible metal oxide-supported noble-metal catalysts (e.g., Pt/TiO2) are promising materials for HDO chemistry, with recent work aiming to elucidate the role of various Pt/TiO2 actives sites (i.e., Pt-metal, TiO2-support, Pt-TiO2-interfacial sites) in the competing desired deoxygenation and undesired decarboxylation/decarbonylation reactions for important classes of CFP vapor model compounds. Carboxylic acids are one important class of bio-derived compounds that has not been studied extensively for HDO on Pt/TiO2, particularly at the atomic level. Past experimental work evaluating Pt/C and Pt/TiO2 catalysts for acetic acid HDO (AA-HDO) demonstrated selectivity toward C-C and C-O bond-dissociation products, respectively. This work utilizes atomic-scale modeling to discern the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in AA-HDO. Using Pt(111) and Pt-TiO2-interface surface models, adsorption and reaction energetics calculated by density functional theory provide fundamental insights into the role of interfacial sites and oxygen vacancies in promoting desired deoxygenation pathways over undesired decarboxylation/decarbonylation pathways.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

DOE OSTI2025

Aggressively Dissipative Dark Dwarfs: The Effects of Atomic Dark Matter on the Inner Densities of Isolated Dwarf Galaxies

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.

cold dark matter↗

DOE OSTI2018

Structural features of ISG borosilicate nuclear waste glasses revealed from high-energy X-ray diffraction and molecular dynamics simulations

Vitrification is a widely accepted method to immobilize nuclear waste. Detailed structural information is critical to understand the physical and chemical behaviors, including the long-term chemical durability of these glasses that are to be stored in geological storage sites. High-energy X-ray diffraction studies are used to obtain accurate structural information of the International Simple Glass (ISG), a six-component borosilicate model nuclear waste glass with a composition (mol%) of 60.2SiO 2 -16.0B 2 O 3 -12.6Na 2 O-3.8Al 2 O 3 -5.7CaO-1.7ZrO 2 . Classical molecular dynamics (MD) simulations with recently developed effective potentials are utilized to generate structural models of the ISG glasses to provide interpretation of the high-energy X-ray structural data and additional short and medium range structural details. The atomic structure model generated from MD simulations shows an excellent agreement with high-energy X-ray diffraction, measured for the first time for ISG, with an R x value of 5.5%. Systematic investigations of the medium-range structural features of ISG are performed, built upon recent study of the short-range structure. In particular, the glass former cation oxygen polyhedral analysis show that [SiO 4 ]-[SiO 4 ] and [BO 4 ]-[SiO 4 ] connections are the most abundant while [ZrO 6 ] octohedra having higher probability to connect to another [SiO 4 ], [BO 4 ] and [BO 3 ]. Moreover, we have studied the ISG surface structure using MD simulations. Noticeable differences in terms of both chemical compositions and structural features between the bulk and surface are observed. There exist increased small-sized rings, under coordinated Si species, high 3-coordinated boron content (86.5%) on the surface. Very importantly, sodium ions are found to be enriched on the glass surface, consistent with recent experimental results. Furthermore, both bulk and surface structural information are discussed in terms of the dissolution and corrosion mechanisms of ISG and related nuclear waste glasses.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Temperature and density dependence of Kr L-shell spectrum in hot dense plasmas

Kr L-shell spectroscopy modeling results are discussed in this paper, focusing on the n = 4 to n = 2 line transitions of Be- and Li-like Kr ions. Collisional radiative atomic kinetic and Stark-broadened spectral line shape calculations show electron temperature Te and density ne sensitivity in the spectrum. The combination of the Te dependence due to the relative intensity of Be-like to Li-like line emissions in the range from 1.5 to 3 keV and the ne sensitivity from the Stark broadening effect on the line shapes in the range from 5×1023 to 2×1024/ cc results in a spectrum that can be employed to diagnose Te and ne. Two different collisional radiative atomic kinetic models i.e., Prismspect [J. J. MacFarlane, et al., Int. Fusion Sci. Appl. Conf. Proc. 457 (2003)] and ABAKO [Florido, et al., PRE, 80, 056402 (2009)] produce similar results in level populations and spectra. In x-ray spectroscopy of implosion cores, this Kr L-shell spectrum may prove useful in an intermediate Te range in which Ar is too ionized for its K-shell to be of diagnostic value and Kr is not ionized enough for its K-shell emission to be useful.

Physics↗

DOE OSTI2023

Structure of the Bastnäsite (001) Surface by Crystal Truncation Rod X-ray Diffraction and Ab Initio Molecular Dynamics: Implications for Separations of a Rare Earth Ore Mineral

Bastnäsite ((Ce,La)FCO 3 ) is the primary mineral source of light rare earth elements, but its surface structure is not well understood. This presents a major challenge in improving beneficiation strategies. In this work, a synergistic combination of X-ray scattering and ab initio molecular dynamics (AIMD) was used to gain atomistic insight into the interfacial structure of bastnäsite. Surface X-ray scattering was used to measure crystal truncation rods (CTRs) of the bastnäsite (001) surface, a significant crystal face with a previously unknown termination. The best-fit atomic-scale model of the CTR data features a carbonate layer at the surface, which is stabilized by the relaxation of carbonate groups from their bulk structural positions. AIMD simulations predict similar surface relaxations, which are shown to be influenced by the protonation of oxygen atoms at the surface. Evidence of ordered water at the interface is also observed in the best-fit model and AIMD simulations. Further, the presence of a carbonate layer at this dominant crystal surface is significant for improving separation technologies because most commonly used ligands utilize anionic functional groups to chelate metal cations at particle surfaces. Without modification, anionic ligands are expected to have poor affinity for the carbonate-terminated (001) surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Simple theory of viscosity in liquids

A simplistic model of atomic dynamics in liquid, based on the concept of dynamic fluctuations in local atomic connectivity, is proposed to elucidate the temperature variation of the viscosity of a material in its liquid phase. Within the simplifications and hypothesis made to define the model, it is possible to explain the crossover from a simple liquid to a cooperative liquid and covariance of viscosity near the glass transition temperature. The model is a mechanical analog of the Drude model for the electric transport in metal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1980

An extended soft-cube model for the thermal accommodation of gas atoms on solid surfaces

A numerical soft cube model was developed for calculating thermal accommodation coefficients alpha and trapping fractions f sub t for the interaction of gases incident upon solid surfaces. A semiempirical correction factor c which allows the calculation of alpha and f sub t when the collision times are long compared to the surface oscillator period were introduced. The processes of trapping, evaporation, and detailed balancing were discussed. The numerical method was designed to treat economically and with moderate (+ or - 20 percent) accuracy the dependence of alpha and f sub t on finite and different surface and gas temperatures for a large number of gas/surface combinations. Comparison was made with experiments of rare gases on tungsten and on alkalis, as well as one astrophysical case of H2 on graphite. The dependence of alpha on the soft cube dimensionless parameters is presented graphically.

Burke, J. R.↗

NASA NTRS2017

Influence of Plasma Environment on K-Line Emission in Highly Ionized Iron Atoms Evaluated Using a Debye-Huckel Model

The influence of plasma environment on the atomic parameters associated with the K-vacancy states has been investigated theoretically for several iron ions. To do this, a time-averaged Debye-Huckel potential for both the electron-nucleus and electron-electron interactions has been considered in the framework of relativistic multiconfiguration Dirac-Fock computations. More particularly, the plasma screening effects on ionization potentials, K-thresholds, transition energies, and radiative rates have been estimated in the astrophysical context of accretion disks around black holes. In the present paper, we describe the behavior of those atomic parameters for Ne-, Na-, Ar-, and K-like iron ions.

Deprince, J.↗

DOE OSTI2014

Analytical models of calcium binding in a calcium channel

The anomalous mole fraction effect of L-type calcium channels is analyzed using a Fermi like distribution with the experimental data of Almers and McCleskey [J. Physiol. 353, 585 (1984)] and the atomic resolution model of Lipkind and Fozzard [Biochemistry 40, 6786 (2001)] of the selectivity filter of the channel. Much of the analysis is algebraic, independent of differential equations. The Fermi distribution is derived from the configuration entropy of ions and water molecules with different sizes, different valences, and interstitial voids between particles. It allows us to calculate potentials and distances (between the binding ion and the oxygen ions of the glutamate side chains) directly from the experimental data using algebraic formulas. The spatial resolution of these results is comparable with those of molecular models, but of course the accuracy is no better than that implied by the experimental data. The glutamate side chains in our model are flexible enough to accommodate different types of binding ions in different bath conditions. The binding curves of Na{sup +} and Ca{sup 2+} for [CaCl{sub 2}] ranging from 10{sup −8} to 10{sup −2} M with a fixed 32 mM background [NaCl] are shown to agree with published Monte Carlo simulations. The Poisson-Fermi differential equation—that includes both steric and correlation effects—is then used to obtain the spatial profiles of energy, concentration, and dielectric coefficient from the solvent region to the filter. The energy profiles of ions are shown to depend sensitively on the steric energy that is not taken into account in the classical rate theory. We improve the rate theory by introducing a steric energy that lumps the effects of excluded volumes of all ions and water molecules and empty spaces between particles created by Lennard-Jones type and electrostatic forces. We show that the energy landscape varies significantly with bath concentrations. The energy landscape is not constant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Efficient and selective carbon–carbon coupling on coke-resistant PdAu single-atom alloys

Here, we demonstrate that PdAu single-atom alloy model catalysts offer a heterogeneous route to selective Würtz-type C–C coupling. Specifically, when methyl iodide is exposed to an otherwise unreactive Au(111) surface, single Pd atoms in the surface layer promote C–I dissociation and C–C coupling, leading to the selective formation of ethane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Study on hydrogen isotopes permeation in fluidized state of liquid lithium-lead

Lithium-lead (Li-Pb) is one of the most promising candidate materials for the liquid blanket of fusion reactors. Hydrogen transfer under a fluidized condition of Li-Pb is investigated experimentally to design a Li-Pb blanket system. Li-Pb eutectic alloy flows inside a Ni tube in the experimental system, where H{sub 2} permeates into and out of the forced Li-Pb flow. The overall H{sub 2} permeation rate is analyzed using a mass balance model. Hydrogen atoms diffuse in Ni and Li-Pb. The steady-state H{sub 2} permeation rate obtained by this experiment is smaller than the result of the calculation model. A resistance factor is introduced to the present analysis in order to evaluate the influence of other H{sub 2} transfer mechanisms, such as diffusion in Li-Pb and dissolution reaction between Ni and Li-Pb. The contribution of the resistance to the overall H{sub 2} permeation rate becomes large when the flow rate of Li-Pb is low. This is because the boundary layer thickness between Ni and Li-Pb affects the overall H{sub 2} permeation rate. When the flow velocity of Li-Pb increases, the thickness of the boundary layer becomes thin, and the driving force of H{sub 2} permeation through the Ni wall becomes bigger. (authors)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2015

Hybrid-PIC Modeling of the Transport of Atomic Boron in a Hall Thruster

Computational analysis of the transport of boron eroded from the walls of a Hall thruster is performed by implementing sputter yields of hexagonal boron nitride and velocity distribution functions of boron within the hybrid-PIC model HPHall. The model is applied to simulate NASA's HiVHAc Hall thruster at a discharge voltage of 500V and discharge powers of 1-3 kW. The number densities of ground- and 4P-state boron are computed. The density of ground-state boron is shown to be a factor of about 30 less than the plasma density. The density of the excited state is shown to be about three orders of magnitude less than that of the ground state, indicating that electron impact excitation does not significantly affect the density of ground-state boron in the discharge channel or near-field plume of a Hall thruster. Comparing the rates of excitation and ionization suggests that ionization has a greater influence on the density of ground-state boron, but is still negligible. The ground-state boron density is then integrated and compared to cavity ring-down spectroscopy (CRDS) measurements for each operating point. The simulation results show good agreement with the measurements for all operating points and provide evidence in support of CRDS as a tool for measuring Hall thruster erosion in situ.

Materials Engineering↗