SEARCH · Engineering Papers

Results for “Astronomy software”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗

DOE OSTI2013

DETECTION OF FAST TRANSIENTS WITH RADIO INTERFEROMETRIC ARRAYS

Next-generation radio arrays, including the Square Kilometre Array (SKA) and its pathfinders, will open up new avenues for exciting transient science at radio wavelengths. Their innovative designs, comprising a large number of small elements, pose several challenges in digital processing and optimal observing strategies. The Giant Metre-wave Radio Telescope (GMRT) presents an excellent test-bed for developing and validating suitable observing modes and strategies for transient experiments with future arrays. Here we describe the first phase of the ongoing development of a transient detection system for GMRT that is planned to eventually function in a commensal mode with other observing programs. It capitalizes on the GMRT's interferometric and sub-array capabilities, and the versatility of a new software backend. We outline considerations in the plan and design of transient exploration programs with interferometric arrays, and describe a pilot survey that was undertaken to aid in the development of algorithms and associated analysis software. This survey was conducted at 325 and 610 MHz, and covered 360 deg{sup 2} of the sky with short dwell times. It provides large volumes of real data that can be used to test the efficacies of various algorithms and observing strategies applicable for transient detection. We present examples that illustrate the methodologies of detecting short-duration transients, including the use of sub-arrays for higher resilience to spurious events of terrestrial origin, localization of candidate events via imaging, and the use of a phased array for improved signal detection and confirmation. In addition to demonstrating applications of interferometric arrays for fast transient exploration, our efforts mark important steps in the roadmap toward SKA-era science.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

The Hyper Suprime-Cam software pipeline

Here in this article, we describe the optical imaging data processing pipeline developed for the Subaru Telescope’s Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope’s Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high-level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

JAXtronomy: A JAX port of lenstronomy

Gravitational lensing is a phenomenon where light bends around massive objects, resulting in distorted images seen by an observer. Studying gravitationally lensed systems provides insights into cosmology and astrophysics, including constraints of the expansion rate of the Universe and the distribution of dark matter. Thus, we introduce JAXtronomy, a re-implementation of the gravitational lensing software package lenstronomy (Birrer, 2021; Birrer & Amara, 2018) using JAX (Bradbury et al., 2018). JAX is a Python library that uses an accelerated linear algebra (XLA) compiler to improve the performance of computing software. Our core design principle of JAXtronomy is to maintain an identical API to that of lenstronomy. The main JAX features utilized in JAXtronomy are just-in-time compilation, which can lead to significant reductions in execution time, and automatic differentiation, which allows for the implementation of gradient-based algorithms that were previously impossible. Additionally, JAX allows code to be run on GPUs or parallelized across CPU cores, further boosting the performance of JAXtronomy.

astronomy↗

DOE OSTI2022

PETSc TSAdjoint: A Discrete Adjoint ODE Solver for First-Order and Second-Order Sensitivity Analysis

Here, we present a new software system PETSc TSAdjoint for first-order and second order adjoint sensitivity analysis of time-dependent nonlinear differential equations. The derivative calculation in PETSc TSAdjoint is essentially a high-level algorithmic differentiation process. The adjoint models are derived by differentiating the timestepping algorithms and implementing them based on the parallel infrastructure in PETSc. Full differentiation of the library code, including MPI routines, is avoided, and users do not need to derive their own adjoint models for their specific applications. PETSc TSAdjoint can compute the first-order derivative, that is, the gradient of a scalar functional, and the Hessian-vector product, which carries second-order derivative information, while requiring minimal input (a few callbacks) from the users. The adjoint model employs optimal checkpointing schemes in a manner that is transparent to users. Finally, usability, efficiency, and scalability are demonstrated through examples from a variety of applications.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The Global Dynamical Atlas of the Milky Way Mergers: Constraints from Gaia EDR3–based Orbits of Globular Clusters, Stellar Streams, and Satellite Galaxies

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams, and satellite galaxies in action (J) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy (E) to enhance the contrast between the groupings. For this endeavor, we use Gaia EDR3–based measurements of 170 globular clusters, 41 streams, and 46 satellites to derive their J and E. To detect groups, we use the ENLINK software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of N = 6 groups, including the previously known mergers Sagittarius, Cetus, Gaia–Sausage/Enceladus, LMS-1/Wukong, Arjuna/Sequoia/I'itoi, and one new merger that we call Pontus. All of these mergers, together, comprise 62 objects (≈25% of our sample). We discuss their members, orbital properties, and metallicity distributions. We find that the three most-metal-poor streams of our galaxy—"C-19" ([Fe/H] = –3.4 dex), "Sylgr" ([Fe/H] = –2.9 dex), and "Phoenix" ([Fe/H] = –2.7 dex)—are associated with LMS-1/Wukong, showing it to be the most-metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

