SEARCH · Engineering Papers

Results for “PDF”

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

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Illustrated formalisms for total scattering data: a guide for new practitioners

The total scattering method is the simultaneous study of both the real- and reciprocal-space representations of diffraction data. While conventional Bragg-scattering analysis (employing methods such as Rietveld refinement) provides insight into the average structure of the material, pair distribution function (PDF) analysis allows for a more focused study of the local atomic arrangement of a material. Generically speaking, a PDF is generated by Fourier transforming the total measured reciprocal-space diffraction data (Bragg and diffuse) into a real-space representation. However, the details of the transformation employed and, by consequence, the resultant appearance and weighting of the real-space representation of the system can vary between different research communities. As the worldwide total scattering community continues to grow, these subtle differences in nomenclature and data representation have led to conflicting and confusing descriptions of how the PDF is defined and calculated. This paper provides a consistent derivation of many of these different forms of the PDF and the transformations required to bridge between them. Some general considerations and advice for total scattering practitioners in selecting and defining the appropriate choice of PDF in their own research are presented. This contribution aims to benefit people starting in the field or trying to compare their results with those of other researchers.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Probability Density Function Control of Frequency Fluctuations in Renewable-Rich Power Systems

The stochastic nature of renewable energy sources (RESs) necessitates treating power system frequency response as a random process with a nonstationary probability density function (PDF). Here, based upon the stochastic distribution control theory originated by the second author, this paper proposes a novel stochastic controller to improve the frequency PDF in power grids when integrating a large amount of RESs, thereby minimizing the effects of uncertainties and enhancing overall system stability. The key idea is to manipulate the controllable power generation resources so that the frequency PDF is make to follow a target PDF by using the stochastic distribution control theory originated by the second author. The proposed method can easily be plugged into existing automatic generation controls for multi-area transmission grids. The proposed method is validated via a modified Kundar's two area system and 240-bus Western Electricity Coordinating Council systems. The simulation results show that the proposed control shapes the frequency PDF narrower and sharper, leading to a notable improvement toward minimizing the effects of randomness and uncertainty during grid operation.

frequency↗

DOE OSTI2026

Cation Data for the East River Watershed, Colorado (2014-2025)

