Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Minimizing the diffusivity difference between vacancies and interstitials in multi-principal element alloys

Interstitial atoms usually diffuse much faster than vacancies, which is often the root cause for the ineffective recombination of point defects in metals under irradiation. Here, via ab initio modeling of single-defect diffusion behavior in the equiatomic NiCoCrFe(Pd) alloy, we demonstrate an alloy design strategy that can reduce the diffusivity difference between the two types of point defects. The two diffusivities become almost equal after substituting the NiCoCrFe base alloy with Pd. The underlying mechanism is that Pd, with a much larger atomic size (hence larger compressibility) than the rest of the constituents, not only heightens the activation energy barrier ( E a ) for interstitial motion by narrowing the diffusion channels but simultaneously also reduces E a for vacancies due to less energy penalty required for bond length change between the initial and the saddle states. Our findings have a broad implication that the dynamics of point defects can be manipulated by taking advantage of the atomic size disparity, to facilitate point-defect annihilation that suppresses void formation and swelling, thereby improving radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2022

Deformed spin- 1 2 square lattice in antiferromagnetic NaZnVOPO 4 ( HPO 4 )

In this work, we report the structural and magnetic properties of a new spin-$\frac{1}{2}$ antiferromagnet NaZnVOPO 4 (HPO 4 ) studied via x-ray diffraction, magnetic susceptibility, high-field magnetization, specific heat, and 31 P nuclear magnetic resonance (NMR) measurements, as well as density-functional band-structure calculations. While thermodynamic properties of this compound are well described by the J 1 – J 2 square-lattice model, ab initio calculations suggest a significant deformation of the spin lattice. From fits to the magnetic susceptibility we determine the averaged nearest-neighbor and second-neighbor exchange couplings of $\bar{J}$ 1 ≃ – 1.3 K and $\bar{J}$ 2 ≃ 5.6 K , respectively, resulting in the effective frustration ratio α = $\bar{J}$ 2 – $\bar{J}$ 1 ≃ – 4.3 that implies columnar antiferromagnetic order as the ground state. Experimental saturation field of 15.3 T is consistent with these estimates if 20% spatial anisotropy in J 1 is taken into account. Specific heat data signal the onset of a magnetic long-range order at T N ≃ 2.1 K , which is further supported by a sharp peak in the NMR spin-lattice relaxation rate. The NMR spectra mark the superposition of two P lines due to two nonequivalent P sites where the broad line with the strong hyperfine coupling and short T 1 is identified as the P(1) site located within the magnetic planes, while the narrow line with the weak hyperfine coupling and long T 1 is designated as the P(2) site located between the planes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Particle-hole asymmetric phases in doped twisted bilayer graphene

Twisted bilayer graphene (TBG) has emerged as a paradigmatic platform for exploring the interplay between strong interactions in a multi-band system with nearly flat bands, while offering unprecedented control over the filling fraction of electron/hole carriers. Despite much theoretical work, developing a comprehensive ab initio model for this system has proven challenging due to the inherent trade-off between accurately describing the band structure and incorporating the interactions within the Hamiltonian. Here, we constructed an extended multi-orbital Hubbard model, and performed Hartree-Fock (HF) calculations to explore its phase diagram across commensurate fillings. We found several nearly-degenerate insulating states at charge neutrality, all of which exhibit orbital orders. We also observed that most of the Hartree–Fock ground states exhibit a generalized valley Hund’s-like rule, resulting in valley polarization. Importantly, we show that the incorporation of the intra-valley and inter-valley exchange interactions is crucial to properly stabilize the ordered symmetry-broken states. In agreement with experiments, we find significant particle-hole asymmetry, which underscores the importance of using particle-hole asymmetric models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

High-Precision Spectroscopy of 20 O Benchmarking Ab Initio Calculations in Light Nuclei

