-
Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference
Authors:
Teng Hu,
Vikram Khaire,
Joseph F. Hennawi,
Michael Walther,
Hector Hiss,
Justin Alsing,
Jose Oñorbe,
Zarija Lukic,
Frederick Davies
Abstract:
We present a new approach to measure the power-law temperature density relationship $T=T_0 (ρ/ \barρ)^{γ-1}$ and the UV background photoionization rate $Γ_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$α$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape…
▽ More
We present a new approach to measure the power-law temperature density relationship $T=T_0 (ρ/ \barρ)^{γ-1}$ and the UV background photoionization rate $Γ_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$α$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape of the $b$-$N_{\rm HI}$ distribution is sensitive to the IGM thermal parameters $T_0$ and $γ$, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density d$N$/d$z$, and we demonstrate that precise constraints can also be obtained on $Γ_{\rm HI}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b$-$N_{\rm HI}$ distribution on IGM parameters trained on an ensemble of 624 Nyx hydrodynamical simulations at $z = 0.1$, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS datasets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly$α$ forest spectra at $z\simeq 0.1$, one can measure the thermal and ionization state of the IGM with very high precision ($σ_{\log T_0} \sim 0.08$ dex, $σ_γ\sim 0.06$, and $σ_{\log Γ_{\rm HI}} \sim 0.07$ dex).
△ Less
Submitted 14 July, 2022;
originally announced July 2022.
-
Improving IGM temperature constraints using wavelet analysis on high-redshift quasars
Authors:
Molly Wolfson,
Joseph F. Hennawi,
Frederick B. Davies,
Jose Oñorbe,
Hector Hiss,
Zarija Lukić
Abstract:
The thermal state of the intergalactic medium (IGM) contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young universe. This thermal state is encoded in the small-scale structure of Lyman-$α$ (Ly$α$) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines.…
▽ More
The thermal state of the intergalactic medium (IGM) contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young universe. This thermal state is encoded in the small-scale structure of Lyman-$α$ (Ly$α$) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines. At high redshifts, where the opacity is large, averaging mixes high signal-to-noise ratio transmission spikes with noisy absorption troughs. Wavelet amplitudes are an alternate statistic that maintains spatial information while quantifying fluctuations at the same spatial frequencies as the power spectrum, giving them the potential to more sensitively measure the small-scale structure. Previous Ly$α$ forest studies using wavelet amplitude probability density functions (PDFs) used limited spatial frequencies and neglected strong correlations between PDF bins and across wavelets scales, resulting in sub-optimal and unreliable parameter inference. Here we present a novel method for performing statistical inference using wavelet amplitude PDFs that spans the full range of spatial frequencies probed by the power spectrum and that fully accounts for these correlations. We applied this procedure to realistic mock data drawn from a simple thermal model parameterized by the temperature at mean density, $T_0$, and find that wavelets deliver 1$σ$ constraints on $T_0$ that are on average 7% more sensitive at $z=5$ (12% at $z=6$) than those from the power spectrum. We consider the possibility of combing wavelet PDFs with the power, but find that this does not lead to improved sensitivity.
△ Less
Submitted 6 October, 2021;
originally announced October 2021.
-
A Novel Statistical Method for Measuring the Temperature-Density Relation in the IGM Using the $b$-$N_{\text{HI}}$ Distribution of Absorbers in the Ly$α$ Forest
Authors:
Hector Hiss,
Michael Walther,
José Oñorbe,
Joseph F. Hennawi
Abstract:
We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$α$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0 (ρ/ρ_0)^{γ-1}$, and previous work has inferred $T_0$ and $γ$ by fitting its low-$b$ cutoff. This approach di…
▽ More
We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$α$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0 (ρ/ρ_0)^{γ-1}$, and previous work has inferred $T_0$ and $γ$ by fitting its low-$b$ cutoff. This approach discards the majority of available data, and is susceptible to systematics related to cutoff determination. We present a method that exploits all information encoded in the $b$-$N_{\text{HI}}$ distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, then use principal component decomposition to create an emulator which can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D $(T_0, γ)$ measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at $z=2$. Using 200 absorbers, equivalent in pathlength to a single Ly$α$ forest spectrum, we measure $\log T_0 =4.092^{+0.050}_{-0.055}$ and $γ=1.49^{+0.073}_{-0.074}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of $\log T_0$ and $γ$ with precision on $\log T_0$ ($γ$) nearly two (three) times higher for current dataset sizes.
△ Less
Submitted 28 March, 2019;
originally announced March 2019.
-
A New Measurement of the Temperature Density Relation of the IGM From Voigt Profile Fitting
Authors:
Hector Hiss,
Michael Walther,
Joseph F. Hennawi,
José Oñorbe,
John M. O'Meara,
Alberto Rorai
Abstract:
We decompose the Lyman-α (Lyα) forest of an extensive sample of 74 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-$b$ cutoff in the $b$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff a…
▽ More
We decompose the Lyman-α (Lyα) forest of an extensive sample of 74 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-$b$ cutoff in the $b$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range $2.0\leq z \leq 3.4$, which allows us to measure the evolution of the IGM temperature-density ($T= T_0 (ρ/ ρ_0)^{γ-1}$) relation parameters $T_0$ and $γ$. We calibrate our measurements against Ly$α$ forest simulations, using 21 different thermal models of the IGM at each redshift, also allowing for different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of $T_0$ shows a suggestive peak at $z=2.8$, while our evolution of $γ$ is consistent with $γ\simeq 1.4$ and disfavors inverted temperature-density relations. Our measured evolution of $T_0$ and $γ$ are generally in good agreement with previous determinations in the literature. Both the peak in the evolution of $T_0$ at $z = 2.8$, as well as the high temperatures $T_0\simeq 15000-20000\,$K that we observe at $2.4 < z < 3.4$, strongly suggest that a significant episode of heating occurred after the end of HI reionization, which was most likely the cosmic reionization of HeII.
△ Less
Submitted 16 August, 2018; v1 submitted 2 October, 2017;
originally announced October 2017.
-
A New Precision Measurement of the Small-Scale Line-of-Sight Power Spectrum of the Ly$α$ Forest
Authors:
Michael Walther,
Joseph F. Hennawi,
Hector Hiss,
Jose Oñorbe,
Khee-Gan Lee,
Alberto Rorai,
John O'Meara
Abstract:
We present a new measurement of the Ly$α$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise, and corrects for the bias induced by masking missing data, DLAs, and metal absorption lines. Our measureme…
▽ More
We present a new measurement of the Ly$α$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise, and corrects for the bias induced by masking missing data, DLAs, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes $k > 0.02\,\mathrm{s\,km^{-1}}$, unaccessible to previous SDSS/BOSS analyses. It is well known that these high-$k$ modes are highly sensitive to the thermal state of the intergalactic medium, however contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power, and find a modest effect on modes with $k < 0.1\,\mathrm{s\,km^{-1}} $. As a result, by masking metals and restricting to $k < 0.1\,\mathrm{s\,km^{-1}}$ their impact is completely mitigated. We present an end-to-end Bayesian forward modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Ly$α$ forest simulations. These mocks are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform MCMC fits. Our results agree well with BOSS on scales $k < 0.02\,\mathrm{s\,km^{-1}}$ where the measurements overlap. The combination of BOSS' percent level low-$k$ precision with our $5-15\%$ high-$k$ measurements, results in a powerful new dataset for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables.
△ Less
Submitted 30 November, 2017; v1 submitted 21 September, 2017;
originally announced September 2017.
-
Multi-epoch, high spatial resolution observations of multiple T Tauri systems
Authors:
Gergely Csépány,
Mario van den Ancker,
Péter Ábrahám,
Rainer Köhler,
Wolfgang Brandner,
Felix Hormuth,
Hector Hiss
Abstract:
Context. In multiple pre-main-sequence systems the lifetime of circumstellar disks appears to be shorter than around single stars, and the actual dissipation process may depend on the binary parameters of the systems. Aims. We report high spatial resolution observations of multiple T Tauri systems at optical and infrared wavelengths. We determine if the components are gravitationally bound and orb…
▽ More
Context. In multiple pre-main-sequence systems the lifetime of circumstellar disks appears to be shorter than around single stars, and the actual dissipation process may depend on the binary parameters of the systems. Aims. We report high spatial resolution observations of multiple T Tauri systems at optical and infrared wavelengths. We determine if the components are gravitationally bound and orbital motion is visible, derive orbital parameters and investigate possible correlations between the binary parameters and disk states. Methods. We selected 18 T Tau multiple systems (16 binary and two triple systems, yielding $16 + 2\times2=20$ binary pairs) in the Taurus-Auriga star forming region from the survey by Leinert et al. (1993), with spectral types from K1 to M5 and separations from 0.22" (31 AU) to 5.8" (814 AU). We analysed data acquired in 2006-07 at Calar Alto using the AstraLux lucky imaging system, along with data from SPHERE and NACO at the VLT, and from the literature. Results. We found ten pairs to orbit each other, five pairs that may show orbital motion and five likely common proper motion pairs. We found no obvious correlation between the stellar parameters and binary configuration. The 10 $μ$m infra-red excess varies between 0.1 and 7.2 magnitudes (similar to the distribution in single stars, where it is between 1.7 and 9.1), implying that the presence of the binary star does not greatly influence the emission from the inner disk. Conclusions. We have detected orbital motion in young T Tauri systems over a timescale of $\approx20$ years. Further observations with even longer temporal baseline will provide crucial information on the dynamics of these young stellar systems.
△ Less
Submitted 27 February, 2017;
originally announced February 2017.