This data package contains mean values for cation concentration for water samples taken from the East River Watershed in Colorado. Inductively coupled plasma mass spectrometry (ICP-MS) has been used to measure the concentrations of elements of interest simultaneously for the East River Watershed, Colorado groundwater and surface water samples to inform insights on the biogeochemistry processes within the watershed. The East River is part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States. For samples collected prior to 06-16-2021, the instrumentation, Elan DRC II, PerkinElmer SCIEX, automatically switches among the three models necessary to analyze all 37 elements. These 37 elements include: (1) Lithium (Li), Beryllium (Be), Boron (B), Sodium (Na), Magnesium (Mg), Aluminium (Al), Silicon (Si), Phosphorus (P), Titanium (Ti), Cobalt (Co), Nickel (Ni), Copper (Cu), Zinc (Zn), Germanium (Ge), Arsenic (As), Rubidium (Rb), Strontium (Sr), Zirconium (Zr), Molybdenum (Mo), Silver (Ag), Cadmium (Cd), Tin (Sn), Antimony (Sb), Caesium (Cs), Barium (Ba), Europium (Eu), Lead (Pb), Thorium (Th), Uranium (U) using standard model, argon Ar as reaction gas, (2) Potassium (K), Calcium (Ca), Vanadium (V), Chromium (Cr), Manganese (Mn), Iron (Fe) using dynamic reaction cell (DRC) model, ammonia NH3 as reaction gas, and (3) Phosphorus (P) and Selenium (Se) using DRC model, oxygen O2 as reaction gas. Note for the samples with higher concentrations of chloride (Cl-), asenic (As) concentrations were analysed with DRC model (oxygen O2 as reaction gas) to avoid the interference of chloride. For samples collected on and after 06-16-2021, an advanced Agilent 8900 triple quadrupole inductively coupled plasma mass spectrometry system (Agilent 8900 QQQ ICP-MS, Agilent Technologies) has been used to measure the concentrations of interested 36 elements simultaneously for environmental samples, including (1) Lithium (Li), Beryllium (Be) and Boron (B) using standard no gas mode, (2) Sodium (Na), Magnesium (Mg), Aluminium (Al) Phosphorus (P), Potassium (K), Chromium (Cr), Manganese (Mn), Iron (Fe), Cobalt (Co), Nickel (Ni), Copper (Cu), Zinc (Zn), Germanium (Ge), Arsenic (As), Rubidium (Rb), Strontium (Sr), Zirconium (Zr), Molybdenum (Mo), Silver (Ag), Cadmium (Cd), Tin (Sn), Antimony (Sb), Cesium (Cs), Barium (Ba), Europium (Eu), Lead (Pb), Thorium (Th) and Uranium (U) using standard helium (He) collision mode, (3) Titanium (Ti) and Vanadium (V) using high Energy (HEHe) helium (He) collision mode, and (4) Silicon (Si), Calcium (Ca) and Selenium (Se) using standard H2 reaction mode. All samples were prepared/diluted with 2% (v/v) ultrapure nitric acid in Milli-Q water (18.2 mega ohm-cm), and analyzed under a rigorous quality assurance and quality control (QA/QC) process. This data package contains (1) a zip file (cation_data_2014_2025.zip) containing a total of 5,849 files: 5.848 data files of cation data from across the Lawrence Berkeley National Laboratory (LBNL) Watershed Function Scientific Focus Area (SFA) which is reported in .csv files per location and a locations.csv (1 file) with latitude and longitude for each location; (2) a file-level metadata (v6_20260901_flmd.csv) file that lists each file contained in the dataset with associated metadata; (3) a data dictionary (v6_20260901_dd.csv) file that contains terms/column_headers used throughout the files along with a definition, units, and data type; (4) PDF and docx files for the detemination of Method Detection Limits (MDLs) for ICP-MS PerkinElmer DRC II instrumentation (Detemination_of_Method_Detection_Limits__MDLs__for_ICP_MS__PerkinElmer_Elan_DRC_II__LBL_Bldg74_Lab214D) for samples before November 2021; (5) PDF and docx files for the determination of MDLs for ICP-MS Agilent 8900 QQQ instrumentation (ICP_MS_Analysis_detection_limits_and_QA_QC_WenmingDong_updated_2026-08-06) for samples November 2021 and onward. Missing values within the anion data files are noted as either "-9999" or "0.0" for not detectable (N.D.) data. There are a total of 113 locations containing cation data. Update on 2021-04-11: Added Detemination of Method Detection Limits (MDLs) for ICP-MS document, which can be accessed as a PDF or with Microsoft Word. Update on 2022-06-10: versioned updates to this dataset was made along with these changes: (1) updated cation data for all locations up to 2021-12-31, (2) removal of units from column headers in datafiles, (3) added row underneath headers to contain units of variables, (4) removed suffix and prefix on two variables (“aqberylliumion_asberyllium” and “aqlithiumion_aslithium”), (5) added -9999 for empty numerical cells, and (6) the addition of the file-level metadata (flmd.csv) and data dictionary (dd.csv) were added to comply with the File-Level Metadata Reporting Format. Update on 2022-09-09: Updates were made to reporting format specific files (file-level metadata and data dictionary) to correct swapped file names, add additional details on metadata descriptions on both files, add a header_row column to enable parsing, and add version number and date to file names (v2_20220909_flmd.csv and v2_20220909_dd.csv). Update on 2023-08-08: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-01-05. The file level metadata and data dictionary files were updated to reflect the additional data added. Update on 2024-03-11: Updates were made to both the data files and reporting format specific files. New available anion data was added, up until 2023-10-16. Further, revisions to the data files were made to remove incorrect data points (from 1970 and 2001). The reporting format specific files were updated to reflect the additional data added. Updated versions of the PDF and docx files for determination of MDLs for ICP-MS data were added to this dataset for samples starting in November 2021. Update on 2025-05-15: Updates were made to both the data files and reporting format specific files. New available cation data was added, up until the end of WY2024 (September 30, 2024). International Generic Sample Numbers (IGSNs), when registered, were added to the data files. The reporting format specific files were updated to reflect the additional data added. Update on 2026-09-01: Updates were made to both the data files and reporting format specific files. New available cation data was added, up until the end of WY2025 (September 30, 2025). Updated versions, as of 2026-08-06, of the PDF and docx files for determination of MDLs for ICP-MS data were added to this dataset for samples starting in November 2021.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

