Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution factor”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Estimated Hourly Electricity Demand Profiles for Each County in the Contiguous United States

In this data descriptor, we present hourly electricity demand estimates for each county in the contiguous United States from 2016 to 2023. The demand profiles represent the sum of two components for each county. The first is the weighted average of reported hourly demand profiles for North American Electric Reliability Corporation balancing authority (BA) regions and subregions, scaled to match annual estimates of county-level retail sales and direct use of electricity and weighted by the estimated percentage of county load served by each BA region or subregion. The second is the weighted average of modeled hourly, county- and sector-level distributed photovoltaic (DPV) capacity factor profiles, scaled to match annual estimates of on-site consumption of DPV-generated electricity for each county and weighted by the percentage of consumption attributable to each sector.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Review of Current Software for Analyzing Total X-ray Scattering Data from Liquids

The popularity of the pair distribution function (PDF) analysis of X-ray total scattering data has steadily grown as access to ex situ synchrotron data has expanded. Due to the broadening of the PDF user community, there is a growing demand for software that can be used to extract PDFs and is accessible to non-expert users. While user-friendly options have been developed over the past decade for fast, streamlined data analysis, care must be taken in both processing the data and understanding any limitations, especially in the case of liquids. In this review, the same scattering data are analyzed using different total X-ray scattering software, in order to compare the accuracy of the extracted structure factors and associated pair distribution functions. The goal is to assess the best practices for extracting the most accurate liquid data for each software package. The importance of absolute normalization and the application of the most appropriate corrections are emphasized via quantitative comparisons between liquid sulfur and water. Additionally, an awareness of the competing conventions used to define the PDF in crystallography and liquids/glasses is crucial for both the downstream analyses of the data and a comparison with the previous results in the literature.

Gallington, Leighanne C. (ORCID:0000000203837522)↗
DOE OSTI2025

Analysis of Two Models for the Angular Structure of the Outflows Producing the Swift/XRT “Larger-angle Emission” of Gamma-Ray Bursts

The quasi-instantaneous emission from a relativistic surface endowed with a Lorentz factor that decreases away from the outflow symmetry axis can naturally explain the three phases observed by Swift X-Ray Telescope (XRT) in gamma-ray bursts (GRBs) and their afterglows (GRB tail, afterglow plateau, and postplateau) based only on the angular change of the relativistic Doppler boost across the outflow surface. We develop further the analytical formalism of the “larger-angle emission” model for the case of “n-exponential” outflows (where the Lorentz factor Γ dependence of the angular location θ is Γ ~ exp{-(θ/θ c ) n }), and compare its ability to account for the X-ray emission of XRT afterglows relative to that of “power-law” outflows (Γ ∼ θ −g ). Power-law outflows yield longer afterglow plateaus, followed by slower postplateau flux decays than n-exponential outflows, features which may be used in identifying which type of angular structure is at work in a given afterglow. Identifying the Γ(θ) angular structure that accommodates XRT light curves is slightly complicated by the fact that the afterglow X-ray light curve is also determined by how two characteristics of the comoving-frame emission spectrum (peak energy $E'_p$ and peak intensity $i'_p$) change with the angular location or, equivalently, with the Lorentz factor. Here, we assume power-law Γ dependences of those spectral characteristics and find that, unlike power-law outflows, n-exponential outflows cannot account for plateaus with a temporal dynamical range larger than 100 (2 dex in logarithmic space). To capture all the information contained in XRT afterglow measurements (0.3–10 keV unabsorbed flux and effective spectral slope), we calculate 0.3 and 10 keV light curves using a broken-power-law emission spectrum of peak energy and low- and high-energy slopes that are derived from the effective slope measured by XRT. This economical peak energy determination is found to be consistent with the results of more expensive spectral fits. The angular distributions of the Lorentz factor, comoving frame peak energy, and peak intensity (Γ(θ), $E'_p$(θ), $i'_p$(θ)) constrain the (yet-to-be determined) convolution of various features of the production of relativistic jets by solar-mass black holes and of their propagation through the progenitor/circumburst medium, while the $E'_p$(Γ) and $i'_p$(Γ) dependences may constrain the GRB dissipation mechanism and the GRB emission process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Transverse momentum distributions at large- x

