Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data 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 847 records · Page 47

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗

Nuclear suppression of coherent J / ψ photoproduction in heavy-ion ultraperipheral collisions and leading twist nuclear shadowing

We determine the nuclear suppression factor S Pb ( x ) , where x = M J / ψ 2 / W γ p 2 with M J / ψ the J / ψ mass and W γ p the photon-nucleon energy, for the cross section of coherent J / ψ photoproduction in heavy-ion ultraperipheral collisions at the Large Hadron Collider and Relativistic Heavy Ion Collider by performing the χ 2 fit to all available data on the cross section d σ A A → J / ψ A A / d y as a function of the J / ψ rapidity y and the photoproduction cross section σ γ A → J / ψ A ( W γ p ) as a function of W γ p . We find that while the d σ A A → J / ψ A A / d y data alone constrain S Pb ( x ) for x ≥ 10 − 3 , the combined d σ A A → J / ψ A A / d y and σ γ A → J / ψ A ( W γ p ) data allow us to determine S Pb ( x ) in the wide interval 10 − 5 < x < 0.05 . In particular, the data favor S Pb ( x ) , which decreases with a decrease of x in the 10 − 4 < x < 0.01 interval and can be both decreasing or constant for x < 10 − 4 . Identifying S Pb ( x ) with the ratio of the gluon distributions in Pb and the proton R g ( x , Q 0 2 ) = g A ( x , Q 0 2 ) / [ A g p ( x , Q 0 2 ) ] , we demonstrate that the leading twist approximation for nuclear shadowing provides a good description of all the data on d σ A A → J / ψ A A / d y and σ γ A → J / ψ A ( W γ p ) as well as on the experimental values for S Pb ( x ) derived from σ γ A → J / ψ A ( W γ p ) . We also show that modern nuclear parton distributions reasonably reproduce S Pb ( x ) as well. Published by the American Physical Society 2024

Guzey, V. (ORCID:0000000223938507)↗

Blinded challenge for precision cosmology with large-scale structure: Results from effective field theory for the redshift-space galaxy power spectrum

An accurate theoretical template for the galaxy power spectrum is key for the success of ongoing and future spectroscopic surveys. We examine to what extent the effective field theory (EFT) of large-scale structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a 100 times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge is to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and to pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within subpercent statistical errors corresponding to the total simulation volume.

79 ASTRONOMY AND ASTROPHYSICS↗

Solar neutrinos and ν 2 visible decays to ν 1

Experimental bounds on the neutrino lifetime depend on the nature of the neutrinos and the details of the potentially new physics responsible for neutrino decay. In the case where the decays involve active neutrinos in the final state, the neutrino masses also qualitatively impact how these manifest themselves experimentally. In order to further understand the impact of nonzero neutrino masses, we explore how observations of solar neutrinos constrain a very simple toy model. We assume that neutrinos are Dirac fermions and there is a new massless scalar that couples to neutrinos such that a heavy neutrino— ν 2 with mass m 2 —can decay into a lighter neutrino— ν 1 with mass m 1 —and a massless scalar. We find that the constraints on the new physics coupling depend, sometimes significantly, on the ratio of the daughter-to-parent neutrino masses and that, for large-enough values of the new physics coupling, the “dark side” of the solar neutrino parameter space— sin 2 θ 12 ∼ 0.7 —provides a reasonable fit to solar neutrino data, if only B 8 or Be 7 neutrino data alone are considered, but no allowed region is found in the combined analysis. Our results generalize to other neutrino-decay scenarios, including those that mediate ν 2 → ν 1 ν ¯ 3 ν 3 when the neutrino mass ordering is inverted mass and m 2 > m 1 ≫ m 3 , the mass of ν 3 . Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size L box = 2 Gpc/h and particle mass of 1.01 × 10 10 M ⊙ /h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe" method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, b φ . As a science application, here we investigate the link between halo assembly bias and b φ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between b φ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the b φ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on f NL constraints for the allowed galaxy-halo models (note that σ[f NL ]∝ $\frac{σ[b_φ]}{b_φ}$). We find that the error on f NL is 15, 8, 7 for the LRGs at z = 0.5 and z = 0.8 and QSOs at z = 1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the neutrino-oxygen neutral-current quasielastic cross section using atmospheric neutrinos in the SK-Gd experiment

