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 721 records · Page 40

Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$

Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Plasma diagnostics and modeling of lithium-containing plasmas

Thin-film deposition from chemically reactive multi-component plasmas is complex, and the lack of electron collision cross-sections for even the most common metalorganic precursors and their fragments complicates their modeling based on fundamental plasma physics. This study focuses on understanding the plasma physics and chemistry in argon (Ar) plasmas containing lithium bis (trimethylsilyl) amide used to deposit Li x Si y thin films. These films are emerging as potential solid electrolytes for lithium-ion batteries, and the Li-to-Si ratio is a crucial parameter to enhance their ionic conductivity. We deposited Li x Si y films in an axial flow-through plasma reactor and studied the factors that determine the variation of the Li-to-Si ratio in films deposited at various points on a substrate spanning the entire reactor axis. While the Li-to-Si ratio is 1:2 in the precursor, the Li-to-Si ratio is as high as 3:1 in films deposited near the plasma entrance and decreases to 1:1 for films deposited downstream. Optical emission from the plasma is dominated by Li emission near the entrance, but Li emission disappears downstream, which we attribute to the complete consumption of the precursor. We hypothesized that the axially decreasing precursor concentration affects the electron energy distribution function in a way that causes different dissociation efficiencies for the production of Li and Si. We used Li line intensities to estimate the local precursor concentration and Ar line ratios to estimate the local reduced electric field to test this hypothesis. This analysis suggests that the mean electron energy increases along the reactor axis with decreasing precursor concentration. As a result, the decreasing Li-to-Si ratio with axially decreasing precursor concentration may be explained by Li release from the precursor having lower threshold energy than Si release.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

bnl/pub-Maffettone_2020_08

The application of the XCA package as first demonstrated in aXiv:2008.00283. ABSTRACT: The discovery of new structural and functional materials is driven by phase identification, often using X-ray diffraction (XRD). Automation has accelerated the rate of XRD measurements, greatly outpacing XRD analysis techniques that remain manual, time consuming, error prone, and impossible to scale. With the advent of autonomous robotic scientists or self-driving labs, contemporary techniques prohibit the integration of XRD. Here, we describe a computer program for the autonomous characterization of XRD data, driven by artificial intelligence (AI), for the discovery of new materials. Starting from structural databases, we train an ensemble model using a physically accurate synthetic dataset, which output probabilistic classifications --- rather than absolutes --- to overcome the overconfidence in traditional neural networks. This AI agent behaves as a companion to the researcher, improving accuracy and offering unprecedented time savings, and is demonstrated on a diverse set of organic and inorganic materials challenges. This innovation is directly applicable to inverse design approaches, robotic discovery systems, and can be immediately considered for other forms of characterization such as spectroscopy and the pair distribution function.