The excited states of unstable 20 O were investigated via γ-ray spectroscopy following the 19 O(d, p) ⁢20 O reaction at 8 A⁢ MeV. By exploiting the Doppler shift attenuation method, the lifetimes of the $2^+_2$ and $3^+_1$ states were firmly established. From the γ-ray branching and E⁢2/M⁢1 mixing ratios for transitions deexciting the $2^+_2$ and $3^+_1$ states, the B⁡(E⁢2) and B⁡(M⁢1) were determined. Various chiral effective field theory Hamiltonians, describing the nuclear properties beyond ground states, along with a standard USDB interaction, were compared with the experimentally obtained data. Such a comparison for a large set of γ-ray transition probabilities with the valence space in medium similarity renormalization group ab initio calculations was performed for the first time in a nucleus far from stability. It was shown that the ab initio approaches using chiral effective field theory forces are challenged by detailed high-precision spectroscopic properties of nuclei. Finally, the reduced transition probabilities were found to be a very constraining test of the performance of the ab initio models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Molecular beam epitaxy of PdO on MgO (001)

PdO, widely used in catalysis in powder form for decades, has been predicted recently to be a Dirac semimetal. Synthesis of high-quality single crystals of this material is thus of great interest. Here, by using ozone-assisted molecular beam epitaxy, PdO thin films were grown on a MgO (001) substrate. X-ray diffraction and transmission electron microscopy indicate the film is a / b axis oriented, with the c axis lying in plane. Fully oxygenated PdO films have a low density of holelike carriers, and are insulating at the lowest temperatures. Our density functional theory calculations using the Heyd-Scuseria-Ernzerhof exchange-correlation functional suggest a ~1.0 eV band gap at the M point, where the gap can be reduced by tensile strain along the c axis and the bands begin to invert when the tensile strain is as high as 18%. Although tensile strain of this magnitude is experimentally not viable using epitaxy, electrons can be doped by oxygen reduction. Finally, our results emphasize the need for careful consideration of electron correlation effects and stoichiometry in ab initio modeling of topological semimetals involving transition metal oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tuning band alignment at a semiconductor-crystalline oxide heterojunction via electrostatic modulation of the interfacial dipole

In this paper, we demonstrate that the interfacial dipole associated with bonding across the SrTiO 3 /Si heterojunction can be tuned through space charge, thereby enabling the band alignment to be altered via doping. Oxygen impurities in Si act as donors that create space charge by transferring electrons across the interface into SrTiO 3 . The space charge induces an electric field that modifies the interfacial dipole, thereby tuning the band alignment from type II to III. The transferred charge, accompanying built-in electric fields, and change in band alignment are manifested in electrical transport and hard x-ray photoelectron spectroscopy measurements. Ab initio models reveal the interplay between polarization and band offsets. We find that band offsets can be tuned by modulating the density of space charge across the interface. Modulating the interface dipole to enable electrostatic altering of band alignment opens additional pathways to realize functional behavior in semiconducting hybrid heterojunctions.

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

Transient electron scavengers modulate carrier density at a polar/nonpolar perovskite oxide heterojunction

Here, we show how transient electron scavengers can be utilized to control the carrier concentration at polar/non-polar perovskite interfaces. By combining quantitative synchrotron x-ray-based interface structure determination with ab initio modeling, we demonstrate that Nd vacancy formation and the resulting formation of Nd adatoms, stabilized by oxygen scavengers at the growth front, can quantitatively account for the decreased carrier concentration at the SrTiO 3 /n NdTiO 3 /SrTiO 3 (001) heterojunction for n = 1 unit cell. This study yields new insights into growth mechanisms and the effect of transient species and defects on the electronic properties of oxide heterojunctions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface termination control of charge transfer and band alignment across a semiconductor–crystalline-oxide heterojunction