On the Statistics of Elsasser Increments in Solar Wind and Magnetohydrodynamic Turbulence

In this Letter we investigate the dependency with scale of the empirical probability distribution functions (PDF) of Elsasser increments using large sets of WIND data (collected between 1995 and 2017) near 1 au. The empirical PDF are compared to the ones obtained from high-resolution numerical simulations of steadily driven, homogeneous reduced MHD turbulence on a 2048 3 rectangular mesh. A large statistical sample of Alfvénic increments is obtained by using conditional analysis based on the solar wind average properties. The PDF tails obtained from observations and numerical simulations are found to have exponential behavior in the inertial range, with an exponential decrement that satisfies power laws of the form α l ∝ l –μ , where l is the scale size, with μ between 0.17 and 0.25 for observations and 0.43 for simulations. PDF tails were extrapolated assuming their exponential behavior extends to arbitrarily large increments in order to determine structure function scaling laws at very high orders. Our results point to potentially universal scaling laws governing the PDF of Elsasser increments and to an alternative approach to investigate high-order statistics in solar wind observations.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

A comprehensive model to determine the effects of temperature and species fluctuations on reactions in turbulent reacting flows

The joint probability distribution function (pdf), which is a modification of the bivariate Gaussian pdf, is discussed and results are presented for a global reaction model using the joint pdf. An alternative joint pdf is discussed. A criterion which permits the selection of temperature pdf's in different regions of turbulent, reacting flow fields is developed. Two principal approaches to the determination of reaction rates in computer programs containing detailed chemical kinetics are outlined. These models represent a practical solution to the modeling of species reaction rates in turbulent, reacting flows.

Antaki, P. J.↗

NASA NTRS1985

Effect of reactive intermediates on the propagation of turbulent premixed flames

The most simple version of the Bray and Libby model of turbulent flame propagation has been extended to include a two-step reaction mechanism in which reactants decompose into reactive intermediates which then form products. A transport equation for the joint pdf of intermediate and product concentrations was derived, the pdf parameterized, and transport equations for the parameters obtained. The one-dimensional concentration equations were then solved for the parameters of the pdf and the turbulent kinetic energy. Flame speed and angle are recovered as eigenvalues of the problem. The results show the essential features observed in experiments. The product pdf is composed of peaks at zero and unity product concentration, with a continuous structure representing partially reacted gas. The intermediates pdf has a peak at zero concentration and a broad bell-shaped distribution at nonzero concentration, representing partially reacted gas. The results also show the effect of partitioning of energy release on the flame thickness and speed.

Azzazy, M.↗

NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

NASA NTRSDate not supplied

The Impact of Pixel Size on the Characterization of Deep Convective Clouds for Calibration

The NASA CERES project provides the scientific community the observed TOA SW and LW fluxes for climate monitoring and climate model validation. CERES utilizes hourly geostationary imager derived broadband fluxes, which rely on the channel radiances and associated cloud retrievals, are used to estimate the broadband fluxes between CERES observations. This requires stable and consistent cross-platform imager visible channel calibration. The CERES project utilizes deep convective clouds (DCC) as an invariant Earth target to both monitor the stability of sensors and for radiometric scaling. GSICS, an international collaboration, is also evaluating and implementing the DCC invariant target calibration methodology to provide consistent calibration coefficients across geostationary imagers anchored to the AquaMODIS or the NOAA-20 VIIRS calibration reference. Tropical DCC are the brightest, coldest, most Lambertian, top of the atmosphere Earth targets. The DCC invariant target calibration methodology relies on a large ensemble of tropical D CC-identified pixel-level reflectances, which are aggregated as probability density functions (PDF). By assuming the monthly PDF shape is otherwise consistent in time excepting shifts in reflectance caused by changes in the sensor calibration, the imager stability is monitored. Radiometric scaling is accomplished by ratioing the sensor pair DCC PDF reflectance values. The success of the DCC methodology relies on consistent PDF distributions. The goal of this study is to determine the impact of pixel resolution on the DCC reflectance distribution. Single SNPP-VIIRS 750-m and Landsat 8 OLI 30-m granules are aggregated to degrade the pixel resolution from the native level. The DCC pixels are identified using a BT threshold. Most of the brightest DCC pixels are also the coldest, although there are exceptions. It was found that increasing the BT threshold exponentially increased the number of darker pixels. The pixel resolution did not seem to impact the DCC reflectance PDF distribution for pixel resolutions less than 3 km, which suggests that imagers of varying pixel resolutions may be radiometrically scaled to each other using DCC targets.