We investigate the collinear matching of transverse momentum dependent (TMD) distributions at large values of x, computing and resumming the leading large-x asymptotics for matching coefficients. The large-x resummation is done directly within TMD distributions, ensuring the process-independence of the result. The derived resummation formulas are valid for all TMD distributions (except the pretzelosity). Their application improves perturbative convergence, provides practical estimation for unknown higher-order contributions, and sets restrictions for the nonperturbative part of models. Using the known anomalous dimensions, resummation can reach N 3 LL, often exceeding the accuracy of known coefficient functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Whole-Genome SNP Analysis Identifies Putative Mycobacterium bovis Transmission Clusters in Livestock and Wildlife in Catalonia, Spain

The high-resolution WGS analyses of MTBC strains have provided useful insight for determining sources of infection for animal tuberculosis. In Spain, tuberculosis in livestock is caused by Mycobacterium bovis and Mycobacterium caprae, where wildlife reservoirs play an important role. We analyzed a set of 125 M. bovis isolates obtained from livestock and wildlife from Catalonia to investigate strain diversity and identify possible sources and/or causes of infection. Whole-genome SNP profiles were used for phylogenetic reconstruction and pairwise SNP distance analysis. Additionally, SNPs were investigated to identify virulence and antimicrobial resistance factors to investigate clade-specific associations. Putative transmission clusters (≤12 SNPs) were identified, and associated epidemiological metadata were used to determine possible explanatory factors for transmission. M. bovis distribution was heterogeneous, with 7 major clades and 21 putative transmission clusters. In order of importance, the explanatory factors associated were proximity and neighborhood, residual infection, livestock-wildlife interaction, shared pasture, and movement. Genes related to lipid transport and metabolism showed the highest number of SNPs. All isolates were pyrazinamide resistant, and five were additionally resistant to isoniazid, but no clade-specific associations could be determined. Our findings highlight the importance of high-resolution molecular surveillance to monitor bovine tuberculosis dynamics in a low-prevalence setting.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Investigating the weak charge of 48 Ca using a dispersive optical model

A new nonlocal dispersive-optical-model analysis has been carried out for neutrons and protons in 48 Ca that reproduces the weak-form-factor measurement of CREX. In addition to elastic-scattering angular distributions, total and reaction cross sections, single-particle energies, the neutron and proton numbers, and the charge distribution, the CREX-measured weak form factor has been fit to extract the neutron and proton self-energies both above and below the Fermi energy. The resulting single-particle propagators yield a weak form factor of F w = 0.125 ± 0.05 and a neutron skin of R skin = 0.152 ± 0.05 fm, in good agreement with CREX. The rearrangement of the neutron distribution to accommodate such a thin neutron skin results in the high-momentum content of the neutrons exceeding that of the protons, in contrast to what is expected from high-energy two-nucleon knockout measurements by the CLAS collaboration and ab initio asymmetric matter calculations. The present analysis also emphasizes the importance of neutron experimental data in constraining weak charge observables necessary for a precise description of neutron densities. Notably, the neutron reaction cross section and further parity-violating experiments weak form factor measurements are essential to generate a unique way to determine the 48 Ca neutron distribution in this framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hourly Electricity Demand Profiles for Each County in the Contiguous United States

This dataset provides estimated hourly electricity demand for each county in the contiguous United States from 2016-2023. The demand profiles represent the sum of two components: (1) Weighted averages of reported hourly demand profiles for North American Electric Reliability Corporation balancing authority (BA) regions and subregions, scaled to match annual estimates of county-level retail sales and direct use of electricity and weighted by the estimated percentage of county load served by each BA region or subregion. (2) Weighted averages of modeled hourly, county- and sector-level distributed photovoltaic (DPV) capacity factor profiles, scaled to match annual estimates of on-site consumption of DPV-generated electricity for each county and weighted by the percentage of consumption attributable to each sector Annual county-level retail sales are estimated by aggregating utility-reported sales to the state level and allocating the results to counties according to each county's share of state population. Annual county-level direct use is calculated by aggregating power plant-reported direct use values. Annual county-level on-site consumption of DPV-generated electricity is estimated by aggregating utility-reported net metering data to determine the amount of DPV-generated electricity sold back to the grid for each state, subtracting those values from modeled state-level DPV generation estimates, and allocating the results to counties according to each county's share of statewide modeled DPV generation. The open-source Python code used to develop this dataset is available at "Historical Load Data Repository" link below.