Charge redistribution across heterojunctions has long been utilized to induce functional response in materials systems. Here we examine how the composition of the terminating surface affects charge transfer across a heterojunction consisting of Si and the crystalline complex oxide SrTiO 3 . Itinerant electrons in Si migrate across the interface toward the surface of SrTiO 3 due to surface depletion. The electron transfer in turn creates an electric field across the interface that modifies the interfacial dipole associated with bonding between SrTiO 3 and Si. The modification in the dipole leads to a change in band alignment, in which the conduction band of SrTiO 3 moves from being above the valence band of Si in energy, to below it. By capping the SrTiO 3 surface with ultrathin (≤ 1 nm) layers of BaO, SrO or TiO 2 , charge transfer across the interface can be weakened or inhibited. Ab initio modeling implicates the adsorption of oxygen associated with exposure to ambient conditions as driving the surface depletion in SrTiO 3 . In conclusion, the electronic coupling between the surface and buried interface expands the functionality of semiconductor-crystalline oxide heterojunctions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bulk properties of honeycomb lattices of superconducting microwave resonators

We have realized different honeycomb lattices for microwave photons in the 4 to 8 GHz band using superconducting spiral resonators. Each lattice comprises a few hundred sites. Two designs have been studied, one leading to two bands touching at the Dirac points and one where a gap opens at the Dirac points. Using a scanning laser technique to image the eigenmodes of this new type of photonic lattices, we are able to reconstruct their band structure. The measured bands are in excellent agreement with ab initio models that combine numerical simulations of the electromagnetic properties of the spiral resonator and analytical calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Chemically tuned intermediate band states in atomically thin Cu x GeSe/SnS quantum material for photovoltaic applications

A new generation of quantum material derived from intercalating zerovalent atoms such as Cu into the intrinsic van der Waals gap at the interface of atomically thin two-dimensional GeSe/SnS heterostructure is designed, and their optoelectronic features are explored for next-generation photovoltaic applications. Advanced ab initio modeling reveals that many-body effects induce intermediate band (IB) states, with subband gaps (~0.78 and 1.26 electron volts) ideal for next-generation solar devices, which promise efficiency greater than the Shockley-Queisser limit of ~32%. The charge carriers across the heterojunction are both energetically and spontaneously spatially confined, reducing nonradiative recombination and boosting quantum efficiency. Using this IB material in a solar cell prototype enhances absorption and carrier generation in the near-infrared to visible light range. Tuning the active layer’s thickness increases optical activity at wavelengths greater than 600 nm, achieving ~190% external quantum efficiency over a broad solar wavelength range, underscoring its potential in advanced photovoltaic technology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

FY23 Progress Report on Viscosity and Thermal Conductivity Measurements of Molten Salts

As presented in this report, thermal conductivity and viscosity measurements were performed on key chloride pseudo-binary molten salt systems of relevance to molten salt reactor developers. Thermal conductivity measurements were conducted with a variable gap technique, in which a known heat flux is driven across a molten salt specimen and the temperature difference is measured, allowing calculation of the thermal conductivity. This is achieved by establishing a small gap between the bottom of a cylindrical inner containment, which houses electrical heating elements, and an outer containment, which houses cooling channels; the gap size can be varied by compression of a formed bellows. A new calibration scheme was developed herein, involving a correction factor to the heat flux based on He measurements at various temperatures. Furthermore, the data processing methodology was improved to minimize the impact of radiative heat transfer in the calculation of salt specimen thermal conductivity from the temperature difference measurements. Viscosity measurements were conducted with a rolling ball viscometer, in which a ball rolls some known distance in an angled tubular crucible, and the terminal velocity can be used to calculate the viscosity of the salt. The measurement can be performed in a quartz crucible, with which a standard camera can be used to track the ball, or in a metal crucible, with which x-ray radiography is required to track the ball. A new custom x-ray system was made and dedicated to the rolling ball viscometer to enable high throughput automated measurements with salts which require containment with metal. Both systems have been integrated with a ventilation stack which allows for off-gassing of radioactive material, enabling future measurements with U-bearing salts. The thermal conductivity measurements performed herein were with NaCl-KCl (44 mol% NaCl). This salt system was measured in the previous fiscal year, however the thermal conductivity values obtained were comparatively low, and so the system was remeasured with the aforementioned calibration scheme and improved post-processing techniques. The newly obtained thermal conductivity values for NaCl-KCl (44 mol% NaCl) indicate good agreement with kinetic theory and ab-initio models (within 5–10 %). The viscosity measurements performed herein were with three different compositions of NaCl-KCl: 75, 50, and 25 mol% NaCl; a new calibration scheme was applied to account for variable flow effects in the laminar regime. The results show reasonable agreement with literature (5–20 %, depending on the temperature, composition, and study); however, literature values are likely higher than true values based on pure end-member measurements performed in the comparative studies. The measurements conducted herein do show a trend such that viscosity increases with increasing NaCl concentration, which agrees well with one of two comparative studies. In general, the measurements conducted herein gives confidence in the capability to use these systems to accurately measure thermal conductivity and viscosity of actinide-bearing salts within the next fiscal year.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Computational and Experimental Investigation of Thermal-Mechanical- Chemical Mechanisms of High-burnup Spent Nuclear Fuel (SNF) Processes at Elevated Temperatures and Degradation Behavior in Geologic Repositories

