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 1,657 records · Page 92

The persistence of large scale structures. Part I. Primordial non-Gaussianity

Abstract We develop an analysis pipeline for characterizing the topology of large scale structure and extracting cosmological constraints based on persistent homology . Persistent homology is a technique from topological data analysis that quantifies the multiscale topology of a data set, in our context unifying the contributions of clusters, filament loops, and cosmic voids to cosmological constraints. We describe how this method captures the imprint of primordial local non-Gaussianity on the late-time distribution of dark matter halos, using a set of N-body simulations as a proxy for real data analysis. For our best single statistic, running the pipeline on several cubic volumes of size 40 (Gpc/h) 3 , we detect f NL loc =10 at 97.5% confidence on ~ 85% of the volumes. Additionally wetest our ability to resolve degeneracies betweenthe topological signature of f NL loc and variation of σ 8 and argue that correctly identifying nonzero f NL loc in this case is possible via an optimal template method. Our method relies on information living at $\mathcal{O}$(10) Mpc/h, a complementary scale with respect to commonly used methods such as the scale-dependent bias in the halo/galaxy power spectrum. Therefore, while still requiring a large volume, our method does not require sampling long-wavelength modes to constrain primordial non-Gaussianity. Moreover, our statistics are interpretable: we are able to reproduce previous results in certain limits and we make new predictions for unexplored observables, such as filament loops formed by dark matter halos in a simulation box.

Astronomy & Astrophysics↗

Constraining the H i–Halo Mass Relation from Galaxy Clustering

We study the dependence of galaxy clustering on H i mass using ∼16,000 galaxies with redshift in the range of 0.0025 10{sup 8} M{sub ⊙}, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with M{sub HI} above different thresholds and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak H i mass dependence, we find both the clustering amplitudes on scales above a few megaparsecs and the bias factors to increase significantly with increasing H i mass for M{sub HI}>10{sup 9} M{sub ⊙}. For H i mass thresholds below ∼10{sup 9} M{sub ⊙}, the inferred galaxy bias factors are systematically lower than the minimum halo bias from mass-selected halo samples. We extend the simple halo model, in which the galaxy content is only determined by halo mass, by including the halo formation time as an additional parameter. A model that puts H i-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of H i mass with halo mass and formation time, as well as the halo occupation distributions and H i mass functions for central and satellite galaxies. These results are compared with the predictions from semianalytic galaxy formation models and hydrodynamic galaxy formation simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

RElative QUantitation Inferred by Evaluating Mixtures (REQUIEM)

Motivated by the lack of easily implementable and generally applicable strategies to increase and assess data accuracy, we devised a novel label-free approach, termed REQUIEM, to address challenges in relative quantitation. For comparing the relative amounts of analytes in two samples, a mixture is prepared from aliquots of the samples, and the samples and the mixture are analyzed in parallel according to the intended workflow. Processing of the resulting data using the REQUIEM algorithm yields unbiased analyte fold-changes and associated statistics, allowing several types of errors to be diagnosed or eliminated. Extensive simulations and analysis of carefully prepared standard samples demonstrated the rigorous foundations of REQUIEM. We applied REQUIEM to several real-world analytical techniques and workflows, notably to tandem mass spectrometry analysis by using isomeric oligosaccharides as test analytes. Finally, we conclude that REQUIEM can reveal inaccuracies in the data that are difficult to identify by using traditional approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimal structure and parameter learning of Ising models

Reconstruction of the structure and parameters of an Ising model from binary samples is a problem of practical importance in a variety of disciplines, ranging from statistical physics and computational biology to image processing and machine learning. The focus of the research community shifted toward developing universal reconstruction algorithms that are both computationally efficient and require the minimal amount of expensive data. Here, we introduce a new method, interaction screening, which accurately estimates model parameters using local optimization problems. The algorithm provably achieves perfect graph structure recovery with an information-theoretically optimal number of samples, notably in the low-temperature regime, which is known to be the hardest for learning. Here, the efficacy of interaction screening is assessed through extensive numerical tests on synthetic Ising models of various topologies with different types of interactions, as well as on real data produced by a D-Wave quantum computer. Finally, this study shows that the interaction screening method is an exact, tractable, and optimal technique that universally solves the inverse Ising problem.

97 MATHEMATICS AND COMPUTING↗

A Population Study of Gaseous Exoplanets

