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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Light mesons within the basis light-front quantization framework

In this work, we study the light-unflavored mesons as relativistic bound states in the nonperturbative Hamiltonian formalism of the basis light-front quantization approach. The dynamics for the valence quarks of these mesons is specified by an effective Hamiltonian containing the one-gluon exchange interaction and the confining potentials both introduced in our previous work on heavy quarkonia, supplemented additionally by a pseudoscalar contact interaction. We diagonalize this Hamiltonian in our basis function representation to obtain the mass spectrum and the light-front wave functions (LFWFs). Based on these LFWFs, we then study the structure of these mesons by computing the electromagnetic form factors, the decay constants, the parton distribution amplitudes, and the parton distribution functions. Our results are comparable to those from experiments and other theoretical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Measurement of differential t-channel single top (anti)quark production cross-sections at 13 TeV with the ATLAS detector

The production of single top quarks and top antiquarks via the t-channel exchange of a virtual W boson is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider. The full Run 2 data sample recorded with the ATLAS detector in the years 2015–2018 is used, corresponding to an integrated luminosity of 140 fb−1. The absolute and normalised production cross-sections are measured differentially as a function of the transverse momentum and absolute rapidity of the top quark and top antiquark. In addition, the ratio of top quark to top antiquark production cross-sections is measured. The measured distributions are compared with next-to-leading-order quantum chromodynamics predictions obtained with different combinations of matrix-element generators, parton-shower programs and proton parton distribution functions, as well as to next-to-next-to-leading-order calculations. Overall, good agreement is observed between the measurements and the theoretical predictions. For most measured distributions, the sensitivity to differences between the predictions is limited by the systematic uncertainties in the measurement. The measured differential distributions are also interpreted in an effective field theory approach to constrain the Wilson-Coefficient CQq3,1$$ {C}_{Qq}^{3,1} $$ associated with a four-quark operator. The interpretation accounts for the effect of the selection efficiency, which is altered significantly by non-zero contributions from CQq3,1$$ {C}_{Qq}^{3,1} $$.

Aad, G↗
DOE OSTI2020

Mechanistic insight of KBiQ 2 (Q = S, Se) using panoramic synthesis towards synthesis-by-design

Solid-state synthesis has historically focused on reactants and end products; however, knowledge of reaction pathways, intermediate phases and their formation may provide mechanistic insight of solid-state reactions. With an increased understanding of reaction progressions, design principles can be deduced, affording more predictive power in materials synthesis. In pursuit of this goal, in situ powder X-ray diffraction is employed to observe crystalline phase evolution over the course of the reaction, thereby constructing a “panoramic” view of the reaction from beginning to end. We conducted in situ diffraction studies in the K–Bi–Q (Q = S, Se) system to understand the formation of phases occurring in this system in the course of their reactions. Powder mixtures of K 2 Q to Bi 2 Q 3 in 1:1 and 1.5:1 ratios were heated to 800 °C or 650 °C, while simultaneously collecting diffraction data. Three new phases, K 3 BiS 3 , β-KBiS 2 , and β-KBiSe 2 , were discovered. Panoramic synthesis showed that K 3 BiQ 3 serves an important mechanistic role as a structural intermediate in both chalcogen systems (Q = S, Se) in the path to form the KBiQ 2 structure. Thermal analysis and calculations at the density functional theory (DFT) level show that the cation-ordered β-KBiQ 2 polymorphs are the thermodynamically stable phase in this compositional space, while Pair Distribution Function (PDF) analysis shows that all α-KBiQ 2 structures have local disorder due to stereochemically active lone pair expression of the bismuth atoms. The formation of the β-KBiQ 2 structures, both of which crystallize in the α-NaFeO 2 structure type, show a boundary where the structure can be disordered or ordered with regards to the alkali metal and bismuth. A cation radius tolerance for six-coordinate cation site sharing of ($\frac{r^+}{r^{3+}}$) ~ 1.3 is proposed. The mechanistic insight the panoramic synthesis technique provides in the K–Bi–Q system is progress towards the overarching goal of synthesis-by-design.

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

Significance of radiative corrections on measurements of the EMC effect

