Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of parton distributions in a pion with Bézier parametrizations

We explore the role of parametrizations for nonperturbative QCD functions in global analyses, with a specific application to extending a phenomenological analysis of the parton distribution functions (PDFs) in the charged pion realized in the xFitter fitting framework. The parametrization dependence of PDFs in our pion fits substantially enlarges the uncertainties from the experimental sources estimated in the previous analyses. We systematically explore the parametrization dependence by employing a novel technique to automate generation of polynomial parametrizations for PDFs that makes use of Bézier curves. This technique is implemented in a ++ module that is included in the xFitter program. Our analysis reveals that the sea and gluon distributions in the pion are not well disentangled, even when considering measurements in leading-neutron deep inelastic scattering. For example, the pion PDF solutions with a vanishing gluon and large quark sea are still experimentally allowed, which elevates the importance of ongoing lattice and nonperturbative QCD calculations, together with the planned pion scattering experiments, for conclusive studies of the pion structure. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Solar neutrino measurements using the full data period of Super-Kamiokande-IV

An analysis of solar neutrino data from the fourth phase of Super-Kamiokande (SK-IV) from October 2008 to May 2018 is performed and the results are presented. The observation time of the dataset of SK-IV corresponds to 2970 days and the total live time for all four phases is 5805 days. For more precise solar neutrino measurements, several improvements are applied in this analysis: lowering the data acquisition threshold in May 2015, further reduction of the spallation background using neutron clustering events, precise energy reconstruction considering the time variation of the PMT gain. The observed number of solar neutrino events in 3.49–19.49 MeV electron kinetic energy region during SK-IV is 65,443 − 388 + 390 ( stat . ) ± 925 ( syst . ) events. Corresponding B 8 solar neutrino flux is ( 2.314 ± 0.014 ( stat . ) ± 0.040 ( syst . ) ) × 10 6 cm − 2 s − 1 , assuming a pure electron-neutrino flavor component without neutrino oscillations. The flux combined with all SK phases up to SK-IV is ( 2.336 ± 0.011 ( stat . ) ± 0.043 ( syst . ) ) × 10 6 cm − 2 s − 1 . Based on the neutrino oscillation analysis from all solar experiments, including the SK 5805 days dataset, the best-fit neutrino oscillation parameters are sin 2 θ 12 , solar = 0.306 ± 0.013 and Δ m 21 , solar 2 = ( 6.1 0 − 0.81 + 0.95 ) × 10 − 5 eV 2 , with a deviation of about 1.5 σ from the Δ m 21 2 parameter obtained by KamLAND. The best-fit neutrino oscillation parameters obtained from all solar experiments and KamLAND are sin 2 θ 12 , global = 0.307 ± 0.012 and Δ m 21 , global 2 = ( 7.5 0 − 0.18 + 0.19 ) × 10 − 5 eV 2 . Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploring solar neutrino oscillation parameters with the liquid scintillator counter at Yemilab with a comparison to JUNO

We investigate the sensitivities of the liquid scintillator counter (LSC) at Yemilab and JUNO to solar neutrino oscillation parameters, focusing on θ 12 and Δ m 21 2 . We compare the potential of JUNO with LSC at Yemilab utilizing both reactor and solar data in determining those parameters. We find that the solar neutrino data of LSC at Yemilab is highly sensitive to θ 12 enabling its determination with a higher precision compared to reactor experiments. Our study also reveals that if Δ m 21 2 is larger, with a value close to the best-fit value of KamLAND, JUNO reactor data will have about two times better precision than the reactor LSC at Yemilab. On the other hand, if Δ m 21 2 is smaller and closer to the best-fit value of solar neutrino experiments, the precision of the reactor LSC at Yemilab will be better than JUNO. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Testing hadronic-model predictions of depth of maximum of air-shower profiles and ground-particle signals using hybrid data of the Pierre Auger Observatory

We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, X max , and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, S ( 1000 ) , using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists of fitting the measured two-dimensional ( S ( 1000 ) , X max ) distributions using templates for simulated air showers produced with hadronic interaction models pos-, et--04, 2.3d and leaving the scales of predicted X max and the signals from hadronic component at ground as free-fit parameters. The method relies on the assumption that the mass composition remains the same at all zenith angles, while the longitudinal shower development and attenuation of ground signal depend on the mass composition in a correlated way. The analysis was applied to 2239 events detected by both the fluorescence and surface detectors of the Pierre Auger Observatory with energies between 10 18.5 eV to 10 19.0 eV and zenith angles below 60°. We found, that within the assumptions of the method, the best description of the data is achieved if the predictions of the hadronic interaction models are shifted to deeper X max values and larger hadronic signals at all zenith angles. Given the magnitude of the shifts and the data sample size, the statistical significance of the improvement of data description using the modifications considered in the paper is larger than 5 σ even for any linear combination of experimental systematic uncertainties. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Decoding the B → K ν ν excess at Belle II: Kinematics, operators, and masses

