Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling statistics”

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 163 records · Page 9

Forward modeling fluctuations in the DESI LRGs target sample using image simulations

We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Rare Event Sampling Improves Mercury Instability Statistics

Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare-event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately 10 –4 in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event.

79 ASTRONOMY AND ASTROPHYSICS↗

Boosted dark matter and its implications for the features in IceCube HESE data

Here we study the implications of the premise that any new, relativistic, highly energetic neutral particle that interacts with quarks and gluons would create cascade-like events in the IceCube (IC) detector. Such events would be observationally indistinguishable from neutral current deep-inelastic (DIS) scattering events due to neutrinos. Consequently, one reason for deviations, breaks or excesses in the expected astrophysical power-law neutrino spectrum could be the flux of such a particle. Motivated by features in the recent 1347-day IceCube high energy starting event (HESE) data, we focus on particular boosted dark matter (χ) related realizations of this premise. Here, χ is assumed to be much lighter than, and the result of, the slow decay of a massive scalar (φ ) which constitutes a major fraction of the Universe's dark matter (DM) . We show that this hypothesis, coupled with a standard power-law astrophysical neutrino flux is capable of providing very good fits to the present data, along with a possible explanation of other features in the HESE sample. These features include a) the paucity of events beyond ~ 2 PeV b) a spectral feature resembling a dip or a spectral change in the 400 TeV–1 PeV region and c) an excess in the 50-100 TeV region. We consider two different boosted DM scenarios, and determine the allowed mass ranges and couplings for four different types of mediators (scalar, pseudoscalar, vector and axial-vector) which could connect the standard and dark sectors. We consider constraints from gamma-ray observations and collider searches. We find that the gamma-ray observations provide the most restrictive constraints, disfavouring the 1σ allowed parameter space from IC fits, while still being consistent with the 3σ allowed region. We also test our proposal and its implications against the (statistically independent) sample of six year through-going muon track data recently released by IceCube.

79 ASTRONOMY AND ASTROPHYSICS↗

THE DATA REDUCTION PIPELINE FOR THE SDSS-IV MaNGA IFU GALAXY SURVEY

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622–10354 Å and an average footprint of ∼500 arcsec{sup 2} per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low-redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ∼100 million raw-frame spectra and ∼10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline algorithms and centralized metadata framework that produce sky-subtracted spectrophotometrically calibrated spectra and rectified three-dimensional data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13, we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ∼8500 Å and reach a typical 10 σ limiting continuum surface brightness μ = 23.5 AB arcsec{sup −2} in a five-arcsecond-diameter aperture in the g -band. The wavelength calibration of the MaNGA data is accurate to 5 km s{sup −1} rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of σ = 72 km s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗

Evidence for non-exponential elastic proton-proton differential cross-section at low |t| and $\sqrt{s}$ = 8 TeV by TOTEM

The TOTEM experiment has made a precise measurement of the elastic proton–proton differential cross-section at the centre-of-mass energy √s = 8 TeV based on a high-statistics data sample obtained with the β* = 90 m optics. Both the statistical and systematic uncertainties remain below 1%, except for the t -independent contribution from the overall normalisation. This unprecedented precision allows to exclude a purely exponential differential cross-section in the range of four-momentum transfer squared 0.027 < |t| < 0.2 GeV 2 with a significance greater than 7 σ. Two extended parametrisations, with quadratic and cubic polynomials in the exponent, are shown to be well compatible with the data. As a result, using them for the differential cross-section extrapolation to t=0, and further applying the optical theorem, yields total cross-section estimates of (101.5 ± 2.1) mb and (101.9 ± 2.1) mb, respectively, in agreement with previous TOTEM measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enhancing Model Predictability for a Scramjet Using Probabilistic Learning on Manifolds

In this study, the computational burden of a large-eddy simulation for reactive flows is exacerbated in the presence of uncertainty in flow conditions or kinetic variables. A comprehensive statistical analysis, with a sufficiently large number of samples, remains elusive. Statistical learning is an approach that allows for extracting more information using fewer samples. Such procedures, if successful, will greatly enhance the predictability of models in the sense of improving exploration and characterization of uncertainty due to model error and input dependencies, all while being constrained by the size of the associated statistical samples. In this paper, it is shown how a recently developed procedure for probabilistic learning on manifolds can serve to improve the predictability in a probabilistic framework of a scramjet simulation. The estimates of the probability density functions of the quantities of interest are improved together with estimates of the statistics of their maxima. Lastly, it is also demonstrated how the improved statistical model adds critical insight to the performance of the model.

42 ENGINEERING↗

Earth mover’s distance as a measure of CP violation