We report the first measurement of the atmospheric neutrino-oxygen neutral-current quasielastic (NCQE) cross section in the gadolinium-loaded Super-Kamiokande (SK) water Cherenkov detector. In June 2020, SK began a new experimental phase, named SK-Gd, by loading 0.011% by mass of gadolinium into the ultrapure water of the SK detector. The introduction of gadolinium to ultrapure water has the effect of improving the neutron-tagging efficiency. Using a 552.2 day dataset from August 2020 to June 2022, we measure the NCQE cross section to be 0.74 ± 0.22 ( stat ) − 0.15 + 0.85 ( syst ) × 10 − 38 cm 2 / oxygen in the energy range from 160 MeV to 10 GeV, which is consistent with the atmospheric neutrino-flux-averaged theoretical NCQE cross section and the measurement in the SK pure-water phase within the uncertainties. Furthermore, we compare the models of the nucleon-nucleus interactions in water and find that the binary cascade model and the Liège intranuclear cascade model provide a somewhat better fit to the observed data than the Bertini cascade model. Since the atmospheric neutrino-oxygen NCQE reactions are one of the main backgrounds in the search for diffuse supernova neutrino background (DSNB), these new results will contribute to future studies—and the potential discovery—of the DSNB in SK. Published by the American Physical Society 2024

Sakai, S. (ORCID:0000000221900062)↗

Measurement of photonuclear jet production in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector

In ultrarelativistic heavy ion collisions at the LHC, each nucleus acts a sources of high-energy real photons that can scatter off the opposing nucleus in ultraperipheral photonuclear (𝛾 + 𝐴) collisions. Hard scattering processes initiated by the photons in such collisions provide a novel method for probing nuclear parton distributions in a kinematic region not easily accessible to other measurements. ATLAS has measured production of dijet and multijet final states in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV using a dataset recorded in 2018 with an integrated luminosity of 1.72 nb −1 . Photonuclear final states are selected by requiring a rapidity gap in the photon direction; this selects events where one of the outgoing nuclei remains intact. Jets are reconstructed using the anti-𝑘 t algorithm with radius parameter, 𝑅 = 0.4. Triple-differential cross sections, unfolded for detector response, are measured and presented using two sets of kinematic variables. The first set consists of the total transverse momentum (𝐻 T ), rapidity, and mass of the jet system. The second set uses 𝐻 T and particle-level nuclear and photon parton momentum fractions, 𝑥 A and 𝑧 𝛾 , respectively. The results are compared with leading-order perturbative QCD calculations of photonuclear jet production cross sections, where all leading order predictions using existing fits fall below the data in the shadowing region. More detailed theoretical comparisons will allow these results to strongly constrain nuclear parton distributions, and these data provide results from the LHC directly comparable to early physics results at the planned Electron-Ion Collider.

Parton distribution functions↗

Study of High-Transverse-Momentum Higgs Boson Production in Association with a Vector Boson in the qqbb Final State with the ATLAS Detector

This Letter presents the first study of Higgs boson production in association with a vector boson (V=W or Z) in the fully hadronic qqbb final state using data recorded by the ATLAS detector at the LHC in proton-proton collisions at $\sqrt{s}$=13 TeV and corresponding to an integrated luminosity of 137 fb –1 . The vector bosons and Higgs bosons are each reconstructed as large-radius jets and tagged using jet substructure techniques. Dedicated tagging algorithms exploiting b-tagging properties are used to identify jets consistent with Higgs bosons decaying into $b\bar{b}$. Dominant backgrounds from multijet production are determined directly from the data, and a likelihood fit to the jet mass distribution of Higgs boson candidates is used to extract the number of signal events. The VH production cross section is measured inclusively and differentially in several ranges of Higgs boson transverse momentum: 250–450, 450–650, and greater than 650 GeV. The inclusive signal yield relative to the standard model expectation is observed to be μ=1.4$^{+1.0}_{–0.9}$ and the corresponding cross section is 3.1±1.3⁢(stat)$^{+1.8}_{–1.4}$(syst.) pb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Shot-noise-induced lower temperature limit of the nonneutral plasma parallel temperature diagnostic