DCC↗

DOE OSTI2021

A cloud platform for atomic pair distribution function analysis: PDFitc

A cloud web platform for analysis and interpretation of atomic pair distribution function (PDF) data ( PDFitc ) is described. The platform is able to host applications for PDF analysis to help researchers study the local and nanoscale structure of nanostructured materials. The applications are designed to be powerful and easy to use and can, and will, be extended over time through community adoption and development. The currently available PDF analysis applications, structureMining, spacegroupMining and similarityMapping , are described. In the first and second the user uploads a single PDF and the application returns a list of best-fit candidate structures, and the most likely space group of the underlying structure, respectively. In the third, the user can upload a set of measured or calculated PDFs and the application returns a matrix of Pearson correlations, allowing assessment of the similarity between different data sets. structureMining is presented here as an example to show the easy-to-use workflow on PDFitc . In the future, as well as using the PDFitc applications for data analysis, it is hoped that the community will contribute their own codes and software to the platform.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

A High Energy X-ray Diffraction Study of Amorphous Indomethacin

Amorphous pharmaceuticals often possess a wide range of molecular conformations and bonding arrangements. The x-ray pair distribution function (PDF) method is a powerful technique for the characterization of variations in both intra-molecular and inter-molecular packing arrangements. Here, the x-ray PDF of amorphous Indomethacin is shown to be particularly sensitive to the preferred orientations of the chlorobenzyl ring found in isomers in the crystalline state. In some cases, the chlorobenzyl ring has no preferred torsional angle in the amorphous form, while in others evidence of distinct isomer orientations are observed. Amorphous samples with no preferred torsion angles of the chlorobenzyl ring are found to favor enhanced inter-molecular hydrogen bonding, and this is reflected in the intensity of the first sharp diffraction peak. These significant variations in structure rule out amorphous Indomethacin as a possible standard for x-ray PDF measurements. At high humidity, time resolved PDF's for >40 h reveal water molecules forming hydrogen bonds with Indomethacin molecules. Finally, a simple linear hydrogen bond model indicates that water molecules in the wet amorphous form have similar hydrogen bond strengths to those found between Indomethacin dimers or chains in the dry amorphous form.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

High Pressure–Temperature Study of MgF 2 , CaF 2 , and BaF 2 by Raman Spectroscopy: Phase Transitions and Vibrational Properties of AF 2 Difluorides

The phase transition of AF 2 difluorides strongly depends on pressure, temperature, and cationic radius. Here, we have investigated the phase transition of three difluorides, including MgF 2 , CaF 2 , and BaF 2 , at simultaneously high pressures and temperatures using Raman spectroscopy and X-ray diffraction in externally heated diamond anvil cells up to 55 GPa at 300–700 K. Rutile-type difluoride MgF 2 with a small cationic radius undergoes a transition to the CaCl 2 -type phase at 9.9(1) GPa and 300 K, to the HP-PdF 2 -type phase at 21.0(2) GPa, and to the cotunnite-type phase at 44.2(2) GPa. The phase transition pressure to the HP-PdF 2 and cotunnite structure at 300 K for our single crystal was found to be higher than that in previous studies using polycrystalline samples. Elevating the temperature increases the transition pressure from rutile- to the CaCl 2 -type phase but has a negative influence on the transition pressure when MgF 2 transforms from the HP-PdF 2 - to cotunnite-type phase. Meanwhile, the transition pressure from the CaCl 2 - to HP-PdF 2 -type phase for MgF 2 was identified to be independent of the temperature. Raman peaks suspected to belong to the α-PbO 2 -type phase were observed at 14.6–21.0(1) GPa and 400–700 K. At 300 K, difluorides CaF 2 and BaF 2 in the fluorite structure with larger cationic radii transform to the cotunnite-type phase at 9.6(3) and 3.0(3) GPa at 300 K, respectively, and BaF 2 further undergoes a transition to the Ni 2 In-type phase at 15.5(4) GPa. For both CaF 2 and BaF 2 , elevating the temperature leads to a lower transition pressure from fluorite- to the cotunnite-type phase but has little influence on the transition to the Ni 2 In structure. Raman data provide valuable insights for mode Grüneisen parameters. We note that the mode Grüneisen parameters for both difluorides and dioxides vary linearly with the cation radius. Further calculations for the mode Grüneisen parameters at high pressures for MgF 2 , CaF 2 , and BaF 2 yield a deeper understanding of the thermodynamic properties of the difluorides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Influence of Aerosol Embedded in Shallow Cumulus Cloud Fields on the Surface Solar Irradiance

