SEARCH · Engineering Papers

Results for “density evolution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.

Johnson, Praise Noah↗

DOE OSTI2024

Regulation of Electron Cloud Density by Electronic Effects of Substituents to Optimize Photocatalytic H 2 Evolution of Carbon Dots

Carbon dots (CDs) have a wide light absorption range, which are promising candidates for photocatalytic water splitting H 2 evolution, but they show low activity for hydrogen evolution reaction (HER) due to the slow transfer efficiencies of photogenerated carriers. The electron cloud density of photocatalyst is essential for the separation and migration of photogenerated carriers. In this study, the effect mechanism of electron cloud density regulated by the substituents with different electronic effects on the photocatalytic HER of glucose-based CDs is clarified. CDs-SO 3 H and CDs-OH are first obtained by introducing electron-drawing group (–SO 3 H) and electron-donating group (–OH) using a tailoring post-processing strategy. Experimental results show that the H 2 yield catalyzed by CDs-SO 3 H is 89.95 µmol•g –1 in 4 h, which is about four times that of CDs, and 7.4 times that of CDs-OH. The –SO 3 H induces a much negative energy band and high electron cloud density on CDs edges, promoting the separation and transfer of photogenerated carriers; while the CDs-OH exhibits a high positive charge density and an upward energy band, hindering the surface complexation of electron-hole pairs and HER. Furthermore, this study will provide an insight into the design of CDs catalysts with efficient photocatalytic HER at a molecular level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Unraveling the electronic structure and magnetic transition evolution across monolayer, bilayer, and multilayer ferromagnetic Fe 3 GeTe 2

Two-dimensional (2D) van der Waals (vdW) magnets have sparked widespread attention due to their potential in spintronic applications as well as in fundamental physics. Ferromagnetic vdW compound Fe 3 GeTe 2 (FGT) and its Ga variants have garnered significant interest due to their itinerant magnetism, correlated states, and high magnetic transition temperature. Experimental studies have demonstrated the tunability of FGT’s Curie temperature, T C , through adjustments in quintuple layer numbers (QL) and carrier concentrations, n. However, the underlying mechanism remains elusive. In this study, we employ molecular beam epitaxy (MBE) to synthesize 2D FGT films down to 1 QL with precise layer control, facilitating an exploration of the band structure and the evolution of itinerant carrier density. Angle-resolved photoemission spectroscopy (ARPES) reveals significant band structure changes at the ultra-thin limit, while first-principles calculations elucidate the band evolution from 1 QL to bulk, largely governed by interlayer coupling. Additionally, we find that n is intrinsically linked to the number of QL and temperature, with a critical value triggering the magnetic phase transition. Our findings underscore the pivotal role of band structure and itinerant electrons in governing magnetic phase transitions in such 2D vdW magnetic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Visualizing the impact of quenched disorder on 2D electron Wigner solids

Electron Wigner solids (WSs) provide an ideal system for understanding the competing effects of electron–electron and electron–disorder interactions, a central unsolved problem in condensed matter physics. Progress in this topic has been limited by a lack of single-defect-resolved experimental measurements as well as accurate theoretical tools to enable realistic experiment/theory comparison. Here, in this study, we overcome these limitations by combining atomically resolved scanning tunnelling microscopy (STM) with neural-quantum-state quantum Monte Carlo (NQS-QMC) simulation of disordered 2D electron WSs to discover new disorder-induced physical regimes of correlated electron behaviour. STM was used to image the electron density (n e )-dependent evolution of electron WSs in gate-tunable bilayer MoSe 2 (BL-MoSe 2 ) devices with varying long-range (n LR ) and short-range (n SR ) disorder densities. These images were compared with NQS-QMC simulations using realistic disorder maps extracted from experiment, thus allowing the roles of different disorder types to be disentangled. We identify two distinct physical regimes for disordered electron WSs that depend on n SR . For n SR ≲ n e , the WS behaviour is dominated by long-range disorder and features extensive mixed solid–liquid phases, a new type of local re-entrant melting/crystallization and prominent Friedel oscillations. By contrast, when n SR ≫ n e , these features are suppressed and a more robust amorphous WS phase emerges that persists to higher n e , highlighting the importance of short-range disorder in this regime. Our work establishes a powerful framework for studying disordered quantum solids through a combined experimental–theoretical approach.

