Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling size”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

A 2-D Piston Effect Solution for the Relaxation in the MISTE Cell

When a large density stratification is no longer a problem in a microgravity environment, one would like to increase the sample size in order to increase the signal-to-noise ratio for a specific heat measurement. To reduce the equilibration time associated with the large sample size, we designed a cylindrical cell containing a stack of plates that separate the bulk fluid into 60 equally thin layers. To understand the thermal behavior of the whole cell, we analyzed the thermal behavior of a 2-D composite system of a cylindrical near-critical fluid layer in contact with a cylindrical copper plate. In this 2-D analysis, the circumference boundary of the two cylindrical layers is subjected to a step temperature change. The solution of this 2-D composite system includes the piston effect that speeds up the equilibration in the near-critical fluid layer and the pure diffusion in the copper plate. The results of this analysis indicate that the characteristic length for the equilibration of the stacked cell is determined by an effective thickness of a single fluid layer instead of the total height of the cylindrical cell.

Weilert, Mark↗
DOE OSTI2022

Classification of animal sounds in a hyperdiverse rainforest using convolutional neural networks with data augmentation

To protect tropical forest biodiversity, we need to be able to detect it reliably, cheaply, and at scale. Automated detection of sound producing animals from passively recorded soundscapes via machine-learning approaches is a promising technique towards this goal, but it is constrained by the necessity of large training data sets. Using soundscapes from a tropical forest in Borneo and a Convolutional Neural Network model (CNN), we investigate i) the minimum viable training data set size for accurate prediction of call types (‘sonotypes’), and ii) the extent to which data augmentation and transfer learning can overcome the issue of small and imbalanced training data sets. We found that even relatively high sample sizes (>80 per sonotype) lead to mediocre accuracy, which however improved significantly with data augmentation and transfer learning, including at extremely small sample sizes (3 per sonotype), regardless of taxonomic group or call characteristics. Neither transfer learning nor data augmentation alone achieved high accuracy. Our results suggest that transfer learning and data augmentation could make the use of CNNs to classify species’ vocalizations feasible even for small soundscape-based projects with many rare species. Retraining our open-source model requires only basic programming skills which makes it possible for individual conservation initiatives to match their local context, in order to enable more evidence-informed management of biodiversity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Plasma GFAP for populational enrichment of clinical trials in preclinical Alzheimer's disease

Abstract INTRODUCTION Cognitively unimpaired (CU) amyloid beta (Aβ)+ individuals with elevated plasma glial fibrillary acidic protein (GFAP) have an increased risk of Alzheimer's disease (AD)‐related progression. We tested the utility of plasma GFAP for population enrichment CU populations in clinical trials. METHODS We estimated longitudinal progression, effect size, and costs of hypothetical clinical trials designed to test an estimated 25% drug effect on reducing tau positron emission tomography (PET) accumulation in the medial temporal lobe (MTL) and temporal neocortical region (NEO‐T). RESULTS CU GFAP+/Aβ+ individuals present an increased annual rate of change and effect size in tau PET MTL and tau PET NEO‐T compared to the other groups. An enrichment strategy selecting CU GFAP+/Aβ+ individuals would require a smaller sample size (≈ 57% reduction) and fewer Aβ PET scans (≈ 74% reduction) than trials enriched with Aβ PET alone, reducing total clinical trial costs by up to 64%. DISCUSSION Our results suggest that clinical trials focusing on preclinical AD recruiting Aβ+ individuals with elevated GFAP levels would improve cost effectiveness. Highlights Cognitively unimpaired (CU) glial fibrillary acidic protein (GFAP)+/amyloid beta (Aβ)+ shows increased changes in tau positron emission tomography (PET) . CU GFAP+/Aβ+ enriched clinical trials require a reduced sample size compared to Aβ+ only. CU GFAP+/Aβ+ enrichment reduces Aβ PET scans required and costs. CU GFAP+/Aβ+ enrichment allows the selection of individuals at early stages of the Alzheimer's disease continuum.

