Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “likelihood function”

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 OSTI2022

Dark Energy Survey Year 3 results: A 2.7% measurement of baryon acoustic oscillation distance scale at redshift 0.835

We present angular diameter measurements obtained by measuring the position of baryon acoustic oscillations (BAO) in an optimized sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg 2 with 0.6 < z photo < 1.1 and a typical redshift uncertainty of 0.03 ( 1 + z ) . The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift, and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive preunblinding protocol. We used two different statistics; angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance D M at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain D M ( z eff = 0.835 ) / r d = 18.92 ± 0.51 , which is consistent with, but smaller than, the Planck prediction assuming flat Λ CDM , at the level of 2.3 σ . The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at z > 0.75 . When combined with DES 3 x 2 pt + SNIa , they lead to improvements in H 0 and Ω m constraints by ∼ 20 % .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Importance of Convective Stage on Out-of-Cloud Convectively Induced Turbulence from High-Resolution Simulations

Turbulence (clear-air, mountain wave, convectively induced) is an aviation hazard that is a challenge to forecast due to the coarse resolution ultilized in operational weather models. Turbulence indices are commonly used to aid pilots in avoiding turbulence, but these indices have been designed and calibrated for midlatitude clear-air turbulence prediction (e.g., the Ellrod index). A significant limitation with current convectively induced turbulence (CIT) prediction is the lack of storm stage dependency. In this study, six high-resolution simulations of tropical oceanic and midlatitude continental convection are performed to characterize the turbulent environment near various convective types during the developing and mature stages. Second-order structure functions, a diagnostic commonly used to identify turbulence in turbulence prediction systems, are used to characterize the probability of turbulence for various convective types. Turbulence likelihood was found to be independent of region (i.e., tropical vs midlatitude) but dependent on convective stage. The probability of turbulence increased near developing convection for the majority of cases. Additional analysis of static stability and vertical wind shear, indicators of turbulence potential, showed that the convective environment near developing convection was more favorable for turbulence production than mature convection. Near developing convection, static stability decreased and vertical wind shear increased. Vertical wind shear near mature and developing convection was found to be weakly correlated to turbulence intensity in both the tropics and the midlatitudes. Finally, this study emphasizes the need for turbulence avoidance guidelines for the aviation community that are dependent on convective stage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Factors linked to participant attrition in a longitudinal occupational health surveillance program

For occupational medical screening programs focused on long-term health surveillance, participant attrition is a significant barrier to success. Here we investigate demographic, medical history, and clinical data from National Supplemental Screening Program (NSSP) examinees for association with likelihood of return for a second exam (rescreening). A total of 15,733 individuals completed at least one NSSP exam before December 31, 2016; of those, 4832 also completed a second exam on or before December 31, 2019. Stepwise logistic regression models were used to identify variables associated with whether a participant was rescreened in the NSSP. Individuals were less likely to return for rescreening if they had a history of any cancer; cardiovascular problems; diabetes or kidney disease; or if they used insulin. Age at time of first exam and job site category significantly influenced likelihood of return. Workers categorized as “guests” were more likely to return. Participants were less likely to return if they had an abnormal urinalysis, abnormal pulmonary function, pneumoconiosis, aortic atherosclerosis, or hearing loss at their initial exam. Participants who received a chest X-ray at their initial screening were more likely to return. The presence of health problems is strongly linked to screening program attrition. Participants who are older at the time of their initial screening exam are less likely to return. The discovery of several strong demographic, medical, and job associations reveals the importance for medical screening programs to understand and address factors that influence participant retention and, consequently, the effectiveness of long-term health surveillance activities.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Output-weighted sampling for multi-armed bandits with extreme payoffs

We present a new type of acquisition function for online decision-making in multi-armed and contextual bandit problems with extreme payoffs. Specifically, we model the payoff function as a Gaussian process and formulate a novel type of upper confidence bound acquisition function that guides exploration towards the bandits that are deemed most relevant according to the variability of the observed rewards. This is achieved by computing a tractable likelihood ratio that quantifies the importance of the output relative to the inputs and essentially acts as an attention mechanism that promotes exploration of extreme rewards. Our formulation is supported by asymptotic zero-regret guarantees, and its performance is demonstrated across several synthetic benchmarks, as well as two realistic examples involving noisy sensor network data. Finally, we provide a JAX library for efficient bandit optimization using Gaussian processes.

Science & Technology - Other Topics↗
DOE OSTI2026

BICEP/Keck XXI: Constraints on early-Universe parity violation from multipole-dependent birefringence

