Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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 721 records · Page 40

Energy calibration of germanium detectors for the Majorana Demonstrator

The Majorana Demonstrator was a search for neutrinoless double-beta decay (0νββ) in the 76 Ge isotope. It was staged at the 4850-foot level of the Sanford Underground Research Facility (SURF) in Lead, SD. The experiment consisted of 58 germanium detectors housed in a low background shield and was calibrated once per week by deploying a 228 Th line source for 1 to 2 hours. The energy scale calibration determination for the detector array was automated using custom analysis tools. We describe the offline procedure for calibration of the Demonstrator germanium detectors, including the simultaneous fitting of multiple spectral peaks, estimation of energy scale uncertainties, and the automation of the calibration procedure.

47 OTHER INSTRUMENTATION↗

ALICE luminosity determination for Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Luminosity determination within the ALICE experiment is based on the measurement, in van der Meer scans, of the cross sections for visible processes involving one or more detectors (visible cross sections). In 2015 and 2018, the Large Hadron Collider provided Pb–Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{NN}}$ = 5.02 TeV. Two visible cross sections, associated with particle detection in the Zero Degree Calorimeter (ZDC) and in the V0 detector, were measured in a van der Meer scan. This article describes the experimental set-up and the analysis procedure, and presents the measurement results. The analysis involves a comprehensive study of beam-related effects and an improved fitting procedure, compared to previous ALICE studies, for the extraction of the visible cross section. The resulting uncertainty of both the ZDC-based and the V0-based luminosity measurement for the full sample is 2.5%. The inelastic cross section for hadronic interactions in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, obtained by efficiency correction of the V0-based visible cross section, was measured to be 7.67 ± 0.25 b, in agreement with predictions using the Glauber model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparison of point cloud and image-based models for calorimeter fast simulation

Score based generative models are a new class of generative models that have been shown to accurately generate high dimensional calorimeter datasets. Recent advances in generative models have used images with 3D voxels to represent and model complex calorimeter showers. Point clouds, however, are likely a more natural representation of calorimeter showers, particularly in calorimeters with high granularity. Furthermore, point clouds preserve all of the information of the original simulation, more naturally deal with sparse datasets, and can be implemented with more compact models and data files. In this work, two state-of-the-art score based models are trained on the same set of calorimeter simulation and directly compared.

47 OTHER INSTRUMENTATION↗

Photon classification with Gradient Boosted Trees at CLAS12

Dihadron semi-inclusive deep inelastic scattering (SIDIS) of 10.6 GeV longitudinally polarized electrons off the proton has been measured using the CLAS12 detector at Jefferson Lab. Two separate channels, π + π 0 and π - π 0 , were analyzed, requiring the reconstruction of diphoton pairs. Here, in this analysis, we addressed the problem of false neutral particles being reconstructed by CLAS12's event builder, polluting the otherwise physical combinatorial background underneath the π 0 peak. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained with Monte Carlo simulations to reduce the amount of background π 0 's. We show that the nearest-neighbor features learned by the model lead to a substantial increase in signal vs. background discrimination compared to previous CLAS12 π^0 analyses. The machine learning approach recovers several times more dihadron statistics for the dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Batch VUV4 characterization for the SBC-LAr10 scintillating bubble chamber

The Scintillating Bubble Chamber (SBC) collaboration purchased 32 Hamamatsu VUV4 silicon photomultipliers (SiPMs) for use in SBC-LAr10, a bubble chamber containing 10 kg of liquid argon. A dark-count characterization technique, which avoids the use of a single-photon source, was used at two temperatures to measure the VUV4 SiPMs breakdown voltage (V BD ), the SiPM gain (g SiPM ), the rate of change of g SiPM with respect to voltage (m), the dark count rate (DCR), and the probability of a correlated avalanche (P CA ) as well as the temperature coefficients of these parameters. A Peltier-based chilled vacuum chamber was developed at Queen's University to cool down the Quads to 233.15 ± 0.2 K and 255.15 ± 0.2 K with average stability of ±20 mK. An analysis framework was developed to estimate V BD to tens of mV precision and DCR close to Poissonian error. The temperature dependence of V BD was found to be 56 ± 2 mV K -1 , and m on average across all Quads was found to be (459 ± 3(stat.)±23(sys.))× 10 3 e- PE -1 V -1 . The average DCR temperature coefficient was estimated to be 0.099 ± 0.008 K -1 corresponding to a reduction factor of 7 for every 20 K drop in temperature. The average temperature dependence of P CA was estimated to be 4000 ± 1000 ppm K -1 . P CA estimated from the average across all SiPMs is a better estimator than the P CA calculated from individual SiPMs, for all of the other parameters, the opposite is true. All the estimated parameters were measured to the precision required for SBC-LAr10, and the Quads will be used in conditions to optimize the signal-to-noise ratio.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Beam-beam backgrounds for the Cool Copper Collider