We introduce a new unbinned two sample test statistic sensitive to CP violation utilizing the optimal transport plan associated with the Wasserstein (earth mover’s) distance. The efficacy of the test statistic is shown via two examples of CP asymmetric distributions with varying sample sizes: the Dalitz distributions of B 0 → K + π – π 0 and of D 0 → π + π – π 0 decays. The windowed version of the Wasserstein distance test statistic is shown to have comparable sensitivity to CP violation as the commonly used energy test statistic, but also retains information about the localized distributions of CP asymmetry over the Dalitz plot. For large statistic datasets we introduce two modified Wasserstein distance based test statistics — the binned and the sliced Wasserstein distance statistics, which show comparable sensitivity to CP violation, but improved computing time and memory scalings. Finally, general extensions and applications of the introduced statistics are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Minimum local distance density estimation

We present a local density estimator based on first-order statistics. To estimate the density at a point, $x$, the original sample is divided into subsets and the average minimum sample distance to $x$ over all such subsets is used to define the density estimate at $x$. The tuning parameter is thus the number of subsets instead of the typical bandwidth of kernel or histogram-based density estimators. The proposed method is similar to nearest-neighbor density estimators but it provides smoother estimates. We derive the asymptotic distribution of this minimum sample distance statistic to study globally optimal values for the number and size of the subsets. Simulations are used to illustrate and compare the convergence properties of the estimator. The results show that the method provides good estimates of a wide variety of densities without changes of the tuning parameter, and that it offers competitive convergence performance.

97 MATHEMATICS AND COMPUTING↗

Microstructure Development of 308L Stainless Steel During Additive Manufacturing

$In$ $situ$ high-energy X-ray diffraction measurements were completed during deposition of 308L stainless steel wire onto a 304L stainless steel substrate. Attempts were made to extract microstructural features such as phase fraction and internal stress, as well as temperature evolution immediately following the deposition. The limited data that could be collected during deposition and rapid solidification are critically examined. High-energy X-rays coupled with relatively slow detectors were utilized to enable determination of orientation-dependent lattice parameters accurately enough to comment on phase strain evolution between austenite and ferrite. Here, information about the hydrostatic and deviatoric stress states of the constituent phases was determined on time scales that are relevant to their development. However, the time resolution of the technique was insufficient to monitor phase evolution during the solid–solid phase transformation and, more so, during solidification. Moreover, the accurate and absolute determination of inherently statistical parameters, such as phase fraction, depends critically on the ability to sample a statistically significant numbers of grains in the microstructure.

36 MATERIALS SCIENCE↗

Quantification of Organic Porosity and Water Accessibility in Marcellus Shale Using Neutron Scattering

Pores within organic matter (OM) are a significant contributor to the total pore system in gas shales. These pores contribute most of the storage capacity in gas shales. Here we present a novel approach to characterize the OM pore structure (including the porosity, specific surface area, pore size distribution, and water accessibility) in Marcellus shale. By using ultrasmall and small-angle neutron scattering, and by exploiting the contrast matching of the shale matrix with suitable mixtures of deuterated and protonated water, both total and water-accessible porosity were measured on centimeter-sized samples from two boreholes from the nanometer to micrometer scale with good statistical coverage. Samples were also measured after combustion at 450 °C. Analysis of scattering data from these procedures allowed quantification of OM porosity and water accessibility. OM hosts 24–47% of the total porosity for both organic-rich and -poor samples. This porosity occupies as much as 29% of the OM volume. In contrast to the current paradigm in the literature that OM porosity is organophilic and therefore not likely to contain water, our results demonstrate that OM pores with widths >20 nm exhibit the characteristics of water accessibility. In conclusion, our approach reveals the complex structure and wetting behavior of the OM porosity at scales that are hard to interrogate using other techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of MARSSIM and MARSAME Guidance to Footprint Reduction Projects at a DOE Site (LA-UR-15-28730) - 16216

Department of Energy regulations require that potentially contaminated items and materials being released to the public are monitored and that any residual radioactive contamination does not contribute more than 10 μSv/y (1 mrem/y). In footprint reduction and demolition projects, buildings with potential residual contamination must be carefully and systematically surveyed prior to recycling materials or releasing debris for burial in commercial landfills. Multi-agency guidance is available to assure statistically representative sampling of the items and materials, appropriate measurement quality objectives, and valid statistical approaches for release decisions. Additionally, implementing a cost-effective evaluation program for releases requires working closely with regulators and waste generators and developing a technically sound process to account for naturally occurring radioactive material (NORM) in building materials. Los Alamos National Laboratory has successfully implemented the available guidance for footprint reduction projects to help divert clean, recyclable materials from landfills and limit the possibility of mistakenly disposing of NORM as low level waste. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development and Comparison of Purgeable Mercury Values in SRR Samples Measured by SRNL and Eurofins FGS