The overarching goal of the combined computational and experimental R&D activities proposed in this project is to enhance understanding of the mechanisms and thermal-mechanical-chemical (TMC) parameters controlling the instant release fraction (IRF) and matrix dissolution of high-burnup (HB; burnup) spent nuclear fuels (SNFs) and the subsequent formation, stability, and phase transformations of SNF alteration products under long-term storage and geological disposal conditions. Uranium dioxide may undergo oxidative corrosion/alteration, and the IRF may be increased for HB SNF, both of which may affect environmental systems associated with SNF long-term storage and disposal. The oxidative matrix dissolution may form various complex uranyl-based phases, including a rich variety of oxides, silicates, carbonates and other secondary minerals in varied geological environments (e.g., studtite, metastudtite, amorphous uranyl peroxide, uranium trioxide, triuranium octoxide, schoepite, dehydrated schoepite, metaschoepite, becquerelite, soddyite, rutherfordine,...). These uranyl phases generally have higher mobility UO 2 +2 species than less soluble U 4+ phases. However, limited information on the thermodynamic properties and formation kinetics of these uranyl-bearing phases is available to predict explicitly paragenesis under the conditions relevant to long-term storage or disposal. The proposed project draws on complementary expertise and research backgrounds from the team members: (i) to apply a combined ab initio modeling (UNLV/UTEP and SNL) and experimental (UNLV) strategy investigating the high-temperature TMC mechanisms of alteration of SNF under α-radiolysis conditions; (ii) to investigate the mechanistic of phase transformations in UNF degradation products under various conditions expected in long-term storage systems (e.g. (UO 2 )O 2 (H 2 O) 4 → (UO 2 )O 2 (H 2 O) 2 → U 2 O 7 → UO 3 → U 3 O 8 ); (iii) to determine high-accuracy TMC parameters for complex uranyl-based phases formed in storage or geological disposal environments (e.g. UO 3 (H 2 O) 2 , Ca[(UO 2 ) 6 O 4 (OH) 8 ] 8 H 2 O, (UO 2 ) 2 (SiO 4 ) 3 2H 2 O,…). The unforeseen COVID-19 pandemic led to the laboratory/campus closure since March 2020, that resulted in a significant delay in reaching milestones in a satisfactory manner, due to (i) the statewide recommendation from stop-working to later limited work in the lab and work-from-home (WFH), (ii) no in-person interactions, and (iii) a hiring freeze at UNLV. Therefore, a no cost extension (10/01/2021- 9/30/2022) was requested to help make up the time we lost during the global pandemic in 2020-2021, leading to paradigm shifts in the focus of the project in the following three main tasks: Task 1 (Computational), Task 2 (Experimental), and Task 3 (Final report, due on 12/29/2022).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Understanding and Tailoring Diffusion and Co-Adsorption Inside the Confined Pores of Metal-Organic Frameworks (Final Scientific/Technical Report for Award DE-SC0019902)

