Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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

Resolving length-scale-dependent transient disorder through an ultrafast phase transition

Abstract Material functionality can be strongly determined by structure extending only over nanoscale distances. The pair distribution function presents an opportunity for structural studies beyond idealized crystal models and to investigate structure over varying length scales. Applying this method with ultrafast time resolution has the potential to similarly disrupt the study of structural dynamics and phase transitions. Here we demonstrate such a measurement of CuIr 2 S 4 optically pumped from its low-temperature Ir-dimerized phase. Dimers are optically suppressed without spatial correlation, generating a structure whose level of disorder strongly depends on the length scale. The redevelopment of structural ordering over tens of picoseconds is directly tracked over both space and time as a transient state is approached. This measurement demonstrates the crucial role of local structure and disorder in non-equilibrium processes as well as the feasibility of accessing this information with state-of-the-art XFEL facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Impact of LHC precision measurements of inclusive jet and dijet production on the CTEQ-TEA global PDF fit

In this study, we investigate the impact of new Large Hadron Collider (LHC) inclusive jet and dijet measurements on parton distribution functions (PDFs) that describe the proton structure, with a particular focus on the gluon distribution at large momentum fraction, x , and the corresponding partonic luminosities. We assess constraints from these datasets using next-to-next-to-leading-order (NNLO) theoretical predictions, accounting for a range of uncertainties from scale dependence and numerical integration. From the scale choices available for the calculations, our analysis shows that the central predictions for inclusive jet production show a smaller scale dependence than dijet production. We examine the relative constraints on the gluon distribution provided by the inclusive jet and dijet distributions, and also explore the phenomenological implications for inclusive H , t t ¯ , and t t ¯ H production at the LHC at 14 TeV. Published by the American Physical Society 2025

Ablat, Alim (ORCID:0009000519877755)↗
DOE OSTI2026

A Review of Advanced Characterization Methods and Techniques for Amorphous and Poorly Crystalline Materials Applicable to Cementitious Waste Forms

This review covers advanced characterization methods of cementitious materials. The particular focus is to understand cementitious compositions with amorphous content because conventional techniques such as X-ray diffraction are more applicable to materials with long range order. Also of interest is advanced characterization that provides deeper insights into properties such as porosity. There are five groups of characterization methods discussed in this review which are: electron beam/microscopy, X-ray, neutron diffraction and scattering, nuclear magnetic resonance, and optical techniques. The techniques can be further grouped by the categories of nanoscale structural understanding (transmission electron microscopy and pair distribution function), microscale structural understanding (focused ion beam and x-ray tomography), microscopic elemental distributions (energy dispersive spectroscopy and X-ray fluorescence), and chemical binding information (electron energy loss spectroscopy, x-ray absorption spectroscopy, nuclear magnetic resonance, Raman spectroscopy, and infrared spectroscopy), with advantages and disadvantages of each method discussed. The recommendation for further research is to begin with techniques that are easier to prepare samples for, measure, and analyze, such as the optical techniques of Raman and infrared spectroscopy. If more information is needed, pursue collaborations, other techniques, or proposals for beam time at the synchrotron sources.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

A functional global sensitivity measure and efficient reliability sensitivity analysis with respect to statistical parameters

Sensitivity analysis and reliability assessment are two important aspects of structural and system safety. Epistemic uncertainty with respect to probabilistic model of input parameters due to lack of knowledge is present in many scarce-data applications and complicates the characterization of uncertainty in model response. In this article, we present two importance measures to evaluate the impact of distribution parameters on the probability distribution function (PDF) of the output and the failure probability. The epistemic uncertainty associated with the distribution parameters is modeled as random variables. Additionally, a modified extended polynomial chaos expansion (MEPCE) approach is introduced in which aleatory and epistemic random variables are modeled and propagated simultaneously while allowing the separate assessment for any single epistemic variable. A MEPCE-based kernel density estimation (KDE) construction provides a composite map from each epistemic variable to the response PDF. The functional global sensitivity index of the PDF with respect to the distribution parameters is thus derived, as a function of output, which is both more informative and more efficient than standard scalar sensitivity measures. Reliability sensitivity indices can be readily evaluated by integrating the global sensitivity index function over the failure zone. Three illustrative examples are used to demonstrate the proposed methodology.