Savannah River Remediation (SRR) requested the development of mercury speciation capabilities at the Savannah River National Laboratory (SRNL) to support the Liquid Waste Operations at SRS. As part of that method development, SRR requested that SRNL Analytical Development (AD) compare their results with those obtained from their outside contract laboratory, Eurofins Frontier Global Sciences (FGS). This document reports on this method development work performed at SRNL as well as the comparative analyses conducted between the two laboratories. Development, optimization, and validation were undertaken at SRNL to produce a method for the species-specific analysis of dissolved gaseous mercury (DGM) in samples from SRR. It was determined that DGM present in SRR samples represented a combination of various soluble and volatile mercury species- and is thus referred to in this document as purgeable mercury. A method has been developed and analytically validated for the analysis of purgeable mercury in caustic nuclear waste. Using multivariate optimization, an efficient and sensitive method was produced. Total recovery of spiked purgeable mercury from deionized water was 99.4 ±10.6% and 91.6 ±14.9% from spiked high-activity tank waste samples. Calibration, by external calibration and standard addition calibration, maintained a mean accuracy of 101 ±5% and 99.9 ±1.7%, respectively. An interlaboratory comparison was performed between SRNL and Eurofins FGS in the co-analysis of three quarterly SRR Tank 50 samples. The statistical evaluations presented in this document support the following conclusions: for internal validation data, the two laboratories demonstrate biases of -4.3% for SRNL and 4.4% for Eurofins FGS, though neither are statistically significant at the 5% level. While the percent relative standard deviation for each of the laboratories is rather large (21.5% for SRNL and 18.4% for Eurofins FGS), there is no indication of a statistically significant difference between the two laboratories in the precision of these results. For the three sets of SRR quarterly Tank 50 samples, there is no indication of a statistically significant bias between the two laboratories at the 5% significance level. These results indicate that SRNL may provide data to the customer of statistically similar quality to Eurofins FGS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Gamma Ray Spectrometry Assay Results of 242 Pu Samples at the Chemistry and Metallurgy Research Facility

Thirty individual items at the Chemistry and Metallurgy Research (CMR) facility with accountable amounts of 242 Pu required assay by gamma ray spectrometry in order to establish their masses with reasonable precision. There were two phases to this measurement campaign: 1) a preliminary screening phase in which each item was measured for a short 600 second count time to establish the relative masses between similar or like items, and 2) a second phase in which groups of items were counted for a considerably longer time in order to establish the group mass of 242 Pu in a reasonable time frame. The data from the two phases were assessed in a manner that allowed the project team to assign a 242 Pu mass to each individual item in the inventory, with the exception of three column tubes that could not be measured individually because of radiological contamination issues. The project team performed technical analyses with the SNAP™ analysis routine, PC-FRAM plutonium isotopic routine, and the NEN-2 PeakEasy routine. When total measurement uncertainties are factored in, the individual mass quantities calculated for all samples are statistically equivalent to the expected mass of 0.88 grams of total plutonium per sample.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

D-Guaranteed Discrimination of Statistical Hypotheses: a Review of Results and Unsolved Problems

We compare two sequential d-guaranteed tests and an optimal d-guaranteed test based on a fixed number of observations with respect to the average number of observations within the most accepted practical applications of Bayesian probabilistic models. We consider the sequential “first skipping” test, the sequential locally efficient test based on the score statistic, and the test based on a fixed number of observations that minimizes the necessary sample size. We study various characteristics of these tests connected with the number of observations within three probabilistic models, namely, the normal (ϑ, σ{sup 2}) distribution of the observed random variable and the normal a priori distribution of ϑ with fixed σ{sup 2}; the exponential distribution with the intensity parameter ϑ and the a priori gamma distribution of ϑ; and the Bernoulli sampling with the success probability ϑ and the a priori beta distribution of ϑ. We discuss the connection of the d-posterior approach with the compound decision problem as applied to the analysis of data provided by microchips (when the false discovery rate, or FDR, for short, is treated as the d-risk of the first kind). We present the vast data on characteristics of the mentioned tests obtained by the method of mathematical modeling in several tables. We discuss unsolved problems of the d-guaranteed discrimination of hypotheses with the minimal number of observations and approaches to their solution.

97 MATHEMATICS AND COMPUTING↗

Statistical Methods for Improved Evaluation of Environmental Sample Data quarterly FY18Q2