14 SOLAR ENERGY↗
DOE OSTI2024

EPCAPE-PT-LANL Measurements: Wideband Integrated Bioaerosol Sensor

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Wideband Integrated Bioaerosol Sensor (Droplet Measurements Technology) Data Notes: The WIBS is an online single-particle measurement that detects FBAPs (within a size range of 0.5 - 30 microns in diameter) based on the excitation and emission wavelengths of the individual particles. Using two xenon lamps, the WIBS excites FBAPs at 280 nm and 370 nm. Their emission is detected across two wavebands of 310-400 nm and 420-650 nm. We classified the FBAPs into seven different categories (A, B, C, AB, BC, AC, and ABC) using the classification scheme in Perring et. al. (2015) [1]. Averaged number concentration of FBAPs (total and by category) and particles that non-fluorescent bioaerosols particles (NFBAPs). In separate files, we also present one-minute-averaged size distributions and the asymmetry factor (AF, a surrogate for shape) of all FBAPs and NFBAPs. The logarithmic bin width of the size bins are the same as the average bin width of the AOS's optical particle counter (OPC, Grimm) for the range of sizes in which they overlap (26 bins from 0.5 - 30 microns). AF of the particles ranges from 0-100 and is divided into five bins with a linear spacing at increments of 20. The smallest AF bin represents more spherical particles while the largest bin represents more rod-shaped particles. [1] Perring, A. E., et al. (2015), Airborne observations of regional variation in fluorescent aerosol across the United States, J. Geophys. Res. Atmos., 120, 1153–1170, doi:10.1002/2014JD022495. Abstract and description of the campaign can be found here : https://www.arm.gov/research/campaigns/amf2023epcape-pt-lanl. Files data_10min_WIBS_AFDist.csv Header: - FBAP_AFDist[/cm3]_Bin_1 to Bin_5: Concentration of fluorescent bioaerosol particles in the each of 5 AF bins, measured in particles per cubic centimeter. Each bin represents a specific range of particle AF, capturing the shapes of FBAPs detected during the measurement. - NFBAP_AFDist[/cm3]_Bin_1 to Bin_5: Concentration of non-fluorescent bioaerosol particles in the each of 5 AF bins, measured in particles per cubic centimeter. Each bin represents a specific range of particle AF, capturing the shapes of NFBAPs detected during the measurement. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement. AF Bins: • Bin 1: 0 – 20 [unitless] • Bin 2: 21 – 40 [unitless] • Bin 3: 41 – 60 [unitless] • Bin 4: 61 – 80 [unitless] • Bin 5: 81 – 100 [unitless] Files data_10min_WIBS_Conc.csv Header: - NumberConcentrationA[/cm3]: Number concentration of bioaerosol particles detected by fluorescence channel A, measured in particles per cubic centimeter. - NumberConcentrationB[/cm3]: Number concentration of bioaerosol particles detected by fluorescence channel B, measured in particles per cubic centimeter. - NumberConcentrationC[/cm3]: Number concentration of bioaerosol particles detected by fluorescence channel C, measured in particles per cubic centimeter. - NumberConcentrationAB[/cm3]: Combined number concentration of bioaerosol particles detected by both fluorescence channels A and B, measured in particles per cubic centimeter. - NumberConcentrationBC[/cm3]: Combined number concentration of bioaerosol particles detected by both fluorescence channels B and C, measured in particles per cubic centimeter. - NumberConcentrationAC[/cm3]: Combined number concentration of bioaerosol particles detected by both fluorescence channels A and C, measured in particles per cubic centimeter. - NumberConcentrationABC[/cm3]: Combined number concentration of bioaerosol particles detected by all three fluorescence channels A, B, and C, measured in particles per cubic centimeter. - NumberConcentrationNFBAP[/cm3]: Number concentration of non-fluorescent bioaerosol particles, measured in particles per cubic centimeter. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement. Files data_10min_WIBS_SizeDist.csv Header: - FBAP_SizeDist[/cm3]_Bin_1 to FBAP_SizeDist[/cm3]_Bin_26: Number concentrations of FBAP in each of 26 size bins, measured in particles per cubic centimeter. Each bin represents a specific range of particle sizes, capturing the size distribution of FBAPs detected during the measurement. - NFBAP_SizeDist[/cm3]_Bin_1 to NFBAP_SizeDist[/cm3]_Bin_26: Number concentrations of NFBAP in each of 26 size bins, measured in particles per cubic centimeter. Similar to FBAP, each bin covers a specific range of particle sizes, detailing the size distribution of NFBAPs detected. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement. Size Bins: • Bin 1: 0.48 to 0.57 μm • Bin 2: 0.57 to 0.67 μm • Bin 3: 0.67 to 0.79 μm • Bin 4: 0.79 to 0.93 μm • Bin 5: 0.93 to 1.1 μm • Bin 6: 1.1 to 1.29 μm • Bin 7: 1.29 to 1.52 μm • Bin 8: 1.52 to 1.8 μm • Bin 9: 1.8 to 2.11 μm • Bin 10: 2.11 to 2.5 μm • Bin 11: 2.5 to 2.94 μm • Bin 12: 2.94 to 3.46 μm • Bin 13: 3.46 to 4.08 μm • Bin 14: 4.08 to 4.81 μm • Bin 15: 4.81 to 5.67 μm • Bin 16: 5.67 to 6.68 μm • Bin 17: 6.68 to 7.88 μm • Bin 18: 7.88 to 9.29 μm • Bin 19: 9.29 to 10.96 μm • Bin 20: 10.96 to 12.92 μm • Bin 21: 12.92 to 15.23 μm • Bin 22: 15.23 to 17.96 μm • Bin 23: 17.96 to 21.17 μm • Bin 24: 21.17 to 24.96 μm • Bin 25: 24.96 to 29.43 μm • Bin 26: 29.43 to 34.70 μm

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Distributed-Memory Parallel JointNMF