We present here the analysis of 30 gaseous extrasolar planets, with temperatures between 600 and 2400 K and radii between 0.35 and 1.9 R {sub Jup}. The quality of the HST/WFC3 spatially scanned data combined with our specialized analysis tools allow us to study the largest and most self-consistent sample of exoplanetary transmission spectra to date and examine the collective behavior of warm and hot gaseous planets rather than isolated case studies. We define a new metric, the Atmospheric Detectability Index (ADI) to evaluate the statistical significance of an atmospheric detection and find statistically significant atmospheres in around 16 planets out of the 30 analyzed. For most of the Jupiters in our sample, we find the detectability of their atmospheres to be dependent on the planetary radius but not on the planetary mass. This indicates that planetary gravity plays a secondary role in the state of gaseous planetary atmospheres. We detect the presence of water vapour in all of the statistically detectable atmospheres, and we cannot rule out its presence in the atmospheres of the others. In addition, TiO and/or VO signatures are detected with 4σ confidence in WASP-76 b, and they are most likely present in WASP-121 b. We find no correlation between expected signal-to-noise and atmospheric detectability for most targets. This has important implications for future large-scale surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of Aluminum Titanate (Al 2 TiO 5 ) Doping on the Mechanical Performance of Solid Oxide Fuel Cell Ni-YSZ Anode

The mechanical behavior of un-doped and 1–10 wt.% aluminum titanate doped NiO-YSZ anodes was evaluated in the oxidized and reduced state of the material. Sample bars 25 × 5 × 2 mm were fabricated and tested in a three-point bending apparatus and statistical analyzes of the collected data were performed. A remarkable enhancement of the flexural strength was found both for the reduced and oxidized state of the material when compared to the un-doped samples. In both cases, the development of a secondary phase was observed proportional to the doping amount of aluminum titanate. Morphological analysis along with preliminary phase identification and suspected mechanisms for strength enhancement are presented and discussed.

25 ENERGY STORAGE↗

Experimental Evidence of Mechanical Isotropy in Porcine Lung Parenchyma

Pulmonary injuries are a major source of morbidity and mortality associated with trauma. Trauma includes injuries associated with accidents and falls as well as blast injuries caused by explosives. The prevalence and mortality of these injuries has made research of pulmonary injury a major priority. Lungs have a complex structure, with multiple types of tissues necessary to allow successful respiration. The soft, porous parenchyma is the component of the lung which contains the alveoli responsible for gas exchange. Parenchyma is also the portion which is most susceptible to traumatic injury. Finite element simulations are an important tool for studying traumatic injury to the human body. These simulations rely on material properties to accurately recreate real world mechanical behaviors. Previous studies have explored the mechanical properties of lung tissues, specifically parenchyma. These studies have assumed material isotropy but, to our knowledge, no study has thoroughly tested and quantified this assumption. This study presents a novel methodology for assessing isotropy in a tissue, and applies these methods to porcine lung parenchyma. Briefly, lung parenchyma samples were dissected so as to be aligned with one of the three anatomical planes, sagittal, frontal, and transverse, and then subjected to compressive mechanical testing. Stress-strain curves from these tests were statistically compared by a novel method for differences in stresses and strains at percentages of the curve. Histological samples aligned with the anatomical planes were also examined by qualitative and quantitative methods to determine any differences in the microstructural morphology. Our study showed significant evidence to support the hypothesis that lung parenchyma behaves isotropically.

Materials Science↗

Characterizing Natural Attenuation of Cr(VI) in a Chromium Plume at Los Alamos National Laboratory - 18594

Characterization of natural attenuation of hexavalent (Cr(VI)) on aquifer sediments within the Cr plume at Los Alamos National Laboratory has been conducted to (1) refine conceptual and numerical models of chromium fate and transport, (2) evaluate the long-term stability of Cr(III)/Cr(VI) in the complex groundwater system along flow paths from the source to the regional aquifer, and (3) evaluate and develop options for successful remediation. The naturally high chromium abundance in aquifer sediments at Los Alamos pose an unique challenge for characterization of chromium reduction, making conventional approaches such as uptake tests and advanced spectroscopic measurements ineffective in identifying relatively small amounts of anthropogenic Cr(VI) in the sediments. To address this challenge and as complimentary to these conventional methods, sequential leaching tests were conducted on core hole sediments obtained from sonic drilled cores taken from both contaminated and non-contaminated horizons among two mineralogically distinct formations (Puye formation and Miocene Pumiceous Unit) within the Cr plume. A series of three and coherent sequential leaching steps were conducted: 1) an alkaline leach to determine the amount of Cr(VI) sorbed on the sediments. 2) a weak acid leach (using 0.01 M and 0.1 M HNO{sub 3}) to determine the amount of anthropogenic chromium (Cr(III) and Cr(VI)) on aquifer sediment surfaces, and 3) a strong acid leach (using 2 M HNO{sub 3}) to characterize the amount of releasable Cr(III) from aquifer minerals. Additionally we tested the efficiency and detection limit of sequential leaching on bentonite, which consists mainly of smectite, used in drill muds and laboratory-grade biotite spiked with known concentrations of Cr(VI). Smectite and biotite are two phases considered to be potentially important host minerals in aquifer sediments within the plume. Thus, by comparing contaminated and uncontaminated sediment leachate results of chromium and other major and trace elements, the signature and quantity of anthropogenic Cr(VI) can be derived. Our results from the weak acid leaching test indicate that the method detection limit for Cr 'uptake' was about 1 μg/g for biotite and about 0.5 μg/g for bentonite. Our results also indicate that all samples taken from known contaminated horizons in comparison to non-contaminated ones within the Puye formation showed no statistical difference in the amount of Cr leached indicative of its insignificant Cr(VI) attenuation capacity. The amount of excess Cr leached from the magnetic fraction of the contaminated Miocene Pumiceous samples was about 2 μg/g of sediment. Considering that the magnetic fraction was about 11% by mass of the bulk sample, the excess amount of Cr leached from the bulk material would be about 0.2 μg/g of sediment, suggesting the potential for some Cr(VI) attenuation capacity in the Miocene pumice unit. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