We present the first constraints on multipole-dependent cosmic birefringence using CMB polarization data from the BK18 dataset, which combines observations from BICEP2, Array, and BICEP3 at frequencies of 95, 150, and 220 GHz. Photon coupling to an axionlike field leads to the rotation of CMB polarization, inducing nonzero cross-correlations. We show that a multipole-dependent rotation imprints a distinct signature in the polarization spectra that can be constrained. Specifically, we consider an early dark energy (EDE) scenario in which a pseudoscalar field couples to photons through a Chern-Simons interaction, generating a polarization rotation with multipole dependence. We introduce a phenomenological as a step function, obtaining constraints on the step function size consistent with zero, with uncertainties less than (68% CL). In addition, using multi-frequency , , and cross spectra, along with robust BICEP/ foreground treatment and likelihood framework, we derive constraints on the axion-photon coupling amplitude for several choices of EDE parameters. For the baseline best-fit value from the 2018 analysis, we obtain (68% CL), consistent with previous limits.

Ade, P A R↗
DOE OSTI2024

Radiation image reconstruction and uncertainty quantification using a Gaussian process prior

We propose a complete framework for Bayesian image reconstruction and uncertainty quantification based on a Gaussian process prior (GPP) to overcome limitations of maximum likelihood expectation maximization (ML-EM) image reconstruction algorithm. The prior distribution is constructed with a zero-mean Gaussian process (GP) with a choice of a covariance function, and a link function is used to map the Gaussian process to an image. Unlike many other maximum a posteriori approaches, our method offers highly interpretable hyperparamters that are selected automatically with the empirical Bayes method. Furthermore, the GP covariance function can be modified to incorporate a priori structural priors, enabling multi-modality imaging or contextual data fusion. Lastly, we illustrate that our approach lends itself to Bayesian uncertainty quantification techniques, such as the preconditioned Crank–Nicolson method and the Laplace approximation. The proposed framework is general and can be employed in most radiation image reconstruction problems, and we demonstrate it with simulated free-moving single detector radiation source imaging scenarios. We compare the reconstruction results from GPP and ML-EM, and show that the proposed method can significantly improve the image quality over ML-EM, all the while providing greater understanding of the source distribution via the uncertainty quantification capability. Furthermore, significant improvement of the image quality by incorporating a structural prior is illustrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A knock-on deuteron imager for measurements of fuel and hotspot asymmetry in direct-drive inertial confinement fusion implosions (invited)

A knock-on deuteron imager (KoDI) has been implemented to measure the fuel and hotspot asymmetry of cryogenic inertial confinement fusion implosions on OMEGA. Energetic neutrons produced by D–T fusion elastically scatter (“knock on”) deuterons from the fuel layer with a probability that depends on ρR. Deuterons above 10 MeV are produced by near-forward scattering, and imaging them is equivalent to time-integrated neutron imaging of the hotspot. Deuterons below 6 MeV are produced by a combination of side scattering and ranging in the fuel, and encode information about the spatial distribution of the dense fuel. Here, the KoDI instrument consists of a multi-penumbral aperture positioned 10–20 cm from the implosion using a ten-inch manipulator and a detector pack at 350 cm from the implosion to record penumbral images with magnification of up to 35×. Range filters and the intrinsic properties of CR-39 are used to distinguish different charged-particle images by energy along the same line of sight. Image plates fielded behind the CR-39 record a 10 keV x-ray image using the same aperture. A maximum-likelihood reconstruction algorithm has been implemented to infer the source from the projected penumbral images. The effects of scattering and aperture charging on the instrument point-spread function are assessed. Synthetic data are used to validate the reconstruction algorithm and assess an appropriate termination criterion. Significant aperture charging has been observed in the initial experimental dataset, and increases with aperture distance from the implosion, consistent with a simple model of charging by laser-driven EMP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energetics of silicon in the bulk and near surfaces of tungsten: a first-principles study