An excess in the branching fraction for B + → K + ν ν recently measured at Belle II may be a hint of new physics. We perform thorough likelihood analyses for different new physics scenarios such as B → K X with a new invisible particle X , or B → K χ χ through a scalar, vector, or tensor current with χ being a new invisible particle or a neutrino. We find that vector-current three-body decay with m X ≃ 0.6 GeV—which may be dark matter—is most favored, while two-body decay with m X ≃ 2 GeV is also competitive. The best-fit branching fractions for the scalar and tensor cases are a few times larger than for the two-body and vector cases. Past measurements provide further discrimination, although the best-fit parameters stay similar. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Robustness of the Galactic Center excess morphology against masking

The Galactic Center excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the gamma-ray emission from point sources and the galactic disk are excluded. Using different masks constructed from Fermi point source catalogs and a wavelet method, we find that the GCE morphology, particularly its ellipticity and cuspiness, is relatively independent of the choice of mask for energies above 2–3 GeV. The GCE morphology systematically favors an approximately spherical shape, as expected for dark matter annihilation. Compared to various stellar bulge profiles, a spherical dark matter annihilation profile better fits the data across different masks and galactic diffuse emission backgrounds, except for the stellar bulge profile which provides a similar fit to the data. Modeling the GCE with two components, one from dark matter annihilation and one tracing the Coleman bulge, we find this two-component model outperforms any single component or combinations of dark matter annihilation and other stellar bulge profiles. Uncertainty remains about the exact fraction contributed by each component across different background models and masks. Furthermore, when the Coleman bulge dominates, its corresponding spectrum lacks characteristics typically associated with millisecond pulsars, suggesting that it mostly models the emission from other sources instead of the GCE that is still present and spherically symmetric.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dark Energy Survey: Galaxy sample for the baryonic acoustic oscillation measurement from the final dataset

In this paper, we present and validate the galaxy sample used for the analysis of the baryon acoustic oscillation (BAO) signal in the Dark Energy Survey (DES) Y6 data. The definition is based on a color and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.6, while ensuring a high-quality photo- z determination. The optimization is performed using a Fisher forecast algorithm, finding the optimal i -magnitude cut to be given by i < 19.64 + 2.894 z ph . For the optimal sample, we forecast an increase in precision in the BAO measurement of ∼ 25 % with respect to the Y3 analysis. Our BAO sample has a total of 15,937,556 galaxies in the redshift range 0.6 < z ph < 1.2 , and its angular mask covers 4 , 273.42 deg 2 to a depth of i = 22.5 . We validate its redshift distributions with three different methods: directional neighborhood fitting algorithm (DNF), which is our primary photo- z estimation; direct calibration with spectroscopic redshifts from VIPERS, which is a spectroscopic galaxy sample that overlaps with our BAO sample and is complete within our selection cuts; and clustering redshift using SDSS galaxies. The fiducial redshift distribution is a combination of these three techniques performed by modifying the mean and width of the DNF distributions to match those of VIPERS and clustering redshift. In this paper, we also describe the methodology used to mitigate the effect of observational systematics, which is analogous to the one used in the Y3 analysis. This paper is one of the two dedicated to the analysis of the BAO signal in DES Y6. In its companion paper, we present the angular diameter distance constraints obtained through the fitting to the BAO scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Improving resbos for the precision needs of the LHC

The resummation calculation (esos) is a widely used tool for the simulation of single vector boson production at colliders. In this work, we develop a significant improvement over the esos code by increasing the accuracy from NNLL + NLO to N 3 LL + NNLO and release the esos v2.0 code. Furthermore, we propose a new nonperturbative function that includes information about the rapidity of the system (IFY). The IFY functional form was fitted to data from fixed target experiments, the Tevatron, and the LHC. We find that the nonperturbative function has mild rapidity dependence based on the results of the fit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Superhorizon isocurvature fluctuations relax tensions