Deep inelastic scattering (DIS) from nuclear targets probes the parton distribution functions (PDFs) in nuclei. Comparisons of the PDFs from heavy nuclei and the deuteron show deviations that demonstrate a non-trivial nuclear dependence to these distributions, referred to as the EMC effect. A global analysis of the worlds data on the EMC effect reveals tensions between different extractions. Precise measurements at Jefferson Lab, studying the dependence on both the quark momentum fraction, x, and nuclear mass, show systematic discrepancies among experiments, making the extraction of the A dependence of the EMC effect sensitive to the selection of datasets. Further, by comparing various methods and assumptions used to calculate radiative corrections, we have identified differences that, while not large, significantly impact the EMC ratios and show that using a consistent radiative correction procedure resolves this discrepancy, leading to a more coherent global picture, and allowing for a more robust extraction of the EMC effect for infinite nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Predicting the Galactic population of free-floating planets from realistic initial conditions

ABSTRACT We present the first prediction for the mass distribution function of Galactic free-floating planets (FFPs) that aims to accurately include the relative contributions of multiple formation pathways and stellar populations. We derive our predicted distribution from dedicated simulations of planet birth, growth, migration, and ejection around circumbinary systems and extend these results to also include the contributions from single and wide binary systems. Our resulting FFP mass distribution shows several distinct features, including a strong peak at ${\sim} 8\, {\rm M}_{\oplus }$ arising from the transition between pebble and gas accretion regimes and a trough at ${\sim} 1\, {\rm M}_{\oplus }$ due to the shift in the dominant ejection process from planet–planet scattering to ejection through interactions with stars in circumbinary systems. We find that interactions with the central binary in close circumbinary systems are likely the dominant progenitor for FFPs more massive than Earth, leading to a steep power-law dependence in mass that agrees well with existing observations. In contrast, we find planet–planet scattering events in single and wide binary systems likely produce the majority of planets at Mars mass and below, resulting in a shallower power-law dependence. Our results suggest that existing extrapolations into the sub-terrestrial mass range may significantly overestimate the true FFP abundance. The features we predict in the mass distribution of FFPs will be detectable by upcoming space-based microlensing surveys and, if observed, will provide key insight into the origins of FFPs and the environments in which they form.

Coleman, Gavin A. L. (ORCID:0000000151118963)↗
DOE OSTI2021

Deep-inelastic scattering with positron beams

Here we discuss several advantages of using positron beams in addition to electron beams for the study of parton distribution functions (PDFs) of the nucleon in deep-inelastic scattering. Specific topics include the flavor separation of the unpolarized PDFs as well as the spin-dependent distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗
DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Theory and applications of parton pseudodistributions

We review the basic theory of the parton pseudodistributions approach and its applications to lattice extractions of parton distribution functions. The crucial idea of the approach is the realization that the correlator $M(z,p)$ of the parton fields is a function $M(ν,-z^2)$ of Lorentz invariants $ν=-(zp)$, the Ioffe time, and the invariant interval $z^2$. This observation allows to extract the Ioffe-time distribution $M(ν,-z^2)$ from Euclidean separations z accessible on the lattice. Another basic feature is the use of the ratio $M(ν,-z^2)≡M(ν,-z^2)/M(0,-z^2)$, that allows to eliminate artificial ultraviolet divergence generated by the gauge link for spacelike intervals. The remaining $z^2$-dependence of the reduced Ioffe-time distribution $M(ν,-z^2)$ corresponds to perturbative evolution, and can be converted into the scale-dependence of parton distributions $f(x,μ^2)$ using matching relations. The ν-dependence of $M(ν,-z^2)$ governs the x-dependence of parton densities $f(x,μ^2)$. The perturbative evolution was successfully observed in exploratory quenched lattice calculation. The analysis of its precise data provides a framework for extraction of parton densities using the pseudodistributions approach. It was used in the recently performed calculations of the nucleon and pion valence quark distributions. We also discuss matching conditions for the pion distribution amplitude and generalized parton distributions, the lattice studies of which are now in progress.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Magnetic Energy Dissipation and γ-Ray Emission in Energetic Pulsars