Ge, Zhehao [University of California, Berkeley, CA↗

DOE OSTI2025

Quantum entanglement correlations in double quark PDFs

Methods from Quantum Information Theory are used to scrutinize quantum correlations encoded in the two-quark density matrix over light-cone momentum fractions x 1 and x 2 . A nonperturbative three quark model light-cone wavefunction predicts significant nonclassical correlations associated with the “entanglement negativity” measure for asymmetric and small quark momentum fractions. We perform one step of QCD scale evolution of the entire density matrix, not just its diagonal (dPDF), by computing collinearly divergent corrections due to the emission of a gluon. Finally, we present first qualitative numerical results for single-step scale evolution of quantum entanglement correlations in double quark PDFs. At a higher Q 2 scale, the nonclassical correlations manifest in the dPDF for nearly symmetric momentum fractions.

Entanglement in field theory↗

DOE OSTI2026

Mechanism-Resolved PFM of Ferroionic and Ferroelectric Responses in Thickness-Gradient Hf 0.5 Zr 0.5 O 2 Libraries

Resolving growth mechanisms and thickness evolution of functional properties is one of the key tasks in materials discovery and optimization involving thin-film materials, traditionally requiring significant experimental budgets. Here we introduce the combination of thickness-gradient libraries and automated scanning probe microscopy as a systematic pathway to elucidate growth modes and disentangle ferroelectric and electrochemical contributions in ferroelectric thin films. As a model system, we explore the Hf 0.5 Zr 0.5 O 2 (HZO) gradient thin films grown on La x Sr 1-x MnO 3 (LSMO) bottom electrode thin films. Automated piezoresponse force microscopy, spectroscopy, and lithography reveals that irreversible topographic deformation arises from electrochemical activity at the LSMO surface, whereas reversible phase inversion in HZO reflects ferroelectric switching. Automated topography height-map scans are used to further quantify nucleation density, particle-size evolution, and roughness correlations across the thickness-gradient, demonstrating that improved plume stabilization during growth suppresses interfacial reactions and promotes dense, fine-grained HZO conducive to ferroelectric phase formation. This combined materials-engineering and automated-SPM framework establishes a platform for high-throughput, mechanism-resolved characterization of ferroionic and ferroelectric responses in complex oxide films.

FOS: Physical sciences↗

DOE OSTI2026

Structural and Electrocatalytic Studies of Pulsed Laser Deposited Epitaxial RuO 2 Thin Films

Two sets of high-quality epitaxial ruthenium oxide (RuO 2 ) thin films with different thicknesses were synthesized in situ on cost-competitive sapphire substrates by using a pulsed laser deposition technique. The first set of films, with a thickness of 40 nm and a sheet resistance of 15.7 Ω/□, was prepared using 2100 laser pulses, while the second set of films, with a thickness of 87 nm and a sheet resistance of 6.7 Ω/□, was prepared using 4800 laser pulses. All other deposition parameters were kept the same. The post-deposition structural and morphological measurements showed that both sets of films grew at the same growth rate, had the same crystallinity, similar grain boundary density, and slightly different surface roughness. The thicker RuO 2 films achieved an overpotential of 280 mV for the oxygen evolution reaction at a current density of 100 μA/cm 2 , comparable to or exceeding the performance of films grown on more expensive substrates. Comparatively, the thinner RuO 2 films, which have a significantly higher charge transfer resistance (250 Ω versus 100 Ω for thick films), display a higher overpotential of 320 mV. These results indicate that the lower electrical resistance of thicker films promotes charge transfer through the film body, leading to superior electrocatalytic properties.

charge-transfer↗

DOE OSTI2024

Bottomonium suppression from the three-loop QCD potential

We compute the suppression of bottomonium in the quark-gluon plasma using the three-loop QCD static potential. The potential describes the spin-averaged bottomonium spectrum below threshold with a less than 1% error. Within potential nonrelativistic quantum chromodynamics and an open quantum systems framework, we compute the evolution of the bottomonium density matrix. The values of the quarkonium transport coefficients are obtained from lattice QCD measurements of the bottomonium in-medium width and thermal mass shift; we additionally include for the first time a vacuum contribution to the dispersive coefficient γ . Using the three-loop potential and the values of the heavy quarkonium transport coefficients, we find that the resulting bottomonium nuclear modification factor is consistent with experimental observations, while at the same time reproducing the lattice measurements of the in-medium width. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Tracking surface charge dynamics on single nanoparticles

Surface charges play a fundamental role in physics and chemistry, in particular in shaping the catalytic properties of nanomaterials. However, tracking nanoscale surface charge dynamics remains challenging due to the involved length and time scales. Here, we demonstrate time-resolved access to the nanoscale charge dynamics on dielectric nanoparticles using reaction nanoscopy. We present a four-dimensional visualization of the spatiotemporal evolution of the charge density on individual SiO 2 nanoparticles under strong-field irradiation with femtosecond-nanometer resolution. The initially localized surface charges exhibit a biexponential redistribution over time. Our findings reveal the influence of surface charges on surface molecular bonding through quantum dynamical simulations. We performed semi-classical simulations to uncover the roles of diffusion and charge loss in the surface charge redistribution process. Understanding nanoscale surface charge dynamics and its influence on chemical bonding on a single-nanoparticle level unlocks an increased ability to address global needs in renewable energy and advanced health care.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

A Multifunctional Isostructural Bilayer Oxygen Evolution Electrode for Durable Intermediate-Temperature Electrochemical Water Splitting

The overarching goal of the proposed research is to address SOEC’s degradation problem by advancing a new isostructural highly electrocatalytically active bilayer oxygen evolution reaction (OER) electrode, consisting of a LSCF (La 1-x Sr x Co 1-y Fe y O 3-δ ) core and a SCT (SrCo 0.9 Ta 0.1 O 3-δ ) shell, to achieve high and sustainable rate of oxygen evolution matching operating current densities without encountering delamination. To realize this goal, the project has adopted a combined experimental and theoretical approach to conduct research in the following six areas closely associated with SOPO tasks: 1) Development of electrocatalytically active bilayer oxygen electrodes (SOPO task-1) 2) Development of new symmetric three electrode cell (STEC) methodology to extract electrokinetic data of oxygen electrodes (SOPO task-2) 3) Quantification of electrokinetics of bilayer oxygen electrodes and correlation with degradation and delamination (SOPO task-2) 4) Performances of bilayer oxygen electrodes under fuel cells and electrolyzers modes (SOPO task-3) 5) Microscale modeling of oxygen electrode/electrolyte interface in solid oxide electrolysis cells (SOPO task-4) 6) Prediction of crack growth rate at oxygen electrode/electrolyte interface in solid oxide electrolysis cells (SOPO task-4)