Abstract Siliconization of the tokamak walls is a candidate method to improve plasma confinement in fusion tokamaks containing tungsten plasma facing components (W PFCs). To understand the interactions of silicon (Si) with W, the Si behavior in bulk W, and near three low-index W surfaces ((100), (110) and (111)) has been investigated using first-principles density functional theory. In bulk W, Si interstitial atoms have a low solution ability and high mobility, and Si atoms can be strongly trapped by W vacancies. The interaction between two Si adatoms is responsible for the stability of adatom superstructures on W surfaces, consistent with previous experimental observation (Tsong and Casanova 1981 Phys. Rev. Lett. 47 113). Although the coverage dependence of Si adsorption and diffusion energetics on surfaces is related to surface orientation, the W(110) surface has lower Si adsorption affinity and higher Si diffusivity than either the W(111) or W(100) surfaces. The most stable Si adatom superstructure on W surfaces is: square c(2 × 2) pattern on W(100) covered with 0.5 ML Si; rectangular c(4 × 2) pattern on W(110) with 0.25 ML Si; and rhombus p(1 × 1) pattern on W(111) with 1 ML Si. The coverage dependence of Si mobility on/toward W surfaces is generally related to the stability of the Si superstructures as a function of coverage on each surface. Interestingly, Si adatoms prefer to transport below the surface and into W subsurface by an exchange mechanism with W atoms, indicating the likelihood of epitaxial growth of W silicide layers on W surfaces during the operation of W PFCs.

Yang, L. (ORCID:0000000322166071)↗
DOE OSTI2024

Evolving Competitive Markets in SAPP: Leveraging Competitive Wholesale Electricity Markets to Drive Renewable Generation Capacity in the Southern African Power Pool (SAPP)

The SADC region has significant natural resource potential to increase renewable energy generation, improve electricity reliability, and support economic development. This research finds an apparent lack of confidence from electricity infrastructure investors in SAPP wholesale electricity markets, which increases risk perception and lowers the likelihood of capital deployment. With respect to free market fundamentals, competitive market obstacles and renewable energy development obstacles are characterized. Stakeholders identified the top obstacles to well-functioning competitive markets as insufficient transmission infrastructure for interconnection and regional movement of electricity, dominance of national single-buyer markets, and lack of or weak nation-state regulatory frameworks. Stakeholders prioritized the top three obstacles for renewable energy development as a lack of viable commercial arrangements for variable renewable energy (VRE) balancing, lack of functional and consistent nation-level regulations, and higher project costs related to reliance on imported equipment. With respect to potential solution options, stakeholders prioritized the development of new cost allocation and finance methods to facilitate new transmission expansion, training to educate new or potential new market entrants on SAPP processes, as well as modeling and analysis of regional SAPP participation benefits disaggregated to the nation-state level. From these perspectives, this research identified strategy options for consideration including transitioning SAPP to a regional transmission operator (RTO) for operation and planning of cross-border transmission facilities and market administration, shifting operations of SAPP member transmission systems to Independent System Operators (ISOs), establishing a regional regulatory authority and enhancing market data transparency. Implementing these reforms is expected to be challenging, but not insurmountable, given the domestic political, legal, and jurisdictional complexities of the SADC region.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Interpreting High-resolution Spectroscopy of Exoplanets using Cross-correlations and Supervised Machine Learning

We present a new method for performing atmospheric retrieval on ground-based, high-resolution data of exoplanets. Our method combines cross-correlation functions with a random forest, a supervised machine-learning technique, to overcome challenges associated with high-resolution data. A series of cross-correlation functions are concatenated to give a “CCF-sequence” for each model atmosphere, which reduces the dimensionality by a factor of ∼100. The random forest, trained on our grid of ∼65,000 models, provides a likelihood-free method of retrieval. The precomputed grid spans 31 values of both temperature and metallicity, and incorporates a realistic noise model. We apply our method to HARPS-N observations of the ultra-hot Jupiter KELT-9b and obtain a metallicity consistent with solar (logM = − 0.2 ± 0.2). Our retrieved transit chord temperature (T=6000{sub −200}{sup +0}K) is unreliable as strong ion lines lie outside of the extent of the training set, which we interpret as being indicative of missing physics in our atmospheric model. We compare our method to traditional nested sampling, as well as other machine-learning techniques, such as Bayesian neural networks. We demonstrate that the likelihood-free aspect of the random forest makes it more robust than nested sampling to different error distributions, and that the Bayesian neural network we tested is unable to reproduce complex posteriors. We also address the claim in Cobb et al. 2019 that our random forest retrieval technique can be overconfident but incorrect. We show that this is an artifact of the training set, rather than of the machine-learning method, and that the posteriors agree with those obtained using nested sampling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Pulse-shape discrimination against low-energy Ar-39 beta decays in liquid argon with 4.5 tonne-years of DEAP-3600 data

The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from 39 Ar beta decays and is suppressed using pulse-shape discrimination (PSD). We use two types of PSD estimator: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the mean single-photoelectron charge, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulse shape and for afterpulsing in the light detectors. The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

