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

Measurement of the Higgs boson production in association with top quarks in multilepton final states in pp collisions at s=13 TeV with the ATLAS detector

A measurement of the associated production of a top-quark pair with the Higgs boson (tt¯H$$ t\overline{t}H $$) in multilepton final states is presented. The analysis is based on a data sample of proton-proton collisions at s=13$$ \sqrt{s}=13 $$ TeV recorded with the ATLAS detector at the CERN Large Hadron Collider and corresponding to an integrated luminosity of 140 fb−1. Six final states defined by the number and flavour of reconstructed charged leptons are combined in a simultaneous likelihood fit to extract the tt¯H$$ t\overline{t}H $$ signal and constrain the most relevant backgrounds. The measured tt¯H$$ t\overline{t}H $$ cross-section normalised to Standard Model (SM) prediction is σtt¯H/σSM=0.63−0.19+0.20$$ {\sigma}_{t\overline{t}H}/{\sigma}^{\mathrm{SM}}=0.{63}_{-0.19}^{+0.20} $$. This result corresponds to an observed (expected) significance of 3.3σ (5.3σ). Additionally, two other fits are used to measure the tt¯H$$ t\overline{t}H $$ cross-section differentially in bins of the Higgs boson transverse momentum in the simplified template cross-section framework, and to extract the associated production cross-section of a single top-quark with the Higgs boson (tH) together with the tt¯H$$ t\overline{t}H $$ one. The CP structure of the top quark-Higgs boson Yukawa coupling is probed through an analysis of tt¯H$$ t\overline{t}H $$ and tH events. The results are compatible with the SM hypothesis, and values of the mixing angle between CP-even and CP-odd top-Higgs Yukawa couplings of |α| > 62° are excluded at 68% confidence level.

Aad, G↗
DOE OSTI2026

QCD+QED PDF implications for the Higgs sector

In this work, we examine the implications of electroweak corrections beyond leading order for processes of special interest in the Higgs sector. We especially explore the role of these corrections given the introduction of an explicit parton distribution function (PDF) for the photon in the proton, an object which emerges necessarily in global PDF fits which include QED effects (i.e., ‘QCD+QED PDFs’). We concentrate on several representative cases, including total Higgs-production cross sections through gluon fusion, gg → H, vector-boson fusion (VBFH), and associated production, pp → V H; we also examine differential distributions, taking a representative Higgs-strahlung process, pp → W + H. We find that the recently developed LUX formalism for the photon PDF significantly stabilizes the PDF dependence of both QED-PDF and electroweak corrections in the Higgs sector, while leaving overall ~ 3–4% cross-section-level variations, depending on the chosen QCD+QED PDF. We illustrate this QCD+QED PDF dependence by exploring predictions based upon recent analyses of the CTEQ-TEA, MSHT, and NNPDF analysis groups, fitted either at NNLO or approximate N3LO in QCD.

Higgs production↗
DOE OSTI2021

The benefits of diligence: how precise are predicted gravitational wave spectra in models with phase transitions?

Models of particle physics that feature phase transitions typically provide predictions for stochastic gravitational wave signals at future detectors and such predictions are used to delineate portions of the model parameter space that can be constrained. The question is: how precise are such predictions? Uncertainties enter in the calculation of the macroscopic thermal parameters and the dynamics of the phase transition itself. We calculate such uncertainties with increasing levels of sophistication in treating the phase transition dynamics. Currently, the highest level of diligence corresponds to careful treatments of the source lifetime; mean bubble separation; going beyond the bag model approximation in solving the hydrodynamics equations and explicitly calculating the fraction of energy in the fluid from these equations rather than using a fit; and including fits for the energy lost to vorticity modes and reheating effects. The lowest level of diligence incorporates none of these effects. We compute the percolation and nucleation temperatures, the mean bubble separation, the fluid velocity, and ultimately the gravitational wave spectrum corresponding to the level of highest diligence for three explicit examples: SMEFT, a dark sector Higgs model, and the real singlet-extended Standard Model (xSM). In each model, we contrast different levels of diligence in the calculation and find that the difference in the final predicted signal can be several orders of magnitude. Our results indicate that calculating the gravitational wave spectrum for particle physics models and deducing precise constraints on the parameter space of such models continues to remain very much a work in progress and warrants care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On the decaying-sterile-neutrino solution to the electron (anti)neutrino appearance anomalies