Abstract Some of the most energetic pulsars exhibit rotation-modulated γ -ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power ( γ -ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower γ -ray efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Neutron Structure Functions at Large Bjorken X

The E12-10-002 (F2) experiment ran in Hall C at JLab to extract the F2 structure functions of proton and neutron from inclusive H(e, e0) and D(e, e0) reactions in the high Bjorken x region. Those extractions put constraints on the Parton Distribution Functions PDFs (especially at large Bjorken x), and facilitate the study of Quark Hadron Duality. The F2 experiment measurements covers a large kinematic range in x from 0.2 to 1.0, and in Q2 from 4 to 16 GeV2 . While the proton F2 structure function has been studied extensively through inelastic electron-proton scattering, much less is known about neutron structure due to the unavailability of high density, free neutron targets. The BONuS12 experiment was proposed to measure the neutron F2 on a nearly free neutron within a weakly bound deuteron target via the spectator tagging method. Tagging the slow backward moving spectator protons minimizes both off-shell and final-state interaction effects, and the measured proton momentum is used to correct for the initial-state momentum of the struck neutron. The recoil detector is used to detect spectator protons with momenta 70 < ps < 150 MeV/c. This dissertation outlines the design, construction and testing of the recoil detector, Radial Time Projection Chamber (RTPC) and discusses the results of the H(e, e0) and D(e, e0) cross sections and the F2 structure functions obtained from the Hall C F2 data analysis.

Nadeeshani, Sooriyaarachchilage↗
DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Exploring the non-Gaussianity of the cosmic infrared background and its weak gravitational lensing

ABSTRACT Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the cosmic infrared background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussianity of the CIB, revealing additional detail on top of its well-measured power spectrum. To achieve this, we use needlet-like multipole-band filters to calculate the variance and higher-point correlations. Using our simulations, we show the two-, three- and four-point spectra, and compare our calculated power spectra and bispectra to Planck values. We then lens the CIB, shell-by-shell with corresponding convergence maps, to capture the broad redshift extent of both the CIB and its lensing convergence. The lensing of the CIB changes the three- and four-point functions by a few tens of per cent at large scales, unlike with the power spectrum, which changes by less than two per cent. We expand our analyses to encompass the full intensity probability distribution functions (PDFs) involving all n-point correlations as a function of scale. In particular, we use the relative entropy between lensed and unlensed PDFs to create a spectrum of templates that can allow estimation of lensing. The underlying CIB model is missing the important role of star bursting, which we test by adding a stochastic lognormal term to the intensity distributions. The novel aspects of our filtering and lensing pipeline should prove useful for any radiant background, including line intensity maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Forward production of prompt neutrinos from charm in the atmosphere and at high energy colliders

The high-energy atmospheric neutrino flux is dominated by neutrinos from the decays of charmed hadrons produced in the forward direction by cosmic ray interactions with air nuclei. We evaluate the charm contributions to the prompt atmospheric neutrino flux as a function of the center-of-mass energy √s of the hadronic collision and of the center-of-mass rapidity y of the produced charm hadron. Uncertainties associated with parton distribution functions are also evaluated as a function of y. We find that the y coverage of LHCb for forward heavy-flavour production, complemented by the angular coverage of present and future forward neutrino experiments at the LHC, bracket the most interesting y regions for the prompt atmospheric neutrino flux. At √s = 14 TeV foreseen for the HL-LHC phase, nucleon collisions in air contribute to the prompt neutrino flux prominently below E ν ~ 10 7 GeV. Measurements of forward charm and/or forward neutrinos produced in hadron collisions up to √s = 100 TeV, which might become possible at the FCC, are relevant for the prompt atmospheric neutrino flux up to E ν = 10 8 GeV and beyond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Machine Learning Approach to Modeling of Neutral Particles Transport in Plasma

A propagator‐based approach is investigated for Monte‐Carlo (MC) modeling of neutral particle transport in fusion boundary plasmas. The propagator is based on a Green's function for the neutral kinetic equation, which depends on the plasma profiles. A neural network (NN)‐based model for the propagator provides a fast and accurate solution for the neutral distribution function in plasma. Preliminary results from a small 1D test problem look encouraging. The proposed approach, a propagator‐based NN model for neutral transport in plasma, has potential for generalization to higher dimensions and efficient coupling with plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurement of the top quark pole mass using $ \textrm{t}\overline{\textrm{t}} $+jet events in the dilepton final state in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A measurement of the top quark pole mass $\mathcal{m}^\text{pole}_\text{t}$ in events where a top quark-antiquark pair ($\text{t}\bar{\text{t}}$) is produced in association with at least one additional jet ($\text{t}\bar{\text{t}}$ +jet) is presented. This analysis is performed using proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the CMS experiment at the CERN LHC, corresponding to a total integrated luminosity of 36.3 fb -1 . Events with two opposite-sign leptons in the final state (e + e – , μ + μ – , e ± μ ∓ ) are analyzed. The reconstruction of the main observable and the event classification are optimized using multivariate analysis techniques based on machine learning. The production cross section is measured as a function of the inverse of the invariant mass of the $\text{t}\bar{\text{t}}$ +jet system at the parton level using a maximum likelihood unfolding. Given a reference parton distribution function (PDF), the top quark pole mass is extracted using the theoretical predictions at next-to-leading order. For the ABMP16NLO PDF, this results in $\mathcal{m}^\text{pole}_\text{t}$ = 172.93 ± 1.36 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Z-boson production in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 8.16 TeV and Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5.02 TeV

Measurement of Z-boson production in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 8 . 16 TeV and Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5 . 02 TeV is reported. It is performed in the dimuon decay channel, through the detection of muons with pseudorapidity - 4 < η μ < - 2 . 5 and transverse momentum $ {p}_{\mathrm{T}}^{\mu } $ > 20 GeV/ c in the laboratory frame. The invariant yield and nuclear modification factor are measured for opposite-sign dimuons with invariant mass 60 < m μμ < 120 GeV/ c 2 and rapidity 2 . 5 < $ {y}_{\mathrm{cms}}^{\mu \mu} $ < 4. They are presented as a function of rapidity and, for the Pb-Pb collisions, of centrality as well. The results are compared with theoretical calculations, both with and without nuclear modifications to the Parton Distribution Functions (PDFs). In p-Pb collisions the center-of-mass frame is boosted with respect to the laboratory frame, and the measurements cover the backward ( - 4 . 46 < $ {y}_{\mathrm{cms}}^{\mu \mu} $ < - 2 . 96) and forward (2 . 03 < $ {y}_{\mathrm{cms}}^{\mu \mu} $ < 3 . 53) rapidity regions. For the p-Pb collisions, the results are consistent within experimental and theoretical uncertainties with calculations that include both free-nucleon and nuclear-modified PDFs. For the Pb-Pb collisions, a 3 . 4 σ deviation is seen in the integrated yield between the data and calculations based on the free-nucleon PDFs, while good agreement is found once nuclear modifications are considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Determination of proton PDF uncertainties with Markov chain Monte Carlo

We present an analysis of parton distribution functions (PDFs) of the proton using Markov chain Monte Carlo (MCMC) methods. The MCMC approach naturally implements Bayes’ theorem and, thus, provides a means to directly sample the underlying probability distribution—in this case, the probability distribution of the PDF parameters. This allows for a straightforward propagation of the resulting uncertainties into any PDF-dependent observable, preserving their simple probabilistic interpretation. In our analysis we include a broad set of deep inelastic scattering data from HERA, BCDMS and NMC experiments along with the Drell-Yan, 𝑊 and 𝑍 boson data from LHC and Tevatron experiments, which combined with theoretical calculations at next-to-next-to-leading order in QCD allow for realistic determination of PDFs. The main focus of this analysis is to explore alternative methods for PDF uncertainty estimation that are more firmly grounded in statistical principles. We show that the flexibility of the Bayes framework, allowing one, e.g., to account for non-Gaussianity or inconsistencies of datasets, is crucial to extract realistic uncertainties when such assumptions are not fulfilled. We also demonstrate that MCMC allows one to determine the Δ⁢𝜒 2 value corresponding to a given confidence level in the sample, which can, in turn, be used as a statistically well-founded tolerance criterion used in the Hessian method, thus addressing one of its main long-standing drawbacks.

Risse, Peter Clemens [Universität Münster (Germany↗