Neurosciences & Neurology↗
DOE OSTI2021

The Initial–Final Mass Relation for Hydrogen-deficient White Dwarfs

The initial–final mass relation represents the total mass lost by a star during the entirety of its evolution from the zero age main sequence to the white-dwarf cooling track. The semiempirical initial–final mass relation (IFMR) is largely based on observations of DA white dwarfs, the most common spectral type of white dwarf and the simplest atmosphere to model. We present a first derivation of the semiempirical IFMR for hydrogen-deficient (non-DA) white dwarfs in open star clusters. We identify a possible discrepancy between the DA and non-DA IFMRs, with non-DA white dwarfs ≈0.07 M {sub ⊙} less massive at a given initial mass. Such a discrepancy is unexpected based on theoretical models of non-DA formation and observations of field white dwarf mass distributions. If real, the discrepancy is likely due to enhanced mass loss during the final thermal pulse and renewed post-AGB evolution of the star. However, we are dubious that the mass discrepancy is physical and instead is due to the small sample size, to systematic issues in model atmospheres of non-DAs, and to the uncertain evolutionary history of Procyon B (spectral type DQZ). A significantly larger sample size is needed to test these assertions. In addition, we also present Monte Carlo models of the correlated errors for DA and non-DA white dwarfs in the initial–final mass plane. We find the uncertainties in initial–final mass determinations for individual white dwarfs can be significantly asymmetric, but the recovered functional form of the IFMR is grossly unaffected by the correlated errors.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

Salvaging Data Records with Missing Data: Data Imputation using the Multivariate t Distribution

When doing multivariate data analysis, one commonobstacle is the presence of incomplete observations, i.e., observationsfor which one or more key fields are blank. Missing datais often countered by deleting entire observations that containmissing data. The negative effects of deleting entire observationsare multiple: deleting observations reduces sample size andcan also result in biased inferences even if data is missing atrandom. In addition, knowledge contained within incompleteobservations is knowledge lost when they are deleted– and theeffort spent collecting that knowledge is effort wasted. Data imputationmethods, or methods of statistically “filling-in” missingdata, can help combat small sample sizes by using the existinginformation in partially complete observations with the end goalof producing less biased and higher confidence inferences. Whena sample from a multivariate normal population is only partiallycomplete, and the missing data meets appropriate assumptions(missing at random), robust data imputation of the missing datacan be implemented with monotone data augmentation (MDA)using the multivariate t distribution.Missing data imputation is applied to data from the NASA InstrumentCost Model (NICM) using the MDA algorithm underthe assumption of having a multivariate t distribution with fixeddegrees of freedom. A sensitivity analysis to the degrees offreedom parameter is presented to demonstrate robustness ofthe multivariate t distribution when dealing with small samplesas compared to the multivariate normal distribution.

DiNicola, Michael↗
DOE OSTI2025

Creating Accurate Methane Emission Inventories through Data-Driven Airborne Survey Strategies