Joint Nonnegative Matrix Factorization (JointNMF) is a hybrid method for mining information from datasets that contain both feature and connection information. We propose distributed-memory parallelizations of three algorithms for solving the JointNMF problem based on Alternating Nonnegative Least Squares, Projected Gradient Descent, and Projected Gauss-Newton. We extend well-known communication-avoiding algorithms using a single processor grid case to our coupled case on two processor grids. We demonstrate the scalability of the algorithms on up to 960 cores (40 nodes) with 60% parallel efficiency. The more sophisticated Alternating Nonnegative Least Squares (ANLS) and Gauss-Newton variants outperform the first-order gradient descent method in reducing the objective on large-scale problems. We perform a topic modelling task on a large corpus of academic papers that consists of over 37 million paper abstracts and nearly a billion citation relationships, demonstrating the utility and scalability of the methods.

Eswar, Srinivas↗
DOE OSTI2021

Correlating conductivity and Seebeck coefficient to doping within crystalline and amorphous domains in poly(3-(methoxyethoxyethoxy)thiophene)

Molecular doping of conjugated polymers (CPs) plays a vital role in optimizing organic electronic and energy applications. For the case of organic thermoelectrics, it is commonly believed that doping CPs with a strong dopant could result in higher conductivity (σ) and thus better power factor (PF). Herein, by investigating thermoelectric performance of a polar side-chain bearing CP, poly(3- (methoxyethoxyethoxy)thiophene) (P3MEET), vapor doped with fluorinatedderivative of tetracyanoquinodimethane FnTCNQ (n = 1, 2, 4), we show that using strong dopants can in fact have detrimental effects on the thermoelectric performance of CPs. Despite possessing higher electron affinity, doping P3MEET with F4TCNQ only results in a σ (27.0 S/cm) comparable to samples doped with other two weaker dopants F2TCNQ and F1TCNQ (26.4 and 20.1 S/cm). Interestingly, F4TCNQ-doped samples display a marked reduction in the Seebeck coefficient (α) compared to F1TCNQ- and F2TCNQ-doped samples from 42 to 13 μV/ K, leading to an undesirable suppression of the PF. Additionally, structural characterizations coupled with Kang-Snyder modeling of the α–σ relation show that the reduction of α in F4TCNQ-doped P3MEET samples originates from the generation of low mobility carrier within P3MEET's amorphous domain. Our results demonstrate that factors such as dopant distribution and doping efficiency within the crystalline and amorphous domains of CPs should play a crucial role in advancing rational design for organic thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gluon double-spin asymmetry in the longitudinally polarized p + p collisions

