Strongly lensed neutral hydrogen emission: detection predictions with current and future radio interferometers
Strong gravitational lensing provides some of the deepest views of the Universe, enabling studies of high-redshift galaxies only possible with next-generation facilities without the lensing phenomenon. To date, 21 cm radio emission from neutral hydrogen has only been detected directly out to , limited by the sensitivity and instantaneous bandwidth of current radio telescopes. We discuss how current and future radio interferometers such as the Square Kilometre Array (SKA) will detect lensed H i emission in individual galaxies at high redshift. Our calculations rely on a semi-analytic galaxy simulation with realistic H i disks (by size, density profile and rotation), in a cosmological context, combined with general relativistic ray tracing. Wide-field, blind H i surveys with the SKA are predicted to be efficient at discovering lensed H i systems, increasingly so at . This will be enabled by the combination of the magnification boosts, the steepness of the H i luminosity function at the high-mass end, and the fact that the H i spectral line is relatively isolated in frequency. These surveys will simultaneously provide a new technique for foreground lens selection and yield the highest redshift H i emission detections. More near term (and existing) cm-wave facilities will push the high redshift H i envelope through targeted surveys of known lenses.
keywords:gravitational lensing: strong, galaxies: evolution, galaxies: ISM
The evolution of neutral cold gas ( K) is fundamental to our understanding of galaxy evolution. Our current cosmological picture is that following recombination at , the Universe was predominately neutral until reionisation occurred between , driven (in poorly constrained proportion) by the onset of star formation and supermassive black hole accretion in the earliest galactic haloes. The intergalactic medium is rapidly ionized about 1 Gyr after the Big Bang, with the vast majority of remaining neutral hydrogen found in galaxies where it is sufficiently shielded and/or replenished (e.g. Barkana & Loeb, 2001; Fan et al. 2006).
Within galactic haloes, cold gas is a significant mass component and evolutionary ingredient: it is the pristine fuel for star formation, as well as the exhaust of supernova explosions. The biggest obstacle in uncovering the precise evolutionary role of H i is that direct H i observations are currently restricted to the modern Universe (). Many important questions remain observationally unanswered as a result: does the cold gas continually condense onto galactic disks, resupplying that which is driven out by AGN and supernovae (e.g. Davé et al. 2013, Morganti et al. 2013)? How do H i morphology and dispersion evolve out to cosmic noon () and beyond? What is the link between the evolution of H i and cosmic star formation history of the Universe (e.g. Obreschkow & Rawlings 2009)? Addressing these questions demands H i observations at look-back times much greater than dynamical times (i.e., at ). Since H i extends to large galactocentric radii, high- H i observations may also enable studies of the halo structure evolution and allow an important measurement of the baryon angular momentum (Obreschkow et al. 2015).
Over the next decade we are set to revolutionize our view of the cosmic evolution of neutral hydrogen. A broad range of radio telescopes are being developed to tune into the rest-frame 21 cm spectral line out to large redshifts. There are a number of upcoming radio facilities, including the Australian Square Kilometre Array Pathfinder (ASKAP; Johnston et al., 2007), MeerKAT (Jonas, 2009), SKA (e.g. Dewdney et al. 2013), as well the upgraded Karl G. Jansky Very Large Array (VLA), Westerbork Synthesis Radio Telescope (WSRT) with the APERture Tile In Focus project (APERTIF, Verheijen et al. 2008), and the Giant Metrewave Radio Telescope (GMRT, Ananthakrishnan, 2005). These telescopes will all perform deep H i surveys, some of which will detect individual galaxies out to redshifts of and beyond. These telescopes will enable the study of H i in galaxies with sample sizes and at distances not practically possible with current facilities. This will renew our understanding of neutral hydrogen and its role in galaxy evolution. While H i absorption studies have made important contributions to this field (e.g. Gupta 2006, Allison et al. 2015), there remain significant sources of uncertainty in absorption systems that limit extrapolation to the cosmic evolution of H i in galaxies (e.g. spin temperature, gas morphology).
One of the byproducts of the significantly larger survey speed and extended redshift range of next-generation instruments will be the dramatically increased probability that an H i-detected source is strongly lensed by a foreground galaxy. Here we use the canonical definition of strongly lensed as a total magnification . This corresponds to a factor saving in integration time to detect a given source, provided it is unresolved in the image plane. This will enable the detection of H i at even greater cosmological distances than otherwise expected with future instruments.
In this letter, we use an N-body simulation with semi-analytic prescriptions for the spatially-extended H i sources to predict the magnification statistics and resultant number counts of strongly lensed H i emission. Our simulation uses a 150 deg field of view mock observing cone out to a redshift limit of , enabling the prediction of detection rates of future H i surveys. We assume cosmological parameter values of = 0.25, = 0.75, and = 73 km s Mpc (Millennium run cosmology, Springel et al. 2005). The best estimate of the cosmological parameters today (e.g. the latest Planck cosmological parameters, Planck collaboration 2015) are slightly different, but this has a marginal difference on both the semi-analytic results (Guo et al., 2013) and magnification estimates, and thus on our predictions.
2.1 Mock sky
The traditional approach in estimating the lensing probability (also known as the lensing opacity) is to use analytic volume densities of both the foreground and background source populations and integrate for a given redshift range. This approach has been readily applied in lensing predictions for point-like sources (e.g. optical quasars, high redshift supernovae, Oguri & Marshall 2010; compact radio AGN, McKean et al. 2015), however, the large spatial extent and non-trivial density-velocity structure of H i disks pose a challenge to estimating total magnification for H i sources. This challenge is further complicated by the requirement of modelling large cosmological volumes ( galaxies) in order to obtain sufficiently converged statistics to support actual galaxy surveys. We address this challenge using the S-SAX model (Obreschkow et al. 2009), a semi-analytic model (SAM) extending the SAMs by Croton et al. (2006) and De Lucia & Blaizot (2007), relying on the Millennium -body simulation of cosmic structure (Springel et al., 2005). S-SAX is the only model to date that provides a mock sky of millions of H i disks with realistic spatial distributions and velocity profiles (Obreschkow et al., 2009b). This model is well-suited to generate H i lensing predictions since the mock galaxies have all the properties needed for both the source and lens galaxies. Explicitly, (i) the mock galaxies are correctly clustered due to the underlying Millennium simulation; (ii) the H i (needed in source-galaxies) has realistic spatial distributions, accounting for the large H i disk size ( stellar effective radius) and central troughs; (iii) the H i has realistic velocity profiles that also satisfy the Tully-Fisher relation and local H i velocity function (Obreschkow et al. 2013); (iv) the galaxies have multiple components (halo+disk+bulge) with realistic mass distributions (needed for lens galaxies), (v) the properties of the mock galaxies evolve with redshift according to the prescriptions of Croton et al. (2006) and Obreschkow et al. (2009). A velocity dispersion of km s is added to the SAM-generated H i velocity profile, consistent with observed properties of nearby H i galaxies (e.g. de Blok et al. 2008). The typical HI velocity dispersions and scale heights potentially increase with redshift, similarly to the observed increase in the molecular velocity dispersion (e.g. Förster Schreiber et al. 2006). We decided to neglect this effect, since the rest-frame channel widths and pixel sizes of our lensing calculations are generally a few times larger than the dispersion scale and thickness of the disk. Crucially, the use of an N-body simulation and full source-to-lens plane ray-tracing enables the H i flux of all lensed galaxies to be split into individual frequency channels, resulting in a realistic prediction of the observables (e.g. peak channel flux, velocity line widths) and hence optimal detection strategies.
2.2 Lens candidate selection
Based on the S-SAX model we produce five statistically independent mock observing cones (as described in Obreschkow et al., 2009b), each covering a field-of-view of 30 deg and extending from to (30 deg corresponds to the cross-section of the Millennium box at ). Together these five cones represent a mock sky covering 150 deg with negligible cosmic variance. From this sky we extract all lens-source pair candidates with a background source integrated H i flux Jy km s and foreground halo virial mass . In addition, we only consider lens candidates with a source redshift and an intrinsic lens-source impact factor arcsec) AND (), where is the foreground halo virial radius. These cuts dramatically decrease the computational demand and have a negligible effect on the inferred lensing statistics (following a trial and error process). For example, the majority of lensing galaxies have Einstein radii (and hence deflection angles) of order a few arcseconds (e.g. Treu 2010, and references therein). The above selection criteria result in lens-source pair candidates, i.e. approximately 2 per arcmin.
2.3 H i magnification factor
The lens models are split into three components: (a) a dark matter halo, (b) bulge, and (c) a disk component. The dark matter halo is assumed to have a spherically-symmetric NFW density profile with a concentration parameter derived from a redshift and halo mass dependent fit (Hennawi et al. 2007), the distribution of which is consistent with empirical results from X-ray and strong lensing observations of galaxy clusters (Comerford & Natarajan 2007). The latter consistency is particularly important since magnification bias will preferentially select high surface mass density foreground dark matter haloes.
The bulge and disk effective radii are derived directly from the self-consistent dynamics of the semi-analytic prescription to all galaxies in the simulation. The bulge is set to a de Vaucoleurs density profile with an ellipticity drawn from a Gaussian distribution with mean and 1 standard deviation of (truncated between ), consistent with observational properties of bulges and previous approaches (Jorgensen & Franx 1994, Möller et al. 2007, Oguri & Marshall 2010). The halo+bulge density profiles used in SAX-Lens are realistic not only in that the distribution of their concentration parameters match available observations; but also in that global H i linewidths (outer rotation) and CO linewidths (inner rotation) are both consistent with observations (Obreschkow et al. 2009b). The disk is modeled as a spherically-symmetric component with an exponential density profile. It is sub-dominant at all radii and adds a marginal fraction of total convergence to the system as seen in Möller et al. (2007). We do not consider the additional convergence contributed by satellite or neighbouring galaxies associated with the foreground lens, which may add up to several tens of percent to some of the largest image separations (and hence magnification, Oguri 2006), however this effect decreases for lower mass, lower redshift lenses as discussed by Möller et al. (2007).
With both source and lens profiles defined, we perform the general-relativistic ray-tracing using the glafic software package (Oguri 2010). To this end, the H i emission of all sources is split into 50 kHz channels (intrinsic velocity widths of km s for ), each of which is individually lensed. We thus have noise-free image plane (lensed) and source plane (unlensed) observing cubes (two spatial dimensions, one frequency dimension) of all H i sources. Fig. 1 shows example source and image plane channel maps of a source (at 200 kHz channel resolution for clarity). The associated lensed and unlensed global profiles are shown in the bottom panel. Optimal detection strategies for such a distorted profile are discussed in Sec. 3.3. The global magnification factor is defined as the ratio between the integrated H i fluxes of the magnified image and the unmagnified source. All systems with are written to a catalogue, referred to as the Semi-Analytic eXtragalalactic Lens (SAX-Lens) catalogue. This catalog contains lens-source pairs, i.e. about 300 per deg.
3.1 Survey-independent lens statistics
The primary output of the methodology laid out in Sec. 2 is the SAX-Lens catalogue of simulated H i sources that are magnified by a factor of 2 or more. In Fig. 2 (top panel) we show the redshift distribution of both the foreground lenses and the background H i sources. This demonstrates the lensing probability only becomes appreciable above , well above the highest redshift H i emission detection to date (, Catinella et al. 2008). The distribution of strongly lensed H i disks with Jy km s peaks at , which corresponds to an observed frequency of MHz. The bottom panel in Fig. 2 demonstrates that higher redshift sources are, on average, more magnified due to higher alignment probability and smaller apparent source size. This means that surveys designed to detect these sources will be highly biased towards strongly lensed systems. It is likely that future radio telescopes able to detect lensed H i sources, will, due to the wide () instantaneous band widths of modern receivers, simultaneously detect the H i emission or absorption of the foreground lens. This will greatly expedite the followup observations of lensed sources discovered with an entirely different selection function to other lens catalogues (e.g. CLASS, Myers et al. 2003; SLACS, Bolton et al. 2006).
In Fig. 3 we show some of the statistical properties of the lens-source pairs, including the foreground lens halo (total) mass, impact parameter (angular distance between lens and source centroids), magnification (colour scale), as well as the H i disk angular size. The ellipse axis ratios correspond to that of the background H i source (with min=0.01). The ellipse position angle represents the angle between the impact vector and H i semi-major axis (see caption). Fig. 3 shows a rough envelope above which virtually no sources are strongly lensed (marginalized over a number of parameters, e.g. redshifts, halo concentration parameter). An approximate fit to this envelope is: impact arcsec.
The highest magnification sources () in the SAX-Lens catalogue have an Einstein ring (see Fig. 1) in at least one of their channel maps and thus have the largest apparent solid angles. However, current and future radio facilities will not have sufficient angular resolution (for the vast majority of their collecting area) at the appropriate frequency to spatially resolve the H i Einstein rings. This maximises their probability of detection since the surveys will benefit from the total H i magnification boost. This also implies that macro lens models derived from optical/NIR/mm observations will be required to constrain the intrinsic H i disk mass and size. Of all the lensed H i disks in the SAX-Lens catalogue, just 6 percent have an extent larger than 0.5 arcsec. This implies that even SKA1-MID, with a arcsec PSF for GHz, will not resolve the vast majority of lensed H i disks. Of course, the sources that are resolved will be of great interest for individual followup to perform more detailed modelling of the astrophysical properties. In this paper we limit the scope to the detection statistics, assuming all sources are unresolved.
3.2 Survey-specific predictions
COSMOS HI Large Extragalactic Survey (Fernandez et al., 2015)
Looking at the Distant Universe with the MeerKAT Array (Holwerda, Blyth & Baker, 2012)
Deep Investigations of Neutral Gas Origins (Meyer, 2009)
Having presented the basic statistical properties of the lensed H i systems, we now turn to predicted detection rates with planned surveys. There are a wide range of H i surveys planned over the course of the next decade which will dramatically expand our view and understanding of H i in galaxies. In Tab. 1 we list the predicted number of H i lenses (and mean redshift) that will be detectable for some of the major future surveys. We define the detection threshold as , where the mean H i profile flux in units of Jy and the noise corresponds to that expected in a channel width equal to the H i profile FWHM for each individual galaxy. The latter are determined from the 100 kHz channel sensitivities listed in Tab. 1. A major caveat in any such prediction is the ability to distinguish lensed from non-lensed sources, given that they will be unresolved (in H i ) and the vast majority will only be lensed by a factor of two (see inset, Fig. 3). Multi-wavelength data will therefore be essential to unambiguously identify low () magnification H i lenses, but may also be used to lower the detection threshold and discover many more lensed systems, as suggested by the number of detections indicated in parentheses.
Table 1 shows that none of the pre-SKA surveys are particularly optimized for detecting large numbers of lensed H i systems (e.g. detections for DINGO and LADUMA). However, lensed H i systems will increase the number counts and enhance these already planned H i surveys beyond their original science cases. The predicted mean redshifts of these lensed H i sources are much higher than the current record of (except for DINGO-D). This combined with the magnification boosts implies that it will be possible to detect high redshift systems with much lower H i mass than that possible for non-lensed systems. Furthermore, the relative isolation of the H i line, the steep high-end H i mass function and the systematically larger magnification factors at higher redshift suggest that some of the most distant H i sources will be detected and identified within a fraction of the total survey duration. See Serjeant (2014) for a discussion on using the high-mass end of the H i mass function to select lenses. The above factors, combined with the increasing H i lens opacity at higher redshift (see Fig. 2, bottom panel) will make a significant impact on measurements at . Indeed, even values with pre-SKA facilities may require strong-lensing correction factors. Tab. 1 also suggests that targeted surveys for lensed H i may be the optimal strategy in the pre-SKA era.
The SKA1 surveys will revolutionize our view of H i in galaxies at high redshift. The predictions suggest lensed detections, well over an order of magnitude larger than the presently known lens population at any wavelength. The mean and maximum redshifts of sources detected in these surveys will provide important insights on the evolution of cold gas towards cosmic noon (), a perspective not possible with any other facility.
The full SKA is expected to detect up to non-lensed H i galaxies out to through surveys aimed primarily at high precision measurements of baryonic acoustic oscillations (Rawlings et al. 2004, Maartens et al. 2015, Bull et al. 2015). However, even with the SKA, strong gravitational lensing will be required to characterize Milky Way like galaxies to the epoch of peak cosmic star formation and AGN activity at . Stacking methods have been proposed to pursue this goal, however, the uncertainties and associated contamination using these methods are unclear. Furthermore, calibration and source modelling limitations may place upper limits on the sensitivity of stacking methods (Makhathini et al., in prep.). Well-modeled strong gravitational lenses, with the benefit of excellent multi-wavelength coverage from ALMA, LSST, JWST and the E-ELT will clearly provide an excellent opportunity to push back this observational frontier. While these surveys will dramatically enhance our knowledge of high redshift cold gas and disk formation, a further valuable contribution will be a new probe of the foreground lens population using a unique selection function that is not biased by dust attenuation (in both the source and lens galaxies).
3.3 Optimal detection strategy
In all of the lensed statistics presented above, the detection thresholds assume a channel width matched to the FWHM of the H i line width. This optimizes the S/N ratio and hence maximizes detection probability for the majority of lensed systems. However, there are cases when a narrow channel width has a significantly larger magnification than that of the total intensity map. This is simply due to the channel flux solid angle and position relative to a caustic. Such differential magnification is seen in the relative magnification of the CO (1-0) channel maps of a lensed disk galaxy (Deane et al. 2013a), as well as the order of magnitude difference in magnification of the radio core to CO emission (Deane et al. 2013b). Differential magnification can be exploited to increase the detection probability by running spectral line source finders with narrow channel widths as well as those optimized to the full galaxy line width. An illustration of this effect is seen in Fig. 1 which shows a H i source with total magnification , however one particular 50 kHz channel is magnified by a factor . An appropriate channel selection yields a higher probability of detection than matching to the line width. Clearly, there is a great deal of optimization to be explored in both targeted and blind lensed H i surveys, which we investigate in an upcoming paper.
The use of an N-body simulation with applied semi-analytics enables the following predictions relevant to lensed Hi.
A full ray-tracing analysis of lens candidates in the 530 deg simulation reveal 46,000 strong lenses () with unlensed velocity-integrated H i fluxes Jy km s.
The predicted number counts imply that the dramatic increase in sensitivity and frequency coverage of upcoming radio facilities will enable H i detection at redshifts where strong lensing has been highly beneficial to galaxy evolution studies that utilise other spectral lines (e.g. CO, [C ii], [O iii]) as well as continuum observations.
Planned H i surveys (LADUMA, DINGO) on future telescopes will detect some () lenses, which in the case of LADUMA may require consideration when determining . However, none of these surveys are particularly optimal for detecting lensed systems. Lensing does however enable the detection much lower H i mass systems at high redshift, enhancing the scientific output of planned surveys.
Current telescopes (e.g. GMRT, VLA) fitted with their new broadband receivers, and most definitely future telescopes (e.g. MeerKAT), are very well placed to detect lensed H i in targeted surveys and expand our view of cold neutral gas to higher redshift.
The depth and area of the SKA1-MID surveys will lead to efficient lens selection at high redshift. Overall, these surveys are predicted to detect over an order of magnitude more H i lenses () than the presently known lens population at any wavelength.
The large solid angles covered by H i disks leads to significant differential magnification in some lensed sources. Appropriately optimized three-dimensional source finders may therefore enhance the detection rates presented here.
We thank Masamune Oguri for making the glafic software package publicly available and for providing very helpful advice. We also thank the anonymous referee for their comments. The financial assistance of the South African SKA Project (SKA SA) towards this research is hereby acknowledged. Opinions expressed and conclusions arrived at are those of the author and are not necessarily to be attributed to the SKA SA. www.ska.ac.za.
-  Allison, J. R., et al., 2015, arXiv:1503.01265
-  Ananthakrishnan, S., 2005, International Cosmic Ray Conference, 10, 125
-  Barkana, R., & Loeb, A., 2001, Physics Reports, 349, 125
-  Bolton, A. S., et al., 2006, ApJ, 638, 703
-  Bull, P., et al., 2015, arXiv:1501.04088
-  Catinella, B., et al., 2008, ApJL, 685, L13
-  Comerford, J. M., & Natarajan, P., 2007, MNRAS, 379, 190
-  Croton, D. J., et al., 2006, MNRAS, 365, 11
-  Davé, R., et al., 2013, MNRAS, 434, 2645
-  de Blok, W. J. G., et al., 2008, AJ, 136, 2648
-  De Lucia, G., & Blaizot, J., 2007, MNRAS, 375, 2
-  Deane, R. P., et al., 2013, MNRAS, 434, 23
-  Deane, R. P., et al., 2013, MNRAS, 434, 3322
-  Fan, X., Carilli, C. L., & Keating, B., 2006, ARA&A, 44, 415
-  Fernandez, X., van Gorkom, J. H., Momjian, E., & Chiles Team 2015, American Astronomical Society Meeting Abstracts, 225, #427.03
-  Förster Schreiber, N. M., et al., 2006, ApJ, 645, 1062
-  Guo, Q., et al., 2013, MNRAS, 428, 1351
-  Hennawi, J. F., et al., 2007, ApJ, 654, 714
-  Holwerda, B. W., Blyth, S.-L., & Baker, A. J., 2012, IAU Symposium, 284, 496
-  Johnston, S., Bailes, M., Bartel, N., et al., 2007, PASA, 24, 174
-  Jonas, J. L. 2009, IEEE Proceedings, 97, 1522
-  Jorgensen, I., & Franx, M., 1994, ApJ, 433, 553
-  Möller, O., Kitzbichler, M., & Natarajan, P., 2007, MNRAS, 379, 1195
-  Maartens, R., et al., 2015, arXiv:1501.04076
-  McKean, J. P., et al., 2015, arXiv:1502.03362
-  Meyer, M. 2009, Panoramic Radio Astronomy: Wide-field 1-2 GHz Research on Galaxy Evolution, 15
-  Morganti, R., et al., 2013, Science, 341, 1082
-  Myers, S. T., et al., 2003, MNRAS, 341, 1
-  Obreschkow, D., & Rawlings, S., 2009, ApJL, 696, L129
-  Obreschkow, D., et al., 2009, ApJ, 698, 1467
-  Obreschkow, D., et al., 2009, ApJ, 703, 1890
-  Obreschkow, D., et al., 2015, arXiv:1501.01048
-  Obreschkow, D., et al., 2013, ApJ, 766, 137
-  Oguri, M., 2006, MNRAS, 367, 1241
-  Oguri, M., 2010, PASJ, 62, 1017
-  Oguri, M., & Marshall, P. J., 2010, MNRAS, 405, 2579
-  Planck Collaboration, et al., 2014, A&A, 571, AA16
-  Rawlings, S., et al., 2004, New Astronomy Reviews, 48, 1013
-  Serjeant, S., 2014, ApJL, 793, L10
-  Springel, V., et al., 2005, Nature, 435, 629
-  Treu, T., 2010, ARA&A, 48, 87
-  Verheijen, M. A. W., et al., 2008, The Evolution of Galaxies Through the Neutral Hydrogen Window, 1035, 265