In this paper, we present a comprehensive characterization of beam-beam backgrounds for the Cool Copper Collider (C 3 ), a proposed linear e + e - collider designed for precision Higgs studies at center-of-mass energies of 250 and 550 GeV. Using a simulation pipeline based on the Key4hep framework, we evaluate incoherent pair production and hadron photoproduction backgrounds through the SiD detector for baseline, power-efficiency, and high-luminosity C 3 operating scenarios. The occupancy induced by the beam-beam background is evaluated for each scenario, validating the compatibility of the existing SiD detector design with operations at C 3 without substantial modifications. Furthermore, at the same time, the modular simulation framework and analysis methodology presented in this paper offer a versatile toolkit for background studies in future collider proposals, contributing to a common platform for different machine designs.

Analysis and statistical methods↗

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

Bayesian characterization of explosive sources using infrasonic signals

In this paper, a Bayesian approach is applied to the estimation of source characteristics for above ground explosive events using infrasonic observations at propagation distances of hundreds of kilometers. Propagation-based, stochastic path geometry models have previously been investigated for infrasonic source localization with promising results; here, statistical transmission loss models are investigated using archived atmospheric specifications and propagation simulation tools to demonstrate that such an approach is applicable to source characterization. The resulting transmission loss models models are utilized in estimating near-source spectral characteristics for explosive events from infrasonic observations at regional distances of several hundred kilometers or more using a Bayesian framework to compute a combined estimate. The near-source spectral characterization can be used along with a simplistic explosive source model to compute an estimate of explosive yield with quantifiable uncertainty. Finally, the yield estimation framework is evaluated using above-ground explosive events from the western US during the fall of 2012 with promising results.

58 GEOSCIENCES↗

Anisotropy as an indicator for reservoir changes: example from the Rotokawa and Ngatamariki geothermal fields, New Zealand

We investigate the relation between geothermal field production and fracture density and orientation in the Ngatamariki and Rotokawa geothermal fields, located in the Taupo Volcanic Zone, New Zealand using shear wave splitting (SWS). We determine the SWS parameters for 17 702 microseismic events across 38 stations spanning close to 4 yr from 2012 to 2015. We compare the strength of anisotropy to changes in field production and injection. We also compare the orientation of the anisotropy to in situ and regional measurements of maximum horizontal stress orientation. ($S_{\mathrm{ H}_{\mathrm{ max}}}$). Due to the volume of unique events (approximately 160 000), shear wave phases are picked automatically. We carry out automatic SWS measurements using the Multiple Filter Automatic Splitting Technique (MFAST). The SWS measurements are interpreted in the context of stress aligned microcracks. Outside both fields and within Ngatamariki, fast polarizations align with the NE–SW regional orientation of $S_{\mathrm{ H}_{\max}}$. Within Rotokawa a greater complexity is observed, with polarizations tending toward N–S. We observe increases in per cent anisotropy coinciding with the start of production/injection in Ngatamariki and then a later correlated drop in per cent anisotropy and v P /v S ratios in southern Ngatamariki as injection is shifted to the north. This relationship is consistent with pore fluid pressure within the reservoir being affected by local changes in production and injection of geothermal fluids causing cracks to open and close in response.

58 GEOSCIENCES↗

A framework for measuring weak-lensing magnification using the Fundamental Plane

ABSTRACT Galaxy–galaxy lensing is an essential tool for probing dark matter haloes and constraining cosmological parameters. While galaxy–galaxy lensing measurements usually rely on shear, weak-lensing magnification contains additional constraining information. Using the Fundamental Plane (FP) of elliptical galaxies to anchor the size distribution of a background population is one method that has been proposed for performing a magnification measurement. We present a formalism for using the FP residuals of elliptical galaxies to jointly estimate the foreground mass and background redshift errors for a stacked lens scenario. The FP residuals include information about weak-lensing magnification κ, and therefore foreground mass, since to first order, non-zero κ affects galaxy size but not other FP properties. We also present a modular, extensible code that implements the formalism using emulated galaxy catalogues of a photometric galaxy survey. We find that combining FP information with observed number counts of the source galaxies constrains mass and photo-z error parameters significantly better than an estimator that includes number counts only. In particular, the constraint on the mass is 17.0 per cent if FP residuals are included, as opposed to 27.7 per cent when only number counts are included. The effective size noise for a foreground lens of mass $M_\mathrm{ H}=10^{14}\, \mathrm{M}_\odot$, with a conservative selection function in size and surface brightness applied to the source population, is σκ, eff = 0.250. We discuss the improvements to our FP model necessary to make this formalism a practical companion to shear analyses in weak-lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Implications of the search for optical counterparts during the second part of the Advanced LIGO’s and Advanced Virgo’s third observing run: lessons learned for future follow-up observations

ABSTRACT Joint multimessenger observations with gravitational waves and electromagnetic (EM) data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on 2019 April 1; during the 11 months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no EM counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from 2019 October to 2020 March. We highlight two neutron star–black hole candidates (S191205ah and S200105ae), two binary neutron star candidates (S191213g and S200213t), and a binary merger with a possible neutron star and a ‘MassGap’ component, S200115j. Assuming that the gravitational-wave (GW) candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential GW candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage versus observational depth to understand the likelihood of counterpart detection.