MAJOR-MERGER GALAXY PAIRS AT Z = 0: DUST PROPERTIES AND COMPANION MORPHOLOGY

We present an analysis of dust properties of a sample of close major-merger galaxy pairs selected by K {sub s} magnitude and redshift. The pairs represent the two populations of spiral–spiral (S+S) and mixed morphology spiral–elliptical (S+E). The Code Investigating GALaxy Emission software is used to fit dust models to the Two Micron All Sky Survey, Wide-Field Infrared Survey Explorer , and Herschel flux density measurements, and to derive the parameters describing the polycyclic aromatic hydrocarbons contribution, interstellar radiation field, and photodissociation regions. Model fits verify our previous Spitzer Space Telescope analysis that S+S and S+E pairs do not have the same level of enhancement of star formation and differ in dust composition. The spirals of mixed-morphology galaxy pairs do not exhibit the enhancements in interstellar radiation field and therefore dust temperature for spirals in S+S pairs in contrast to what would be expected according to standard models of gas redistribution due to encounter torques. This suggests the importance of the companion environment/morphology in determining the dust properties of a spiral galaxy in a close major-merger pair.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The gMOSS: the galaxy survey and galaxy populations of the large homogeneous field

ABSTRACT We present the gMOSS (Galaxies of Medium-band One-meter Schmidt telescope Survey) catalogue of ∼19 000 galaxies in 20 filters (4 broad-band SDSS and 16 medium-band filters). We observed 2.386 deg2 on the central part of the HS47.5-22 field with the 1-m Schmidt telescope of the Byurakan Astrophysical Observatory. The gMOSS is a complete flux-limited sample of galaxies with a threshold magnitude of r SDSS ≤ 22.5 AB. From photometric measurements with 16 medium-band filters and u SDSS, we get spectral energy distributions for each object in the field, which are used for further analysis. Galaxy classification and photometric redshift estimation based on spectral template matching with zebra software. The obtained redshift accuracy is σNMAD < 0.0043. Using the SED-fitting cigale code, we obtained the main properties of the stellar population of galaxies, such as rest-frame (u − r)res colour, stellar mass, extinction, and mass-weighted age with a precision of 0.16 ± 0.07 mag, 0.14 ± 0.04 dex, 0.27 ± 0.1 mag, and 0.08 ± 0.04 dex, respectively. Using a dust-corrected colour–mass diagram, we divided the full sample into populations of red and blue galaxies and considered the dependencies between stellar mass and age. Throughout cosmic time, red sequence galaxies remain older and more massive than blue cloud galaxies. The star formation history of a complete subsample of galaxies selected in the redshift range 0.05 ≤ z ≤ 0.015 with <$\mathrm{log} M \mathrm{\gt }_\mathrm{[M_\odot ]}$>8.3 shows an increase in the SFRD up to z ∼ 3, under the results obtained in earlier studies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Phase-curve Pollution of Exoplanet Transit Depths

The next generation of space telescopes will enable transformative science to understand the nature and origin of exoplanets. In particular, transit spectroscopy will reveal the chemical composition of the exoplanet atmospheres with unprecedented detail thanks to precise measurements of the visible-to-infrared transit depths down to 10 parts per million. Such a level of instrumental precision raises the challenge to obtain even more precise astrophysical models so as not to significantly influence the interpretation of the observed data. We must therefore critically revisit some of the commonly accepted assumptions that were adequate for analyzing past and current observations. A common approximation in the analysis of exoplanetary primary transits is that the planet does not contribute to the recorded flux, so-called dark planet hypothesis. In this paper, we investigate the impact of the dark planet hypothesis on the parameters obtained from the analysis of transits with particular attention to the transit depth. We develop mathematical formulae and release new software to estimate the magnitude of the potential bias. These tools will be useful in the preparation of observing proposals, as well as within the scientific consortia of the James Webb Space Telescope (JWST) and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (ARIEL) missions. We probe the accuracy of the mathematical formulae through the analysis of synthetic observations with the JWST Mid-InfraRed Instrument. We find that self-blending from nightside emission attenuates the transit depth by >3σ for some of the known exoplanet systems, in agreement with previous work. An additional unreported effect caused by the nightside rotating into view can also impart a significant effect, but in the opposite direction (increasing the transit depth); this effect can largely be removed with conventional detrending practices, at the expense of a slight increase in noise, and mixing astrophysical variations and instrumental drifts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

VARIABILITY OF DISK EMISSION IN PRE-MAIN SEQUENCE AND RELATED STARS. III. EXPLORING STRUCTURAL CHANGES IN THE PRE-TRANSITIONAL DISK IN HD 169142