We explore the hypothesis that the unexplained data from LSND and MiniBooNE are evidence for a new, heavy neutrino mass-eigenstate that mixes with the muon-type neutrino and decays into an electron-type neutrino and a new, very light scalar particle. We consider two different decay scenarios, one with Majorana neutrinos, one with Dirac neutrinos; both fit the data equally well. We find a reasonable, albeit not excellent, fit to the data of MiniBooNE and LSND. The decaying-sterile-neutrino hypothesis, however, cleanly evades constraints from disappearance searches and precision measurements of leptonic meson decays, as long as $1~{\rm MeV}\gtrsim m_4\gtrsim 10$~keV. The SBN program at Fermilab should be able to definitively test the decaying-sterile-neutrino hypothesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Global QCD analysis and dark photons

We perform a global QCD analysis of high energy scattering data within the JAM Monte Carlo framework, including a coupling to a dark photon that augments the Standard Model (SM) electroweak coupling via kinetic mixing with the hypercharge B boson. We first set limits on the dark photon mass and mixing parameter assuming that the SM is the true theory of Nature, taking into account also the effect on g – 2 of the muon. If instead we entertain the possibility that the dark photon may play a role in deep-inelastic scattering (DIS), we find that the best fit is preferred over the SM at 6.5σ, even after accounting for missing higher order uncertainties. The improvement in χ 2 with the dark photon is stable against all the tests we have applied, with the improvements in the theoretical predictions spread across a wide range of x and Q 2 . The largest improvement corresponds to the fixed target and HERA DIS data, while the best fit yields a value of g – 2 which significantly reduces the disagreement with the latest experimental determination.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A low-energy perspective on the minimal left-right symmetric model

We perform a global analysis of the low-energy phenomenology of the minimal left-right symmetric model (mLRSM) with parity symmetry. We match the mLRSM to the Standard Model Effective Field Theory Lagrangian at the left-right-symmetry breaking scale and perform a comprehensive fit to low-energy data including mesonic, neutron, and nuclear β -decay processes, Δ F = 1 and Δ F = 2 CP-even and -odd processes in the bottom and strange sectors, and electric dipole moments (EDMs) of nucleons, nuclei, and atoms. We fit the Cabibbo-Kobayashi-Maskawa and mLRSM parameters simultaneously and determine a lower bound on the mass of the right-handed W R boson. In models where a Peccei-Quinn mechanism provides a solution to the strong CP problem, we obtain $ {M}_{W_R} $ ≳ 5 . 5 TeV at 95% C.L. which can be significantly improved with next-generation EDM experiments. In the P -symmetric mLRSM without a Peccei-Quinn mechanism we obtain a more stringent constraint $ {M}_{W_R} $ ≳ 17 TeV at 95% C.L., which is difficult to improve with low-energy measurements alone. In all cases, the additional scalar fields of the mLRSM are required to be a few times heavier than the right-handed gauge bosons. We consider a recent discrepancy in tests of first-row unitarity of the CKM matrix. We find that, while TeV-scale W R bosons can alleviate some of the tension found in the V ud,us determinations, a solution to the discrepancy is disfavored when taking into account other low-energy observables within the mLRSM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A flavor of SO(10) unification with a spinor Higgs

We investigate Higgs Parity unification — a realization of SO(10) grand unification based on the Higgs Parity mechanism in which the Standard Model (SM) Higgs resides in a spinor representation. The theory has an intermediate left-right symmetric stage where the SU(2)R symmetry breaking scale is fixed by the vanishing of the SM Higgs quartic coupling. The strong CP problem is solved by parity. Gauge coupling unification successfully predicts αs(MZ) to within 1%. The spinor Higgs naturally leads to a seesaw origin for SM flavor observables. We identify a novel mechanism where large mixing of third generation fermions with additional heavy vector-like fermions accounts for the anarchical nature of the PMNS matrix and the lack of hierarchy in the neutrino mass spectrum, relative to the up-quarks. A fit to quark and lepton masses and mixings, with a minimal parameter set, predicts 1) A testable relation between the top quark mass and αs(MZ) which is about (1 – 2)σ from current best fit values, 2) The order of magnitude of the baryon asymmetry of the universe, via leptogenesis from second-generation right-handed neutrino decays. 3) The proton decay and the neutron EDM are likely observable in next generation experiments, and 4) A normal ordered neutrino mass spectrum where 0νββ decay and the mass of the lightest neutrino are out of reach of next generation experiments.

Baryo-and Leptogenesis↗
DOE OSTI2024

Measurement of the effective leptonic weak mixing angle