Here, we present a new class of models that have potential to alleviate tensions present in the cosmological data today. We postulate about variation in the sound horizon scale on superhorizon scales, i.e., on scales that are larger than that of the present observable low-redshift Universe (≳ 1Gpc) while at the same time smaller than the largest scales probed by the cosmic microwave background (CMB) (≲ 10Gpc). In this scenario, CMB peaks are naturally smoothed as preferred by the Planck data, while at the same time the low-redshift baryon acoustic oscillation calibration is partially decoupled from the CMB. Taking superhorizon variations in baryon fraction as an example and using approximate modeling, we find improvement in the best fit Planck power spectrum model Δχ 2 ~ 6 for 1 extra degree of freedom with the relevant extension parameter 10 3 σ b = 2.10 ± 0.60, implying about 10% variations in baryon fraction across the Universe. At the same time, S 8 drops by about 1 sigma, easing tension with weak lensing surveys. We find that the combination of Planck 2018 data, eBOSS BAO data, and Riess et al. distance ladder Hubble parameter determination produce a satisfactory fit in our model if we allow for a phantom dark energy equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DESI 2024: Constraints on physics-focused aspects of dark energy using DESI DR1 BAO data

Baryon acoustic oscillation data from the first year of the Dark Energy Spectroscopic Instrument (DESI) provide near percent-level precision of cosmic distances in seven bins over the redshift range z=0.1–4.2. Here, this paper is the follow-up to the original DESI BAO cosmology paper [A. G. Adame et al. (DESI Collaboration), arXiv:2404.03002], which considered the conventional w 0 w a cold dark matter (CDM) model. We use the novel DESI data, together with other cosmic probes, to constrain the background expansion history using some well-motivated physical classes of dark energy. In particular, we explore three physics-focused behaviors of dark energy from the equation of state and energy density perspectives: the thawing class (matching many simple quintessence potentials), emergent class (where dark energy comes into being recently, as in phase transition models), and mirage class [where phenomenologically the distance to cosmic microwave background (CMB) last scattering is close to that from a cosmological constant Λ despite dark energy dynamics]. All three classes fit the data at least as well as Λ ⁢CDM, and indeed can improve on it by Δ⁢χ 2 ≈ –5 to –17 for the combination of DESI BAO with CMB and supernova data while having one more parameter. The mirage class does essentially as well as w 0 ⁢w a CDM and exhibits moderate to strong Bayesian evidence preference with respect to Λ⁢ CDM. These classes of dynamical behaviors highlight worthwhile avenues for further exploration into the nature of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Appraising constrained second-order power corrections in HQET with Λ b →Λ c Ι⁢ν

We derive the Λ b →Λ c form factors for the Standard Model and beyond at second order in heavy quark effective theory, applying the recently proposed residual chiral expansion to reduce the set of unknown subsubleading hadronic functions to a single, highly constrained function, that is fully determined by hadron mass parameters at zero recoil. We fit a form factor parametrization based on these results to all available lattice QCD predictions and experimental data. We find that the constrained and predictive structure of the form factors under the residual chiral expansion is in excellent agreement with lattice QCD predictions and experimental data, as well as prior heavy-quark-effective-theory-based fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Finite width effects in nonleptonic D -meson decays

Many analyses of two-body nonleptonic decays of D -mesons rely on flavor S U ( 3 ) symmetry relations and fits of experimental data of decays rates to extract the universal transition amplitudes. Such fits assume that the final state mesons are well-defined asymptotic states of QCD. We develop a technique to take into account the finite width effects of the final state mesons and study their effects on the extracted values of transition amplitudes. Published by the American Physical Society 2025

Kumar, Girish (ORCID:0000000160512495)↗
DOE OSTI2025

Precision measurements of EFT parameters and BAO peak shifts for the Lyman- α forest

We present precision measurements of the bias parameters of the one-loop power spectrum model of the Lyman- α (Ly- α ) forest, derived within the effective field theory (EFT) of large-scale structure. We fit our model to the three-dimensional flux power spectrum measured from the ACCEL 2 hydrodynamic simulations. The EFT model fits the data with an accuracy of below 2% up to k = 2 h Mpc − 1 . Further, we analytically derive how nonlinearities in the three-dimensional clustering of the Ly- α forest introduce biases in measurements of the baryon acoustic oscillations (BAOs) scaling parameters in radial and transverse directions. From our EFT parameter measurements, we obtain a theoretical error budget of Δ α ∥ = − 0.2 % ( Δ α ⊥ = − 0.3 % ) for the radial (transverse) parameters at redshift z = 2.0 . This corresponds to a shift of − 0.3 % (0.1%) for the isotropic (anisotropic) distance measurements. We provide an estimate for the shift of the BAO peak for Ly- α -quasar cross-correlation measurements assuming analytical and simulation-based scaling relations for the nonlinear quasar bias parameters resulting in a shift of − 0.2 % ( − 0.1 % ) for the radial (transverse) dilation parameters, respectively. This analysis emphasizes the robustness of Ly- α forest BAO measurements to the theory modeling. We provide informative priors and an error budget for measuring the BAO feature—a key science driver of the currently observing Dark Energy Spectroscopic Instrument (DESI). Our work paves the way for full-shape cosmological analyses of Ly- α forest data from DESI and upcoming surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST. Published by the American Physical Society 2025