From binaries to multiples. II. Hierarchical multiplicity of F and G dwarfs

Statistics of hierarchical multiplicity among solar-type dwarfs are studied using the distance-limited sample of 4847 targets presented in the accompanying Paper I. Known facts about binaries (multiplicity fraction 0.46, lognormal period distribution with median period 100 yr and logarithmic dispersion 2.4, and nearly uniform mass-ratio distribution independent of the period) are confirmed with a high statistical significance. The fraction of hierarchies with three or more components is 0.13 ± 0.01, and the fractions of targets with n = 1, 2, 3, ... components are 54:33:8:4:1. Subsystems in the secondary components are almost as frequent as in the primary components, but in half of such cases both inner pairs are present. The high frequency of those 2+2 hierarchies (4%) suggests that both inner pairs were formed by a common process. The statistics of hierarchies can be reproduced by simulations, assuming that the field is a mixture coming from binary-rich and binary-poor environments. Periods of the outer and inner binaries are selected recursively from the same lognormal distribution, subject to the stability constraint and accounting for the correlation between inner subsystems. The simulator can be used to evaluate the frequency of multiple systems with specified parameters. However, it does not reproduce the observed excess of inner periods shorter than 10 days, caused by tidal evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigations of element spatial correlation in Mn-promoted Co-based Fischer–Tropsch synthesis catalysts

Making connections between performance and structure in bimetallic catalysts requires knowledge of how the two elements are spatially associated. Elemental maps obtained by analytical TEM methods are an invaluable tool for identifying the location of different elements, but for many samples, visual inspection of elemental maps is insufficient for assessing the degree of element spatial correlation. This is particularly true for beam-sensitive materials where short mapping acquisition times lead to images with high noise and low color depth. In these situations, statistical analysis of elemental maps can be used to identify spatial correlations among the elements in a sample. In this work, the relationship between catalyst performance and bimetallic spatial association was explored using Mn-promoted Co-based Fischer-Tropsch synthesis catalysts prepared by different pretreatment methods. Mn was used as a catalyst additive to suppress methane formation. Catalysts that underwent calcination before reduction produced more methane and fewer long-chain hydrocarbons than catalysts that were directly reduced. The extent to which Co and Mn were spatially associated was assessed using correlation metrics, colocation plots, and histograms generated using data from STEM-EDS maps. Although both catalysts yielded visually similar elemental maps, the results of statistical analysis suggested that the calcined catalyst exhibited greater spatial segregation between the Co and Mn. These findings support the hypothesis that having Mn in close proximity to the Co is essential for the manifestation of Mn promotion effects in Co-based FTS catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reconsidering barriers to wind power projects: community engagement, developer transparency and place

In 2016, we undertook a nationally representative wind power perceptions survey of individuals living within 8 km of over 600 projects in the United States, generating 1705 telephone, web, and mail responses. We sought information on a variety of topics, including procedural fairness and its relationship to project attitude, the foci of the present analysis. Here, we present a series of descriptive statistics and regression results, emphasizing those residents who were aware of their local project prior to construction. Sample weighting is employed to account for stratification and non-response. We find that a developer being open and transparent, a community being able to influence the outcome, and having a say in the planning process are all statistically significant predictors of a process perceived as being ‘fair,’ with an open and transparent developer having the largest effect. We also find developer transparency and ability to influence outcomes to have statistically significant relationships to a more positive attitude, with those findings holding when aesthetics, landscape, and wind turbine sound considerations are controlled for. The results indicate that jurisdictions might consider developing procedures, which ensure citizens are consulted and heard, and benchmarks or best practices for developer interaction with communities and citizens.

17 WIND ENERGY↗

Mitigating the impact of the DESI fiber assignment on galaxy clustering