Using pp collision data at $\sqrt{s}$ = 13 TeV, recorded by the LHCb experiment between 2016 and 2018 and corresponding to an integrated luminosity of 5.4 fb −1 , the forward-backward asymmetry in the pp → Z/γ* → μ + μ − process is measured. The measurement is carried out in ten intervals of the difference between the muon pseudorapidities, within a fiducial region covering dimuon masses between 66 and 116 GeV, muon pseudorapidities between 2.0 and 4.5 and muon transverse momenta above 20 GeV. These forward-backward asymmetries are compared with predictions, at next-to-leading order in the strong and electroweak couplings. The measured effective leptonic weak mixing angle is ${\sin}^2{\theta}_{\textrm{eff}}^{\ell }=0.23147\pm 0.00044\pm 0.00005\pm 0.00023,$ where the first uncertainty is statistical, the second arises from systematic uncertainties associated with the asymmetry measurement, and the third arises from uncertainties in the fit model used to extract ${\sin}^2{\theta}_{\textrm{eff}}^{\ell }$ from the asymmetry measurement. This result is based on an arithmetic average of results using the CT18, MSHT20, and NNPDF31 parameterisations of the proton internal structure, and is consistent with previous measurements and with predictions from the global electroweak fit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A comparison of explosively driven shock wave radius versus time scaling approaches

Abstract Explosively driven shock wave radius versus time profiles are frequently used to document energy release and relative explosive performance. Recently, two universal shock wave radius versus time profiles have been presented in the literature, which demonstrate the ability to represent explosively driven shock wave profiles for all explosive sources in any fluid environment. These two universal shock wave profiles are examined here relative to each other and relative to a commonly used nonlinear shock wave profile, which is fit to experimental data for individual explosive materials. The nonlinear profile, originally developed by Dewey, is examined here, and a universal non-dimensional form of the equation is proposed. The universal shock wave profiles are all found to be relatively similar, but with slight variations in a transition region of non-dimensional radii $$0.15\lesssim R^*\lesssim 2$$ 0.15 ≲ R ∗ ≲ 2 . The variations in this region result in different estimations of energy release or blast strength between the curve fits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Continuous and discrete modeling of HIV-1 decline on therapy

Mathematical models have shed light on the dynamics of HIV- 1 infection in vivo. In this paper, we generalize continuous mathematical models of drug therapy for HIV-1 by Perelson et al. (Science 271:1582–1586, 1996) and Perelson and Nelson (SIAM Rev 41:3–44, 1999) on time scales, i.e., a nonempty closed subset of real numbers in order to derive new discrete models that predict the total concentration of plasma virus as a function of time. Here, one of our main goals is to compare discrete mathematical models with the continuous model in Perelson et al. (1996) where HIV infected patients were given protease inhibitors and sampled frequently thereafter. For the comparison, we use experimental data collected in Perelson et al. (1996) and estimate the parameters such as the virion clearance rate and the rate of loss of infected cells by fitting the total concentration of plasma virus to this data set. Our results show that discrete systems describe the best fit. In the previous models of this study, the efficacy of protease inhibitor is assumed to be perfect. Motivated by Perelson and Nelson (1999), we end the paper with a mathematical model of imperfect protease inhibitor and reverse transcriptase (RT) inhibitor combination therapy of HIV-1 infection on time scales with its stability analysis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Yb:Lu 2 O 3 hydrothermally grown single-crystal high-resolution absorption spectra obtained between 8 and 300 K

Here we present new high-resolution absorption data for the important sesquioxide laser material Yb:Lu 2 O 3 for the spectral range of 880–1020 nm, at various temperatures between 8 and 300 K, and for the zero-phonon region from 960 to 990 nm, at temperatures from 8 to 300 K. We have experimentally observed the C 3i (0,1)–(1,3) transition for the first time, located at 880.7 nm at 8 K. Based on high confidence fitting functions to the experimental data, we provide the first complete compilation of all observed electronic and electronic–vibrational transitions. Detailed fitting and plots of the C 2 and C 3i zero-line data show an evolution of the linewidth from being predominantly electronic below about 100 K to being dominated by thermal processes above 100 K. We have also found evidence for a “soft” phase transition between 80 and 100 K that changes the local coordination environment.

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

Reynolds stress tensor measurements using magnetic resonance velocimetry: expansion of the dynamic measurement range and analysis of systematic measurement errors