Abstract We develop a new algorithm to estimate the temperature of a nonneutral plasma in a Penning-Malmberg trap. The algorithm analyzes data obtained by slowly lowering a voltage that confines one end of the plasma and collecting escaping charges, and is a maximum likelihood estimator based on a physically-motivated model of the escape protocol presented in (Beck in Measurement of the magnetic and temperature dependence of the electron-electron anisotropic temperature relaxation rate. PhD thesis, 1990). Significantly, our algorithm may be used on single-count data, allowing for improved fits with low numbers of escaping electrons. This is important for low-temperature plasmas such as those used in antihydrogen trapping. We perform a Monte Carlo simulation of our algorithm, and assess its robustness to intrinsic shot noise and external noise. The assumptions in this paper allow for a lower bound for measurable plasma temperatures of approximately $3\,\mathrm{K}$ 3 K for plasmas of length $1\,\mathrm{cm}$ 1 cm , with approximately 100 particle counts needed for an accuracy of $\pm 10 \%$ ± 10 % .

Zhong, Adrianne (ORCID:0000000162618736)↗

Parametric Study of the Influence of Support Type, Presence of Platinum on Support, and Ionomer Content on the Microstructure of Polymer Electrolyte Fuel Cell Catalyst Layers

We report Ultra-small angle X-ray scattering (USAXS) was employed to investigate the effects of carbon support type, the presence of platinum on carbon, and ionomer loading on the microstructure of polymer electrolyte fuel cell (PEFC) catalyst layers (CLs). Particle size distributions (PSDs), obtained from fitting the measured scattering data were used to interpret the size of carbon aggregates (40–300 nm) and agglomerates (>400 nm) from two-component carbon/ionomer and three-component platinum/carbon/ionomer CLs. Two types of carbon supports were investigated: high surface area carbon (HSC) and Vulcan XC-72. CLs with a range of perfluorosulfonic acid (PFSA) ionomer to carbon (I/C) ratios (0.2–1.2) and also with perfluoroimide acid (PFIA) ionomer were studied to evaluate the effect of ionomer on CL microstructure. The carbon type, the presence of platinum, and ionomer loading were all found to significantly impact carbon agglomeration. The extent of Pt/C agglomeration in the CL was found to increase with increasing ionomer and platinum concentration and to decrease with increasing carbon surface area. Platinum electrochemically-active surface area (ECSA) and local oxygen transport resistance (R nF ) were correlated to the CL microstructure to yield relationships affecting electrode performance.

25 ENERGY STORAGE↗

Modeling Sorption of Water and Vanadium Cations by Ion-Exchange Membranes

We develop a macroscopic model for sorption of ions and solvent by ion-exchange membranes and use it to investigate how different membranes behave when exposed to electrolytes for vanadium flow batteries. The model combines a classical expression for the Gibbs free energy of a moderately concentrated electrolyte containing several ions with a term that describes elastic swelling of a polymer membrane. Parameters describing pairwise interactions between mobile ions and fixed membrane ions were fit to published experimental data and discussed in the context of the competition for sites in ion-exchange membranes. The model provides insights into why sorption of vanadium cations by Nafion is lower than predicted by Donnan theory and why switching to an anion-exchange membrane does not dramatically reduce vanadium sorption.

Electrochemistry↗

Modeling and Analysis of Polarization Losses in Solid Oxide Fuel Cells with Siloxane Contamination

In this study, the degradation of the solid-oxide fuel cell (SOFC) nickel-yttria stabilized zirconia anode under decamethyltetrasiloxane (L4) contamination is examined with experiments and modeling. A model is developed for the polarization losses based on the charge transfer coefficient,α, and diffusion layer thickness,δ, and fitted to the experimental data to understand how the siloxane degrades the SOFC performance with time. The results of the model indicate that the total polarization losses increase approximately 44% over the course of the 180 min experiment at 350 mA cm −2 . Activation losses dominate the polarization losses initially but decrease in their total contribution while concentration losses increase. Scanning electron microscopy (SEM) with wavelength dispersive X-ray spectroscopy (WDS) elemental mapping indicates that silicon deposition is highest at the outer edge of the anode and forms a barrier layer to fuel diffusion, increasing concentration losses. When the model is applied to other previous D4 and L4 siloxane experiments conducted over a period of 40 h, similar trends in polarization losses are observed. Polarization losses increase more rapidly with D4 compared to L4 siloxane contamination, with concentration losses increasing the fastest with both types of siloxane.