We derive the first-ever small-x expression for the inclusive gluon production cross section in the central rapidity region of the longitudinally polarized proton-proton collisions. The cross section depends on the polarizations of both protons, therefore comprising the numerator of the longitudinal double-spin asymmetry ALL for the produced gluons. The cross section is calculated in the shock wave formalism and is expressed in terms of the polarized dipole scattering amplitudes on the projectile and target protons. We show that the small-x evolution corrections are included into our cross section expression if one evolves these polarized dipole amplitudes using the double-logarithmic helicity evolution derived in [1–4]. Our calculation is performed for the gluon sector only, with the quark contribution left for future work. When that work is complete, the resulting formula will be applicable to longitudinally polarized proton-proton and proton-nucleus collisions, as well as to polarized semi-inclusive deep inelastic scattering (SIDIS) on a proton or a nucleus. Our results should allow one to extend the small-x helicity phenomenology analysis of [5] to the jet/hadron production data reported for the longitudinally polarized proton-proton collisions at RHIC and to polarized SIDIS measurements at central rapidities to be performed at the EIC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

All order factorization for virtual Compton scattering at next-to-leading power

We discuss all-order factorization for the virtual Compton process at next-to leading power (NLP) in the Λ QCD /Q and $\sqrt{-t}$/Q expansion (twist-3), both in the double deeply-virtual case and the single-deeply-virtual case. We use the soft-collinear efective theory (SCET) as the main theoretical tool. We conclude that collinear factorization holds in the double-deeply virtual case, where both photons are far of-shell. The agreement is found with the known results for the hard matching coefcients at leading order $α^0_s$, and we can therefore connect the traditional approach with SCET. In the single-deeply-virtual case, commonly called deeply virtual Compton scattering (DVCS), the contribution of non-target collinear regions complicates the factorization. These include momentum modes collinear to the real photon and (ultra)soft interactions between the photon-collinear and target-collinear modes. However, such contributions appear only for the transversely polarized virtual photon at the NLP accuracy and in fact it is the only NLP ~ (Λ QCD /Q) 1 ~ ( $\sqrt{-t}$/Q) 1 contribution in that case. We therefore conclude that the DVCS amplitude for a longitudinally polarized virtual photon, where the leading power ~ (Λ QCD /Q) 0 ~ ($\sqrt{-t}$/Q) 0 contribution vanishes, is free of non-target collinear contributions and the collinear factorization in terms of twist-3 GPDs holds in that case as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Opportunities for installed combined heat and power (CHP) to increase grid flexibility in the U.S.

Increasing use of renewable energy requires sufficient grid flexibility to address uncertainty and variability in electricity generation. Previous studies suggest that combined heat and power (CHP) systems may support grid flexibility but they do not consider operating hours. In this paper, we used CHP operating data and determined annual and monthly availability of the installed CHP capacity from various sectors (e.g., utility, independent power producer, commercial, and industrial) in all seven U.S. independent system operators (ISOs) and regional transmission organizations (RTOs). Also, we estimated hourly CHP availability installed in five facility types (i.e., hospitals, universities, hotels, offices, and manufacturing) in the state of New York. The results show that regardless of ISO/RTO, sector, or season, more than 40% of the installed CHP capacity (0.7–8.7 GW) was not fully utilized in 2019; the results are similar for 2018. This available CHP capacity accounted for up to 9% of the ISO/RTO's peak electric demand, which may yield cost savings up to $16 billion by avoiding installation costs of new natural gas combustion or combined-cycle turbines. To exploit the available CHP capacity to enhance grid flexibility, we recommend different policy implications including flexible contract lengths between CHP owners and grid operators, improved market designs, and simplified interconnection standards.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Multiscale electric-field imaging of polarization vortex structures in PbTiO 3 /SrTiO 3 superlattices

In ferroelectric heterostructures, the interaction between intrinsic polarization and the electric field generates a rich set of localized electrical properties. The local electric field is determined by several connected factors, including the charge distribution of individual unit cells, the interfacial electromechanical boundary conditions, and chemical composition of the interfaces. However, especially in ferroelectric perovskites, a complete description of the local electric field across micro-, nano-, and atomic-length scales is missing. Here, by applying four-dimensional scanning transmission electron microscopy (4D STEM) with multiple probe sizes matching the size of structural features, we directly image the electric field of polarization vortices in (PbTiO 3 ) 16 /(SrTiO 3 ) 16 superlattices and reveal different electric field configurations corresponding to the atomic scale electronic ordering and the nanoscale boundary conditions. The separability of two different fields probed by 4D STEM offers the possibility to reveal how each contributes to the electronic properties of the film.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deterministic High-Fidelity Neutronics Simulation of Pebble Bed Reactors Using Pebble Tracking Transport

