Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density of states”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sound Velocity and Grüneisen Parameter in Shock‐Melted Silica at Deep Earth Pressures

Silica is a primary component of rocky planet interiors and its melt properties are important for understanding planetary formation and differentiation, magma oceans, and the deep mantle. Although well understood in the solid state, the high-pressure behavior of liquid silica is poorly constrained at lower mantle pressures. Using laser interferometry to measure shock wave profiles, we report measured stress-density states and longitudinal sound speeds in shock-synthesized stishovite, from fused silica staring material, and across the solid-liquid phase boundary up to 154 GPa. Our results constrain completion of melt at 80 GPa and show that at pressures relevant to the deep mantles of Earth-sized, rocky planets, the Grüneisen parameter for liquid silica increases with compression. This finding is consistent with a continuous increase in Si-O coordination above six for liquid silica at core-mantle boundary relevant pressures and temperatures.

Ocampo, I. K. [Princeton University, NJ (United St↗
NASA NTRS1974

Equation of state at ultrahigh densities. I

Particular attention is given to the physics involved in the derivation of an equation of state for a density range from 10 g/ccm to 200,000 billion g/ccm and for the density region from 200,000 billion g/ccm to 2,000,000 billion g/ccm. Conditions in a high density regime with densities exceeding 500,000 billion g/ccm are examined. Equations of state for a gas of pure neutrons and for a gas of nucleons, electrons, and hyperons are considered.

Canuto, V.↗
DOE OSTI2025

Rb 4 CuSb 2 Cl 11 and Rb 2 In 0.91(0.2) Sb 0.09 Cl 5 ·H 2 O: Wide Band Gap 0D Metal Halide Semiconductors

Herein, we report the discovery, structural and photophysical characterization of a new zero-dimensional (0D) lead-free all-inorganic halide, Rb 4 CuSb 2 Cl 11 , which adopts a new structure type. Single-crystal X-ray diffraction (SCXRD) shows that the structure consists of isolated, distorted seesaw [SbCl 4 ] − and trigonal planar [CuCl 3 ] 2− units, separated by Rb + cations that provide charge balance. Optoelectronic measurements and density functional theory (DFT) calculations indicate an indirect band gap of 2.89 eV, making it a candidate for wide-bandgap optoelectronic applications. Electrical resistivity was measured at 1.29 × 10 10 Ω·cm, and the trap-state density (n trap ) was found to be 7.44 × 10 10 cm −3 . Attempts to synthesize substitution analogs of Rb 4 CuSb 2 Cl 11 led to the synthesis of Rb 2 In 0.91(0.2) Sb 0.09 Cl 5 ·H 2 O, which was erroneously reported as Rb 2 SbCl 5 O in a previous study. Rb 2 In 0.91(0.2) Sb 0.09 Cl 5 ·H 2 O adopts a vacancyordered perovskite structure and exhibits broad-band yellow emission under UV excitation. The measured photoluminescence quantum yield (PLQY) for Rb 2 In 0.91 (0.2)Sb 0.09 Cl 5 ·H 2 O is 18.2%. These findings add to the growing class of quaternary metal halides with multiple cation and anion compositions, expanding the chemical phase space for the discovery of new materials with functional properties.

Crystal structure↗
NASA NTRS1977

The effect of oxidation-expanded defects upon MOS parameters

Implantation of Ne ions was chosen as a means of introducing a uniformly distributed quantitatively controllable and reproducible amount of microdefects at the silicon surface. Each wafer was ion-implanted on one half of its surface, the other half remaining as a control. Seven implant fluences from 1 trillion to 5 by 10 to the 14th power Ne ions/sq cm were used. The fluences, as well as the concentrations of defects expanded by the subsequent gate oxidation, were correlated with seven different MOS parameters. For implant fluences of 10 to the 14th power per sq cm and above, increasing stacking-fault densities were found. The oxide defect density associated with time-dependent breakdown increased significantly with Ne dose. The presence of oxidation-expanded defects drastically increases the surface-generation velocity and decreases the bulk lifetime. Slight increases in flatband voltage and surface-state density occurred only for the wafers having the greatest defect concentrations. Oxidation of a P-containing damaged silicon surface results in a greater pileup of P at the surface than occurs for the oxidation of a similar but undamaged surface. It appears that the expanding defects restrain the diffusion of P inward from the surface.

Prussin, S.↗
DOE OSTI2026

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Identification of In Situ Generated High-Spin Fe(II)-N 4 Active Sites for Acidic Oxygen Reduction Reaction via Operando 57 Fe Mössbauer Spectroscopy

Iron-nitrogen-carbon (FeNC) catalysts are considered among the most promising alternative to Pt catalysts in acidic oxygen reduction reaction (ORR), yet the geometric and electronic structures of the true active site under reaction conditions have not been clearly elucidated. Herein, we synthesized a representative FeNC catalyst by pyrolyzing Fe 3+ absorbed on ZIF-8-derived N-doped carbon at a mild temperature under the H 2 /Ar atmosphere, where a formation mechanism of FeN 4 sites through a Zn-mediated Fe nanoparticle atomization process was proposed. The resulting FeNC-750 catalyst shows high acidic ORR activity with a half-wave potential of 0.838 V and a peak power density of 1 W cm -2 in proton-exchange membrane fuel cell (PEMFC). By using operando 57 Fe Mössbauer spectroscopy on FeNC-750, it was revealed that the pyrrolic N-coordinated high-spin Fe 2+ N 4 sites, which are in situ generated from high-spin Fe 3+ N 4 during ORR, are identified as catalytic active states. Density functional theory calculations further verified that, compared to the pyridinic N-coordinated low-spin Fe 2+ N 4 , the pyrrolic N-coordinated high-spin Fe 2+ exhibits optimized adsorption energy for reaction intermediates, thereby lowering the energy barrier of the rate-determining step (RDS) and facilitating OH* desorption. In conclusion, this work provides both experimental and theoretical evidence of the true active site in the FeNC during acidic ORR, offering significant insights for the rational design of high-performance materials for acidic fuel cells.

FeNC catalyst↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2025

Investigation of the Effect of Gate Oxide Screening with Adjustment Pulse on Commercial SiC Power MOSFETs

This paper presents a method to recover the negative threshold voltage shift during high field gate oxide screening of 1.2 kV 4H-SiC MOSFETs with an additional adjustment gate voltage pulse. To reduce field failure rates of the MOSFETs in operation, manufacturers perform a screening treatment to remove devices with extrinsic defects in the oxide. Current gate oxide screening procedures are limited to oxide fields at or below ~9 MV/cm for short durations (<1 s), which is not enough to remove all the devices with extrinsic defects. The results show that by implementing a lower field gate pulse, the threshold voltage shift can be partially recovered, and therefore the maximum screening field and time can be increased. However, both the initial screening pulse and the adjustment pulse require careful calibration to prevent significant degradation of the device threshold voltage, on-resistance, interface state density, or intrinsic lifetime. With a well calibrated set of pulses, higher screening fields can be utilized without significantly damaging the devices. This leads to an improvement in the overall screening efficiency of the process, reducing the number of devices with extrinsic oxide defects entering the field, and improving the reliability of the SiC MOSFETs in operation.

42 ENGINEERING↗
DOE OSTI2026

Buried Interfaces in Organic Photocathodes for H 2 Evolution: Fermi-Level Pinning and Recombination

Herein, we demonstrate how Fermi-level pinning at buried contacts impacts solar fuel generation in all-polymer photocathodes by systematically comparing the effects of work function, hydroxyl coverage, and hydrogen evolution using chemically modified indium tin oxide (ITO) supports. Photovoltages and net photocathode performance are improved when the ITO is passivated using phosphonic acids, independent of work function, suggesting that the passivation reduces Fermi-level pinning at the buried interface arising from blended heterojunction interactions with surface metal hydroxyls. Transient photovoltage decay reveals differences in recombination mechanisms, supported by light intensity-dependent measurements. Briefly, nonpassivated, hydrophilic contacts exhibit trap-assisted recombination, while passivated, hydrophobic contacts follow bimolecular recombination. We then investigate changes in electroactivity of hole-transfer processes as a function of scan rate and repetitive cycling using a diffusion-controlled molecular redox probe, analogous to a hole-only device achieved via the electrolyte. The nonpassivated buried contacts exhibit higher overpotentials for oxidation, indicative of hole injection/extraction barriers. We observe irreversible electron transfer via the hole-transport level of the blended heterojunction and a strong cycle dependence, consistent with changes in the hole trap state density. Passivation results in more reversible redox behaviors, consistent with more Ohmic-like contacts. Collectively, these results provide context toward the realization of durable organic photoelectrodes with optimized photovoltages and net solar-to-hydrogen conversion efficiencies via fundamental understanding of the rates of carrier generation, recombination, and transport in high-dielectric aqueous environments and opportunities to characterize buried interfaces under device-relevant electric fields.

buried interfaces↗
DOE OSTI2026

Pre-Equilibrium De-Excitations in Neutrino-Nucleus Interactions

Pre-equilibrium de-excitation is a well-established stage of nucleon-induced nuclear reactions but has not previously been incorporated into neutrino event generators. This work implements the Koning-Duijvestijn exciton model within the MARLEY neutrino event generator to simulate energy redistribution and particle emission before compound nucleus equilibrium. The implementation includes modular calculations of particle-hole state densities, internal transition rates, emission rates, and Monte Carlo cascade sampling, and is validated against the TALYS-2.2 nuclear reaction code. A new event record is also proposed to enable future uncertainty reweighting studies. This work establishes the first framework for studying pre-equilibrium effects in neutrino-induced reactions, with the goal of improving predictions of particle multiplicities, γ-ray production, and MeV-scale detector signals relevant to experiments such as DUNE.

Visser, Erin [Michigan State U., East Lansing (mai↗
DOE OSTI2026

Microscopy X-ray imaging enriched with small angle X-ray scattering for few nanometer resolution reveals shock waves and compression in intense short pulse laser irradiation of solids

Understanding how laser pulses compress solids into high-energy-density states requires diagnostics that simultaneously resolve macroscopic geometry and nanometer-scale structure. Here we present a combined X-ray imaging (XRM) and small-angle X-ray scattering (SAXS) approach that bridges this diagnostic gap. Using the Matter in Extreme Conditions end station at LCLS, we irradiated 25 μm copper wires with 45 fs, 0.9 J, 800 nm pulses at 3.5 × 10 19 W/cm 2 while probing with 8.2 keV XFEL pulses. XRM visualizes the evolution of ablation, compression, and inward-propagating fronts with ∼ 200 nm resolution, while SAXS quantifies their nanometer-scale sharpness via the time-resolved evolution of scattering streaks. The joint analysis reveals that an initially smooth compression steepens into a nanometer-sharp shock front after t sh ≈ (18 ± 3) ps, consistent with an analytical steepening model and hydrodynamic simulations. The front reaches a velocity of c sh ≈ 25 k m / s and a lateral width of several tens of microns, demonstrating direct observation of shock formation and decay at solid density for the first time with few-nanometer precision. This integrated XRM–SAXS method establishes a quantitative, multi-scale diagnostic of laser-driven shocks in dense plasmas relevant to inertial confinement fusion, warm dense matter, and planetary physics.

Kluge, Thomas [Helmholtz-Zentrum Dresden-Rossendor↗
DOE OSTI2026

Pre-Equilibrium De-Excitations for Neutrino-Nucleus Interactions

The Deep Underground Neutrino Experiment (DUNE) is sensitive to MeV-scale energy depositions from low-energy astrophysical neutrinos, including those from core-collapse supernovae. Interpreting these detector signals requires accurate modeling of the nuclear de-excitation that follows from the neutrino-nucleus interaction. The MARLEY (Model of Argon Reaction Low-Energy Yields) event generator specializes in the low-energy regime. MARLEY currently assumes the residual nucleus equilibrates immediately after the primary interaction. This omits the intermediate pre-equilibrium stage in which energy redistributes among nucleons until statistical equilibrium is reached. While pre-equilibrium effects are well established for nucleon-induced reactions, they have not previously been studied for neutrino-nucleus interactions. This work addresses that gap by implementing a two-component exciton model, which is the first dedicated treatment of pre-equilibrium de-excitation for neutrino-nucleus interactions, restructured around MARLEY's existing class hierarchy to prepare for direct integration, including particle-hole state densities, internal transition rates, and pre-equilibrium particle emission. The calculations show encouraging agreement with the TALYS-2.2 nuclear reaction code for neutron-nucleus interactions. We further propose a concrete integration path into the full MARLEY event generator, including derived class structure and an extended event record for pre-equilibrium vertices in support of future reweighting. Remaining work focuses on refining the emission width calculation, adding $\gamma$-ray emission, and completing this integration to quantify the impact of pre-equilibrium effects on the expected low-energy neutrino signals in DUNE and similar experiments.

Visser, Erin [Michigan State U., East Lansing (mai↗
DOE OSTI2025

Demonstrating the impact of ro-vibrationally excited H 2 on divertor detachment via population modelling

EDGE2D-EIRENE profiles characteristic for low recycling, high recycling and detached plasmas are postprocessed applying a vibrationally resolved collisional radiative (CR) model based on the Yacora solver to predict excited state densities and thus the molecular Lyman-, Werner- and Fulcher band emission intensities. These emission intensities are compared to predictions from the AMJUEL database standardly used in EIRENE. Lower Fulcher band and higher Lyman- and Werner band emission is predicted by the Yacora model in comparison to AMJUEL suggesting that inconsistent results from model based evaluation of measured spectra are expected depending on which model and molecular emission band is considered. The implications on detachment rate predictions are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Strong Proton‐Phonon Coupling Drives Fast Ion Transport in Perovskites

Conduction of protons in solids is a cooperative process propelled by phonons, with molecular details obscured by the irregular movements in the thermal bath. It is shown that substitution with Y forms an imaginary phonon mode, instrumental for the function as proton conductor and effectively lowering the activation barrier for proton transport. To untangle the interplay in the exemplary proton conductor BaSn 0.9 Y 0.1 O 3 , its crystallographic structure is determined with high resolution neutron diffractometry and its phonon density of states with density functional theory calculations, experimentally validated by element specific nuclear resonant vibration spectroscopy. Based on phonon analysis, a quantitative transport model is present, which predicts the activation energy and performance by the ratio of ionic radii. Rather than individual vibrational modes, it is the oxygen sub-lattice which exerts its momentum on the protons. The extent of this momentum transfer is governed by the ratio of ionic radii. This model extends the transition state theory by the phonon-phonon interaction and complements the previously proposed idea that lattice dynamics is decisive for proton transport and specifies which properties of the material exactly define the vibration properties.

activation barriers↗
DOE OSTI2025

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗
DOE OSTI2025

Repartitioning the Hamiltonian in many-body second-order Brillouin–Wigner perturbation theory: Uncovering new size-consistent models

Second-order Møller-Plesset perturbation theory is well-known as a computationally inexpensive approach to the electron correlation problem that is size-consistent with a size-consistent reference but fails to be regular. On the other hand, the less well-known many-body version of Brillouin-Wigner perturbation theory has the reverse properties: it is regular but fails to be size-consistent when used with the standard MP partitioning. Consequently, its widespread use remains limited. In this work, we analyze the ways in which it is possible to use alternative non-MP partitions of the Hamiltonian to yield variants of BW2 that are size-consistent as well as regular. We show that there is a vast space of such BW2 theories and also show that it is possible to define a repartitioned BW2 theory from the ground state density alone, which regenerates the exact correlation energy. We also provide a general recipe for deriving regular, size-consistent, and size-extensive partitions from physically meaningful components, and we apply the result to small model systems. The scope of these results appears to further set the stage for a revival of BW2 in quantum chemistry.

Ab initio perturbation↗