Search NASA⌕ Search

SEARCH · Search NASA

Results for “First principles calculations”

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

NiGa 2 O 4 interfacial layers in NiO/Ga 2 O 3 heterojunction diodes at high temperature

NiO/Ga 2 O 3 heterojunction diodes have attracted attention for high-power applications, but their high temperature performance and reliability remain underexplored. Here, we report the time evolution of the electrical properties in the widely studied p-NiO/n-Ga 2 O 3 heterojunction diodes and formation of NiGa 2 O 4 interfacial layers at high temperatures. Results of our thermal cycling experiment show an initial leakage current increase which stabilizes after sustained thermal load, due to reactions at the NiO–Ga 2 O 3 interface. High-resolution TEM microstructure analysis of the devices after thermal cycling indicates that the NiO–Ga 2 O 3 interface forms a ternary compound at high temperatures, and thermodynamic calculations suggest the formation of the spinel NiGa 2 O 4 layer between NiO and Ga 2 O 3 . First-principles defect calculations find that NiGa 2 O 4 shows low p-type intrinsic doping and hence can serve to limit electric field crowding at the interface. Vertical NiO/ Ga 2 O 3 diodes with intentionally grown 5 nm thin spinel-type NiGa 2 O 4 interfacial layers show an excellent device ON/OFF ratio of >10 10 (± 3 V), V ON of ~1.9 V, and increased breakdown voltage of ~1.2 kV for an initial unoptimized 300 lm diameter device. These p–n heterojunction diodes are promising for high-voltage, high temperature applications.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1975

Elastic deformation of the gyro rotor in the relativity gyro experiment

The elastic deformation of a spinning sphere is calculated from first principles. Numerical results are quoted for a quartz sphere spinning at 200 Hz, the most likely configuration for the relativity gyro experiment. The difference in principal moments of inertia is also calculated for the spinning gyro.

Eby, P. B.↗
DOE OSTI2025

Derivation of at-birth neutron energy distribution of americium-beryllium source based on ISO8529 spectral data

In many practical applications, the accuracy of neutron transport calculations depends on an estimate of the energy distribution provided in the source definition. While the energy distribution of neutrons emitted by spontaneous fission sources is well-defined, the spectra of neutrons produced in (α,n) reactions in compound mixtures, such as americium-beryllium (AmBe), depend on the microstructural properties of the source material, which may vary. Although the energy of neutrons produced in nuclear reactions can be calculated from first principles, such computations require detailed knowledge of the key characteristics of the source material, which are often unknown. Here, this study focuses on deriving the at-birth neutron energy distribution for an AmBe source by applying reverse transport calculations to the externally measured high-precision spectrum recommended in the latest revision of the international standard. The derived at-birth spectrum was validated through indirect energy-sensitive neutron emission rate measurements of AmBe sources calibrated at national primary metrology institutions. The resulting at-birth neutron energy distribution serves as essential input data for radiation transport modeling, providing a validated spectrum for broader experimental and computational applications.

AmBe↗
DOE OSTI2025

Fermi surface reconstruction due to orthorhombic distortion in the Dirac semimetal YbMnSb 2

A Dirac semimetal with magnetic atoms as constituents delivers an interesting platform to investigate the interplay of Fermi surface (FS) topology, electron correlation, and magnetism. One such family of semimetals is YbMn⁢ 𝑃⁢𝑛 2 (𝑃⁢𝑛 = Sb,Bi), which is being actively studied due to the intertwined spin and charge degrees of freedom. Here, in this study, we investigate the relationship between the magnetic/crystal structures and FS topology of YbMnSb 2 using single crystal x-ray diffraction, neutron scattering, magnetic susceptibility, magnetotransport measurements, and complimentary density functional theory (DFT) calculations. Contrary to previous reports, the x-ray and neutron diffraction reveal that YbMnSb 2 crystallizes in an orthorhombic 𝑃⁢𝑛⁢𝑚⁢𝑎 structure with a notable antiphase displacement of the magnetic Mn ions that increases in magnitude upon cooling. As a result of this orthorhombicity, first-principles DFT calculation reveals a reduced Brillouin zone and more anisotropic FS of YbMnSb 2 compared to YbMnBi 2 . Moreover, the hole-type carrier density drops by two orders of magnitude as YbMnSb 2 orders antiferromagnetically, indicating band folding in a magnetic ordered state. In addition, the Landau level fan diagram yields a nontrivial nature of the Shubnikov–de Haas quantum oscillation frequency arising from the Dirac-like Fermi pocket. These results imply that YbMnSb 2 is an ideal platform to explore the interplay of subtle lattice distortion, magnetic order, and topological transport arising from relativistic quasiparticles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1995