Because natural gas emits less carbon than other fossil fuels, it holds promise as a green energy transition fuel. However, the overall carbon footprint of natural gas is significantly elevated by methane emissions that occur during its production and transmission (Cusworth et al. 2022). Methane “super-emitters,” while comprising only about 1% of sites, are responsible for the majority of oil- and gas-sourced methane emissions, making their detection and mitigation critical in reducing the climate impact of natural gas and in meeting national and global sustainability goals (Sherwin et al. 2024). Yet, despite advancements in detection, significant uncertainties remain regarding the size, frequency, and duration distributions of methane emissions (e.g., Frankenberg et al. 2016, Cusworth et al. 2022, Chen, Sherwin et al. 2022, Conrad et al. 2023, Johnson et al. 2023, Sherwin et al. 2024) underscoring the need for comprehensive emissions inventories segmented by basin across the US. Airborne surveys are well-suited for collecting data to build these comprehensive, basin-level inventories because they allow for extensive spatial coverage, and have the spatial resolution, and the sensitivity to pinpoint individual methane sources. As remote sensing technologies enable rapid basin-scale surveys, it is imperative to establish scientifically and statistically robust standards to generate reliable and actionable emissions inventories. Recent work has shown that differences in airborne sampling strategies, detection technologies, and analysis can lead to large differences between survey conclusions if not correctly accounted for (Chen et al. 2024). This elevates the importance of incorporating proper sampling and analysis techniques when designing a methane emissions monitoring campaign to produce accurate results and facilitate cross-study comparisons. In this paper, we describe a survey strategy designed using the latest conclusions from the literature to align results from different aerial surveys. We identify several sampling and analysis principles, including large sample sizes, balanced sampling across oil and gas production, careful survey area definition, and a unified protocol for analysis, to be vital to producing an unbiased estimate of basin-scale emissions. We present results from a Department of Energy-funded project that deployed this survey strategy in two understudied oil and gas- producing regions in the United States: the Haynesville Basin in Texas and Louisiana, and the Woodford Shale in the Anadarko Basin in Oklahoma.

03 NATURAL GAS↗
DOE OSTI2023

A Multi-Armed Bayesian Ordinal Outcome Utility-Based Sequential Trial with a Pairwise Null Clustering Prior

A multi-armed trial based on ordinal outcomes is proposed that leverages a flexible non-proportional odds cumulative logit model and numerical utility scores for each outcome to determine treatment optimality. This trial design uses a Bayesian clustering prior on the treatment effects that encourages the pairwise null hypothesis of no differences between treatments. A group sequential design is proposed to determine which treatments are clinically different with an adaptive decision boundary that becomes more aggressive as the sample size or clinical significance grows, or the number of active treatments decreases. A simulation study is conducted for 3 and 5 treatment arms, which shows that the design has superior operating characteristics (family wise error rate, generalized power, average sample size) compared to utility designs that do not allow clustering, a frequentist proportional odds model, or a permutation test based on empirical mean utilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hollow fiber mid-ID spectrometer with UV laser ablation sampling for fine spatial reoslution of isotope ratios in solids

We describe a system that combines isotope ratio analysis via mid-infrared (Mid-IR) laser absorption spectroscopy with fine spatial resolution sampling using a UV pulsed laser. The UV laser ablates a pit in a solid on the order of 10 microns in diameter. The sample-derived particulates resulting from laser ablation pass through a micro-combustor, and the resulting gas is analyzed using Mid-IR laser absorption spectroscopy in a capillary absorption spectrometer (CAS). The CAS uses a hollow fiber optic waveguide with a reflective inner coating and a small internal volume on the order of 1ml. The hollow fiber both guides the laser light from source to detector and contains the gas sample at reduced pressure. Near unity overlap between the laser beam and sample enables sensitive analysis with ultra-small sample size. A prototype system has been demonstrated to enable stable carbon isotopic analysis (d13C) with 1 per mil precision using <1 picomole of carbon and is currently being used to study nutrient exchange in soil/root/microbial rhizosphere studies. The smaller sample size of this system is enabling fine spatial resolution analysis (on the order of 10 microns), which is roughly an order of magnitude smaller than was possible with an isotope ratio mass spectrometer (IRMS). In addition to rhizosphere studies, the system can provide a useful tool for fine scale isotope analysis with applications in biological, forensic, and environmental science.

capillary absorption spectroscopy, laser ablation,↗
NASA NTRSDate not supplied

Exploring the Relationship Between Science Funding and Science Productivity on Extended Missions

