Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maximum Likelihood”

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 OSTI2024

Joint inference of multiplicative and additive systematics in galaxy density fluctuations and clustering measurements

Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical uncertainties get smaller, correcting for these systematic contaminants becomes increasingly important for unbiased cosmological analysis. We present and validate a new method for understanding and mitigating both additive and multiplicative systematics in galaxy clustering measurements (two-point function) by joint inference of contaminants in the galaxy overdensity field (one-point function) using a maximum-likelihood estimator (MLE). We test this methodology with Kilo-Degree Survey-like mock galaxy catalogues and synthetic systematic template maps. We estimate the cosmological impact of such mitigation by quantifying uncertainties and possible biases in the inferred relationship between the observed and the true galaxy clustering signal. Our method robustly corrects the clustering signal to the sub-percent level and reduces numerous additive and multiplicative systematics from 1.5σ to less than 0.1σ for the scenarios we tested. In addition, we provide an empirical approach to identifying the functional form (additive, multiplicative, or other) by which specific systematics contaminate the galaxy number density. Even though this approach is tested and geared towards systematics contaminating the galaxy number density, the methods can be extended to systematics mitigation for other two-point correlation measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗
DOE OSTI2024

Optimal 1D Ly α forest power spectrum estimation – III. DESI early data

ABSTRACT The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and discuss the key sources of systematic error that we need to address with the upcoming DESI Year 1 analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DESI peculiar velocity survey – Fundamental Plane

The Dark Energy Spectroscopic Instrument (DESI) peculiar velocity survey aims to measure the peculiar velocities of early- and late-type galaxies within the DESI footprint using both the Fundamental Plane and optical Tully–Fisher relations. Direct measurements of peculiar velocities can significantly improve constraints on the growth rate of structure, reducing uncertainty by a factor of approximately 2.5 at redshift 0.1 compared to the DESI Bright Galaxy Survey’s redshift space distortion measurements alone. We assess the quality of stellar velocity dispersion measurements from DESI spectroscopic data. These measurements, along with photometric data from the Legacy Survey, establish the Fundamental Plane relation and determine distances and peculiar velocities of early-type galaxies. During survey validation, we obtain spectra for 6698 unique early-type galaxies, up to a photometric redshift of 0.15. 64 per cent of observed galaxies (4267) have relative velocity dispersion errors below 10 per cent. This percentage increases to 75 per cent if we restrict our sample to galaxies with spectroscopic redshifts below 0.1. We use the measured central velocity dispersion, along with photometry from the DESI Legacy Imaging Surveys, to fit the Fundamental Plane parameters using a 3D Gaussian maximum likelihood algorithm that accounts for measurement uncertainties and selection cuts. In addition, we conduct zero-point calibration using the absolute distance measurements to the Coma cluster, leading to a value of the Hubble constant, H 0 = 76.05 ± 0.35 (statistical) ±0.49 (systematic Fundamental Plane) ±4.86 (statistical due to calibration) km s –1 Mpc –1 ⁠. This H 0 value is within 2σ of Planck cosmic microwave background results and within 1σ of other low-redshift distance indicator-based measurements.

cosmological parameters↗
DOE OSTI2022

The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Phylogenetic Diversity and Mycotoxin Potential of Emergent Phytopathogens Within the Fusarium tricinctum Species Complex