Improved Capacitive Liquid-Level Probe

Improved capacitive probe designed to measure level of dielectric liquid in tank. In original application, liquid is hydrogen at low temperature, but design principles applicable to measurement of levels of other dielectric liquids in various temperature ranges. Outstanding features of design those making probe insensitive to temperature-induced and other spurious capacitive effects that otherwise distort capacitive liquid-level readings, and those making it possible to calibrate probe mostly by calculation from first principles.

Ziemke, Robert A.↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

An energy-modified quantum defect method for the analysis of Rydberg spectra: Application to 2-butyne

The high resolution Rydberg absorption spectrum of 2-butyne C 4 H 6 recorded previously at the SOLEIL synchrotron facility has been interpreted using multichannel quantum defect theory (MQDT). The calculations are based on the continuum scattering calculations of Xu et al., J. Chem. Phys. 136, 154303 (2012) and of Jacovella et al., J. Phys. Chem. A 119, 12339 (2015) pertaining to the dipole-allowed excited state symmetries in absorption from the ground state. In contrast to the traditional approach of calculating low-lying electronic states first and then attempting to extend the calculations to ever higher energy, here the analysis proceeds through the extension of these detailed calculations of the electronic continuum scattering down into the discrete region of the spectrum. The continuum reaction matrices and dipole transition moments are adapted to the discrete Rydberg region via the use of an energy-modified formulation of MQDT theory and associated energy dependences of the quantum defects. The analysis reproduces more than 40 Rydberg states from n ≈ 10 down to the 3 d and 4 s levels with an rms error of better than 20 cm −1 . Further, these belong to five Rydberg series with three different molecular symmetries. While the approach predicts many additional series, most of these are calculated and observed to carry only little oscillator strength. The analysis shows that the Rydberg spectrum is dominated by the excitation of an e″ symmetry electron of fδ and gπ type, in line with what previous studies of the above-threshold shape resonance of 2-butyne have shown. The present study is intended to serve as an example showing how first principles continuum calculations may be useful for the interpretation of highly bound discrete states in a range that poses problems for the standard ab initio techniques. The quantitative treatment of the dipole absorption cross sections is deferred to a future paper.

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

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Rearrangement collision theory of phonon-driven exciton dissociation

Understanding the processes governing the dissociation of excitons to free charge carriers in semiconductors and insulators is of central importance for photovoltaic applications. Dyson's S-matrix formalism provides a framework for computing scattering rates between quasiparticle states derived from the same underlying Hamiltonian, often reducing to familiar Fermi's “golden rule” like expressions at first order. By presenting a rigorous formalism for multichannel scattering, we extend this approach to describe scattering between composite quasiparticles and, in particular, the process of exciton dissociation mediated by the electron-phonon interaction. Subsequently, we derive rigorous expressions for the exciton dissociation rate, a key quantity of interest in optoelectronic materials, which enforce correct energy conservation and may be readily used in ab initio calculations. We apply our formalism to a three-dimensional model system to compare temperature-dependent exciton rates obtained for different scattering channels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dynamical correlations leading to site and orbital selective Mott insulator transition in hydrogen doped SmNiO 3

Electron doping induces metal -to -insulator transition (MIT) in SmNiO 3 as realized by experiments. While earlier density functional theory (DFT) studies with static correlations fell short of explaining the recent MIT observations at lower hydrogen concentrations, we present a comprehensive computational investigation employing an advanced approach. We combine DFT with dynamical mean field theory (DFT + DMFT) to efficiently analyze the insulating behavior of hydrogen -doped SmNiO 3 . In contrast to previous theoretical works, our calculations predict an insulator transition occurring at a reduced doping level of H:Ni = 0.5:1. Specifically, while the DFT + U method reveals a gap opening between p-to-d orbitals, the DMFT approach highlights a gap opening between d-to-d orbitals. Our findings uncover a selective Mott transition in site and orbital characteristics, with the Ni ions proximate to the doped hydrogen exhibiting Mott -like traits. Notably, DMFT calculations highlight a pronounced dependence on Hund's parameter J, implying the presence of Hundness in the Mott insulator. Finally, this study underscores the necessity of accounting for dynamical correlations to accurately describe the electronic structure of strongly correlated electron-doped rare-earth nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Theoretical Calculation Study of Nanomaterials