The aim of this program was to gain a fundamental understanding of the behavior of various guest molecules in nano-confined environments, such as metal organic frameworks (MOFs), using a combination of novel synthesis, ab initio modeling, and in situ characterization. Through this project, we developed a concise understanding of the mechanisms that control adsorption/desorption of gaseous molecules and their mixtures, leading to design/synthesis guidelines for MOFs with desired functionality. We further developed methods to disentangle kinetic from thermodynamic effects during adsorption, as well as to characterize the interactions at play. In the first funding cycle, the focus was on the unambiguously characterization of co-adsorption and diffusion of gasses/vapors and their mixtures. In the second funding cycle, the focus was on characterizing the effects of the nano-confinement on the kinetics and thermodynamics of adsorption processes inside MOFs, again with an emphasis on mixtures of gasses and vapors. The nano-confinement can tip the thermodynamic vs. kinetic balance, and current understanding and theory based on single-component analysis can lead to incorrect predictions for mixtures. This is of particular interest in real-world applications, where gasses/vapors are typically mixed, contain impurities, or are often exposed to humid conditions. Our main findings were: (i) within confined environments the adsorption behavior of mixed gasses/vapors can be drastically different from the “sum” of the corresponding single phases; (ii) co-adsorption is often competitive and detrimental to performance, but it can also be cooperative and beneficial; (iii) in some co-adsorbed gasses/vapors, molecules that are strongly bound in the single-component phase can be replaced by molecules that are nominally weaker bound (molecular exchange) due to guest-guest interactions that lower the kinetic barriers and favor the final adsorption state; (iv) kinetic and thermodynamic effects can be precisely controlled through pore-size engineering and synthesis; and, (v) kinetic effects can be identified and disentangled from thermodynamic effects during adsorption through a series of sequential and simultaneous gas loading measurements. The short-term goal of this program was the controlling and understanding of common MOF systems in real-world situations where gasses/vapors are mixed, which will have an important impact on industrial processes and applications from gas storage and sequestration to catalysis and sensors. The long-term goals include the development of theoretical and experimental methods for gaining a fundamental understanding of adsorption/reaction processes within MOFs, as well as new guidelines for synthesizing MOFs with tailored physical and chemical properties.

36 MATERIALS SCIENCE↗
NASA NTRS1993

A generalized version of the Rankine-Hugoniot relations including ionization, dissociation and related phenomena

For purposes of computing shocks in stellars atmospheres and winds we have developed a generalized version of the Rankine-Hugoniot relations including ionization, dissociation, radiation and related phenomena such as excitation, rotation and vibration of molecules. The new equations are given in analytical form. They are valid as long as the internal energy E, the total pressure P, and the first adiabatic coefficient gamma(sub 1) can be evaluated. However, we have not treated shock structures. In the case of non-LTE we have to employ an approximation for gamma(sub 1) because in that case no definition exists. Our new version of the Rankine-Hugoniot relations can easily be used for many purposes including ab-initio modeling. In our derivation we introduce a parameter gamma(sub H), which is definded as the ratio of the enthalpy H (sometimes called heat function w) to the internal energy E (sometimes called U). Using this parameter we solve the equations for changing mu and (d(natural log P)/d(natural log rho))(sub ad) identically equal to gamma(sub 1) on both sides of the shock. Both gamma(sub H) and gamma(sub 1), and also mu are functions of pressure P and temperature T. We present: (1) the derivation, (2) examples of gamma(sub 1) (P,T) and gamma(sub H) (P,T) which include/exclude ionization and radiation, and (3) as an example the differences in post-shock parameters as function of the pre-shock temperature for the case with ionization and without radiation.

Nieuwenhuijzen, H.↗