42 ENGINEERING↗
DOE OSTI2024

Algebraic model to study the internal structure of pseudoscalar mesons with heavy-light quark content

The internal structure of all lowest-lying pseudoscalar mesons with heavy-light quark content is studied in detail using an algebraic model that has been applied recently, and successfully, to the same physical observables of pseudoscalar and vector mesons with hidden-flavor quark content, from light to heavy quark sectors. The algebraic model consists on constructing simple and evidence-based of the meson’s Bethe-Salpeter amplitude (BSA) and quark’s propagator in such a way that the Bethe-Salpeter wave function (BSWF) can then be readily computed algebraically. Its subsequent projection onto the light front yields the light front wave function (LFWF) whose form allows us a simple access to the valence-quark parton distribution amplitude (PDA) by integrating over the transverse momentum squared. We exploit our current knowledge of the PDAs of lowest-lying pseudoscalar heavy-light mesons to compute their generalized parton distributions (GPDs) through the overlap representation of LFWFs. From these three dimensional knowledge, different limits/projections lead us to deduce the related parton distribution functions (PDFs), electromagnetic form factors (EFFs), and impact parameter space GPDs (IPS-GPDs). When possible, we make explicit comparisons with available experimental results and earlier theoretical predictions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Non-singlet quark helicity $\mathrm{PDF}$s of the nucleon from pseudo-distributions

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional z 2 contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional z 2 contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale μ 2 . The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned z 2 contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Node Distortion as a Tunable Mechanism for Negative Thermal Expansion in Metal–Organic Frameworks

Chemically functionalized series of metal–organic frameworks (MOFs), with subtle differences in local structure but divergent properties, provide a valuable opportunity to explore how local chemistry can be coupled to long-range structure and functionality. Using in situ synchrotron X-ray total scattering, with powder diffraction and pair distribution function (PDF) analysis, we investigate the temperature dependence of the local- and long-range structure of MOFs based on NU-1000, in which Zr 6 O 8 nodes are coordinated by different capping ligands (H 2 O/OH, Cl – ions, formate, acetylacetonate, and hexafluoroacetylacetonate). We show that the local distortion of the Zr 6 nodes depends on the lability of the ligand and contributes to a negative thermal expansion (NTE) of the extended framework. Using multivariate data analyses, involving non-negative matrix factorization (NMF), we demonstrate a new mechanism for NTE: progressive increase in the population of a smaller, distorted node state with increasing temperature leads to global contraction of the framework. The transformation between discrete node states is noncooperative and not ordered within the lattice, i.e., a solid solution of regular and distorted nodes. Density functional theory calculations show that removal of ligands from the node can lead to distortions consistent with the Zr···Zr distances observed in the experiment PDF data. Control of the node distortion imparted by the nonlinker ligand in turn controls the NTE behavior. Furthermore, these results reveal a mechanism to control the dynamic structure of MOFs based on local chemistry.

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

Small-𝑥 asymptotics of the leading-twist flavor-singlet quark TMDs

In this paper, we investigate the small-𝑥 behavior of the flavor-singlet, leading-twist quark transverse-momentum-dependent parton distribution functions (TMDs) using the light-cone operator treatment. This formalism allows us to express TMD operators at small 𝑥 in terms of polarized dipole amplitudes, enabling a systematic approach to their small-𝑥 evolution. We derive the evolution equations for these TMDs and solve them within the large-𝑁 𝑐 approximation under the linearized, double-logarithmic approximation, where 𝑁 𝑐 represents the number of quark colors. Expanding on previous work on unpolarized and helicity TMDs, we present the small-𝑥 asymptotics for a comprehensive set of TMDs, including the Sivers function, helicity worm-gear, transversity, pretzelosity, Boer-Mulders, and transversity worm-gear distributions. Our results provide a complete picture of the small-𝑥 asymptotic behavior for all leading-twist flavor-singlet quark TMDs. We also discuss the implications of our findings for phenomenological applications and outline potential avenues for further research, particularly in understanding nonlinear effects and extending beyond the double-logarithmic approximation and large-𝑁 𝑐 approximations.