Ubiquitous shallow cumulus clouds are associated with complex variability in surface solar irradiance (SSI). Aerosol embedded in the cloud field typically has a much smaller overall radiative effect, but can significantly perturb the shape of the SSI probability density function (PDF). These perturbations have important implications for several applications that utilize SSI, but are poorly quantified and are the subject of this study. Multiple cases of shallow cumulus cloud fields with embedded aerosol are simulated at the Southern Great Plains Atmospheric Observatory using large eddy simulation (LES). The LES-derived cloud and aerosol fields are then ingested into Monte Carlo three-dimensional (3D) radiative transfer to simulate SSI. We find a variety of perturbations to the SSI PDF that depend on aerosol presence and optical properties. The processes leading to these perturbations include extinction of the direct beam that often increases from the clear-sky region toward cloud edge due to aerosol hygroscopic growth, and scattering of radiation by aerosol into cloud shadows. The ability to predict the SSI PDF in the presence of aerosol is assessed by adding three representative aerosol optical properties into an existing machine learning framework. We show that machine learning accurately predicts the SSI PDF across a wide range of conditions with negligible computational expense. Importance metrics reveal the relatively high influence of aerosol optical properties in making the predictions. These new findings highlight the important role that aerosol plays in SSI variability for highly 3D cloud-aerosol environments and provides a computationally efficient route forward for its simulation.

3D radiative transfer↗

DOE OSTI2023

Evaluating Shallow Convection Parameterization Assumptions With a q t – w Quadrant Analysis

Abstract Uncertainties associated with the response of shallow clouds to global warming remain a great challenge for climate projection. Due to the small size of these clouds, parameterizations are required to represent them in both current‐ and next‐generation climate models. We present a quantitative evaluation of several important assumptions used in both mass‐flux and Assumed‐PDF Higher‐Order Closure parameterizations of shallow convection. We use large‐eddy simulations (LESs) of four different shallow convection regimes as benchmarks and apply a q t (total moisture)– w (vertical velocity) quadrant analysis technique to identify the “coherent structures” of moist and dry up/down‐drafts. The statistics of these coherent structures and the environment are then used to evaluate commonly used assumptions in mass‐flux parameterizations related to (a) the convective downdrafts and (b) the subplume variability of convective updrafts. For quantitative evaluation of the PDF closure used in AP‐HOCs, we perform offline calculations with the PDF closure used in the Cloud Layers Unified By Bi‐normals and the Simplified Higher‐Order Closure schemes by supplying the closure with statistical quantities directly calculated from the LESs. Then, the impact of parameterization assumptions embedded in the PDF closure in the same two categories is analyzed and understood through comparisons of the output from the closure with the statistics of the moist and dry up/down‐drafts from the LESs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Effect of solvothermal synthesis parameters on the crystallite size and atomic structure of cobalt iron oxide nanoparticles