08 HYDROGEN↗

DOE OSTI2026

Green Methanol via an Integrated Direct Air Capture, CO 2 Electrolyzer, and Hydrogenation Reactor

This project pioneered a groundbreaking reactor design to produce green methanol by harnessing the electrochemical CO 2 reduction reaction (eCO 2 RR), a cornerstone of power-to-fuels technology. The effort integrated three innovative technologies to achieve carbon-neutral methanol production at a target cost of under $\$$800/ton: 1. Direct Air Capture (DAC): Using a cutting-edge sorbent material developed at Holocene, scalable models were developed to integrate captured atmospheric CO₂ into the reactor system. 2. Intermediate-Temperature CO 2 Electrolyzer: Developed by the University of Tennessee (UTK), this electrolyzer utilizes a cost-effective, proton-conducting solid acid electrolyte (CsH 2 PO 4 , CDP) and a mixed-metal oxide cathode. It achieves high faradaic efficiencies (>98%) by effectively suppressing hydrogen evolution at high current densities, converting CO 2 to CO with remarkable selectivity. 3. Catalysis and Reactor Engineering: Oak Ridge National Laboratory (ORNL) contributed world-class expertise in heterogeneous catalysis and reactor design. Their advanced ASPEN modeling drove systems integration and supported techno-economic and life cycle analyses. This effort was further bolstered by partnerships with industry leaders Air Company and Plug Power, who provided critical guidance on scaling, systems engineering, and the integration of water electrolyzers into large-scale operations. During Phase 1, the team focused on modeling and validating a lab-scale reactor demonstrating the feasibility of the integrated approach. Key accomplishments include a 52% increase in current density at 0.8 V while maintaining >98% CO faradaic efficiency, successful 10× scale-up of the electrolyzer with performance within 5% of coin-cell results, best-in-class durability (168-hour test at 0.6 V with 0.14 mA/cm 2 -h degradation), validated TEA confirming the $\$$800/ton methanol target, and completed preliminary LCA showing potential for net-negative GHG emissions under renewable energy scenarios.