In this paper, we present a simple strategy to mitigate the impact of an incomplete spectroscopic redshift galaxy sample as a result of fiber assignment and survey tiling. The method has been designed for the Dark Energy Spectroscopic Instrument (DESI) galaxy survey but may have applications beyond this. We propose a modification to the usual correlation function that nulls the almost purely angular modes affected by survey incompleteness due to fiber assignment. Predictions of this modified statistic can be calculated given a model of the two point correlation function. The new statistic can be computed with a slight modification to the data catalogues input to the standard correlation function code and does not incur any additional computational time. Lastly, we show that the spherically averaged baryon acoustic oscillation signal is not biased by the new statistic.

79 ASTRONOMY AND ASTROPHYSICS↗

Ozone data and mission sampling analysis

A methodology was developed to analyze discrete data obtained from the global distribution of ozone. Statistical analysis techniques were applied to describe the distribution of data variance in terms of empirical orthogonal functions and components of spherical harmonic models. The effects of uneven data distribution and missing data were considered. Data fill based on the autocorrelation structure of the data is described. Computer coding of the analysis techniques is included.

Robbins, J. L.↗

Convergence Rates for Empirical Estimation of Binary Classification Bounds

Bounding the best achievable error probability for binary classification problems is relevant to many applications including machine learning, signal processing, and information theory. Many bounds on the Bayes binary classification error rate depend on information divergences between the pair of class distributions. Recently, the Henze–Penrose (HP) divergence has been proposed for bounding classification error probability. We consider the problem of empirically estimating the HP-divergence from random samples. We derive a bound on the convergence rate for the Friedman–Rafsky (FR) estimator of the HP-divergence, which is related to a multivariate runs statistic for testing between two distributions. The FR estimator is derived from a multicolored Euclidean minimal spanning tree (MST) that spans the merged samples. We obtain a concentration inequality for the Friedman–Rafsky estimator of the Henze–Penrose divergence. We validate our results experimentally and illustrate their application to real datasets.

97 MATHEMATICS AND COMPUTING↗

A Magellan M2FS Spectroscopic Survey of Galaxies at 5.5 < z < 6.8: Program Overview and a Sample of the Brightest Lyα Emitters

We present a spectroscopic survey of high-redshift, luminous galaxies over four square degrees on the sky, aiming to build a large and homogeneous sample of Lyα emitters (LAEs) at z≈5.7 and 6.5, and Lyman-break galaxies (LBGs) at 5.5<z<6.8. The fields that we choose to observe are well studied, such as by the Subaru XMM-Newton Deep Survey and COSMOS. They have deep optical imaging data in a series of broad and narrow bands, allowing for the efficient selection of galaxy candidates. Spectroscopic observations are being carried out using the multi-object spectrograph M2FS on the Magellan Clay telescope. M2FS is efficient enough to identify high-redshift galaxies, owing to its 256 optical fibers deployed over a circular field of view 30{sup ′} in diameter. We have observed ∼2.5 square degrees. When the program is completed, we expect to identify more than 400 bright LAEs at z≈5.7 and 6.5, and a substantial number of LBGs at z⩾6. This unique sample will be used to study a variety of galaxy properties and to search for large protoclusters. Furthermore, the statistical properties of these galaxies will be used to probe cosmic reionization. We describe the motivation, program design, target selection, and M2FS observations. We also outline our science goals, and present a sample of the brightest LAEs at z≈5.7 and 6.5. This sample contains 32 LAEs with Lyα luminosities higher than 10{sup 43} erg s{sup −1}. A few of them reach ≥3 × 10{sup 43} erg s{sup −1}, comparable to the two most luminous LAEs known at z⩾6, “CR7” and “COLA1.” These LAEs provide ideal targets to study extreme galaxies in the distant universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of AVHRR, CZCS and historical in situ data off the Oregon Coast

The original scientific objectives of this grant were to: (1) characterize the seasonal cycles and interannual variability for phytoplankton concentrations and sea surface temperature (SST) in the California Current using satellite data; and (2) to explore the spatial and temporal relationship between these variables and surface wind forcing. An additional methodological objective was to develop statistical methods for forming mean fields, which minimize the effects of random data gaps and errors in the irregularly sampled CZCS (Coastal Zone Color Scanner) and AVHRR (Advanced Very High Resolution Radiometer) satellite data. A final task was to evaluate the level of uncertainty in the wind fields used for the statistical analysis. Funding in the first year included part of the cost of an image processing system to enable this and other projects to process and analyze satellite data. This report consists of summaries of the major projects carried out with all or partial support from this grant. The appendices include a list of papers and professional presentations supported by the grant, as well as reprints of the major papers and reports.

Strub, P. Ted↗

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