We present the XFaster analysis package, XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-C ℓ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data (Netterfield et al. 2002; Jones et al. 2006), and also used as part of the Planck analysis (Rocha et al. 2011). Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument (SPIDER Collaboration 2021). The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Measurement of the $c$-jet mistagging efficiency in $t\bar{t}$ events using $pp$ collision data at $\sqrt{s}=13$ $\text {TeV}$ collected with the ATLAS detector

A technique is presented to measure the efficiency with which c-jets are mistagged as b-jets (mistagging efficiency) using tt¯ events, where one of the W bosons decays into an electron or muon and a neutrino and the other decays into a quark–antiquark pair. The measurement utilises the relatively large and known W → cs branching ratio, which allows a measurement to be made in an inclusive c-jet sample. The data sample used was collected by the ATLAS detector at √s = 13 TeV and corresponds to an integrated luminosity of 139 fb –1 . Events are reconstructed using a kinematic likelihood technique which selects the mapping between jets and tt¯ decay products that yields the highest likelihood value. The distribution of the b-tagging discriminant for jets from the hadronic W decays in data is compared with that in simulation to extract the mistagging efficiency as a function of jet transverse momentum. The total uncertainties are in the range 3–17%. The measurements generally agree with those in simulation but there are some differences in the region corresponding to the most stringent b-jet tagging requirement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Apprentice for Event Generator Tuning

Apprentice is a tool developed for event generator tuning. It contains a range of conceptual improvements and extensions over the tuning tool Professor. Its core functionality remains the construction of a multivariate analytic surrogate model to computationally expensive Monte-Carlo event generator predictions. The surrogate model is used for numerical optimization in chi-square minimization and likelihood evaluation. Apprentice also introduces algorithms to automate the selection of observable weights to minimize the effect of mis-modeling in the event generators. We illustrate our improvements for the task of MC-generator tuning and limit setting.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Radiative properties of quantum emitters in boron nitride from excited state calculations and Bayesian analysis

Abstract Point defects in hexagonal boron nitride (hBN) have attracted growing attention as bright single-photon emitters. However, understanding of their atomic structure and radiative properties remains incomplete. Here we study the excited states and radiative lifetimes of over 20 native defects and carbon or oxygen impurities in hBN using ab initio density functional theory and GW plus Bethe-Salpeter equation calculations, generating a large data set of their emission energy, polarization and lifetime. We find a wide variability across quantum emitters, with exciton energies ranging from 0.3 to 4 eV and radiative lifetimes from ns to ms for different defect structures. Through a Bayesian statistical analysis, we identify various high-likelihood charge-neutral defect emitters, among which the native V N N B defect is predicted to possess emission energy and radiative lifetime in agreement with experiments. Our work advances the microscopic understanding of hBN single-photon emitters and introduces a computational framework to characterize and identify quantum emitters in 2D materials.

Chemistry↗
DOE OSTI2021

Measurement of the inclusive and differential Higgs boson production cross sections in the leptonic WW decay mode at $\sqrt{s} =$ 13 TeV

Measurements of the fiducial inclusive and differential production cross sections of the Higgs boson in proton-proton collisions at $ \sqrt{s} $ = 13 TeV are performed using events where the Higgs boson decays into a pair of W bosons that subsequently decay into a final state with an electron, a muon, and a pair of neutrinos. The analysis is based on data collected with the CMS detector at the LHC during 2016–2018, corresponding to an integrated luminosity of 137 fb$^{−1}$. Production cross sections are measured as a function of the transverse momentum of the Higgs boson and the associated jet multiplicity. The Higgs boson signal is extracted and simultaneously unfolded to correct for selection efficiency and resolution effects using maximum-likelihood fits to the observed distributions in data. The integrated fiducial cross section is measured to be 86.5 ± 9.5 fb, consistent with the Standard Model expectation of 82.5 ± 4.2 fb. No significant deviation from the Standard Model expectations is observed in the differential measurements.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

LTAU (Loss Trajectory Analysis for Uncertainty) is a technique for estimating the uncertainty in a model's predictions for any regression task by approximating the cumulative distribution function (CDF) of the model's errors over the course of training. The approximated CDF, combined with a similarity search algorithm in the model's descriptor space, can then be used to estimate the likelihood that the model's predictions on any given test point will fall below a chosen threshold. LTAU-FF (LTAU in atomistic Force Fields) is the application of LTAU specifically for use with atomistic force fields.

Vita, JoshuaA↗