The pebble tracking transport (PTT) algorithm offers a high-fidelity deterministic approach for neutron transport for pebble bed reactors (PBRs). This approach requires the mesh for the active-core region to consist exclusively of tetrahedral elements, where each node in the pebble-packing region represents a pebble centroid. This paper investigates the application of PTT for full-scale PBRs, considering both the isothermal and the temperature-dependent core conditions. Macroscopic cross sections are generated using Serpent 2 full-core eigenvalue simulations where pebbles are grouped into disjoint subsets using machine learning. To minimize the need for individual cross-section sets for each pebble in the core, K-means clustering is used to group pebbles by temperature and neutronic environment parameters. Here, we compare the multiplication factor and power rate distributions between PTT simulations using the Griffin reactor physics software and reference solutions from Serpent 2. Our analysis shows that a full-core, high-fidelity PTT calculation produces accurate results with minimal local (pebblewise) errors. Additionally, timing results indicate that PTT simulations converge rapidly on modern supercomputing platforms.

Griffin↗
DOE OSTI2020

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Non-universal stellar initial mass functions: large uncertainties in star formation rates at z ≈ 2–4 and other astrophysical probes

ABSTRACT We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical broken-power-law IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF that varies as a function of the local star formation rate, and explore the effects on the star formation rate density (SFRD), the extragalactic background light, the supernova (both core-collapse and thermonuclear) rates, and the diffuse supernova neutrino background. Our most interesting result is that our adopted varying IMF leads to much greater uncertainty on the SFRD at $z \approx 2-4$ than is usually assumed. Indeed, we find an SFRD (inferred using observed galaxy luminosity distributions) that is a factor of $\gtrsim 3$ lower than canonical results obtained using a universal IMF. Secondly, the non-universal IMF we explore implies a reduction in the supernova core-collapse rate of a factor of $\sim 2$, compared against a universal IMF. The other potential tracers are only slightly affected by changes to the properties of the IMF. We find that currently available data do not provide a clear preference for universal or non-universal IMF. However, improvements to measurements of the star formation rate and core-collapse supernova rate at redshifts $z \gtrsim 2$ may offer the best prospects for discernment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Mixing effects on spectroscopy and partonic observables of heavy mesons with logarithmic confining potential in a light-front quark model

Using the variational principle, we systematically investigate the mass spectra and wave functions of both 1⁢𝑆 and 2⁢𝑆 state heavy pseudoscalar (𝑃) and vector (𝑉) mesons within the light-front quark model. This approach incorporates a Coulomb plus logarithmic confinement potential to accurately describe the constituent quark and antiquark dynamics. Additionally, spin hyperfine interactions are introduced perturbatively to compute the masses of pseudoscalar and vector mesons. The present analyses of the 1⁢𝑆 and 2⁢𝑆 states require the consideration of mixing between them to account for empirical constraints. These constraints include the mass gap Δ⁢𝑀 𝑃 >Δ⁢𝑀 𝑉 , where Δ⁢𝑀 𝑃⁡(𝑉) =𝑀$^{2⁢𝑆}_{𝑃⁡(𝑉)}$−𝑀$^{1⁢𝑆}_{𝑃⁡(𝑉)}$ and the hierarchy of the decay constants 𝑓 1⁢𝑆 >𝑓 2⁢𝑆 . We find the optimal value of the mixing angle to be 𝜃 =1⁢8°, significantly enhancing the consistency between our spectroscopic predictions and the experimental data compiled by the Particle Data Group. Furthermore, based on the predicted mass, the newly observed resonance 𝐵 𝐽⁡ (5840) could be assigned as a 2 1⁢ 𝑆 0 state in the 𝐵 meson family. The study also reports various pertinent observables, including twist-two distribution amplitudes, electromagnetic form factors, charge radii, 𝜉 moments, and transition form factors that are found to be consistent with both available lattice simulations and experimental data. In addition, our predicted branching ratios for the channels of 𝐵 + →𝜏 + ⁢𝜈 𝜏 as well as rare decays of 𝐵 0 and 𝐵$^0_𝑠$ appear in accordance with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