Engineering PapersSearch

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

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory

Trends in spatial correlations, two-phase coexistence, and criticality in a class of type-2 Schloegl models for autocatalysis

A class of type-2 Schloegl models is considered for particles on a square lattice with variable-range cooperativity. These models involve: (i) spontaneous particle annihilation at rate p; (ii) autocatalytic particle creation at unoccupied sites (i, j) with n ⩾ 2 particles within a specified neighborhood, Ω 𝑁 (i, j), of sites at rate $k_n$ = $\frac{^{(^n_2)}}{_{(^N_2)}}$ = $\frac{n{(n-1)}}{_{N(N-1)}}$; and (iii) possible spontaneous particle creation at unoccupied sites at “small” rate ɛ ⩾ 0. In some cases, Ω 𝑁 just includes all symmetry-equivalent sites at a single specific distance 𝑑 (in units of lattice constants) from the unoccupied site, e.g., 𝑑 = 1 (nearest-neighbor sites) where 𝑁 = 4, or 𝑑 = √5 (or √13 or…) where 𝑁 = 8. In other cases, Ω 𝑁 includes sites multiple distances from the unoccupied site, e.g., 𝑑 = {1,√2}, where 𝑁 = 8. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition between high- and low-density states below a critical point, ɛ < ɛ c , with generic two-phase coexistence (2PC) at least for smaller 𝑁. With some exceptions, there is an approach toward mean-field behavior with increasing 𝑁 (so the regime of generic 2PC shrinks, and ɛ c approaches the mean-field value of 1/27). Additional insight into trends is provided by analysis of the exact master equations for the models via hierarchical truncation. These truncations utilize suitably tailored pair approximations which reflect the dominant nonequilibrium spatial correlations. These correlations in turn are shown to reflect the details of the autocatalytic particle creation process. For spatially heterogeneous states, the truncations produce coupled sets of lattice differential equations (LDE) which can describe orientation-dependent propagation of an interface between high- and low-density steady states for ɛ < ɛ c . Pair approximation values of 𝑝 = 𝑝 eq where the interface is stationary, and its orientation-dependence, are in semiquantitative agreement with KMC results. In conclusion, this comparison accounts for propagation failure in the LDE which complicates interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics

Enhanced quantum efficiency from optical interference in alkali antimonide photocathodes: Modeling and experimental results

We present measurements of enhanced quantum efficiency (QE) in thin film alkali antimonide photocathodes from optical interference in the cathode-substrate multilayer. Modulations in the spectral response are observed over a range of visible wavelengths and are shown to increase the QE by more than a factor of two at specific wavelengths. We present a model describing the QE modulations based on the three step photoemission process incorporating cases of both constant density of states and density functional theory-derived density of states and show that the calculated results are in good agreement with the measurements. Model predictions demonstrate that QE can be enhanced by more than a factor of 5 by optimization of cathode and substrate layer thicknesses. Additionally, these calculations reveal that optical interference can yield higher quantum efficiencies in thin films compared to thick, optically dense films. We model the QE vs excitation wavelength of multiple alkali antimonide compounds at different thicknesses. We then discuss the advantages of this interference effect for electron accelerators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

One Beam, Dual Insights: Simultaneous Chemical and Structural Changes in Nanopatterned Ceria under Reaction Conditions

Ceria’s interaction with hydrogen can proceed through multiple chemical forms (hydride, hydroxyl, and oxyhydroxide-like), with consequences for the oxidation state, density, and morphology that are rarely tracked in the same evolving state. Here, in this work, we show that under mild H 2 (and H 2 and CO 2 ) environments nanopatterned ceria undergoes oxidation-state changes accompanied by hydrogen incorporation that increases the effective electron density, establishing the following order: CeO 2 Hy > CeO 2 > CeO 2–x Hy > CeO 2–x . In parallel, the surface roughens in a chemically specific manner, with the largest changes coinciding with conditions where incorporated hydrogen is driven to react with oxygen supplied either by air exposure between experiments or by added CO 2 . We obtained these insights by using a single X-ray beam to simultaneously perform ambient-pressure X-ray photoelectron spectroscopy and grazing-incidence X-ray scattering on the same sample spot. Single-mode measurements can miss key ceria–H 2 transformations relevant to optimizing ceria-based hydrogenation catalysts and supports.

anions

Spatio-Kinetic Electronic Heterogeneity in MoS 2 Revealed by Ultrafast Electron Microscopy: Implications for Optoelectronics

There is broad interest in layered, two-dimensional (2D) semiconductors owing to their high optical absorption and ability to heterogeneously integrate with various substrates. In such layered semiconductors, mesoscopic imperfections, such as surface step edges, may affect charge carrier kinetics and device performance, but studying their effects on carrier transport with sufficient spatial and temporal resolution is a challenge. Using scanning ultrafast electron microscopy, we image the near-surface charge carrier kinetics with μm/ps spatiotemporal resolution in single-crystal MoS 2 samples with surface steps tens of nanometers in height. We find an approximately 3-fold increase in surface trap density near such surface steps compared to the uniform, flat surface regions. Despite the difference in trap state density, filling of the traps occurs over the same time scale (∼100 ps), but carrier recovery is much slower near the step compared to the uniform regions of the material. Furthermore, these results provide an approach to measure ultrafast dynamics and determine the required uniformity of layered materials depending on the speed limit constraints for specific applications.

SEM

The Early History of the Quark-Gluon Plasma

We present the historical antecedents to the field of relativistic heavy ion physics, beginning with early attempts to model the strong interaction and ending with the endorsement of a relativistic heavy ion collider in the 1983 U.S. Long-Range Plan for Nuclear Science. Particular attention is paid to two major themes: 1) A program to study high density states of nuclear matter emerging from the 1974 Bear Mountain conference and 2) Efforts to understand the predictions of QCD for matter at high densities and/or temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

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

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

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

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

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.

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

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

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

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

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