Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional approximation”

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

Coverage, repulsion, and reactivity of hydrogen on High-Entropy alloys

Modeling hydrogen evolution reaction (HER) activity probability on IrPdPtRhRu(1 1 1) high-entropy alloys. Determining hydrogen coverages based on ligand effects and generalized hydrogen–hydrogen repulsion. The rate of H 2 formation is highly impacted by the level of hydrogen coverage on the catalyst surface. In search of optimal catalytic properties high-entropy alloys (HEA) are promising candidates that utilize the compositional space of multiple elements. Based on simulations of HEA model (1 1 1) surfaces with a range of hydrogen coverages, distributions of binding energies are used to construct a framework that approximates the probability that adsorbed hydrogen may lead to the formation of H 2 as a function of applied potential. By optimizing the alloy compositions for the highest activity probability at given potentials the best and most efficient catalyst candidates for HER can be identified. Treating hydrogen–hydrogen repulsion effects and binding energy separately, we find that the repulsion is larger for HEAs than for pure metals. Differing isotherm slopes in the mean adsorption and desorption energies demonstrate a possible hysteresis for hydrogen adsorption on HEAs.

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

Characterization of chromium impurities in β -Ga2O3

Chromium is a common transition-metal impurity that is easily incorporated during crystal growth. It is perhaps best known for giving rise to the 694.3 nm (1.786 eV) emission in Cr-doped Al2O3, exploited in ruby lasers. Chromium has also been found in monoclinic gallium oxide, a wide-bandgap semiconductor being pursued for power electronics. In this work, we thoroughly characterize the behavior of Cr in Ga2O3 through theoretical and experimental techniques. β-Ga2O3 samples are grown with the floating zone method and show evidence of a sharp photoluminescence signal, reminiscent of ruby. We calculate the energetics of formation of Cr from first principles, demonstrating that Cr preferentially incorporates as a neutral impurity on the octahedral site. Cr possesses a quartet ground-state spin and has an internal transition with a zero-phonon line near 1.8 eV. By comparing the calculated and experimentally measured luminescence lineshape function, we elucidate the role of coupling to phonons and uncover features beyond the Franck–Condon approximation. The combination of strong emission with a small Huang–Rhys factor of 0.05 and a technologically relevant host material renders Cr in Ga2O3 attractive as a quantum defect.

Turiansky, Mark E. (ORCID:0000000291543582)↗
DOE OSTI2026

Determination of nuclear PDFs using Markov chain Monte Carlo methods

Global QCD analyses of nuclear parton distribution functions (nPDFs) have traditionally relied on the Hessian method for uncertainty estimation. However, the inherent Gaussian approximation and reliance on local curvature often prove insufficient for nPDF fits, which are frequently characterized by limited data constraints and non-Gaussian likelihoods. In this paper, we present the first nPDF determination based on Markov Chain Monte Carlo (MCMC) techniques, implemented within the nCTEQ framework using an adaptive Metropolis-Hastings algorithm. The MCMC approach enables a direct mapping of the posterior distribution and reveals a highly nontrivial parameter-space structure, including multiple modes and pronounced non-Gaussian behavior, particularly for the valence PDFs. We perform the first single-nucleus global analysis of lead PDFs using exclusively lead data and compare it to a multi-nuclei fit employing a standard analytic A dependence. The inclusion of lighter nuclei reduces quark uncertainties and modifies the shape of the lead PDFs, while leaving the gluon distribution largely unaffected. A complementary Hessian analysis exposes systematic limitations of the Gaussian approximation. Our results demonstrate that MCMC methods provide a more reliable framework for uncertainty quantification in nPDF determinations.