We present near-IR (NIR) and far-UV observations of the pre-transitional (gapped) disk in HD 169142 using NASA's Infrared Telescope Facility and Hubble Space Telescope. The combination of our data along with existing data sets into the broadband spectral energy distribution reveals variability of up to 45% between ∼1.5-10 μm over a maximum timescale of 10 yr. All observations known to us separate into two distinct states corresponding to a high near-IR state in the pre-2000 epoch and a low state in the post-2000 epoch, indicating activity within the ≲1 AU region of the disk. Through analysis of the Pa β and Br γ lines in our data we derive a mass accretion rate in 2013 May of M-dot ≈ (1.5-2.7) × 10{sup –9} M {sub ☉} yr{sup –1}. We present a theoretical modeling analysis of the disk in HD 169142 using Monte-Carlo radiative transfer simulation software to explore the conditions and perhaps signs of planetary formation in our collection of 24 yr of observations. We find that shifting the outer edge (r ≈ 0.3 AU) of the inner disk by 0.05 AU toward the star (in simulation of accretion and/or sculpting by forming planets) successfully reproduces the shift in NIR flux. We establish that the ∼40-70 AU dark ring imaged in the NIR by Quanz et al. and Momose et al. and at 7 mm by Osorio et al. may be reproduced with a 30% scaled density profile throughout the region, strengthening the link to this structure being dynamically cleared by one or more planetary mass bodies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel’dovich measurements from BOSS CMASS and LOWZ halos

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel’dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope DR5 and Planck in combination with the CMASS (mean redshift $\langle z \rangle$ = 0.55 and host halo mass $\langle M_{vir} \rangle$ = 3 × 10 13 M ⊙ ) and LOWZ ( $\langle z \rangle$ = 0.31 , $\langle M_{vir} \rangle$ = 5 × 10 13 M ⊙ ) galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothesis at 6.5σ, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally > 90 σ for the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e., the electron thermal pressure profile of the same CMASS objects, and reject the no-tSZ hypothesis at 10σ. We combine tSZ and kSZ measurements to estimate the electron temperature to 20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ) signal at 2.9 (13.9) σ , and leave their interpretation to future work. This paper and the companion paper demonstrate that current CMB experiments can detect and resolve gas profiles in low mass halos and at high redshifts, which are the most sensitive to feedback in galaxy formation and the most difficult to measure any other way. They will be a crucial input to cosmological hydrodynamical simulations, thus improving our understanding of galaxy formation. These precise gas profiles are already sufficient to reduce the main limiting theoretical systematic in galaxy-galaxy lensing: baryonic uncertainties. Future such measurements will thus unleash the statistical power of weak lensing from the Rubin, Euclid and Roman observatories. Finally, our stacking software ThumbStackis publicly available and directly applicable to future Simons Observatory and CMB-S4 data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

APOGEE-2S Mg–Al anti-correlation of the metal-poor globular cluster NGC 2298

We present detailed elemental abundances and radial velocities of stars in the metal-poor globular cluster (GC) NGC 2298, based on near-infrared high-resolution ($R$ ~ 22 500) spectra of 12 members obtained during the second phase of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) at Las Campanas Observatory as part of the seventeenth Data Release (DR 17) of the Sloan Digital Sky Survey IV (SDSS-IV). We employed the Brussels Automatic Code for Characterizing High accuracy Spectra (BACCHUS) software to investigate abundances for a variety of species including α elements (Mg, Si, and Ca), the odd-Z element Al, and iron-peak elements (Fe and Ni) located in the innermost regions of NGC 2298. We find a mean and median metallicity [Fe/H] = -1.76 and -1.75, respectively, with a star-to-star spread of 0.14 dex, which is compatible with the internal measurement errors. Thus, we find no evidence for an intrinsic [Fe/H] abundance spread in NGC 2298. The typical $α$-element enrichment in NGC 2298 is overabundant relative to the Sun, and it follows the trend of other metal-poor GCs. We confirm the existence of an Al-enhanced population in this cluster, which is clearly anti-correlated with Mg, indicating the prevalence of the multiple-population phenomenon in NGC 2298.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

A Proper Motions Study of the Globular Cluster NGC 3201