Electrochemistry↗

Ferdinand

A code to read, translate and output parameters of R-matrix fits to nuclear resonance data.

Thompson, IanJ↗

FuelsGen

This toolset generates heterogeneous spatial layouts, aimed at representing forest fuels applicable to wildland fire modeling. Methods are included for calibrating (fitting) this model to data extracted from plot survey information.

Gattiker, James↗

mmbo (Multi-modal Bayesian Optimization) [SWR-25-119]

This package implements Bayesian Neural Network (BNN) based surrogate models for multi-modal data. These models are fit using a custom Variational Bayes procedure leveraging conjugate posteriors in the last layer of the network, offering more accurate predictions with better-calibrated uncertainty quantification than mean-field Variational Bayes.

Taylor, Ian [National Laboratory of the Rockies (N↗

Improved Infrared Optical Constants from Pressed Pellets: II. Ellipsometric n and k Values for Ammonium Sulfate with Variability Analysis

Infrared spectral imaging is facilitated via the comparison of spectra recorded in situ to a databank of actual or synthetic infrared reflectance spectra. It has recently been shown that reference spectra corresponding to the many different morphological forms of the same chemical can be generated synthetically using the imaginary, k, and real, n, components of the complex refractive index, n~ = n + ik. One method to obtain the n and k vectors is infrared ellipsometry, which measures the changes in amplitude, tan ?, and phase, ?, of polarized light reflected from the sample both as a function of wavenumber and angle of incidence. The method requires specularly reflected light so best results are usually obtained with polished planar samples of large surface area. Due to the difficulties of obtaining such samples, however, we investigate the possibility of pressing powders of neat materials and obtaining the corresponding optical constants from the pellets. In this paper, variability in the sample pellet and preparation method is investigated, as is variability in the fitting procedure for the derived optical constants. The n/k vectors are derived from the measured ellipsometric parameters, tan ? and ?, as they are fit by an oscillator model which yield n(?) and k(?) vectors as a function of wavenumber, ?. Construction of the oscillator model is not automatic and depends on significant input from the analyst as well as the sample’s physical characteristics. For the pellet pressing, the experimental variability was found to be minimized for size-selected powdered samples as gauged by the minimal variance in ? and ? for three different pellets; similarly, the analytical precision for multiple measurements of the same pellet was also quite good, suggesting that a pressed pellet is a viable sample preparation method. The greatest variability came in the analytic fitting procedure with differences up to 20% arising from four different fits to the same data set. The final ellipsometric n/k data are compared to literature values obtained from crystalline (NH4)2SO4 samples as well as single-angle reflectance measurements that also used pressed pellets. Comparison with the previous literature values shows generally good agreement although higher k-values are observed for the independent sets of data derived from pressed pellets. These data are suggested as an improved set of optical constants for (NH4)2SO4.

optical constants, ellipsometry, ammonium sulfate,↗

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

60 APPLIED LIFE SCIENCES↗

High-pressure phase transition of olivine-type Mg 2 GeO 4 to a metastable forsterite-III type structure and their equations of state

Germanates are often used as structural analogs of planetary silicates. We have explored the high-pressure phase relations in Mg 2 GeO 4 using diamond-anvil cell experiments combined with synchrotron X-ray diffraction and computations based on density functional theory. Upon room temperature compression, forsterite-type Mg 2 GeO 4 remains stable up to 30 GPa. At higher pressures, a phase transition to a forsterite-III type (Cmc2 1 ) structure was observed, which remained stable to the peak pressure of 105 GPa. Using a third-order Birch Murnaghan fit to the experimental data, we obtained V 0 = 305.1(3) Å3, K 0 = 124.6(14) GPa, and $K'_0$ = 3.86 (fixed) for forsterite-type Mg 2 GeO 4 and V 0 = 263.5(15) Å 3 , K 0 = 175(7) GPa, and $K'_0$ = 4.2 (fixed) for the forsterite-III type phase. The forsterite-III type structure was found to be metastable when compared to the stable assemblage of perovskite/post-perovskite + MgO, as observed during laser-heating experiments. Understanding the phase relations and physical properties of metastable phases is crucial for studying the mineralogy of impact sites, understanding metastable wedges in subducting slabs, and interpreting the results of shock compression experiments.

58 GEOSCIENCES↗