Recent studies on multiple continents indicate members of the Fusarium tricinctum species complex (FTSC) are emerging as prevalent pathogens of small-grain cereals, pulses, and other economically important crops. These understudied fusaria produce structurally diverse mycotoxins, among which enniatins (ENNs) and moniliformin (MON) are the most frequent and of greatest concern to food and feed safety. Herein a large survey of fusaria in the Fusarium Research Center and Agricultural Research Service culture collections was undertaken to assess species diversity and mycotoxin potential within the FTSC. A 151-strain collection originating from diverse hosts and substrates from different agroclimatic regions throughout the world was selected from 460 FTSC strains to represent the breadth of FTSC phylogenetic diversity. Evolutionary relationships inferred from a five-locus dataset, using maximum likelihood and parsimony, resolved the 151 strains as 24 phylogenetically distinct species, including nine that are new to science. Of the five genes analyzed, nearly full-length phosphate permease sequences contained the most phylogenetically informative characters, establishing its suitability for species-level phylogenetics within the FTSC. Fifteen of the species produced ENNs, MON, the sphingosine analog 2-amino-14,16-dimethyloctadecan-3-ol (AOD), and the toxic pigment aurofusarin (AUR) on a cracked corn kernel substrate. Interestingly, the five earliest diverging species in the FTSC phylogeny (i.e., F. iranicum, F. flocciferum, F. torulosum, and Fusarium spp. FTSC 8 and 24) failed to produce AOD and MON, but synthesized ENNs and/or AUR. Moreover, our reassessment of nine published phylogenetic studies on the FTSC identified 11 additional novel taxa, suggesting this complex comprises at least 36 species.

Plant Sciences↗
DOE OSTI2023

Seeing through noise in power laws

Despite widespread claims of power laws across the natural and social sciences, evidence in data is often equivocal. Modern data and statistical methods reject even classic power laws such as Pareto’s law of wealth and the Gutenberg–Richter law for earthquake magnitudes. We show that the maximum-likelihood estimators and Kolmogorov–Smirnov (K-S) statistics in widespread use are unexpectedly sensitive to ubiquitous errors in data such as measurement noise, quantization noise, heaping and censorship of small values. This sensitivity causes spurious rejection of power laws and biases parameter estimates even in arbitrarily large samples, which explains inconsistencies between theory and data. We show that logarithmic binning by powers of λ > 1 attenuates these errors in a manner analogous to noise averaging in normal statistics and that λ thereby tunes a trade-off between accuracy and precision in estimation. Binning also removes potentially misleading within-scale information while preserving information about the shape of a distribution over powers of λ, and we show that some amount of binning can improve sensitivity and specificity of K-S tests without any cost, while more extreme binning tunes a trade-off between sensitivity and specificity. We therefore advocate logarithmic binning as a simple essential step in power-law inference.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Probing the PeV region in the astrophysical neutrino spectrum using 𝜈 𝜇 from the Southern sky

IceCube has observed a diffuse astrophysical neutrino flux over the energy region from a few TeV to a few PeV. At PeV energies, the spectral shape is not yet well measured due to the low statistics of the data. This analysis probes the gap between 1 and 10 PeV by using high-energy downgoing muon neutrinos. Here, to reject the large atmospheric muon background, two complementary techniques are combined. The first technique selects events with high stochasticity to reject atmospheric muon bundles whose stochastic energy losses are smoothed due to high muon multiplicity. The second technique vetoes atmospheric muons with the IceTop surface array. Using 9 yrs of data, we found two neutrino candidate events in the signal region, consistent with expectation from background, each with relatively high signal probabilities. A joint maximum likelihood estimation is performed using this sample and an independent 9.5-yr sample of tracks to measure the neutrino spectrum. A likelihood ratio test is done to compare the single power-law (SPL) vs SPL+cutoff hypothesis; the SPL+cutoff model is not significantly better than the SPL. High-energy astrophysical objects from four source catalogs are also checked around the direction of the two events. No significant coincidence was found.