With a high value of heliocentric radial velocity, a retrograde orbit, and suspected to have an extragalactic origin, NGC 3201 is an interesting globular cluster for kinematical studies. Our purpose is to calculate the relative proper motions (PMs) and membership probability for the stars in the wide region of globular cluster NGC 3201. PM based membership probabilities are used to isolate the cluster sample from the field stars. The membership catalog will help address the question of chemical inhomogeneity in the cluster. Archive CCD data taken with a wide-field imager (WFI) mounted on the ESO 2.2 m telescope are reduced using the high-precision astrometric software developed by Anderson et al. for the WFI images. The epoch gap between the two observational runs is ∼14.3 years. To standardize the BVI photometry, Stetson’s secondary standard stars are used. The CCD data with an epoch gap of ∼14.3 years enables us to decontaminate the cluster stars from field stars efficiently. The median precision of PMs is better than ∼0.8 mas yr{sup −1} for stars having V<18 mag that increases up to ∼1.5 mas yr{sup −1} for stars with 18<V<20 mag. Kinematic membership probabilities are calculated using PMs for stars brighter than V∼20 mag. An electronic catalog of positions, relative PMs, BVI magnitudes, and membership probabilities in the ∼19.7 × 17 arcmin{sup 2} region of NGC 3201 is presented. We use our membership catalog to identify probable cluster members among the known variables and X-ray sources in the direction of NGC 3201.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

PBjam: A Python Package for Automating Asteroseismology of Solar-like Oscillators

Asteroseismology is an exceptional tool for studying stars using the properties of observed modes of oscillation. So far the process of performing an asteroseismic analysis of a star has remained somewhat esoteric and inaccessible to nonexperts. In this software paper we describe PBjam, an open-source Python package for analyzing the frequency spectra of solar-like oscillators in a simple but principled and automated way. The aim of PBjam is to provide a set of easy-to-use tools to extract information about the radial and quadropole oscillations in stars that oscillate like the Sun, which may then be used to infer bulk properties such as stellar mass, radius, age, or even structure. Asteroseismology and its data analysis methods are becoming increasingly important as space-based photometric observatories are producing a wealth of new data, allowing asteroseismology to be applied in a wide range of contexts such as exoplanet, stellar structure and evolution, and Galactic population studies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE THEORETICAL ASTROPHYSICAL OBSERVATORY: CLOUD-BASED MOCK GALAXY CATALOGS

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires expert knowledge of galaxy modeling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalog suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO’s features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

CosTuuM: Polarized Thermal Dust Emission by Magnetically Oriented Spheroidal Grains

We present the new open-source C++-based Python library CosTuuM that can be used to generate infrared absorption and emission coefficients for arbitrary mixtures of spheroidal dust grains that are (partially) aligned with a magnetic field. We outline the algorithms underlying the software, demonstrate the accuracy of our results using benchmarks from literature, and use our tool to investigate some commonly used approximative recipes. We find that the linear polarization fraction for a partially aligned dust grain mixture can be accurately represented by an appropriate linear combination of perfectly aligned grains and grains that are randomly oriented, but that the commonly used picket fence alignment breaks down for short wavelengths. We also find that for a fixed dust grain size, the absorption coefficients and linear polarization fraction for a realistic mixture of grains with various shapes cannot both be accurately represented by a single representative grain with a fixed shape, but that instead an average over an appropriate shape distribution should be used. Insufficient knowledge of an appropriate shape distribution is the main obstacle in obtaining accurate optical properties. CosTuuM is available as a standalone Python library and can be used to generate optical properties to be used in radiative transfer applications.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Atmospheric Transmission Spectroscopy of Hot Jupiter KELT-10b using Synthetic Telluric Correction Software

High-resolution spectroscopic visible data were obtained with the Ultraviolet and Visible Echelle Spectrograph on the Very Large Telescope. Our goal was to analyze the data in an effort to detect the presence of sodium in the atmosphere of hot Jupiter exoplanet KELT-10b, as well as characterize the orbit of the planet via the Rossiter-McLaughlin effect. Eighty spectra were collected during a single transit of KELT-10b. After standard spectroscopic calibration using ESO-Reflex, the synthetic telluric modeling software molecfit was applied to remove terrestrial atmospheric effects, and to refine the wavelength calibration. Sodium is recognized by its characteristic absorption doublet located at 5895.924 and 5889.951 Å, which can be seen in the planet atmosphere transmission spectrum and through excess absorption during the transit. The radial velocity of the host star was analyzed by measuring the average shift of absorption features from spectrum to spectrum. Our results indicate a sodium detection in the planet transmission spectrum with a line contrast of 0.66% and 0.43% ± 0.09% for the sodium DII and DI lines, respectively. Excess absorption measurements agree to within one half combined standard deviation between the planet transmission spectrum (0.143% ± 0.020%, a 7σ detection) and during the time series (0.124% ± 0.034%, a 3.6σ detection) in a band 1.25 Å wide. The wavelength grid corrections provided by molecfit were insufficient to determine radial velocities and measure the Rossiter-McLaughlin effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Measurements of the z > 5 Lyman-α forest flux autocorrelation functions from the extended XQR-30 data set

We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, $\langle{F}\rangle$. This increase may come from decreasing $\langle{F}\rangle$ or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data.

79 ASTRONOMY AND ASTROPHYSICS↗