Adamiak, Daniel [Thomas Jefferson National Acceler↗
DOE OSTI2024

Large-scale Linear Magnetic Holes with Magnetic Mirror Properties in Hybrid Simulations of Solar Wind Turbulence

Magnetic holes (MHs) are coherent magnetic field dips whose size ranges from fluid to kinetic scale, ubiquitously observed in the heliosphere and in planetary environments. Despite the long-standing effort in interpreting the abundance of observations, the origin and properties of MHs are still debated. In this Letter, we investigate the interplay between plasma turbulence and MHs, using a 2D hybrid simulation initialized with solar wind parameters. We show that fully developed turbulence exhibits localized elongated magnetic depressions, whose properties are consistent with linear MHs frequently encountered in space. The observed MHs develop self-consistently from the initial magnetic field perturbations by trapping hot ions with large pitch angles. Ion trapping produces an enhanced perpendicular temperature anisotropy that makes MHs stable for hundreds of ion gyroperiods, despite the surrounding turbulence. We introduce a new quantity, based on local magnetic field and ion temperature values, to measure the efficiency of ion trapping, with potential applications to the detection of MHs in satellite measurements. We complement this method by analyzing the ion velocity distribution functions inside MHs. Our diagnostics reveal the presence of trapped gyrotropic ion populations, whose velocity distribution is consistent with a loss cone, as expected for the motion of particles inside a magnetic mirror. Our results have potential implications for the theoretical and numerical modeling of MHs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Extreme Risk Mitigation in Reinforcement Learning using Extreme Value Theory

Risk-sensitive reinforcement learning (RL) has garnered significant attention in recent years due to the growing interest in deploying RL agents in real-world scenarios. A critical aspect of risk awareness involves modelling highly rare risk events (rewards) that could potentially lead to catastrophic outcomes. These infrequent occurrences present a formidable challenge for data-driven methods aiming to capture such risky events accurately. While risk-aware RL techniques do exist, they suffer from high variance estimation due to the inherent data scarcity. Our work proposes to enhance the resilience of RL agents when faced with very rare and risky events by focusing on refining the predictions of the extreme values predicted by the state-action value distribution. To achieve this, we formulate the extreme values of the state-action value function distribution as parameterized distributions, drawing inspiration from the principles of extreme value theory (EVT). We propose an extreme value theory based actor-critic approach, namely, Extreme Valued Actor-Critic (EVAC) which effectively addresses the issue of infrequent occurrence by leveraging EVT-based parameterization. Importantly, we theoretically demonstrate the advantages of employing these parameterized distributions in contrast to other risk-averse algorithms. Our evaluations show that the proposed method outperforms other risk averse RL algorithms on a diverse range of benchmark tasks, each encompassing distinct risk scenarios.

Wang, Yu↗
DOE OSTI2024

3D Lead-Organoselenide-Halide Perovskites and their Mixed-Chalcogenide and Mixed-Halide Alloys

Here, we incorporate Se into the 3D halide perovskite framework using the zwitterionic ligand: SeCYS ( + NH 3 (CH 2 ) 2 Se - ), which occupies both the X - and A + sites in the prototypical ABX 3 perovskite. The new organoselenide-halide perovskites: (SeCYS)PbX 2 (X=Cl, Br) expand upon the recently discovered organosulfide-halide perovskites. Single-crystal X-ray diffraction and pair distribution function analysis reveal the average structures of the organoselenide-halide perovskites, whereas the local lead coordination environments and their distributions were probed through solid-state 77 Se and 207 Pb NMR, complemented by theoretical simulations. Density functional theory calculations illustrate that the band structures of (SeCYS)PbX 2 largely resemble those of their S analogs, with similar band dispersion patterns, yet with a considerable band gap decrease. Optical absorbance measurements indeed show band gaps of 2.07 and 1.86 eV for (SeCYS)PbX 2 with X=Cl and Br, respectively. We further demonstrate routes to alloying the halides (Cl, Br) and chalcogenides (S, Se) continuously tuning the band gap from 1.86 to 2.31 eV–straddling the ideal range for tandem solar cells or visible-light photocatalysis. The comprehensive description of the average and local structures, and how they can fine-tune the band gap and potential trap states, respectively, establishes the foundation for understanding this new perovskite family, which combines solid-state and organo-main-group chemistry.

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

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Semi-classical approximation description of static properties of nuclei

In this report, we present recent development in semi-classical description of static properties of nuclei. By employing the Wigner transform, we derive simple semiclassical approximations for evaluating properties of nuclear system. We present results of calculations of static properties, using the energy density functional associated with Skyrme effective nucleon–nucleon interaction. In particular, we consider properties of the Wigner distribution function, line of beta stability, Coulomb and symmetry energies, deformation energies, equation of state, nuclear radii, and single-particle level density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

F2 extraction from Inclusive Cross Section Data at Large Bjorken x

While we have made significant progress in probing nuclear structure at low Bjorken x and high Q 2 , we still have gaps in our knowledge regarding the high Bjorken x and intermediate Q 2 kinematic region. This is not an accident, but instead due to the diffculty in assessing the non-perturbative region of nuclear physics. For this reason, the E12-10-002 experiment measured the H(e,e') and D(e,e') inclusive cross section in the resonance region. This will allow us to study both perturbative and non-perturbative physics in a kinematic region that lacks precision measurements. Measurements were made up to a Bjorken x of 0.99 and Q 2 up to 17 GeV 2 . These measurements are made using the independent HMS and SHMS spectrometers at the upgraded beam energy of 11 GeV in Hall C at Jefferson Lab. Further, we extract F 2 structure functions from the measured data and study the effect of their inclusion in a global Parton Distribution Function analysis.

Araiza Gonzalez, Fernando↗
DOE OSTI2025

Global QCD analysis of spin PDFs in the proton with high-𝑥 and lattice constraints

We perform a comprehensive global QCD analysis of spin-dependent parton distribution functions (PDFs), combining all available data on inclusive and semi-inclusive deep-inelastic scattering (DIS), as well as inclusive weak boson and jet production in polarized 𝑝⁢𝑝 collisions, simultaneously extracting spin-averaged PDFs and fragmentation functions. Including recent Jefferson Lab DIS data at high 𝑥, together with subleading power corrections to the leading-twist framework, allows us to verify the stability of the PDFs for 𝑊 2 ≥ 4 GeV 2 and quantify the uncertainties on the spin structure functions more reliably. We explore the use of new lattice QCD data on gluonic pseudo-Ioffe time distributions, which, together with jet production and high-𝑥 DIS data, improve the constraints on the polarized gluon PDF. The expanded kinematic reach afforded by the data into the high-𝑥 region allows us to refine the bounds on higher-twist contributions to the spin structure functions, and test the validity of the Bjorken sum rule.

Perturbative QCD↗
DOE OSTI2020

Particlization in fluid dynamical simulations of heavy-ion collisions: The IS3D module

The IS3D particlization module simulates the emission of hadrons from heavy-ion collisions via Monte-Carlo sampling of the Cooper–Frye formula which converts fluid dynamical information into local phase-space distributions for hadrons. The code package includes multiple choices for the non-equilibrium correction to these distribution functions: the 14-moment approximation, first-order Chapman–Enskog expansion, and two types of modified equilibrium distributions. This makes it possible to explore to what extent heavy-ion experimental data are sensitive to different choices for δf n , presently the main source of theoretical uncertainty in the particlization stage. Here, we validate our particle sampler with a high degree of precision by generating several million hadron emission events from a longitudinally boost-invariant hypersurface and comparing the event-averaged particle spectra and space–time distributions to the Cooper–Frye formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov– Poisson–Lenard–Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as $f$ = $\mathscr{E}$[$ρ_f$] + g , where $\mathscr{E}$ is a local equilibrium distribution, depending on the macroscopic moments $ρ_f$ = $∫_{\mathbb{R}}$ $efdv$ = $\langle$$ef$$\rangle$$_{\mathbb{R}}$, where e = (1, $v$, $\frac{1}{2}$ $v^2$) T , and $g$, the microscopic distribution, is defined such that $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $ρ_f$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. Here, we present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