We here investigate how the synthesis method affects the crystallite size and atomic structure of cobalt iron oxide nanoparticles. By using a simple solvothermal method, we first synthesized cobalt ferrite nanoparticles of ca. 2 and 7 nm, characterized by Transmission Electron Microscopy (TEM), Small Angle X-ray scattering (SAXS), X-ray and neutron total scattering. The smallest particle size corresponds to only a few spinel unit cells. Nevertheless, Pair Distribution Function (PDF) analysis of X-ray and neutron total scattering data shows that the atomic structure, even in the smallest nanoparticles, is well described by the spinel structure, although with significant disorder and a contraction of the unit cell parameter. These effects can be explained by the surface oxidation of the small nanoparticles, which is confirmed by X-ray near edge absorption spectroscopy (XANES). Neutron total scattering data and PDF analysis reveal a higher degree of inversion in the spinel structure of the smallest nanoparticles. Neutron total scattering data also allow magnetic PDF (mPDF) analysis, which shows that the ferrimagnetic domains correspond to ca. 80% of the crystallite size in the larger particles. A similar but less well-defined magnetic ordering was observed for the smallest nanoparticles. Finally, we used a co-precipitation synthesis method at room temperature to synthesize ferrite nanoparticles similar in size to the smallest crystallites synthesized by the solvothermal method. Structural analysis with PDF demonstrates that the ferrite nanoparticles synthesized via this method exhibit a significantly more defective structure compared to those synthesized via a solvothermal method.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Origin and resummation of threshold logarithms in the lattice QCD calculations of PDFs

Many present lattice QCD approaches to calculate the parton distribution functions (PDFs) rely on a factorization formula or effective theory expansion of certain Euclidean matrix elements in boosted hadron states. In the quasi- and pseudo-PDF methods, the matching coefficient in the factorization or expansion formula includes large logarithms near the threshold, which arise from the subtle interplay of collinear and soft divergences of an underlying 3D momentum distribution. We use the standard prescription to resum such logarithms in the Mellin-moment space at next-to-leading logarithmic accuracy, which also accounts for the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution, and we show that it can suppress the PDF at large x . Unlike the deep inelastic scattering and Drell-Yan cross sections, the resummation formula is away from the Landau pole. We then apply our formulation to reanalyze the recent lattice results for the pion valence PDF, and find that within the current data sensitivity, the effect of threshold resummation is marginal for the accessible moments and the PDF at large x.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Position-dependent Voronoi probability distribution functions for matter and halos

Here, we measure the Voronoi density probability distribution function (PDF) for both dark matter and halos in N-body simulations. For the dark matter, Voronoi densities represent the matter density field smoothed on a uniform mass scale, which approximates the Lagrangian density field. For halos, the Voronoi densities contain information about the local environment of each halo. We measure the halo virial masses, the total amount of dark matter within each halo Voronoi cell, and the halo Voronoi cell volumes, and we show how halo abundances depend on these three quantities. We then study the position-dependent Voronoi density PDF, measured within finite subregions of the Universe, using separate universe simulations. We demonstrate that the spatial variation of the position-dependent PDF is due to large-scale density fluctuations, indicating that the position-dependent PDF is a biased tracer of large-scale structure. We measure this bias for the dark matter, and interpret it as the bias of regions of the Lagrangian density field that are selected based on density. For the halos, this bias can be interpreted as a form of assembly bias. We present the mapping from late-time to early-time Voronoi density for each simulation dark matter particle, which is highly stochastic. We compare the median of this stochastic map with spherical collapse calculations and discuss challenges involved in modeling the evolution of the density field on these scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Impact of heavy quark and quarkonium data on nuclear gluon PDFs

A clear understanding of nuclear parton distribution functions (nPDFs) plays a crucial role in the interpretation of collider data taken at the Relativistic Heavy Ion Collider, the Large Hadron Collider (LHC), and in the near future at the Electron-Ion Collider. Even with the recent inclusions of vector boson and light meson production data, the uncertainty of the gluon PDF remains substantial and limits the interpretation of heavy ion collision data. To obtain new constraints on the nuclear gluon PDF, we extend our recent nCTEQ15WZ+SIH analysis to inclusive quarkonium and open heavy-flavor meson production data from the LHC. This vast new data set covers a wide kinematic range and puts strong constraints on the nuclear gluon PDF down to x ≲ 10 -5 . The theoretical predictions for these data sets are obtained from a data-driven approach, where proton-proton data are used to determine effective scattering matrix elements. This approach is validated with detailed comparisons to existing next-to-leading order calculations in nonrelativistic QCD for quarkonia and in the general-mass variable-flavor-number scheme for the open heavy-flavored mesons. In addition, the uncertainties from the data-driven approach are determined using the Hessian method and accounted for in the PDF fits. This extension of our previous analyses represents an important step toward the next generation of PDFs not only by including new data sets, but also by exploring new methods for future analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