Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRedshifts for galaxies in radio continuum surveys from Bayesian model fitting of HI 21-cm lines
65
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce a new Bayesian HI spectral line fitting technique capable of obtaining spectroscopic redshifts for millions of galaxies in radio surveys with the Square Kilometere Array (SKA). This technique is especially well-suited to the low signal-to-noise regime that the redshifted 21-cm HI emission line is expected to be observed in, especially with SKA Phase 1, allowing for robust source detection. After selecting a set of continuum objects relevant to large, cosmological-scale surveys with the first phase of the SKA dish array (SKA1-MID), we simulate data corresponding to their HI line emission as observed by the same telescope. We then use the MultiNest nested sampling code to find the best-fitting parametrised line profile, providing us with a full joint posterior probability distribution for the galaxy properties, including redshift. This provides high quality redshifts, with redshift errors $Delta z / z <10^{-5}$, from radio data alone for some 1.8 million galaxies in a representative 5000 square degree survey with the SKA1-MID instrument with up-to-date sensitivity profiles. Interestingly, we find that the SNR definition commonly used in forecast papers does not correlate well with the actual detectability of an HI line using our method. We further detail how our method could be improved with per-object priors and how it may be also used to give robust constraints on other observables such as the HI mass function. We also make our line fitting code publicly available for application to other data sets.
rate research
Read More
We present results from a search for 21 cm associated HI absorption in a sample of 29 radio sources selected from the Australia Telescope 20 GHz survey. Observations were conducted using the Australia Telescope Compact Array Broadband Backend, with which we can simultaneously look for 21 cm absorption in a redshift range of 0.04 < z < 0.08, with a velocity resolution of 7 km/s . In preparation for future large-scale H I absorption surveys we test a spectral-line finding method based on Bayesian inference. We use this to assign significance to our detections and to determine the best-fitting number of spectral-line components. We find that the automated spectral-line search is limited by residuals in the continuum, both from the band-pass calibration and spectral-ripple subtraction, at spectral-line widths of Deltav_FWHM > 103 km/s . Using this technique we detect two new absorbers and a third, previously known, yielding a 10 per cent detection rate. Of the detections, the spectral-line profiles are consistent with the theory that we are seeing different orientations of the absorbing gas, in both the host galaxy and circumnuclear disc, with respect to our line-of-sight to the source. In order to spatially resolve the spectral-line components in the two new detections, and so verify this conclusion, we require further high-resolution 21 cm observations (~0.01 arcsec) using very long baseline interferometry.
Maximally Smooth Functions (MSFs) are a form of constrained functions in which there are no inflection points or zero crossings in high order derivatives. Consequently, they have applications to signal recovery in experiments where signals of interest are expected to be non-smooth features masked by larger smooth signals or foregrounds. They can also act as a powerful tool for diagnosing the presence of systematics. The constrained nature of MSFs makes fitting these functions a non-trivial task. We introduce maxsmooth, an open source package that uses quadratic programming to rapidly fit MSFs. We demonstrate the efficiency and reliability of maxsmooth by comparison to commonly used fitting routines and show that we can reduce the fitting time by approximately two orders of magnitude. We introduce and implement with maxsmooth Partially Smooth Functions, which are useful for describing elements of non-smooth structure in foregrounds. This work has been motivated by the problem of foreground modelling in 21-cm cosmology. We discuss applications of maxsmooth to 21-cm cosmology and highlight this with examples using data from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) and the Large-aperture Experiment to Detect the Dark Ages (LEDA) experiments. We demonstrate the presence of a sinusoidal systematic in the EDGES data with a log-evidence difference of $86.19pm0.12$ when compared to a pure foreground fit. MSFs are applied to data from LEDA for the first time in this paper and we identify the presence of sinusoidal systematics. maxsmooth is pip installable and available for download at: https://github.com/htjb/maxsmooth
Observations of the EoR with the 21-cm hyperfine emission of neutral hydrogen (HI) promise to open an entirely new window onto the formation of the first stars, galaxies and accreting black holes. In order to characterize the weak 21-cm signal, we need to develop imaging techniques which can reconstruct the extended emission very precisely. Here, we present an inversion technique for LOFAR baselines at NCP, based on a Bayesian formalism with optimal spatial regularization, which is used to reconstruct the diffuse foreground map directly from the simulated visibility data. We notice the spatial regularization de-noises the images to a large extent, allowing one to recover the 21-cm power-spectrum over a considerable $k_{perp}-k_{para}$ space in the range of $0.03,{rm Mpc^{-1}}<k_{perp}<0.19,{rm Mpc^{-1}}$ and $0.14,{rm Mpc^{-1}}<k_{para}<0.35,{rm Mpc^{-1}}$ without subtracting the noise power-spectrum. We find that, in combination with using the GMCA, a non-parametric foreground removal technique, we can mostly recover the spherically average power-spectrum within $2sigma$ statistical fluctuations for an input Gaussian random rms noise level of $60 , {rm mK}$ in the maps after 600 hrs of integration over a $10 , {rm MHz}$ bandwidth.
We use the results of previous work building a halo model formalism for the distribution of neutral hydrogen, along with experimental parameters of future radio facilities, to place forecasts on astrophysical and cosmological parameters from next generation surveys. We consider 21 cm intensity mapping surveys conducted using the BINGO, CHIME, FAST, TianLai, MeerKAT and SKA experimental configurations. We work with the 5-parameter cosmological dataset of {$Omega_m, sigma_8, h, n_s, Omega_b$} assuming a flat $Lambda$CDM model, and the astrophysical parameters {$v_{c,0}, beta$} which represent the cutoff and slope of the HI- halo mass relation. We explore (i) quantifying the effects of the astrophysics on the recovery of the cosmological parameters, (ii) the dependence of the cosmological forecasts on the details of the astrophysical parametrization, and (iii) the improvement of the constraints on probing smaller scales in the HI power spectrum. For an SKA I MID intensity mapping survey alone, probing scales up to $ell_{rm max} = 1000$, we find a factor of $1.1 - 1.3$ broadening in the constraints on $Omega_b$ and $Omega_m$, and of $2.4 - 2.6$ on $h$, $n_s$ and $sigma_8$, if we marginalize over astrophysical parameters without any priors. However, even the prior information coming from the present knowledge of the astrophysics largely alleviates this broadening. These findings do not change significantly on considering an extended HIHM relation, illustrating the robustness of the results to the choice of the astrophysical parametrization. Probing scales up to $ell_{rm max} = 2000$ improves the constraints by factors of 1.5-1.8. The forecasts improve on increasing the number of tomographic redshift bins, saturating, in many cases, with 4 - 5 redshift bins. We also forecast constraints for intensity mapping with other experiments, and draw similar conclusions.
We present the results of a survey for intervening HI 21-cm absorbers at intermediate and low redshift (0<z<1.2). For our total sample of 24 systems, we obtained high quality data for 17 systems, the other seven being severely affected by radio frequency interference (RFI). Five of our targets are low redshift (z<0.17) optical galaxies with small impact parameters (<20 kpc) toward radio-bright background sources. Two of these were detected in 21-cm absorption, showing narrow, high optical depth absorption profiles, the narrowest having a velocity dispersion of only 1.5 km/s, which puts an upper limit on the kinetic temperature of T_k<270 K. Combining our observations with results from the literature, we measure a weak anti-correlation between impact parameter and integral optical depth in local (z<0.5) 21-cm absorbers. Of eleven CaII and MgII systems searched, two were detected in 21-cm absorption, and six were affected by RFI to a level that precludes a detection. For these two systems at z~0.6 we measure spin temperatures of T_s=(65+/-17) K and T_s>180 K. A subset of our systems were also searched for OH absorption, but no detections were made.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