Introduction: The goal of this summer internship research project was to explore the relationship between the level of funding for science team research and science productivity for planetary spaceflight missions. The productivity and funding for the prime and extended missions were compared when possible. Determining the actual rather than anecdotal relationship between the level of science team funding and their productivity will help determine the adequacy of science team funding. Additionally, determining whether productivity increases or decreases over time (within missions, or within the years between different missions) provides insights into ongoing trends. This project was conducted using data from 11 extended missions within the Planetary Science Division. The main goals of this project were to 1) quantify the level of productivity of the science team during the extended portion of each mission, 2) compare the level of productivity between different types of missions, and 3) compare the level of productivity between the prime missions and extended missions. Methodology: The central metric enters around the “publication rate” of each mission. This rate is normalized in a variety of ways, including the number of publications per year, per Co-I, and per instruments. For the most part, this data was attainable through mission websites, team rosters, and data collected from previous projects. The publication data for this project utilized the SAO/NASA Astrophysics Data System, and it consisted of a search for any papers from the extended mission that had a co-investigator as an author. The time period for these searches also included an additional two years after the end of the extended mission to account for additional time it typically takes to finish data analysis and publish papers. Data from prime missions was collected from a previous project that measured similar metrics for the same missions. When comparing the different types of missions, three factors were used: mission size, dates of the prime and extended missions, and the type of the mission (rover, flyby, orbiter, or lander). To determine mission size, missions were compared by the number of Co-I’s on the mission, the number of instruments in the payload, and the mission class. Findings: On average, prime mission publication rates were generally higher than those of extended missions. Out of the nine missions used for this comparison, six had a higher prime mission publication rate and three had a higher extended mission publication rate. This distinction was interesting to observe, but it should be noted that the methodology for determing the number of published papers differed for the prime missions and the extended missions Discovery missions were found to be the most scientifically productive. Additionally, more recent extended missions tend to possess greater productivity. When comparing missions of different types, it was found that orbiters were much more productive than other types of missions. Future Work: There are a number of opportunities to extend and improve on this project in the future. Even though the metrics used were the same, the prime mission data could be recollected from the same sources as the extended mission data (Astrophysics Data System) to insure greater consistency. The impact level of the publications could also be accounted for, and financial data could be used to compare teams with different levels of funding. Additionally, more missions could be examined to increase the sample size of these comparisons, especially the comparison of different mission types which had an incredibly small sample size for some categories.

Cameron Maxwell Thompson↗
DOE OSTI2023

Adaptive sampling quasi-Newton methods for zeroth-order stochastic optimization

Here, we consider unconstrained stochastic optimization problems with no available gradient information. Such problems arise in settings from derivative-free simulation optimization to reinforcement learning. We propose an adaptive sampling quasi-Newton method where we estimate the gradients using finite differences of stochastic function evaluations within a common random number framework. We develop modified versions of a norm test and an inner product quasi-Newton test to control the sample sizes used in the stochastic approximations and provide global convergence results to the neighborhood of a locally optimal solution. We present numerical experiments on simulation optimization problems to illustrate the performance of the proposed algorithm. When compared with classical zeroth-order stochastic gradient methods, we observe that our strategies of adapting the sample sizes significantly improve performance in terms of the number of stochastic function evaluations required.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

The narrow pulse approximation and long length scale determination in xenon gas diffusion NMR studies of model porous media

We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Pade approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Pade interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Pade length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter. c. 2002 Elsevier Sciences (USA).

NASA Discipline Life Sciences Technologies↗
NASA NTRSDate not supplied

Small-Sample Estimator Decisions – Certainly Uncertain in Human Spaceflight at NASA

Within statistics when estimating means, we rely heavily on the Central Limit Theorem (CLT) to aid in our inference, however in cases where the CLT does not apply we begin to be strongly limited in our choices. If a violation to the CLT comes in the form of small samples and unknown population variance, yet we maintain the assumption of normality, the use of t-distributions is perfectly valid. We explore cases where population normality is not assumed, where sample sizes are small, and where some actionable estimate of the central tendency of a distribution is needed such as for spaceflight operational decisions. Here we identify different scenarios based on measures taken from spaceflight, characterizing performance of estimators of central tendency within our set of simulations. We examine bias and variability of standard estimators of central tendency as they apply to varying small sample sizes and varying population distributions. Further work in this area is needed to develop a framework or set of guidelines for individuals set in these situations. Even with limited information, decisions need to be made.