10 SYNTHETIC FUELS↗

DOE OSTI2025

Disentangling the intertwined orders in a magnetic kagome metal

Intertwined orders appear when multiple orders are strongly interacting, and kagome metals have emerged as new platforms to explore exotic phases. FeGe has been found to develop a charge density wave (CDW) order within magnetic phase, suggesting an intricate interplay of the lattice, charge, and spin degrees of freedom. Recently, postgrowth annealing has been proposed to tune the CDW order from long-range to complete suppression, offering a tuning knob for the CDW order. Here, by comparing the electronic structures of FeGe subjected to different annealing conditions and distinct CDW properties, we report spectral evolution associated with the lattice and spin degrees of freedom. We find band evolution linked to a spin density wave (SDW) order present in both samples with and without CDW order, and another evolution connected to the lattice distortions that onset with the long-range CDW order and revert with the SDW order. Our results reveal a rare competitive cooperation of the lattice, spin, and charge in FeGe.

Oh, Ji Seop↗

DOE OSTI2026

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.

Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗

DOE OSTI2025

Computing Nonlinear Power Spectra Across Dynamical Dark Energy Model Space with Neural ODEs

I show how to compute the nonlinear power spectrum across the entire $w(z)$ dynamical dark energy model space. Using synthetic ΛCDM data, I train a neural ordinary differential equation (ODE) to infer the evolution of the nonlinear matter power spectrum as a function of the background expansion and mean matter density across ∼9 Gyr of cosmic evolution. After training, the model generalises to any dynamical dark energy model parameterised by $w(z)$. With little optimisation, the neural ODE is accurate to within 4% up to $k = 5\, h\, {\mathrm Mpc}^{−1}$. Unlike simulation rescaling methods, neural ODEs naturally extend to summary statistics beyond the power spectrum that are sensitive to the growth history.

cosmology↗

DOE OSTI2026

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗

DOE OSTI2025

Close-up view of the interplay between lattice distortions and charge density waves: Case study on rare-earth nickel carbides

Using variable temperature total x-ray scattering and large-scale structure modeling, we study the temperature and composition evolution of lattice distortions and charge density waves (CDWs) in the archetypal strongly correlated system RENi⁢ C 2 (RE=Dy, Tb, Gd, and Sm). Joint analysis of reciprocal and real space data reveals the presence of strong lattice distortions that, depending on the RE species and temperature, appear periodic or remain local, forming a complex CDW phase diagram. Apparently, the CDWs in RENi⁢ C 2 materials arise from a common to the system lattice instability involving diverse distortion modes that persist when magnetic order also sets in at very low temperature. In conclusion, the results support the notion that intrinsic lattice distortions are important to the emergence of exotic electronic phases in strongly correlated systems and call for their further investigation using the advanced approach adopted here.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