Abbasi, R. [Loyola University Chicago] (ORCID:0000↗
DOE OSTI2025

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗
DOE OSTI2020

Calorimetry for low-energy electrons using charge and light in liquid argon

Precise calorimetric reconstruction of 5–50 MeV electrons in liquid argon time projection chambers (LArTPCs) will enable the study of astrophysical neutrinos in DUNE and could enhance the physics reach of oscillation analyses. Liquid argon scintillation light has the potential to improve energy reconstruction for low-energy electrons over charge-based measurements alone. Here we demonstrate light-augmented calorimetry for low-energy electrons in a single-phase LArTPC using a sample of Michel electrons from decays of stopping cosmic muons in the LArIAT experiment at Fermilab. Michel electron energy spectra are reconstructed using both a traditional charge-based approach as well as a more holistic approach that incorporates both charge and light. A maximum-likelihood fitter, using LArIAT’s well-tuned simulation, is developed for combining these quantities to achieve optimal energy resolution. A sample of isolated electrons is simulated to better determine the energy resolution expected for astrophysical electron-neutrino charged-current interaction final states. In LArIAT, which has very low wire noise and an average light yield of 18 pe/MeV, an energy resolution of σ/E≃9.3%/E⊕1.3% is achieved. Samples are then generated with varying wire noise levels and light yields to gauge the impact of light-augmented calorimetry in larger LArTPCs. At a charge-readout signal-to-noise of S/N≃30, for example, the energy resolution for electrons below 40 MeV is improved by ≈10%, ≈20%, and ≈40% over charge-only calorimetry for average light yields of 10 pe/MeV, 20 pe/MeV, and 100 pe/MeV, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Dark matter trigger for early dark energy coincidence

Current cosmological measurements present a persistent tension in the value of the current cosmic expansion rate, the Hubble constant, as inferred from cosmic microwave background (CMB) and large-scale structure (LSS) data compared to that inferred from the classical distance ladder. Early dark energy (EDE), whose cosmological role is localized in time around the epoch of matter-radiation equality just prior to the release of the CMB photons, is designed to resolve this "Hubble tension". However, the model introduces a new coincidence problem: Why should the EDE dynamics occur near matter-radiation equality if EDE is decoupled from both matter and radiation? The resolution of this problem may lie in an early dark sector (EDS), wherein the dark matter mass is dependent on the EDE scalar field. In this work, we construct such an EDS model and show that it naturally resolves the EDE coincidence problem at the background level without any fine-tuning of the coupling to dark matter or of the initial conditions. When fitting to current cosmological data, including that from the local distance ladder, CMB, and LSS, our EDS maximum-likelihood model performs comparably to EDE for resolving the Hubble tension. However, fitting the Planck CMB data requires a specific range of initial field positions to balance the scalar field fluctuations that drive acoustic oscillations, providing testable differences with other EDE models and a platform for future model-building.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Event-by-event direction reconstruction of solar neutrinos in a high light-yield liquid scintillator

The direction of individual B 8 solar neutrinos has been reconstructed using the SNO + liquid scintillator detector. Prompt, directional Cherenkov light was separated from the slower, isotropic scintillation light using time information, and a maximum likelihood method was used to reconstruct the direction of individual scattered electrons. A clear directional signal was observed, correlated with the solar angle. The observation was aided by a period of low primary fluor concentration that resulted in a slower scintillator decay time. This is the first time that event-by-event direction reconstruction in high light-yield liquid scintillator has been demonstrated in a large-scale detector. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Parametrized uncertainties in the spectral function model of neutrino charged-current quasielastic interactions for oscillation analyses

A substantial fraction of systematic uncertainties in neutrino oscillation experiments stem from the lack of precision in modeling the nuclear target in neutrino-nucleus interactions. Whilst this has driven significant progress in the development of improved nuclear models for neutrino scattering, it is crucial that the models used in neutrino data analyses be accompanied by parameters and associated uncertainties that allow the coverage of plausible nuclear physics. Based on constraints from electron scattering data, we propose such a set of parameters, which can be applied to nuclear shell models, and test their application to the Benhar [] spectral function model. The parametrization is validated through a series of maximum likelihood fits to cross section measurements made by the T2K and MINERvA experiments, which also permit an exploration of the power of near-detector data to provide constraints on the parameters in neutrino oscillation analyses. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Measurement of e + e − → η J / ψ cross section from s = 3.808 GeV to 4.951 GeV

Using data samples with an integrated luminosity of 22.42 fb − 1 collected by the BESIII detector operating at the BEPCII storage ring, we measure the cross sections of the e + e − → η J / ψ process at center-of-mass energies from 3.808 to 4.951 GeV. Three structures are observed in the line shape of the measured cross sections. A maximum-likelihood fit with ψ ( 4040 ) , two additional resonances, and a nonresonant component are performed. The mass and width of the first additional state are ( 4219.7 ± 2.5 ± 4.5 ) MeV / c 2 and ( 80.7 ± 4.4 ± 1.4 ) MeV , respectively, consistent with the ψ ( 4230 ) . For the second state, the mass and width are ( 4386 ± 13 ± 17 ) MeV / c 2 and ( 177 ± 32 ± 13 ) MeV , respectively, consistent with the ψ ( 4360 ) . The first uncertainties are statistical, and the second ones are systematic. The statistical significance of ψ ( 4040 ) is 8.0 σ and those for ψ ( 4230 ) and ψ ( 4360 ) are more than 10.0 σ . Published by the American Physical Society 2024

Ablikim, M.↗
DOE OSTI2024

Evidence for B + → K + ν ν ¯ decays

We search for the rare decay B + → K + ν ν ¯ in a 362 fb − 1 sample of electron-positron collisions at the ϒ ( 4 S ) resonance collected with the Belle II detector at the SuperKEKB collider. We use the inclusive properties of the accompanying B meson in ϒ ( 4 S ) → B B ¯ events to suppress background from other decays of the signal B candidate and light-quark pair production. We validate the measurement with an auxiliary analysis based on a conventional hadronic reconstruction of the accompanying B meson. For background suppression, we exploit distinct signal features using machine learning methods tuned with simulated data. The signal-reconstruction efficiency and background suppression are validated through various control channels. The branching fraction is extracted in a maximum likelihood fit. Our inclusive and hadronic analyses yield consistent results for the B + → K + ν ν ¯ branching fraction of [ 2.7 ± 0.5 ( stat ) ± 0.5 ( syst ) ] × 10 − 5 and [ 1.1 − 0.8 + 0.9 ( stat ) − 0.5 + 0.8 ( syst ) ] × 10 − 5 , respectively. Combining the results, we determine the branching fraction of the decay B + → K + ν ν ¯ to be [ 2.3 ± 0.5 ( stat ) − 0.4 + 0.5 ( syst ) ] × 10 − 5 , providing the first evidence for this decay at 3.5 standard deviations. The combined result is 2.7 standard deviations above the standard model expectation. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Calibrating two jets at once

Jet-energy calibration is an important aspect of many measurements and searches at the LHC. Currently, these calibrations are performed on a per-jet basis, i.e., agnostic to the properties of other jets in the same event. In this work, we propose taking advantage of the correlations induced by momentum conservation between jets in order to improve their jet-energy calibration. By fitting the p T asymmetry of dijet events in simulation, while remaining agnostic to the p T spectra themselves, we are able to obtain correlation-improved maximum likelihood estimates. This approach is demonstrated with simulated jets from the CMS detector, yielding a 3%–5% relative improvement in the jet-energy resolution, corresponding to a quadrature improvement of approximately 35%. Published by the American Physical Society 2024

Gambhir, Rikab (ORCID:0000000251080448)↗
DOE OSTI2021

Room-Temperature Mechanical Resonator with a Single Added or Subtracted Phonon

A room-temperature mechanical oscillator undergoes thermal Brownian motion with an amplitude much larger than the amplitude associated with a single phonon of excitation. This motion can be read out and manipulated using laser light using a cavity-optomechanical approach. By performing a strong quantum measurement (i.e., counting single photons in the sidebands imparted on a laser), we herald the addition and subtraction of single phonons on the 300 K thermal motional state of a 4 GHz mechanical oscillator. Further, to understand the resulting mechanical state, we implement a tomography scheme and observe highly non-Gaussian phase-space distributions. Using a maximum likelihood method, we infer the density matrix of the oscillator, and we confirm the counterintuitive doubling of the mean phonon number resulting from phonon addition and subtraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