Central Tendency↗
DOE OSTI2020

Learning Spectral Templates for Photometric Redshift Estimation from Broadband Photometry

Estimating redshifts from broadband photometry is often limited by how accurately we can map the colors of galaxies to an underlying spectral template. Current techniques utilize spectrophotometric samples of galaxies or spectra derived from spectral synthesis models. Both of these approaches have their limitations: either the sample sizes are small and often not representative of the diversity of galaxy colors, or the model colors can be biased (often as a function of wavelength), which introduces systematics in the derived redshifts. In this paper, we learn the underlying spectral energy distributions from an ensemble of ∼100 K galaxies with measured redshifts and colors. We show that we are able to reconstruct emission and absorption lines at a significantly higher resolution than the broadband filters used to measure the photometry for a sample of 20 spectral templates. We find that our training algorithm reduces the fraction of outliers in the derived photometric redshifts by up to 28%, bias up to 91%, and scatter up to 25%, when compared to estimates using a standard set of spectral templates. We discuss the current limitations of this approach and its applicability for recovering the underlying properties of galaxies. Our derived templates and the code used to produce these results are publicly available in a dedicated Github repository: https://github.com/dirac-institute/photoz-template-learning.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Transcriptomics-based Machine Learning Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), has typically limited machine learning (ML) in space studies and further study of radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNAseq) data from 6 mouse liver GeneLab datasets (GLDS) with a total of 113 spaceflight and ground-control samples to determine top features relevant to spaceflight including the effect of radiation exposure. Data was normalized within each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. The top MRMR features were used to predict spaceflight vs. ground-control samples using a Random Forest (RF) classifier with 5-fold cross validation (CV). The ML-based gene sets were further compared against differential gene expression results from individual GLDS. CV training using the top 100 MRMR genes show averages of 86% accuracy and 0.95 AUC value on the validation set over 5 folds (Figure 1A). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 811 or 68 DEGs overlapping between at least 2 or 3 studies, respectively (Figure 1B). Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism. Set analysis between the MRMR features and the DEGs showed 60 or 8 genes overlapping with at least 1 or 2 studies, respectively. MRMR feature selection and ensemble ML methods (e.g. RF) improve performance relative to a Naïve Bayes classifier when NGS data sets are analyzed. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise ratio. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from RNASeq analysis. Non-intersecting sets introduce opportunity to explore spaceflight relevant genes and implementing ML methods across existing NGS datasets may overcome sample size limitations. ML coupled with existing analytical methods enhances understanding of disease by revealing common underlying pathways across datasets.

Machine Learning↗
NASA NTRS2019

Computational Bayesian Methods Applied to Complex Problems in Bio and Astro Statistics

In this dissertation we apply computational Bayesian methods to three distinct problems. In the first chapter, we address the issue of unrealistic covariance matrices used to estimate collision probabilities. We model covariance matrices with a Bayesian Normal-Inverse-Wishart model, which we fit with Gibbs sampling. In the second chapter, we are interested in determining the sample sizes necessary to achieve a particular interval width and establish non-inferiority in the analysis of prevalences using two fallible tests. To this end, we use a third order asymptotic approximation. In the third chapter, we wish to synthesize evidence across multiple domains in measurements taken longitudinally across time, featuring a substantial amount of structurally missing data, and fit the model with Hamiltonian Monte Carlo in a simulation to analyze how estimates of a parameter of interest change across sample sizes.

Elrod, Chris↗
DOE OSTI2021

Subaru Hyper Suprime-Cam revisits the large-scale environmental dependence on galaxy morphology over 360 deg2 at z = 0.3–0.6