In recent years, theoretical simulations using computational physics and chemistry models have been used widely to predict the physical, chemical, and structural properties of nanomaterials. To reflect such rapid progress, two Special Issues on the theoretical study of nanomaterials were published. The 36 articles collected cover broad re-search fields and applications. In this Editorial, we provide an overview of these pub-lications and summarize the main findings. Based on these research outcomes, we lay out some future research directions.

First-principles calculations↗
DOE OSTI2026

Impact of absorption due to zero-field splitting on loss in dielectrics: A case study in sapphire

The coherence times of superconducting qubits are limited by loss mechanisms, whose microscopic origins have remained elusive. We propose a mechanism caused by transitions between zero-field-split states of paramagnetic impurities or defects. We derive the absorption cross section for a magnetic dipole transition and apply it to calculate the bulk loss tangent. For Cr, Fe, and V impurities in sapphire, we find loss tangents at 4.5 GHz in the range of 10 −9 –10 −8 , comparable to the loss measured in experiments. This value suggests that magnetic loss may be a limiting factor in the coherence times of superconducting qubits.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Size-dependent attraction of Cu solutes to clusters formed at Ag grain boundaries

We report a size-dependent solute clustering mechanism at grain boundaries in a sputtered ultrafine-grained Ag-Cu alloy, where large Cu clusters form despite weak individual solute-solute interactions. X-ray diffraction confirms limited Cu solubility in the Ag matrix, while scanning transmission electron microscopy reveals Cu clustering at both ordinary GBs and GB junctions. Density functional theory calculations show that 12-atom Cu clusters are energetically preferred, while smaller three-atom clusters are significantly less stable. Additional calculations demonstrate a marked increase in solute-cluster attraction energy with cluster size. As a result, these findings point to a previously unrecognized pathway for grain-boundary solute clustering in immiscible systems, driven by collective solute-cluster interactions, with implications for segregation behavior and stability in nanocrystalline and ultrafine-grained alloys.

Alloys↗
DOE OSTI2024

Two-Dimensional Moiré Polaronic Electron Crystals

Two-dimensional moiré materials have emerged as the most versatile platform for realizing quantum phases of electrons. Here, we explore the stability origins of correlated states in WSe 2 /WS 2 moiré superlattices. We find that ultrafast electronic excitation leads to partial melting of the Mott states on timescales 5 times longer than predictions from the charge hopping integrals and that the melting rates are thermally activated, with activation energies of 1⁢8±3 and 1⁢3±2 meV for the one- and two-hole Mott states, respectively, suggesting significant electron-phonon coupling. A density functional theory calculation of the one-hole Mott state confirms polaron formation and yields a hole-polaron binding energy of 16 meV. These findings reveal a close interplay of electron-electron and electron-phonon interactions in stabilizing the polaronic Mott insulators at transition metal dichalcogenide moiré interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

Fitting Vibrational Energy Parameters to Falloff Experiments: CH3 + CH3 + He (0.6-2 Torr 200-298 K)

There are many reported experiments and theoretical calculations for the combination of methyl free radicals diluted in a non-reactive 'deactivator. At high pressure the ethane (C2H6) combination product is the sole product; with decreasing pressure the chemically activated C2H6 will decompose. The falloff in the observed rate coefficient is the result of the competition between collisional stabilization of the chemically activated CA and unimolecular decomposition. The dependence of the rate coefficient on pressure, temperature and collision properties is complex and can not be calculated from first principles. The understanding of this system is not only of fundamental importance but is relevant to the recent detection of methyl free radicals in the atmospheres of Saturn - and Neptune. The temperatures of these outer planet atmospheres are in the 140-200 K region with total pressures (predominately H2 and He) less than 0.2 Torr. Experimentally determined rate coefficients have been reported for this reaction at T = 296-906 K and T = 200-408 K mostly with argon as the deactivator. At T = 200 K only the high pressure rate coefficient has been determined. Complete falloff curves over a wide temperature range (200-1600 K) with a variety of weak collider models used to simulate argon as the deactivator have also been reported by Klippenstein and Harding (KH). More recently we have reported the experimental rate coefficients in the falloff region with helium as the deactivator at 200 and 298 K. In this paper we have used the calculated falloff curves reported by KH for argon to determine the average energy transferred per collision for helium in our recently reported experiments. Collision rates were converted using Lennard Jones parameters; the temperature dependence of this conversion factor is noted. The helium experiments were consistent with a down of approximately 100 cm (exp-1); the temperature dependence was slight. The magnitude of down and its temperature dependence will be compared to other systems with similar substrate complexity and reaction energetics.

Tardy, Dwight C.↗