Astro2020 Science White Paper The hidden circumgalactic medium
Thematic Areas: Planetary Systems Star and Planet Formation Formation and Evolution of Compact Objects Cosmology and Fundamental Physics Stars and Stellar Evolution Resolved Stellar Populations and their Environments Galaxy Evolution Multi-Messenger Astronomy and Astrophysics
Name: Claudia Cicone Institution: INAF - Osservatorio astronomico di Brera, Milano (Italy) Email: email@example.com Phone: +39 02 72320381
Co-authors: Carlos De Breuck (ESO); Chian-Chou Chen (ESO); Eelco van Kampen (ESO); Desika Narayanan (U.Florida); Tony Mroczkowski (ESO); Paola Andreani (ESO); Pamela Klaassen (UK ATC); Axel Weiss (MPIfR); Kotaro Kohno (U.Tokyo); Jens Kauffmann (MIT); Jeff Wagg (SKAO); Dominik Riechers (Cornell/MPIA); Bitten Gullberg (Durham); James Geach (U.Hertfordshire); Sijing Shen (U.Oslo); J. Colin Hill (IAS/CCA); Simcha Brownson (U.Cambridge).
Abstract: The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.
1 Scientific context and motivation
Galaxies grow and evolve within a network of dark matter through a finely-tuned exchange of baryons with the cosmic web, a process that is challenging to understand and model. The importance of describing baryons, their interplay and behaviour in and out of galaxies is two-fold. On the one hand, if our research were restricted to the cosmological framework, we would not be able to comprehend how the bright components of galaxies - the stars - formed out of dense gas clouds, hence how planetary systems originated, and eventually how life began. On the other hand, since the only observables against which we can test theoretical predictions are produced by baryons, a proper understanding of the complex mechanisms involving baryonic matter is essential to interpret the observable data and so test our cosmological paradigm. Hence, understanding the baryon cycle is an endeavour at the forefront of modern astrophysics, and it is relevant to many different fields, from Galactic studies to Cosmology.
The baryon cycle problem can be visualised as follows: galaxies are embedded in large gas haloes that we define “circumgalactic medium (CGM)”, extending by hundreds of kpc (up to the virial radius). A small portion of CGM material, the interstellar medium (ISM), has already cooled and flowed towards the central regions of galaxies, and it is directly involved in the processes of star formation and accretion onto supermassive black holes (SMBHs). However, most of the CGM mass lies outside the ISM realm and represents the interface between the ISM and the matter distributed outside galaxies, the intergalactic medium (IGM). The CGM is fed - from the inside - by galactic outflows and fountains (Biernacki+Teyssier18, ; Costa+18, ) - and, from the outside - by cosmic streams (Keres+05, ; Dekel+09, ) and galaxy mergers (Hani+18, ). Furthermore, the CGM gets depleted via accretion onto galaxies, and through expulsion into the IGM by energetic feedback processes (Silk+Rees98, ; DiMatteo+05, ; Hopkins+14, ). The CGM is therefore intrinsically linked to the feedback and accretion processes in galaxies, and carries the imprint of their evolution.
Although traditionally the CGM is thought to be dominated by hot and warm gas (Tumlinson+17, ), with little or no cold gas beyond the ISM scales, this picture has been challenged by recent breakthroughs in this field. On the one hand, the lack of cold gas in most simulated galaxy haloes appears to be a bias of hydrodynamical simulations rather than a theoretical prediction. More specifically, several studies are now finding that increasing the spatial resolution on halo scales (i.e. beyond the star forming ISM) naturally leads to an enhanced cold gas fraction in the CGM, because it allows to properly capture cooling processes and to spatially resolve smaller cold clumps embedded in a warmer medium (Suresh+19, ; Hummels+18, ; vandeVoort+19, ). On the other hand, a massive cold ( K) CGM reservoir may result from cosmic feeding, which at high redshift is dominated by cold filamentary inflows, but also from powerful feedback. Indeed, recent studies have shown that intense episodes of star formation and active galactic nuclei (AGN) can drive massive () outflows of Hi and H gas. These can extend by several kpc (Cicone+14, ; Cicone+18b, ) and embed very dense ( cm) clouds (Aalto+12, ; Aalto+15, ). However, although H outflows can be powerful enough to clear the central regions of galaxies within a few Myr, they are likely incapable of expelling all of the molecular ISM from the halo, both in starbursts (Leroy+15, ) and in luminous quasars (Fluetsch+19, ; Alatalo15, ). Instead, cumulative effects of such outflows may be: (i) the transportation of material to CGM scales, where most of the (metal-enriched) H and Hi gas stalls without escaping the halo (Costa+18, ; Shen+13, ), and (ii) the enhancement of gas cooling on halo scales (Biernacki+Teyssier18, ; Prieto+17, ; Costa+15, ). Therefore, we expect the CGM to be multiphase and entrain not only diffuse ionised gas, but also colder and denser clumps. Such cold and non-pristine gas represents fuel that, once (re)-accreted onto the ISM, would be readily available to feed star formation, and hence it would have a crucial role for galaxy evolution.
As we argue in this paper, only (sub-)millimeter ((sub-)mm) and far-infrared (FIR) observations can probe the multiphase CGM including its coldest component, and so deliver a complete census of baryons in galaxy haloes. Such measurements are however extremely challenging even for ALMA because they require high sensitivity to large-scale emission and wide field of views (FoVs) which only a large (e.g. 50m) single-dish can provide.
2 State-of-the-art and forthcoming experiments
Currently, we have very little observational constraints on the mass and physical properties of cold gas on CGM scales. The CGM of high-z (radio quiet) AGNs has so far been investigated predominantly in its diffuse ionised and atomic phases through absorption features detected in rest-frame UV spectra. Interestingly, some of these works hinted at a higher cold gas fraction in the CGM of luminous quasars at compared to inactive galaxies at similar redshifts and with similar masses, especially at large projected velocities of km s (Johnson+15, ). Recently, the ubiquitous detection of low surface brightness kpc-size Ly nebulae, permitted by instruments such as MUSE at the VLT and KCWI at the Keck telescope, has enabled direct imaging of the CGM at (Arrigoni-Battaia+19, ; Wisotzki+18, ; Borisova+16, ). The properties of Ly haloes around quasars suggest entrainment of cold ( K) and dense (density, cm) clumps up to kpc in these sources (Arrigoni-Battaia+15, ). However, the absorption of Ly photons and the resonant nature of Ly emission prevent a detailed kinematical study of the CGM using this tracer. Furthermore, the presence of a ionising source (such as a quasar) can significantly boost the Ly flux, compared to what is expected, due to the ubiquitous fluorescent signal stimulated by the high UV background (Cantalupo+12, ). Therefore, Ly imaging cannot provide an unbiased census of the high-z CGM.
The ionised and neutral CGM, including any cold and dense H gas, can be directly imaged through FIR fine structure lines (such as [Cii] 158m, [Ci] 609m & 370m, [Nii] 205m & 122m, [Oiii] 88m, etc) and molecular rotational transitions (e.g. CO, HCN, CN, CS). At high-, FIR lines are redshifted into the (sub-)mm atmospheric windows, allowing ground-based observations. Thanks to the bright and multiphase (Hii, Hi, H) nature of [Cii] (Pineda+13, ), sensitive [Cii] line observations hold the potential to directly image faint and diffuse CGM components at K. Cicone et al. (Cicone+15, ) demonstrated the existence of copious amounts of [Cii]-emitting gas in the CGM of a luminous quasar at , which was also found to host the most extended ( kpc) high-velocity ( km s) outflow of cold gas known. By modelling the uv visibilities, 70% of the [Cii] flux at low projected velocities was found to trace a 20 kpc-scale nebula, while only 30% of [Cii] arises from the compact ISM. This result completely overturned the scenario favoured by earlier interferometric [Cii] observations of the same source, which, lacking short antenna spacings, resolved out all of the extended emission.
The cold CGM may be dense and bright enough to be detected also in CO line emission. The so-called ‘Spiderweb’ galaxy, a protocluster at with a central bright radio galaxy, hosts a 70 kpc-size halo of molecular gas (Emonts+16, ; Emonts+18, ), which does not arise from individual galaxies but from the extended CGM. Surprisingly, its carbon abundance and excitation properties are similar to a normal star forming disk, hence supporting a ‘recycled’ origin of the CGM - due to a combination of outflows, mass transfer among galaxies, gas accretion, and mergers. Interestingly, the CO(1-0) emission seen in the Spiderweb galaxy is so diffuse that it was resolved out by the JVLA, while being detected at high significance by the shorter baselines of ATCA (see also Emonts+18_ngVLA, ). Molecular haloes may not be unique to powerful quasars. Ginolfi et al. (Ginolfi+17, ) discovered a 40 kpc-size CO(4-3) structure around a massive galaxy at , where most of the molecular line emission is not associated with the ISM of the main galaxy or its satellites. These pilot studies suggest that FIR and (sub-)mm line observations of the high-z CGM carried out with a high sensitivity to large-scale structures hold a huge discovery potential.
The variety of (sub-)mm and FIR lines would allow us not only to detect any hidden massive CGM reservoir, but also to constrain its physical properties such as density, temperature, ionised-to-neutral gas fraction, ionisation parameter, and metallicity (Fernandez-Ontiveros+16, ; Bethermin+16, ). Moreover, the unique chemical properties of some molecular ions (OH, CH, HO) make them suitable to explore physical processes that can critically affect galaxy growth, such as turbulence and shocks (Rangwala+14, ; Gonzalez-Alfonso+18, ). These mechanisms are strictly associated to the baryon cycle (feedback through outflows and feeding by inflows) and so are particularly relevant to CGM gas. For example, ALMA imaging of CH lines, detectable in absorption from diffuse molecular gas (Godard+12, ), has recently revealed hidden massive reservoirs of highly turbulent cool ( K) gas, extending to kpc, far beyond the kpc-size ISM of (lensed) dusty star forming galaxies at (Falgarone+17, ). The same observations detected extremely broad (FWZI) CH emission tracing very dense gas () entrained in powerful galactic outflows, which are probably feeding the large-scale diluted turbulent H CGM reservoir seen in absorption.
In parallel, high sensitivity, wide FoV, broadband continuum observations will open interesting new avenues to the investigation of the hidden CGM. Bright high-z FIR/(sub-)mm continuum emitters can be used as background sources to detect diffuse, low column density CGM gas in absorption through virtually any molecular tracer (e.g. CO, HCN, HCO, OH, OH, CH, (Muller+06, ; Combes08, )). Such experiments require a large FoV to sample several continuum emitters and so overcome the likely low filling factor of the CGM absorbers. They also require high spectral resolution, since molecular line absorbers usually have line widths of km s.
Finally, sensitive continuum observations will enable the detection of the thermal Sunyaev-Zeldovich (tSZ) effect from the warm/hot phase of the CGM. The tSZ signal, being proportional to the integral of electron pressure, provides calorimetry of the ionised gas (e.g. Mroczkowski2019, ; Mroczkowski19b, ). This, in turn, provides strong constraints on galaxy formation models, whose predictions span a wide range of tSZ signals depending on the AGN feedback prescriptions (see e.g. (Battaglia2017, ), and white paper by Battaglia & Hill et al.). In addition, the integrated tSZ signal can give an estimate of the total mass of a system. Stacking of Planck, ACT, and SPT data (Greco2015, ; Ruan2015, ) has allowed to reach down to galaxy mass scales, with typical fluxes of a few tens of Jy111For example, a total mass of produces a negative continuum signal of Jy at .. Thus tSZ measurements of the CGM require good sensitivity and likely next generation (sub-)mm instruments.
3 Opportunities to make a big step forward
Improving over ALMA: the quest for a large (sub)-millimetre single dish:
The cold neutral CGM, and in particular the molecular CGM, is an almost completely unexplored avenue, due to the lack of (sub-)mm facilities that are able to detect and trace large-scale, low surface brightness structures. The latter are indeed filtered out by interferometers, as explained below, and undetected by current single dish facilities because of their low sensitivity (e.g. the 12-meter APEX telescope or the ALMA Total Power antennas). Larger dishes such as the LMT 50-m or the IRAM 30-m telescopes have a limited frequency coverage, imposed by the atmospheric conditions at the telescope site as well as by the antenna surface accuracy which make observations at GHz almost prohibitive.
ALMA has transformed our understanding of galaxies at both high and low redshift, but there are limitations to its capabilities. Due to spatial scale filtering, an interferometer only detects a fraction of the total flux density for sources with emission on scales larger than its shortest antenna spacing (baseline), even in mosaicked observations. The maximum recoverable scale (MRS) of an interferometer is , where is the shortest baseline (close to the dish diameter), although the exact MRS is best determined through simulations because it depends on both the array configuration and the details of the large-scale structure. For the ALMA 12-m array, at 850 GHz ( for the ACA 7-m array) (ALMATec, ). An MRS of corresponds to a largest detectable angular scale of pc at the distance of M 82 (, distance of 3 Mpc), and of kpc at , hence preventing the study of not only the CGM, but also of any large-scale ISM structure (such as a galactic outflow or a tidal tail) in nearby galaxies. However, at higher-, a corresponds to kpc at , which in principle allows some CGM science with ALMA (see 2) but only limited to much smaller structures compared to those detected in Ly emission (with sizes of kpc, e.g. (Arrigoni-Battaia+19, )).
To further illustrate the advantage of a large single dish (SD) antenna vs ALMA in observing extended sources, we show here a sensitivity comparison. More specifically, we consider the ALMA 50 element array of 12-m diameter antennas, and compare it with a 50-m diameter SD antenna hypothetically placed near the ALMA site. The rms noise for an interferometric observation scales as , where is the noise contribution of the atmosphere and receiver, is the effective collecting area of each element, is the integration time, is the spectral bandwidth, and is the number of baselines, which for array elements equals to (for example, for the ALMA 12-m array and antennas, ). For a SD telescope, equipped with a heterodyne focal plane array (an instrument that delivers spectro-imaging information, i.e. 3D datacubes), the relation is similar, but we redefine to be the number of independent beams on the sky, so that . The effective area of an antenna scales with its geometric area, hence , and so matches for , if one assumes the same , , and . However, one should consider that the FoV of a 50m-antenna is smaller by a factor of than the FoV of ALMA. Hence, the sensitivities per FoV of ALMA and of a 50-m SD become comparable when the source area fills at least beams. At 850 GHz (resolution of 1.45) and for typical beam spacings (i.e. 2f, f=focal ratio), this is true for sources across, hence at larger sizes the mapping speed of the 50-m SD outperforms the full ALMA array. Furthermore, we note that a SD telescope is much easier to keep at the forefront of upgraded technologies compared to a 50-element array. For instance, the new generation heterodyne receivers mounted on the APEX telescope already provide instantaneous spectral bandwidths (per polarization and per pixel) that are twice as large as the current ALMA bandwidth.
Figure 1 compares the capability of such hypothetical 50-m SD telescope (see, e.g. Bertoldi2018 (); deBreuck2018 (); Hargrave2018 (); Klaassen2018 (); Mroczkowski2018 ()) with that of the ACA and ALMA arrays in mapping large-scale emission from the cold CGM at . The input image is drawn from the state-of-the-art SIMBA cosmological simulations with a volume (Dave+19, ). We selected a star-forming () disk galaxy with stellar mass , and total CGM gas mass . The equilibrium molecular line luminosities are calculated simultaneously by employing the thermal-chemical-radiative equilibrium code DESPOTIC (Krumholz14, ). We then projected the particle information onto a uniform 1024 grid, and calculated the line-of-sight velocity-integrated intensities. The latter have been re-normalised to 1 Jy km s, since our purpose is to compare how much flux density is recovered by the different facilities in hours. As shown in Fig. 1, the molecular line flux is spread across a large projected area of , corresponding to 80 kpc at , which is not recovered by the interferometers. Clearly, for such low-z galaxies, ALMA and ACA are not optimised to probe the CGM because of (i) their small beam (even in the most compact configuration) which requires large mosaic observations to probe the full source extent (implying a lower integration time per pointing), and of the (ii) spatial scale filtering of large-scale emission, which affects even mosaic observations. The best results are obtained with the 50-m SD antenna, which we have conservatively assumed to be equipped with a single-pixel ALMA detector (but technology advancements will soon deliver much more sensitive multi-pixel detectors hence significantly shortening the SD observing time, (Mroczkowski2018, )).
Of course, the most important factor to consider for ALMA vs a 50-m SD in the Atacama desert is their complementarity: on the one hand, as shown by Fig. 1, ALMA will never access scales large enough to constrain the full extent - and mass - of the CGM - for this, single dish mapping is required. On the other hand, a 50-m SD will be limited to spatial resolution at 850 GHz. Only joint SD and ALMA mapping can thus provide a complete view of the CGM over a broad range of spatial scales and cosmic times.
Synergies with astronomical facilities planned beyond 2030s:
Capabilities to map low surface brightness Hi 21cm emission from galaxies up to high-z, realised by (near-)future projects including the SKA (SKA, ) and the ngVLA (ngVLA, ), will be highly complementary to the (sub-)mm-based studies of the denser CGM described here. However, until the full development of the SKA, FIR tracers of Hi gas such as the [C ii] line will remain the only way to probe such phase at . Promising avenues to detect the CGM at will be through [Cii] line intensity mapping surveys, to be carried out, e.g., with CCAT-prime ((Stacey+18, ), see also Kovetz et al. A2020 WP), TIME (Crites+17, ), or CONCERTO (Lagache+18, ). These experiments will detect [C ii] both from the host and from the diffuse CGM, without resolving individual structures, hence requiring higher-resolution follow ups. The CGM is also one of the main science goals of Athena (2030s), which will trace the warm/hot phase, of the Simons Observatory (SimonsObs, ) and CMB-S4 (CMB_S4, ) (2020s), which will utilize the thermal and kinematic SZ effects to probe ionized gas, and of LUVOIR (2030s), which will detect the CGM in absorption in the UV band.
- (1) Pawel Biernacki and Romain Teyssier. The combined effect of AGN and supernovae feedback in launching massive molecular outflows in high-redshift galaxies. MNRAS, 475:5688–5703, Apr 2018. arXiv:1712.02794, doi:10.1093/mnras/sty216.
- (2) Tiago Costa, Joakim Rosdahl, Debora Sijacki, and Martin G. Haehnelt. Driving gas shells with radiation pressure on dust in radiation-hydrodynamic simulations. MNRAS, 473:4197–4219, Jan 2018. arXiv:1703.05766, doi:10.1093/mnras/stx2598.
- (3) Dušan Kereš, Neal Katz, David H. Weinberg, and Romeel Davé. How do galaxies get their gas? MNRAS, 363:2–28, Oct 2005. arXiv:astro-ph/0407095, doi:10.1111/j.1365-2966.2005.09451.x.
- (4) A. Dekel, Y. Birnboim, G. Engel, J. Freundlich, T. Goerdt, M. Mumcuoglu, et al. Cold streams in early massive hot haloes as the main mode of galaxy formation. Nature, 457:451–454, Jan 2009. arXiv:0808.0553, doi:10.1038/nature07648.
- (5) Maan H. Hani, Martin Sparre, Sara L. Ellison, Paul Torrey, and Mark Vogelsberger. Galaxy mergers moulding the circum-galactic medium - I. The impact of a major merger. MNRAS, 475:1160–1176, Mar 2018. arXiv:1801.06183, doi:10.1093/mnras/stx3252.
- (6) J. Silk and M. J. Rees. Quasars and galaxy formation. A&A, 331:L1–L4, March 1998.
- (7) T. Di Matteo, V. Springel, and L. Hernquist. Energy input from quasars regulates the growth and activity of black holes and their host galaxies. Nature, 433:604–607, February 2005. arXiv:arXiv:astro-ph/0502199, doi:10.1038/nature03335.
- (8) P. F. Hopkins, D. Kereš, J. Oñorbe, C.-A. Faucher-Giguère, E. Quataert, N. Murray, et al. Galaxies on FIRE (Feedback In Realistic Environments): stellar feedback explains cosmologically inefficient star formation. MNRAS, 445:581–603, November 2014.
- (9) Jason Tumlinson, Molly S. Peeples, and Jessica K. Werk. The Circumgalactic Medium. Annual Review of Astronomy and Astrophysics, 55:389–432, Aug 2017. arXiv:1709.09180, doi:10.1146/annurev-astro-091916-055240.
- (10) J. Suresh, D. Nelson, S. Genel, K. H. R. Rubin, and L. Hernquist. Zooming in on accretion - II. Cold circumgalactic gas simulated with a super-Lagrangian refinement scheme. MNRAS, 483:4040–4059, March 2019. arXiv:1811.01949, doi:10.1093/mnras/sty3402.
- (11) C. B. Hummels, B. D. Smith, P. F. Hopkins, B. W. O’Shea, D. W. Silvia, J. K. Werk, et al. The Impact of Enhanced Halo Resolution on the Simulated Circumgalactic Medium. arXiv e-prints, November 2018. arXiv:1811.12410.
- (12) F. van de Voort, V. Springel, N. Mandelker, F. C. van den Bosch, and R. Pakmor. Cosmological simulations of the circumgalactic medium with 1 kpc resolution: enhanced H I column densities. MNRAS, 482:L85–L89, January 2019. arXiv:1808.04369, doi:10.1093/mnrasl/sly190.
- (13) C. Cicone, R. Maiolino, E. Sturm, J. Graciá-Carpio, C. Feruglio, R. Neri, et al. Massive molecular outflows and evidence for AGN feedback from CO observations. A&A, 562:A21, February 2014.
- (14) Claudia Cicone, Paola Severgnini, Padelis P. Papadopoulos, Roberto Maiolino, Chiara Feruglio, Ezequiel Treister, et al. ALMA [C I] P - P Observations of NGC 6240: A Puzzling Molecular Outflow, and the Role of Outflows in the Global Factor of (U)LIRGs. ApJ, 863:143, Aug 2018. arXiv:1807.06015, doi:10.3847/1538-4357/aad32a.
- (15) S. Aalto, S. Garcia-Burillo, S. Muller, J. M. Winters, P. van der Werf, C. Henkel, et al. Detection of HCN, HCO, and HNC in the Mrk 231 molecular outflow. Dense molecular gas in the AGN wind. A&A, 537:A44, January 2012.
- (16) S. Aalto, S. Garcia-Burillo, S. Muller, J. M. Winters, E. Gonzalez-Alfonso, P. van der Werf, et al. High resolution observations of HCN and HCOJ = 3-2 in the disk and outflow of Mrk 231. Detection of vibrationally excited HCN in the warped nucleus. A&A, 574:A85, February 2015.
- (17) Adam K. Leroy, Fabian Walter, Paul Martini, Hélène Roussel, Karin Sandstrom, Jürgen Ott, et al. The Multi-phase Cold Fountain in M82 Revealed by a Wide, Sensitive Map of the Molecular Interstellar Medium. ApJ, 814:83, Dec 2015. arXiv:1509.02932, doi:10.1088/0004-637X/814/2/83.
- (18) A. Fluetsch, R. Maiolino, S. Carniani, A. Marconi, C. Cicone, M. A. Bourne, et al. Cold molecular outflows in the local Universe and their feedback effect on galaxies. MNRAS, 483:4586–4614, Mar 2019. arXiv:1805.05352, doi:10.1093/mnras/sty3449.
- (19) Katherine Alatalo. Escape, Accretion, or Star Formation? The Competing Depleters of Gas in the Quasar Markarian 231. ApJ, 801:L17, Mar 2015. arXiv:1502.00624, doi:10.1088/2041-8205/801/1/L17.
- (20) Sijing Shen, Piero Madau, Javiera Guedes, Lucio Mayer, J. Xavier Prochaska, and James Wadsley. The Circumgalactic Medium of Massive Galaxies at z ~3: A Test for Stellar Feedback, Galactic Outflows, and Cold Streams. ApJ, 765:89, Mar 2013. arXiv:1205.0270, doi:10.1088/0004-637X/765/2/89.
- (21) Joaquin Prieto, Andrés Escala, Marta Volonteri, and Yohan Dubois. How AGN and SN Feedback Affect Mass Transport and Black Hole Growth in High-redshift Galaxies. ApJ, 836:216, Feb 2017. arXiv:1701.06172, doi:10.3847/1538-4357/aa5be5.
- (22) T. Costa, D. Sijacki, and M. G. Haehnelt. Fast cold gas in hot AGN outflows. MNRAS, 448:L30–L34, March 2015.
- (23) S. D. Johnson, H.-W. Chen, and J. S. Mulchaey. On the origin of excess cool gas in quasar host haloes. MNRAS, 452:2553–2565, September 2015. arXiv:1505.07838, doi:10.1093/mnras/stv1481.
- (24) F. Arrigoni Battaia, J. F. Hennawi, J. X. Prochaska, J. Oñorbe, E. P. Farina, S. Cantalupo, et al. QSO MUSEUM I: a sample of 61 extended Ly -emission nebulae surrounding z ~ 3 quasars. MNRAS, 482:3162–3205, January 2019. arXiv:1808.10857, doi:10.1093/mnras/sty2827.
- (25) L. Wisotzki, R. Bacon, J. Brinchmann, S. Cantalupo, P. Richter, J. Schaye, et al. Nearly all the sky is covered by Lyman- emission around high-redshift galaxies. Nature, 562:229–232, October 2018. arXiv:1810.00843, doi:10.1038/s41586-018-0564-6.
- (26) E. Borisova, S. Cantalupo, S. J. Lilly, R. A. Marino, S. G. Gallego, R. Bacon, et al. Ubiquitous Giant Ly Nebulae around the Brightest Quasars at z ~ 3.5 Revealed with MUSE. ApJ, 831:39, November 2016. arXiv:1605.01422, doi:10.3847/0004-637X/831/1/39.
- (27) F. Arrigoni Battaia, J. F. Hennawi, J. X. Prochaska, and S. Cantalupo. Deep He II and C IV Spectroscopy of a Giant Ly Nebula: Dense Compact Gas Clumps in the Circumgalactic Medium of a z ~ 2 Quasar. ApJ, 809:163, August 2015. arXiv:1504.03688, doi:10.1088/0004-637X/809/2/163.
- (28) S. Cantalupo, S. J. Lilly, and M. G. Haehnelt. Detection of dark galaxies and circum-galactic filaments fluorescently illuminated by a quasar at z = 2.4. MNRAS, 425:1992–2014, September 2012. arXiv:1204.5753, doi:10.1111/j.1365-2966.2012.21529.x.
- (29) J. L. Pineda, W. D. Langer, T. Velusamy, and P. F. Goldsmith. A Herschel [C ii] Galactic plane survey. I. The global distribution of ISM gas components. AA, 554:A103, June 2013. arXiv:1304.7770, doi:10.1051/0004-6361/201321188.
- (30) C. Cicone, R. Maiolino, S. Gallerani, R. Neri, A. Ferrara, E. Sturm, et al. Very extended cold gas, star formation and outflows in the halo of a bright quasar at z6. A&A, 574:A14, February 2015.
- (31) B. H. C. Emonts, M. D. Lehnert, M. Villar-Martín, R. P. Norris, R. D. Ekers, G. A. van Moorsel, et al. Molecular gas in the halo fuels the growth of a massive cluster galaxy at high redshift. Science, 354:1128–1130, December 2016. arXiv:1612.00387, doi:10.1126/science.aag0512.
- (32) B. H. C. Emonts, M. D. Lehnert, H. Dannerbauer, C. De Breuck, M. Villar-Martín, G. K. Miley, et al. Giant galaxy growing from recycled gas: ALMA maps the circumgalactic molecular medium of the Spiderweb in [C I]. MNRAS, 477:L60–L65, June 2018. arXiv:1802.08742, doi:10.1093/mnrasl/sly034.
- (33) B. Emonts, C. Carilli, D. Narayanan, M. Lehnert, and K. Nyland. The Molecular High-z Universe on Large Scales: Low-Surface-Brightness CO and the Strength of the ngVLA Core. In E. Murphy, editor, Science with a Next Generation Very Large Array, volume 517 of Astronomical Society of the Pacific Conference Series, page 587, December 2018. arXiv:1810.06770.
- (34) M. Ginolfi, R. Maiolino, T. Nagao, S. Carniani, F. Belfiore, G. Cresci, et al. Molecular gas on large circumgalactic scales at z = 3.47. MNRAS, 468:3468–3483, July 2017. arXiv:1611.07026, doi:10.1093/mnras/stx712.
- (35) J. A. Fernández-Ontiveros, L. Spinoglio, M. Pereira-Santaella, M. A. Malkan, P. Andreani, and K. M. Dasyra. Far-infrared Line Spectra of Active Galaxies from the Herschel/PACS Spectrometer: The Complete Database. ApJS, 226:19, October 2016. arXiv:1607.02511, doi:10.3847/0067-0049/226/2/19.
- (36) M. Béthermin, C. De Breuck, B. Gullberg, M. Aravena, M. S. Bothwell, S. C. Chapman, et al. An ALMA view of the interstellar medium of the z = 4.77 lensed starburst SPT-S J213242-5802.9. AA, 586:L7, February 2016. arXiv:1601.01682, doi:10.1051/0004-6361/201527739.
- (37) N. Rangwala, P. R. Maloney, J. Glenn, C. D. Wilson, J. Kamenetzky, M. R. P. Schirm, et al. First Extragalactic Detection of Submillimeter CH Rotational Lines from the Herschel Space Observatory. ApJ, 788:147, June 2014. arXiv:1404.7200, doi:10.1088/0004-637X/788/2/147.
- (38) E. González-Alfonso, J. Fischer, S. Bruderer, M. L. N. Ashby, H. A. Smith, S. Veilleux, et al. Outflowing OH in Markarian 231: The Ionization Rate of the Molecular Gas. ApJ, 857:66, April 2018. arXiv:1803.04690, doi:10.3847/1538-4357/aab6b8.
- (39) B. Godard, E. Falgarone, M. Gerin, D. C. Lis, M. De Luca, J. H. Black, et al. Comparative study of CH and SH absorption lines observed towards distant star-forming regions. AA, 540:A87, April 2012. arXiv:1201.5457, doi:10.1051/0004-6361/201117664.
- (40) E. Falgarone, M. A. Zwaan, B. Godard, E. Bergin, R. J. Ivison, P. M. Andreani, et al. Large turbulent reservoirs of cold molecular gas around high-redshift starburst galaxies. Nature, 548:430–433, August 2017. arXiv:1708.08851, doi:10.1038/nature23298.
- (41) S. Muller, M. Guélin, M. Dumke, R. Lucas, and F. Combes. Probing isotopic ratios at z = 0.89: molecular line absorption in front of the quasar PKS 1830-211. AA, 458:417–426, November 2006. arXiv:astro-ph/0608105, doi:10.1051/0004-6361:20065187.
- (42) F. Combes. Molecular absorptions in high-z objects. Ap&SS, 313:321–326, January 2008. arXiv:astro-ph/0701894, doi:10.1007/s10509-007-9632-3.
- (43) T. Mroczkowski, D. Nagai, K. Basu, J. Chluba, J. Sayers, R. Adam, et al. Astrophysics with the Spatially and Spectrally Resolved Sunyaev-Zeldovich Effects. A Millimetre/Submillimetre Probe of the Warm and Hot Universe. Space Sci. Rev., 215:17, February 2019. arXiv:1811.02310, doi:10.1007/s11214-019-0581-2.
- (44) T. Mroczkowski, D. Nagai, P. Andreani, M. Arnaud, J. Bartlett, N. Battaglia, et al. A High-resolution SZ View of the Warm-Hot Universe. arXiv e-prints, March 2019. arXiv:1903.02595.
- (45) N. Battaglia, S. Ferraro, E. Schaan, and D. N. Spergel. Future constraints on halo thermodynamics from combined Sunyaev-Zel’dovich measurements. JCAP, 11:040, November 2017. arXiv:1705.05881, doi:10.1088/1475-7516/2017/11/040.
- (46) J. P. Greco, J. C. Hill, D. N. Spergel, and N. Battaglia. The Stacked Thermal Sunyaev-Zel’dovich Signal of Locally Brightest Galaxies in Planck Full Mission Data: Evidence for Galaxy Feedback? ApJ, 808:151, August 2015. arXiv:1409.6747, doi:10.1088/0004-637X/808/2/151.
- (47) J. J. Ruan, M. McQuinn, and S. F. Anderson. Detection of Quasar Feedback from the Thermal Sunyaev-Zel-dovich Effect in Planck. ApJ, 802:135, April 2015. arXiv:1502.01723, doi:10.1088/0004-637X/802/2/135.
- (48) R. Warmels, A. Biggs, P. Cortes, B. Dent, J. Di Francesco, E. Fomalont, et al. Alma technical handbook, alma doc. 6.3, ver. 1.0, March 2018. URL: www.almascience.org.
- (49) F. Bertoldi. The Atacama Large Aperture Submm/mm Telescope (AtLAST) Project. In Atacama Large-Aperture Submm/mm Telescope (AtLAST), page 3, January 2018. doi:10.5281/zenodo.1158842.
- (50) Carlos De Breuck. Site considerations for atlast, January 2018. URL: https://doi.org/10.5281/zenodo.1158848, doi:10.5281/zenodo.1158848.
- (51) Peter Hargrave. Atlast telescope design working group report, January 2018. URL: https://doi.org/10.5281/zenodo.1159025, doi:10.5281/zenodo.1159025.
- (52) P. Klaassen and J. Geach. Galactic Science Case for AtLAST. In Atacama Large-Aperture Submm/mm Telescope (AtLAST), page 20, January 2018. doi:10.5281/zenodo.1159041.
- (53) T. Mroczkowski and O. Noroozian. AtLAST Instrumentation Considerations and Overview. In Atacama Large-Aperture Submm/mm Telescope (AtLAST), page 26, January 2018. doi:10.5281/zenodo.1159053.
- (54) R. Davé, D. Anglés-Alcázar, D. Narayanan, Q. Li, M. H. Rafieferantsoa, and S. Appleby. Simba: Cosmological Simulations with Black Hole Growth and Feedback. arXiv e-prints, January 2019. arXiv:1901.10203.
- (55) M. R. Krumholz. DESPOTIC - a new software library to Derive the Energetics and SPectra of Optically Thick Interstellar Clouds. MNRAS, 437:1662–1680, January 2014. arXiv:1304.2404, doi:10.1093/mnras/stt2000.
- (56) SKA Organization. Advancing astrophysics with the square kilometre array, September 2015. URL: https://www.skatelescope.org/news/ska-science-book/.
- (57) E. J. Murphy, A. Bolatto, S. Chatterjee, C. M. Casey, L. Chomiuk, D. Dale, et al. The ngVLA Science Case and Associated Science Requirements. In Eric Murphy, editor, Science with a Next Generation Very Large Array, volume 517 of Astronomical Society of the Pacific Conference Series, page 3, Dec 2018. arXiv:1810.07524.
- (58) G. J. Stacey, M. Aravena, K. Basu, N. Battaglia, B. Beringue, F. Bertoldi, et al. CCAT-Prime: science with an ultra-widefield submillimeter observatory on Cerro Chajnantor. In Ground-based and Airborne Telescopes VII, volume 10700 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, page 107001M, July 2018. arXiv:1807.04354, doi:10.1117/12.2314031.
- (59) A. Crites, J. Bock, M. Bradford, B. Bumble, T.-C. Chang, Y.-T. Cheng, et al. Measuring the Epoch of Reionization using [CII] Intensity Mapping with TIME-Pilot. In American Astronomical Society Meeting Abstracts #229, volume 229 of American Astronomical Society Meeting Abstracts, page 125.01, January 2017.
- (60) G. Lagache, M. Cousin, and M. Chatzikos. The [CII] 158 m line emission in high-redshift galaxies. AA, 609:A130, February 2018. arXiv:1711.00798, doi:10.1051/0004-6361/201732019.
- (61) Peter Ade, James Aguirre, Zeeshan Ahmed, Simone Aiola, Aamir Ali, David Alonso, et al. The Simons Observatory: science goals and forecasts. Journal of Cosmology and Astro-Particle Physics, 2019:056, Feb 2019. arXiv:1808.07445, doi:10.1088/1475-7516/2019/02/056.
- (62) Kevork N. Abazajian, Peter Adshead, Zeeshan Ahmed, Steven W. Allen, David Alonso, Kam S. Arnold, et al. CMB-S4 Science Book, First Edition. arXiv e-prints, page arXiv:1610.02743, Oct 2016. arXiv:1610.02743.