We will develop an enhanced statistical framework to evaluate environmental swipe samples collected at nuclear facilities with the goal of identifying anomalies associated with potential undeclared activities. Our approach will be based on Bayesian methodology and time series analysis. This overall framework should result in a more cost effective and operationally effective approach to environmental swipe sample analysis and interpretation.

97 MATHEMATICS AND COMPUTING↗

The Populations of Carina. I. Decoding the Color–Magnitude Diagram

We investigate the color–magnitude diagram (CMD) of the Carina dwarf spheroidal galaxy using data of Stetson et al. and synthetic CMDs based on isochrones of Dotter et al., in terms of the parameters [Fe/H], age, and [ α /Fe], for the cases when (i) [ α /Fe] is held constant and (ii) [ α /Fe] is varied. The data are well described by four basic epochs of star formation, having [Fe/H] = −1.85, −1.5, −1.2, and ∼−1.15 and ages ∼13, 7, ∼3.5, and ∼1.5 Gyr, respectively (for [ α /Fe] = 0.1, constant [ α /Fe], and [ α /Fe] = 0.2, 0.1, −0.2, −0.2, variable [ α /Fe]), with small spreads in [Fe/H] and age of order 0.1 dex and 1–3 Gyr. Within an elliptical radius 13.′1, the mass fractions of the populations, at their times of formation, were (in decreasing age order) 0.34, 0.39, 0.23, and 0.04. This formalism reproduces five observed CMD features (two distinct subgiant branches of old and intermediate-age populations, two younger, main-sequence components, and the small color dispersion on the red giant branch (RGB). The parameters of the youngest population are less certain than those of the others, and given it is less centrally concentrated, it may not be directly related to them. High-resolution spectroscopically analyzed RGB samples appear statistically incomplete compared with those selected using radial velocity, which contain bluer stars comprising ∼5%–10% of the samples. We conjecture these objects may, at least in part, be members of the youngest population. We use the CMD simulations to obtain insight into the population structure of Carina's upper RGB.

79 ASTRONOMY AND ASTROPHYSICS↗

The Discovery of λ Bootis Stars: The Southern Survey I

The λ Boo stars are a class of chemically peculiar Population I A-type stars characterized by under-abundances of the refractory elements, but near-solar abundances of carbon, nitrogen, oxygen, and sulfur. There is some evidence that λ Boo stars have higher frequencies of “bright” debris disks than normal A-type stars. The discovery of four exoplanets orbiting HR 8799, a λ Boo star with a resolved debris disk, suggests that the λ Boo phenomenon may be related to the presence of a dynamic debris disk, perhaps perturbed by migrating planets. However, only 64 λ Boo stars are known, and those stars were discovered with different techniques, making it problematic to use that sample for statistical purposes, including determining the frequency of debris disks. The purpose of this paper is to derive a new sample of λ Boo stars using a technique that does not lead to biases with respect to the presence of infrared excesses. Through spectroscopic observations in the southern hemisphere, we have discovered 33 λ Boo stars and have confirmed 12 others. As a step toward determining the proportion of λ Boo stars with infrared excesses, we have used WISE data to examine the infrared properties of this sample out to 22 μ m. On this basis, we cannot conclude that λ Boo stars have a greater tendency than normal A-type stars to show infrared excesses. However, observing this sample at longer wavelengths may change that conclusion, as many λ Boo debris disks are cool and do not radiate strongly at 22 μ m.

79 ASTRONOMY AND ASTROPHYSICS↗

Huge ensembles – Part 2: Properties of a huge ensemble of hindcasts generated with spherical Fourier neural operators

Abstract. In Part 1, we created an ensemble based on spherical Fourier neural operators. As initial condition perturbations, we used bred vectors, and as model perturbations, we used multiple checkpoints trained independently from scratch. Based on diagnostics that assess the ensemble's physical fidelity, our ensemble has comparable performance to operational weather forecasting systems. However, it requires orders-of-magnitude fewer computational resources. Here in Part 2, we generate a huge ensemble (HENS), with 7424 members initialized each day of summer 2023. We enumerate the technical requirements for running huge ensembles at this scale. HENS precisely samples the tails of the forecast distribution and presents a detailed sampling of internal variability. HENS has two primary applications: (1) as a large dataset with which to study the statistics and drivers of extreme weather and (2) as a weather forecasting system. For extreme climate statistics, HENS samples events 4σ away from the ensemble mean. At each grid cell, HENS increases the skill of the most accurate ensemble member and enhances coverage of possible future trajectories. As a weather forecasting model, HENS issues extreme weather forecasts with better uncertainty quantification. It also reduces the probability of outlier events, in which the verification value lies outside the ensemble forecast distribution.

Mahesh, Ankur↗