This study investigates the role of large-scale environments on the fraction of spiral galaxies at z = 0.3–0.6 sliced to three redshift bins of Δz = 0.1. Here, we sample 276220 massive galaxies in a limited stellar mass of 5 × 1010 solar mass (~M*) over 360 deg2, as obtained from the Second Public Data Release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). By combining projected two-dimensional density information (Shimakawa et al. 2021, MNRAS, 503, 3896) and the CAMIRA cluster catalog (Oguri et al. 2018, PASJ, 70, S20), we investigate the spiral fraction across large-scale overdensities and in the vicinity of red sequence clusters. We adopt transfer learning to reduce the cost of labeling spiral galaxies significantly and then perform stacking analysis across the entire field to overcome the limitations of sample size. Here we employ a morphological classification catalog by the Galaxy Zoo Hubble (Willett et al., 2017, MNRAS, 464, 4176) to train the deep learning model. Based on 74103 sources classified as spirals, we find moderate morphology–density relations on a 10 comoving Mpc scale, thanks to the wide-field coverage of HSC-SSP. Clear deficits of spiral galaxies have also been confirmed, in and around 1136 red sequence clusters. Furthermore, we verify whether there is a large-scale environmental dependence on rest-frame u - r colors of spiral galaxies; such a tendency was not observed in our sample.

Astronomy & Astrophysics↗
DOE OSTI2020

Predicting 3D moisture sorption behavior of materials from 1D investigations

Abstract Moisture in materials can be a source of future outgassing and exacerbate unwanted changes in physical and chemical properties. Here, we investigate the effect of sample size and shape on the moisture transport phenomena through a combined experimental and modeling approach. Several different materials varying in size and shape were investigated over a wide range of relative humidities (0–90%) and temperatures ( $$30{-}70^{\,\circ } \hbox {C}$$ 30 - 70 ∘ C ) using gravimetric type dynamic vapor sorption (DVS). A dynamic triple-mode sorption model, developed previously, was employed to describe the experimental results with good success; the model includes absorption, adsorption, pooling (clustering) of species, and molecular diffusion. Here we show that the full triple-mode sorption model is robust enough to predict the dynamic uptake and outgassing of 3-dimensional (3D) samples using parameters derived from quasi-1D samples. This successful demonstration on three different materials (filled polydimethylsiloxane (PDMS), unfilled PDMS, and ceramic inorganic composite) illustrates that the model is robust at describing the scale-independent physics and chemistry of moisture sorption and diffusion materials. This work demonstrates that while sorption mechanisms manifest in testing of all sample sizes, some of these mechanisms were so subtle that they were overlooked in our initial modeling and assessment, illustrating the importance of multi-scale experiments in the development of robust predictive capabilities. Our study also outlines the challenges and viable solutions for global optimization of a multi-parameter model. The ability to quantify moisture sorption and diffusion, independent of scale, using 1D lab-scale experiments enables prediction of long-term bulk materials behavior in real applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cosmic Reionization on Computers: Statistical Properties of the Distributions of Mean Opacities

Quasar absorption lines provide a unique window to the relationship between galaxies and the intergalactic medium during the Epoch of Reionization. In particular, high redshift quasars enable measurements of the neutral hydrogen content of the universe. However, the limited sample size of observed quasar spectra, particularly at the highest redshifts, hampers our ability to fully characterize the intergalactic medium during this epoch from observations alone. In this work, we characterize the distributions of mean opacities of the intergalactic medium in simulations from the Cosmic Reionization on Computers (CROC) project. We find that the distribution of mean opacities along sightlines follows a non-trivial distribution that cannot be easily approximated by a known distribution. When comparing the cumulative distribution function of mean opacities measurements in subsamples of sample sizes similar to observational measurements from the literature, we find consistency between CROC and observations at redshifts $z\lesssim 5.7$. However, at higher redshifts ($z\gtrsim5.7$), the cumulative distribution function of mean opacities from CROC is notably narrower than those from observed quasar sightlines implying that observations probe a systematically more opaque intergalactic medium at higher redshifts than the intergalactic medium in CROC boxes at these same redshifts. This is consistent with previous analyses that indicate that the universe is reionized too early in CROC simulations.

79 ASTRONOMY AND ASTROPHYSICS↗