This study presents magnetic resonance velocimetry (MRV) Reynolds Stress measurements in a periodic hill channel with a hill Reynolds number of Re = 29,500. The velocity encoding scheme is based on the ICOSA6 method with six icosahedral encoding directions and multiple encoding values are measured to increase the dynamic range. The full Reynolds stress tensor is obtained from a voxel-wise three-dimensional Gaussian fit using the magnitude data of all acquisitions. The MRV results are compared to a wall-resolved large eddy simulation and laser Doppler velocimetry measurements conducted in the same channel. It is shown that the MRV Reynolds stress data have excellent precision and agree qualitatively with the reference data. However, there are apparent systematic deviations. One of the most prominent error contributions is the signal attenuation caused by higher orders of motion, which leads to an overestimation of the turbulence level. Another fundamental error is identified in the assumption that the turbulence is Gaussian distributed. With the presented reconstruction technique, the MRV data are fitted to a statistical model, and depending on the examined flow setup, the Gaussian model can lead to considerable errors. Possible ways of how to reduce all identified errors are presented. In summary, this technique enables Reynolds stress tensor measurements in complex internal flows with high dynamic range and excellent precision. However, several issues need to be resolved to make the turbulence quantification more accurate.

42 ENGINEERING↗
DOE OSTI2020

Calibration of hyperelastic constitutive models: the role of boundary conditions, search algorithms, and experimental variability

The calibration of hyperelastic constitutive models of soft tissue and tissue surrogates is often treated as an exercise in curve-fitting to the average experimental response, and many of the complicating factors such as experimental boundary conditions and data variability are ignored. Here, we focus on three questions that arise in this area: the ramifications of ignoring the experimental boundary conditions, the use of local optimizers, and the role of data variability. Using data from a uniaxial extension experiment on a tissue surrogate, we study how these three factors affect the calibration of isotropic hyperelastic constitutive models. Our results show that even with the simplest of constitutive models, it is necessary to look beyond a "good fit" to the average.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A rapid heating and high magnetic field thermal analysis technique

Abstract A new thermal analysis technique is described that allows measurements to be performed on bulk samples at extreme heating and cooling rates and in high magnetic fields. High heating rates, up to 1000 °C min −1 , are achieved through electromagnetic induction heating of a custom-built apparatus fitted with commercial thermal analysis heads and sensor. Rapid cooling rates, up to 100 °C min −1 , are enabled by gas quenching and the small thermal mass of the induction furnace. The custom apparatus is designed to fit inside a superconducting magnet capable of fields up to 9 Tesla. This study demonstrates that the instrument is capable of collecting accurate thermal analysis data in high magnetic fields and rapidly acquiring data for dynamic processes. While the full potential of the technique is still unrealized, currently, it can provide insight into phenomena at time scales relevant to heat treatment in many industrial processes and into little understood effects of high magnetic field processing.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

On the inference of viscoelastic constants from stress relaxation experiments

Several constitutive theories have been proposed in the literature to model the viscoelastic response of materials, including widely used rheological constitutive models. These models are characterized by certain parameters ("time constants") that define the time scales over which the material relaxes. These parameters are primarily obtained from stress relaxation experiments using curve-fitting techniques. However, the question of how best to estimate these time constants remains open. As a step towards answering this question, we propose an optimal experimental design approach based on ideas from information geometry, namely Fisher information and Kullback-Leibler divergence. The material is modeled as a spring element in parallel with multiple Maxwell elements and described using a one- or two-term Prony series. Treating the time constants as unknowns, we develop expressions for the Fisher information and Kullback-Leibler divergence that allow us to maximize information gain from experimental data. Based on the results of this study, we propose that the largest time constant estimated from a stress relaxation experiment for a linear viscoelastic material should be at most one-fifth of the total time of the experiment in order to maximize information gain. Our results also provide confirmation that the equilibrium modulus of the material cannot be reliably determined from curve-fitting to data from a stress relaxation experiment.

Fisher information↗
DOE OSTI2022

A phase transition for finding needles in nonlinear haystacks with LASSO artificial neural networks

To fit sparse linear associations, a LASSO sparsity inducing penalty with a single hyperparameter provably allows to recover the important features (needles) with high probability in certain regimes even if the sample size is smaller than the dimension of the input vector (haystack). More recently learners known as artificial neural networks (ANN) have shown great successes in many machine learning tasks, in particular fitting nonlinear associations. Small learning rate, stochastic gradient descent algorithm and large training set help to cope with the explosion in the number of parameters present in deep neural networks. Yet few ANN learners have been developed and studied to find needles in nonlinear haystacks. Driven by a single hyperparameter, our ANN learner, like for sparse linear associations, exhibits a phase transition in the probability of retrieving the needles, which we do not observe with other ANN learners. To select our penalty parameter, we generalize the universal threshold of Donoho and Johnstone (Biometrika 81(3):425–455, 1994) which is a better rule than the conservative (too many false detections) and expensive cross-validation. In the spirit of simulated annealing, we propose a warm-start sparsity inducing algorithm to solve the high-dimensional, non-convex and non-differentiable optimization problem. We perform simulated and real data Monte Carlo experiments to quantify the effectiveness of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