79 ASTRONOMY AND ASTROPHYSICS↗

Noise from undetected sources in Dark Energy Survey images

ABSTRACT For ground-based optical imaging with current CCD technology, the Poisson fluctuations in source and sky background photon arrivals dominate the noise budget and are readily estimated. Another component of noise, however, is the signal from the undetected population of stars and galaxies. Using injection of artifical galaxies into images, we demonstrate that the measured variance of galaxy moments (used for weak gravitational lensing measurements) in Dark Energy Survey (DES) images is significantly in excess of the Poisson predictions, by up to 30 per cent, and that the background sky levels are overestimated by current software. By cross-correlating distinct images of ‘empty’ sky regions, we establish that there is a significant image noise contribution from undetected static sources (US), which, on average, are mildly resolved at DES resolution. Treating these US as a stationary noise source, we compute a correction to the moment covariance matrix expected from Poisson noise. The corrected covariance matrix matches the moment variances measured on the injected DES images to within 5 per cent. Thus, we have an empirical method to statistically account for US in weak lensing measurements, rather than requiring extremely deep sky simulations. We also find that local sky determinations can remove most of the bias in flux measurements, at a small penalty in additional, but quantifiable, noise.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Validation of selection function, sample contamination and mass calibration in galaxy cluster samples

ABSTRACT We construct and validate the selection function of the MARD-Y3 galaxy cluster sample. This sample was selected through optical follow-up of the 2nd ROSAT faint source catalogue with Dark Energy Survey year 3 data. The selection function is modelled by combining an empirically constructed X-ray selection function with an incompleteness model for the optical follow-up. We validate the joint selection function by testing the consistency of the constraints on the X-ray flux–mass and richness–mass scaling relation parameters derived from different sources of mass information: (1) cross-calibration using South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) clusters, (2) calibration using number counts in X-ray, in optical and in both X-ray and optical while marginalizing over cosmological parameters, and (3) other published analyses. We find that the constraints on the scaling relation from the number counts and SPT-SZ cross-calibration agree, indicating that our modelling of the selection function is adequate. Furthermore, we apply a largely cosmology independent method to validate selection functions via the computation of the probability of finding each cluster in the SPT-SZ sample in the MARD-Y3 sample and vice versa. This test reveals no clear evidence for MARD-Y3 contamination, SPT-SZ incompleteness or outlier fraction. Finally, we discuss the prospects of the techniques presented here to limit systematic selection effects in future cluster cosmological studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced forLarge Synoptic Survey Telescope The Rubin Observatory Legacy Survey of Space and Time (lsst) Dark Energy Science Collaboration (desc). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. Additionally, we report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-Gaussianity in the weak lensing correlation function likelihood – implications for cosmological parameter biases

ABSTRACT We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in 1D marginal distributions of shear two-point correlation functions in simulated weak lensing data. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in 1D posteriors of Ωm and σ8 due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data (100 deg2). We also present a systematic approach to constructing approximate multivariate likelihoods with 1D parametric functions by assuming independence or more flexible non-parametric multivariate methods after decorrelating the data points using principal component analysis (PCA). While the use of PCA does not modify the non-Gaussianity of the multivariate likelihood, we find empirically that the 1D marginal sampling distributions of the PCA components exhibit less skewness and kurtosis than the original shear correlation functions. Modelling the likelihood with marginal parametric functions based on the assumption of independence between PCA components thus gives a lower limit for the biases. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of ∼5.

79 ASTRONOMY AND ASTROPHYSICS↗

Hemispherical variance anomaly and reionization optical depth

ABSTRACT Cosmic microwave background (CMB) full-sky temperature data show a hemispherical asymmetry in power nearly aligned with the Ecliptic, with the Northern hemisphere displaying an anomalously low variance, while the Southern hemisphere appears consistent with expectations from the best-fitting theory, Lambda Cold Dark Matter (ΛCDM). The low signal-to-noise ratio in current polarization data prevents a similar comparison. Polarization realizations constrained by temperature data show that in ΛCDM the lack of variance is not expected to be present in polarization data. Therefore, a natural way of testing whether the temperature result is a fluke is to measure the variance of CMB polarization components. In anticipation of future CMB experiments that will allow for high-precision large-scale polarization measurements, we study how the variance of polarization depends on ΛCDM-parameter uncertainties by forecasting polarization maps with Planck’s Markov chain Monte Carlo chains. We show that polarization variance is sensitive to present uncertainties in cosmological parameters, mainly due to current poor constraints on the reionization optical depth τ, which drives variance at low multipoles. We demonstrate how the improvement in the τ measurement seen between Planck’s two latest data releases results in a tighter constraint on polarization variance expectations. Finally, we consider even smaller uncertainties on τ and how more precise measurements of τ can drive the expectation for polarization variance in a hemisphere close to that of the cosmic-variance-limited distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

79 ASTRONOMY AND ASTROPHYSICS↗