de Belsunce, Roger (ORCID:0000000336604028)↗
DOE OSTI2025

Bootstrap-determined p values in lattice QCD

We present a general method to determine the probability that stochastic Monte Carlo data, in particular those generated in a lattice QCD calculation, would have been obtained were that data drawn from the distribution predicted by a given theoretical hypothesis. Such a probability, or p -value, is often used as an important heuristic measure of the validity of that hypothesis. The proposed method offers the benefit that it remains usable in cases where the standard Hotelling T 2 methods based on the conventional χ 2 statistic do not apply, such as for uncorrelated fits. Specifically, we analyze q 2 , defined as the correlated χ 2 statistic obtained using an arbitrary covariance matrix estimator, and show how to use the bootstrap as a data-driven method to determine the expected distribution of q 2 for a given hypothesis with minimal assumptions. This distribution can then be used to determine the p -value for a fit to the data. We also describe a bootstrap approach for quantifying the impact upon this p -value of estimating population parameters from a single ensemble of N samples. The overall method is accurate up to a 1 / N bias which we do not attempt to quantify. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Correlation Function Approach for Diffusion in Confined Geometries

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nano channel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

Nanoscale flow, Transport, molecular simulation↗
DOE OSTI2020

Hotspot parameter scaling with velocity and yield for high-adiabat layered implosions at the National Ignition Facility

Here, in this paper, we present a study on hotspot parameters in indirect-drive, inertially confined fusion implosions as they proceed through the self-heating regime. The implosions with increasing nuclear yield reach the burning-plasma regime, hotspot ignition, and finally propagating burn and ignition. These implosions span a wide range of alpha heating from a yield amplification of 1.7–2.5. We show that the hotspot parameters are explicitly dependent on both yield and velocity and that by fitting to both of these quantities the hotspot parameters can be fit with a single power law in velocity. The yield scaling also enables the hotspot parameters extrapolation to higher yields. This is important as various degradation mechanisms can occur on a given implosion at fixed implosion velocity which can have a large impact on both yield and the hotspot parameters. The yield scaling also enables the experimental dependence of the hotspot parameters on yield amplification to be determined. The implosions reported have resulted in the highest yield (1.73×10 16 ±2.6%), yield amplification, pressure, and implosion velocity yet reported at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Characteristics of the Diffuse Astrophysical Electron and Tau Neutrino Flux with Six Years of IceCube High Energy Cascade Data

In this work, we report on the first measurement of the astrophysical neutrino flux using particle showers (cascades) in IceCube data from 2010-2015. Assuming standard oscillations, the astrophysical neutrinos in this dedicated cascade sample are dominated (~90%) by electron and tau flavors. The flux, observed in the sensitive energy range from 16 TeV to 2.6 PeV, is consistent with a single power-law model as expected from Fermi-type acceleration of high energy particles at astrophysical sources. We find the flux spectral index to be γ=2.53±0.07 and a flux normalization for each neutrino flavor of $\phi$ astro =1.66$_{-0.27}^{+0.25}$ at E 0 =100 TeV, in agreement with IceCube's complementary muon neutrino results and with all-neutrino flavor fit results. In the measured energy range we reject spectral indices γ≤2.28 at ≥3σ significance level. Because of high neutrino energy resolution and low atmospheric neutrino backgrounds, this analysis provides the most detailed characterization of the neutrino flux at energies below ~100 TeV compared to previous IceCube results. Results from fits assuming more complex neutrino flux models suggest a flux softening at high energies and a flux hardening at low energies (p value ≥0.06). The sizable and smooth flux measured below ~100 TeV remains a puzzle. In order to not violate the isotropic diffuse gamma-ray background as measured by the Fermi Large Area Telescope, it suggests the existence of astrophysical neutrino sources characterized by dense environments which are opaque to gamma rays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