Derakhshanian, N. [Institute of Nuclear Physics Po↗
DOE OSTI2026

Massive tree-level splitting functions beyond kinematical limits

We present a compact form of the massive 1 → 3 tree-level QCD splitting functions and discuss a decomposition of the results in terms of lower-order expressions, scalar dipole antenna functions and pure higher-order remainders. The two-gluon radiator functions introduced in this context are novel and generalize expressions obtained from the double-soft approximation. Our results are obtained without reference to soft or quasicollinear limits.

Höche, Stefan [Fermi National Accelerator Laborato↗
DOE OSTI2025

Learning local and semi-local density functionals from exact exchange-correlation potentials and energies

Finding accurate exchange-correlation (XC) functionals remains the defining challenge in density functional theory (DFT). Despite 40 years of active development, attaining general purpose chemical accuracy is still elusive with existing functionals. We present a data-driven pathway to learn the XC functional by using the exact density, XC energy, and XC potential. While the exact densities are obtained from accurate configuration interaction (CI), the exact XC energies and XC potentials are obtained via inverse DFT calculations on the CI densities. We demonstrate how simple neural network (NN)–based local density approximation (LDA) and generalized gradient approximation (GGA), trained on just five atoms and two molecules, provide remarkable improvement in total energies and densities. Particularly, the NN-based GGA functional attains similar accuracy as the higher rung SCAN meta-GGA on various thermochemistry datasets. These results underscore the promise of using the XC potential in modeling XC functionals and can pave the way for systematic learning of increasingly accurate XC functionals.

Science & Technology - Other Topics↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Increasing Branching and Functionalization Decreases Crystallinity and Crystallization Rate in Linear EVOH

Here, we report the crystallization behavior of a linear poly(ethylene-co-vinyl alcohol) (LEVOH) under isothermal crystallization as a function of both OH incorporation and undercooling, and compare the results to conventional branched EVOH. This LEVOH is synthesized by post-polymerization functionalization of polycyclooctene, and it exhibits a half-crystallization time and primary crystallization rate nearly an order of magnitude faster than branched EVOH at ~ 6 mol% OH incorporation. While the LEVOH obtains approximately twice the extent of crystallization compared to the branched equivalent, as OH incorporation increases, the crystallization kinetics and crystallinity of LEVOH decrease. The crystal structure of LEVOH is orthorhombic at low functionalization (≤ 11 mol% OH), a mixture of orthorhombic and hexagonal at moderate functionalization (17 - 21 mol%), and hexagonal at high functionalization (23 %). The LEVOH crystallizes into plates, with crystallite widths ~ 20 - 60 times the crystallite thickness. Ultimately, we show LEVOH crystallizes faster and to a greater extent than commercial branched EVOH, a potential advantage of post-polymerization functionalization as part of polymer-to-polymer upcycling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multipole responses in fissioning nuclei and their uncertainties

Electromagnetic multipole responses are key inputs to model the structure, decay, and reactions of atomic nuclei. With the introduction of the finite amplitude method (FAM), large-scale calculations of the nuclear linear response in heavy deformed nuclei have become possible. This work provides a detailed study of multipole responses in actinide nuclei with Skyrme energy density functionals. We quantify both systematic and statistical uncertainties induced by the functional parametrization in FAM calculations. We also extend the FAM formalism to perform blocking calculations with the equal filling approximation for odd-mass and odd-odd nuclei, and analyze the impact of blocking configurations on the response. In conclusion, by examining the entire plutonium isotopic chain from the proton to the neutron drip line, we find a large variability of the response with the neutron number and study how it correlates with the deformation of the nuclear ground state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗
DOE OSTI2025

Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water–methane dimer

Fixed-node diffusion quantum Monte Carlo (FN-DMC) is a widely trusted many-body method for solving the Schrödinger equation, known for its reliable predictions of material and molecular properties. Furthermore, its excellent scalability with system complexity and near-perfect utilization of computational power make FN-DMC ideally positioned to leverage new advances in computing to address increasingly complex scientific problems. Even though the method is widely used as a computational gold standard, reproducibility across the numerous FN-DMC code implementations has yet to be demonstrated. This difficulty stems from the diverse array of DMC algorithms and trial wave functions, compounded by the method’s inherent stochastic nature. Here, this study represents a community-wide effort to assess the reproducibility of the method, affirming that yes, FN-DMC is reproducible (when handled with care). Using the water–methane dimer as the canonical test case, we compare results from eleven different FN-DMC codes and show that the approximations to treat the non-locality of pseudopotentials are the primary source of the discrepancies between them. In particular, we demonstrate that, for the same choice of determinantal component in the trial wave function, reliable and reproducible predictions can be achieved by employing the T-move, the determinant locality approximation, or the determinant T-move schemes, while the older locality approximation leads to considerable variability in results. These findings demonstrate that, with appropriate choices of algorithmic details, fixed-node DMC is reproducible across diverse community codes—highlighting the maturity and robustness of the method as a tool for open and reliable computational science.

Della Pia, Flaviano [Univ. of Cambridge (United Ki↗
DOE OSTI2025

Non-perturbative quarkonium dissociation rates in strongly coupled quark-gluon plasma

Heavy quarks and quarkonia are versatile probes of the transport properties of the hot QCD medium produced in ultra-relativistic heavy-ion collisions (URHICs). A robust description of heavy-flavor transport coefficients requires a microscopic approach that treats the open and hidden heavy-flavor sectors on the same footing. Here, we employ the quantum many-body T -matrix formalism to evaluate the dissociation rates of heavy quarkonia in the quark-gluon plasma (QGP). The basic ingredient is the heavy-light T -matrix, which utilizes a nonperturbative driving kernel constrained by lattice-QCD data. Its resummation in a ladder series provides a much enhanced interaction strength compared to a previously used perturbative coupling to the quasiparticle partons in the QGP. The in-medium quarkonium properties, particularly their temperature-dependent binding energies, are obtained from selfconsistent calculations with the same interaction kernel, including interference effects (also referred to as the imaginary part of the heavy-quark potential) as well as off-shell parton spectral functions. We systematically investigate the interplay of these effects and elaborate on the connections to the dipole approximation used in effective field theory.

effective field theories of QCD↗
DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
DOE OSTI2025

Measuring and unbiasing the BAO shift in the Ly α forest with AbacusSummit

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) places sub- per cent constraints on measurements of the Baryon Acoustic Oscillation (BAO) scaling parameters from the Ly $\alpha$ forest. However, no systematic error budget stemming from non-linearities in the three-dimensional clustering of the Ly $\alpha$ forest is included in the DESI-Ly $\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly $\alpha$ forest mocks produced on the N-body simulation suite AbacusSummit, which adopt the Fluctuating–Gunn–Peterson Approximation (FGPA). Specifically, we measure the Ly $\alpha$ autocorrelation and the Ly $\alpha$-quasar cross-correlation functions. To mitigate the noise, we adopt a linear control variates technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the autocorrelation, we detect a small positive shift in radial direction of $\Delta \alpha _{\parallel }= 0.35~{{\ \rm per\ cent}}$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha _\perp$. From the cross-correlation, we see a similar shift to $\Delta \alpha _\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha _{\perp }=\sim$0.25 per cent, at the 2$\sigma$ level. We also make a connection with the Ly $\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly $\alpha$ auto- and the Ly $\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly $\alpha$ forest analyses based on EFT. This work paves the way for novel full-shape analyses of the currently observing DESI and future surveys such as the PFS, WEAVE-QSO, and 4MOST.

Hadzhiyska, Boryana↗
DOE OSTI2025

Franck-Condon electron emission from polar semiconductor photocathodes

An analytical formulation of (optical-)phonon-mediated and momentum-resonant Franck-Condon emission of photoexcited electrons from polar semiconductors is shown to be very consistent with (i) the observed emission properties of a cesiated Ga⁢As⁢(001) photocathode at 808 nm [J. Phys. D: Appl. Phys. 54, 205301 (2021)] and (ii) the measured spectral emission properties of a Ga⁢N(0001) photocathode from just below its bandgap energy to 5 eV. The theoretical analysis in the parabolic band approximation predicts the form of both the quantum efficiency and mean transverse energy of photoemission as a function of the photocathode’s electron affinity and the electron temperature in the vicinity of its emission face. The good agreement between theory and experimental data also suggests that sub-10-nm rms surface roughness effects are not significant for polar semiconductor photocathodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mapping of the transient electric field causing spark suppression in Micromegas of the ATLAS New Small Wheels

Resistive strips Micromegas are employed in the ATLAS New Small Wheel project. They have already been installed and operate in the experimental cavern of the ATLAS experiment at CERN. This work attempts to describe the mechanism of the surface electric charge spread on a resistive strip of the Micromegas detector and the created transient electric field leading to the suppression of the spark rate in the amplification region. The electric field spatial mapping is calculated versus time by using the accurate 3-layer geometry as a function of the resistive layer thickness around its value used in Micromegas detectors. We also used the approximated method based on 2-layer geometry to investigate its discrepancy from the accurate one.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Global weak solutions for a nonlocal multispecies Fokker–Planck–Landau system

The global-in-time existence of weak solutions to a spatially homogeneous multispecies Fokker–Planck–Landau system for plasmas in the three-dimensional whole space is shown. The Fokker–Planck–Landau system is a simplification of the Landau equations assuming a linearized, velocity-independent, and isotropic kernel. The resulting equations depend nonlocally and nonlinearly on the moments of the distribution functions via the multispecies local Maxwellians. Furthermore, the existence proof is based on a three-level approximation scheme, energy and entropy estimates, as well as compactness results, and it holds for both soft and hard potentials.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Effects of screening and pressure ionization on the electron broadening of spectral lines in dense plasmas

Collisions between electrons and radiating atoms broaden spectral absorption and emission lines in dense plasmas. High densities also introduce screening and pressure ionization effects that distort the wave functions of both bound and free electrons. In order to study how dense plasma effects influence the electron broadening of spectral lines, this paper incorporates electron wave functions from an average-atom (AA) model to calculate the linewidth of the B III 2⁢𝑝−2⁢𝑠 transition at 𝑇 = 10 eV for mass densities ranging from 𝜌 = 10 −4 to 0.4⁢ g⁡/cc. The calculation method uses the impact approximation, allowing the linewidth to be written in terms of electron-collision cross sections and an interference term. Compared to an otherwise identical calculation that uses Coulomb free wave functions, the AA method is found to modify both the cross sections and the resulting linewidth at sufficiently high density by introducing screening and pressure-ionized bound states. Screening lowers the cross sections at low energies and near electron excitation thresholds, while pressure-ionized bound states introduce resonances into the continuum. Thus, as the density increases, the relative linewidth between the AA and Coulomb calculations follows a general decrease because of screening, with sharp increases at various intervals due to pressure ionization. Finally, the AA results are also compared with a common approach to introduce screening through the interaction potential and reduced models that use the Bethe formula for the inelastic electron-collision cross sections.

electronic excitation & ionization↗
DOE OSTI2024

Lattice B -field correlators for heavy quarks

We analyze the color-magnetic (or “ B ”) field two-point function that encodes the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at 1.5 T c , using a range of gradient flow times for noise suppression and operator renormalization. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations to next-to-leading order of this correlation function, matching gradient-flowed correlators to MS ¯ , are used to resolve nontrivial renormalization issues. We perform a spectral reconstruction based on perturbative model fits to estimate the coefficient κ B of the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The approach we present here yields high-precision data for the correlator with all renormalization issues incorporated at next-to-leading order and is also applicable for actions with dynamical fermions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