Maffettone, PhillipM [Brookhaven National Lab. (BN↗

Crystallography companion agent for high-throughput materials discovery

The discovery of new structural and functional materials is driven by phase identification, often using X-ray diffraction (XRD). Automation has accelerated the rate of XRD measurements, greatly outpacing XRD analysis techniques that remain manual, time-consuming, error-prone and impossible to scale. With the advent of autonomous robotic scientists or self-driving laboratories, contemporary techniques prohibit the integration of XRD. Here, we describe a computer program for the autonomous characterization of XRD data, driven by artificial intelligence (AI), for the discovery of new materials. Starting from structural databases, we train an ensemble model using a physically accurate synthetic dataset, which outputs probabilistic classifications—rather than absolutes—to overcome the overconfidence in traditional neural networks. This AI agent behaves as a companion to the researcher, improving accuracy and offering substantial time savings. It is demonstrated on a diverse set of organic and inorganic materials characterization challenges. This method is directly applicable to inverse design approaches and robotic discovery systems, and can be immediately considered for other forms of characterization such as spectroscopy and the pair distribution function.

36 MATERIALS SCIENCE↗

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Presentation of gas-phase-reactant-accessible single-rhodium-atom catalysts for CO oxidation, via MOF confinement of an Anderson polyoxometalate

Geometric or electronic confinement of guests within nanoporous hosts holds promise for imparting catalytic functionality, including single-metal-atom catalytic functionality, to existing materials. When the nanoporous host is a metal–organic framework (MOF), single-metal-ion catalysts have typically been installed by grafting to an open site on an inorganic node, with the node effectively becoming the support for the catalyst. This approach, however, imposes compositional constraints, as the node not only needs to be receptive to grafting, but also must be capable of stabilizing the framework against solvent evacuation, chemical exposure, and heating. Here, we show that disk-like, Anderson polyoxometalate clusters (RhMo 6 O 24 n- and Mo 7 O 24 m- ; POMs) can be confined in pore-specific and orientation-specific fashion within the hierarchically porous, Zr(IV)-based MOF, NU1K. Self-limiting loading of one cluster per pore, and associated nano-confinement, serve to isolate each POM and prevent consolidation caused by sintering. Additionally, the oriented confinement serves to expose individual rhodium atoms to candidate gas-phase reactants, while enabling the rhodium atom to employ a well-defined oxy-molybdenum cluster, rather than a MOF node, as a support. Synchrotron-based difference-electron-density maps and differential pair-distribution-function analyses of scattered X-rays establish cluster siting and orientation and confirm isolation. Nanoconfined (i.e., MOF- and POM-confined) single-rhodium(III)-atoms are catalytically competent for an illustrative gas-phase reaction, CO oxidation by O 2 , with the MOF-isolated POM enormously outperforming nonporous, MOF-free, solid (NH 4 ) 3 [H 6 RhMo 6 O 24 ]·6H 2 O. This paper highlights the value of MOF-based nano-confinement and oriented isolation of planar POMs as a means of uniformly presenting and stabilizing potent single-metal-atom catalysts, in reactant-accessible form, on well-defined supports.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental Verification of Ir 5d Orbital States and Atomic Structures in Highly Active Amorphous Iridium Oxide Catalysts

In iridium oxide catalysts, the electronic states whose energies are in the range of energetics and charge transfer kinetics of the oxygen evolution reaction (OER) originate from the Ir 5d orbital states. However, the understanding of the atomic structures and orbital states underlying catalytic reactivity in amorphous iridium oxide oxygen evolving catalysts (Ir-OECs) is incomplete compared to that of crystalline oxides, owing to a lack of direct experimental verification. Here, we present experimental approaches using resonant inelastic X-ray scattering (RIXS) to directly access Ir 5d orbital excitations at the Ir L 3 edge and atomic pair distribution function (PDF) measurements to characterize electronic and coordination structures at the atomic scale. The so-called iridium blue layer (IrBL) and IrO x were formed from the organometallic precursor complex [Cp*Ir(H 2 O) 3 ]SO 4 and the inorganic precursor IrCl 3 , respectively. Ex situ IrBL and IrO x films for RIXS and PDF measurements were prepared by conditioning electrodeposited films at a low voltage. The incident energy RIXS profile of IrO x exhibited extra weak resonantly enhanced excitation below 2 eV energy loss. The feature was clearly different from a single high-energy excitation above 3 eV of IrBL related to the interband transition between π- and σ-antibonding states. The atomic structure refinement based on PDF measurements revealed the atomic structure domains to have edge- and corner-shared IrO 6 octahedra with trigonal-type distortion. Furthermore, density functional theory calculations guided by the refined atomic structures shed light on the electronic structure corresponding to experimental results, including insulating and metallic phases in ex situ IrBL and IrO x films, respectively. Our study establishes different Ir 5d orbital states and atomic structures in two amorphous Ir oxide OER catalysts in their reduction states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gluon helicity in the nucleon from lattice QCD and machine learning

We present the first lattice QCD determination of the light cone gluon helicity correlation parton distribution function (PDF) with numerical evidence toward disfavoring negative gluon polarization in the nucleon. We present a solution for eliminating an inevitable contamination term that dominates the Euclidean correlations and makes determining gluon helicity PDF unfeasible. The proposed synergy between lattice QCD and artificial intelligence offers a superior platform to alleviate the defining challenge of extracting quark and gluon PDFs from the lattice data that are available in a limited domain due to a finite range of accessible hadron momenta. We suggest a systematically improvable method to extract PDFs from the lattice data, independent of inadequate parametrizations. The result of the gluon helicity will improve our understanding of the role of spin in the strong interaction and the nucleon-spin structure.

Lattice QCD↗

Gluon transverse-momentum-dependent distributions from large-momentum effective theory

We demonstrate that gluon transverse-momentum-dependent parton distribution functions (TMDPDFs) can be extracted from lattice calculations of appropriate Euclidean correlations in large-momentum effective theory (LaMET). Based on perturbative calculations of gluon unpolarized and helicity TMDPDFs, we present a matching formula connecting them and their LaMET counterparts, where the latter are renormalized in a scheme facilitating lattice calculations and converted to the $\overline{MS}$ scheme. The hard matching kernel is given up to one-loop level. We also show that the perturbative result is independent of the prescription used for the pinch-pole singularity in the relevant correlations. Our results offer a guidance for the extraction of gluon TMDPDFs from lattice simulations, and have the potential to greatly facilitate perturbative calculations of the hard matching kernel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fresh look at the nuclear transparency using the generalized parton distributions

Color transparency (CT) is a fundamental phenomenon in QCD in which hadrons produced in high-energy exclusive processes traverse nuclear matter with minimal interactions. Nuclear transparency, which quantifies this attenuation suppression, is a quantity with high sensitivity to CT effects and provides critical insights into QCD dynamics in nuclear environments. In this study, we revisit nuclear transparency using the framework of generalized parton distributions (GPDs). By constructing nuclear GPDs (nGPDs) through the incorporation of nuclear parton distribution functions, we calculate the nuclear transparency 𝑇⁡(𝑄 2 ) for the carbon nucleus as a function of momentum transfer 𝑄 2 considering various definitions and compare the results obtained with available experimental data. Our finding highlights the importance of choosing a physically motivated definition of nuclear transparency. Moreover, we emphasize that a more reliable determination of nGPDs requires a dedicated global analysis incorporating nuclear data. Such an approach is essential for improving the theoretical understanding of CT and for achieving consistency with experimental observations in the high-𝑄 2 regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of threshold resummation for large-x PDF in large momentum effective theory

Parton distribution functions (PDFs) at large x are challenging to extract from experimental data, yet they are essential for understanding hadron structure and searching for new physics beyond the Standard Model. Within the framework of the large momentum P z expansion of lattice quasi-PDFs, we investigate large x PDFs, where the matching coefficient is factorized into the hard kernel, related to the active quark momentum xP z , and the threshold soft function, associated with the spectator momentum (1 − x)P z . The renormalization group equation of the soft function enables the resummation of the threshold double logarithms α k ln 2k (1 − x), which is crucial for a reliable and controllable calculation of large x PDFs. Our analysis with pion valence PDFs indicates that perturbative matching breaks down when the spectator momentum (1 − x)P z approaches Λ QCD , but remains valid when both xP z and (1 − x)P z are much larger than Λ QCD . Additionally, we incorporate leading renormalon resummation within the threshold framework, demonstrating good perturbative convergence in the region where both spectator and active quark momenta are perturbative scales.

factorization↗

Machine learning-based jet and event classification at the Electron-Ion Collider with applications to hadron structure and spin physics

We explore machine learning-based jet and event identification at the future Electron-Ion Collider (EIC). We study the effectiveness of machine learning-based classifiers at relatively low EIC energies, focusing on (i) identifying the flavor of the jet and (ii) identifying the underlying hard process of the event. We propose applications of our machine learning-based jet identification in the key research areas at the future EIC and current Relativistic Heavy Ion Collider program, including enhancing constraints on (transverse momentum dependent) parton distribution functions, improving experimental access to transverse spin asymmetries, studying photon structure, and quantifying the modification of hadrons and jets in the cold nuclear matter environment in electron-nucleus collisions. We establish first benchmarks and contrast the estimated performance of flavor tagging at the EIC with that at the Large Hadron Collider. We perform studies relevant to aspects of detector design including particle identification, charge information, and minimum transverse momentum capabilities. Additionally, we study the impact of using full event information instead of using only information associated with the identified jet. These methods can be deployed either on suitably accurate Monte Carlo event generators, or, for several applications, directly on experimental data. We provide an outlook for ultimately connecting these machine learning-based methods with first principles calculations in quantum chromodynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gaussian processes for inferring parton distributions

The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.

hadronic spectroscopy↗

Transverse momentum structure of proton within the basis light-front quantization framework

We obtain the leading-twist valence quark transverse-momentum-dependent parton distribution functions (TMD PDFs) for the proton within the basis light-front quantization (BLFQ) framework. Our results are consistent with lattice QCD calculations and our previous results for the collinear limit. We also obtain consistency with the Soffer-type bounds. Within our approach, we find that six T-even TMDs in the leading twist are all independent of each other, and previously found model-dependent relations do not hold. This is a promising sign that our results are representative of future, more extensive treatments of QCD. Furthermore, we obtain a non-trivial x-dependence of the <(p ⊥ ) 2 > and some consistency with the Gaussian ansatz but only in the small <(p ⊥ ) 2 > region. Those features suggest our results may be a useful alternative in future experimental extractions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Explainable AI classification for parton density theory

Quantitatively connecting properties of parton distribution functions (PDFs, or parton densities) to the theoretical assumptions made within the QCD analyses which produce them has been a longstanding problem in HEP phenomenology. To confront this challenge, we introduce an ML-based explainability framework, XAI4PDF, to classify PDFs by parton flavor or underlying theoretical model using ResNet-like neural networks (NNs). By leveraging the differentiable nature of ResNet models, this approach deploys guided backpropagation to dissect relevant features of fitted PDFs, identifying x-dependent signatures of PDFs important to the ML model classifications. By applying our framework, we are able to sort PDFs according to the analysis which produced them while constructing quantitative, human-readable maps locating the x regions most affected by the internal theory assumptions going into each analysis. This technique expands the toolkit available to PDF analysis and adjacent particle phenomenology while pointing to promising generalizations.

Artificial Intelligence↗

Prioritizing ICS Beachhead Systems for Cyber Vulnerability Testing

Cyber Testing for Resilient Industrial Control Systems™ (CyTRICS™) is the Department of Energy’s (DOE’s) program for cybersecurity vulnerability testing, digital subcomponent enumeration, and forensic assessment. CyTRICS leverages best-in-class test facilities and analytic capabilities at six DOE National Laboratories and strategic partnerships with key stakeholders including technology developers, manufacturers, asset owners and operators, and interagency partners. During the program’s development, CyTRICS established a unique methodology for prioritizing digital components within operational technology (OT) and industrial control systems (ICS) in the Energy Sector Industrial Base (ESIB) for cyber vulnerability testing. The CyTRICS Prioritization Process leverages multiple characteristics of systems, components, and their contextual deployment to calculate a quantification of individual digital components for CyTRICS testing. The initial version of the CyTRICS Prioritization Process was premised largely upon the impact which could result to an industrial control system if the digital component under testing was compromised, either through malicious means, faulty engineering, or other modes. The worldwide compromise of the SolarWinds Orion platform, first reported in December 2020, through malicious interference with the digital patching cycle was a watershed event in cyber supply chain security. The SolarWinds compromised demonstrated the strategic importance of certain types of ubiquitous software, and the ability to generate widespread cybersecurity effects. To address this challenge and as a part of the Department of Energy’s response to the SolarWinds compromise, DOE’s Office of Cybersecurity, Energy Security, and Emergency Response (CESER) directed the National Laboratories to evolve the CyTRICS Prioritization Process methodology to encompass additional factors related to the strategic importance of digital components. CESER directed CyTRICS researchers to identify, characterize, and append strategic factors to the CyTRICS Prioritization Process to provide additional weight to these characteristics. National Laboratory expert researchers identified functionality, distribution, and platform characteristics for digital components in ICS and OT that they assessed would be likely targeted in strategic initial-access cyber attack. CyTRICS has termed these factors “ICS Beachhead Systems,” leveraging a definition first advanced by Schneider Electric, which is intended as a blanket term to encompass digital components, products, and systems in OT. This paper describes the ICS Beachhead Systems identified and the rationale for inclusion. As a next step in the research and refinement process, the National Laboratories will validate this initial set of characteristics against digital components evaluated by the CyTRICS program and current implementation of the CyTRICS Prioritization Process. After validation, CyTRICS researchers will then develop a scoring methodology to generate a quantitative score to assess the degree to which a digital component is characterized as an ICS Beachhead System. Finally, the National Laboratories will append this scoring to the existing CyTRICS Prioritization Process algorithm.

97 MATHEMATICS AND COMPUTING↗

The photon content of the neutron

In this work, we complete our CT18qed study with the neutron’s photon parton distribution function (PDF), which is essential for the nucleus scattering phenomenology. Two methods, CT18lux and CT18qed, based on the LUXqed formalism and the DGLAP evolution, respectively, to determine the neutron’s photon PDF have been presented. Various low-Q 2 non-perturbative variations have been carefully examined, which are treated as additional uncertainties on top of those induced by quark and gluon PDFs. The impacts of the momentum sum rule as well as isospin symmetry violation have been explored and turned out to be negligible. A detailed comparison with other neutron’s photon PDF sets has been performed, which shows a great improvement in the precision and a reasonable uncertainty estimation. Finally, two phenomenological implications are demonstrated with photon-initiated processes: neutrino-nucleus W-boson production, which is important for the near-future TeV–PeV neutrino observations, and the axion-like particle production at a high-energy muon beam-dump experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