Joint Bayesian Estimation of Quasar Continua and the Lyman-Alpha Forest Flux Probability Distribution Function

Abstract

We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma -1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq 7\%$ and $\simeq 10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma -1$ of the IGM temperature-density relation with a precision of $\pm 8.6\%$ at $z=3$, $\pm 6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.

arXiv comment21 pages (+ Appendix), accepted at ApJ

References in corpus (60)

• emcee: The MCMC Hammer
• Equation of State of the Photoionized Intergalactic Medium
• The Low Redshift Lyman Alpha Forest in Cold Dark Matter Cosmologies
• Gravitational Collapse of Small-Scale Structure as the Origin of the Lyman Alpha Forest
• The Lyman-alpha Forest Power Spectrum from the Sloan Digital Sky Survey
• Recovery of the Power Spectrum of Mass Fluctuations from Observations of the Lyman-alpha Forest
• Probing the Universe with the Lyman-alpha Forest: I. Hydrodynamics of the Low Density IGM
• The thermal history of the intergalactic medium
• The Lyman Alpha Forest from Gravitational Collapse in the CDM+Lambda Model
• The Physics of the Intergalactic Medium
• The Observed Probability Distribution Function, Power Spectrum, and Correlation Function of the Transmitted Flux in the Lyman-alpha Forest
• The Evolution of the Intergalactic Medium
• Evolution of Structure in the Intergalactic Medium and the Nature of the Ly-alpha Forest
• Nyx: A Massively Parallel AMR Code for Computational Cosmology
• The Cosmological Impact of Luminous TeV Blazars I: Implications of Plasma Instabilities for the Intergalactic Magnetic Field and Extragalactic Gamma-Ray Background
• HeII Reionization and its Effect on the IGM
• The Power Spectrum of Mass Fluctuations Measured from the Lyman-alpha Forest at Redshift z=2.5
• A Direct Precision Measurement of the Intergalactic Lyman-alpha Opacity at 2<z<4.2
• A strong redshift dependence of the broad absorption line quasar fraction
• The Evolution of Optical Depth in the Ly-alpha Forest: Evidence Against Reionization at z~6
• The Thermal Memory of Reionization History
• A Refined Measurement of the Mean Transmitted Flux in the Ly-alpha Forest over 2 < z < 5 Using Composite Quasar Spectra
• Measurements of the UV background at 4.6 < z < 6.4 using the quasar proximity effect
• A Measurement of the Temperature-Density Relation in the Intergalactic Medium Using a New Lyman-alpha Absorption Line Fitting Method
• A Measurement of Small Scale Structure in the 2.2 < z < 4.2 Lyman-alpha Forest
• Possible evidence for an inverted temperature-density relation in the intergalactic medium from the flux distribution of the Lyman-alpha forest
• Constraints from the Lyman alpha forest power spectrum
• Reionization and high-redshift galaxies: the view from quasar absorption lines
• An improved measurement of the flux distribution of the Ly-alpha forest in QSO absorption spectra: the effect of continuum fitting, metal contamination and noise properties
• A Principal Component Analysis of quasar UV spectra at z~3
• Galactic and Circumgalactic OVI and its Impact on the Cosmological Metal and Baryon Budgets at 2<z<3.5
• Self-Consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations
• An unbiased measurement of the UV background and its evolution via the proximity effect in quasar spectra
• A consistent determination of the temperature of the intergalactic medium at redshift z=2.4
• The Temperature-Density Relation in the Intergalactic Medium at Redshift <z>=2.4
• On the intergalactic temperature-density relation
• Cosmological and astrophysical constraints from the Lyman-alpha forest flux probability distribution function
• Characterizing the Pressure Smoothing Scale of the Intergalactic Medium
• The First Data Release of the KODIAQ Survey
• Predicting QSO Continua in the Ly Alpha Forest
• IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function
• The intergalactic medium thermal history at redshift z=1.7--3.2 from the Lyman alpha forest: a comparison of measurements using wavelets and the flux distribution
• The Cosmological Impact of Luminous TeV Blazars II: Rewriting the Thermal History of the Intergalactic Medium
• The imprint of inhomogeneous HeII reionization on the HI and HeII Ly-alpha forest
• Mean-flux Regulated PCA Continuum Fitting of SDSS Lyman-alpha Forest Spectra
• The Lyman-alpha forest in a blazar-heated Universe
• The Cosmological Impact of Luminous TeV Blazars III: Implications for Galaxy Clusters and the Formation of Dwarf Galaxies
• Inferring the Linear Power Spectrum from the Lyman-alpha Forest
• A New Method to Directly Measure the Jeans Scale of the Intergalactic Medium Using Close Quasar Pairs
• The Lyman alpha forest flux probability distribution at z>3
• A wavelet analysis of QSO spectra
• Inhomogeneous Helium Reionization and the Equation of State of the Intergalactic Medium
• Sample variance and Lyman-alpha forest transmission statistics
• Exploring the Thermal State of the Low-Density Intergalactic Medium at z=3 with an Ultra-High Signal-to-Noise QSO Spectrum
• The impact of helium reionization on the structure of the intergalactic medium
• Relative Flux Calibration of Keck HIRES Echelle Spectra
• The Probability Distribution of the Lya transmitted flux from a sample of SDSS quasars
• A data-driven model for spectra: Finding double redshifts in the Sloan Digital Sky Survey
• Systematic Continuum Errors in the Lyman-Alpha Forest and The Measured Temperature-Density Relation
• Quasars as probes of cosmological reionization

Cited by in corpus (13)

• The Opacity of the Intergalactic Medium Measured Along Quasar Sightlines at $z\sim 6$
• Detecting and Characterizing Young Quasars I: Systemic Redshifts and Proximity Zones Measurements
• LATIS: The Ly$\alpha$ Tomography IMACS Survey
• Detecting and Characterizing Young Quasars II: Four Quasars at $z\sim 6$ with Lifetimes $<10^4$ years
• A New Method to Measure the Post-Reionization Ionizing Background from the Joint Distribution of Lyman-$\alpha$ and Lyman-$\beta$ Forest Transmission
• Inferring high redshift large-scale structure dynamics from the Lyman-alpha forest
• Anomaly in the Opacity of the Post-Reionization Intergalactic Medium in the Ly$\alpha$ and Ly$\beta$ Forest
• How proper are Bayesian models in the astronomical literature?
• Quasar Factor Analysis -- An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis
• Improving IGM temperature constraints using wavelet analysis on high-redshift quasars
• Forecasting constraints on the mean free path of ionizing photons at $z\ge 5.4$ from the Lyman-$\alpha$ forest flux auto-correlation function
• The Spectroscopic Signature of Variability in High-Redshift Quasars
• Analytic marginalization of absorption line continua