No evidence for enhanced [O III] 88m emission in a z6 quasar compared to its companion starbursting galaxy
We present ALMA band 8 observations of the [O III] 88m line and the underlying thermal infrared continuum emission in the z=6.08 quasar CFHQS J2100–1715 and its dust–obscured starburst companion galaxy (projected distance: 60 kpc). Each galaxy hosts dust–obscured star formation at rates 100 M yr, but only the quasar shows evidence for an accreting 10 M black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, 1 Gyr after the Big Bang. Bright [O III] 88m emission from ionized gas is detected in both systems; the positions and line–widths are consistent with earlier [C II] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [O III] 88m/FIR luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [O III] 5008 Å emission) does not dominantly determine the quasar’s FIR [O III] 88m luminosity. Both the inferred number of photons needed for the creation of O and the typical line ratios can be accounted for without invoking extreme (top–heavy) stellar initial mass functions in the starbursts of both sources.
0000-0003-4793-7880]Fabian Walter \move@AU\move@AF\@affiliationMax Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany \move@AU\move@AF\@affiliationNational Radio Astronomy Observatory, Pete V. Domenici Array Science Center, P.O. Box O, Socorro, NM 87801, USA
0000-0001-9585-1462]Dominik Riechers \move@AU\move@AF\@affiliationCornell University, 220 Space Sciences Building, Ithaca, NY 14853, USA \move@AU\move@AF\@affiliationMax Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
0000-0002-2662-8803]Roberto Decarli \move@AU\move@AF\@affiliationINAFâOsservatorio di Astrofisica e Scienza dello Spazio, via Gobetti 93/3, I-40129, Bologna, Italy
Winona State University, Winona, MN 55987, USA
0000-0001-9024-8322]Bram Venemans \move@AU\move@AF\@affiliationMax Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
0000-0002-2931-7824]Eduardo Bañados \move@AU\move@AF\@affiliationThe Observatories of the Carnegie Institution for Science, 813 Santa Barbara Street, Pasadena, CA 91101, USA
Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
National Radio Astronomy Observatory, Pete V. Domenici Array Science Center, P.O. Box O, Socorro, NM 87801, USA
Steward Observatory, The University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA
University of California Santa Barbara, Santa Barbara, CA, USA
European Southern Observatory, Alonso de Córdova 3107, Vitacura, Región Metropolitana, Chile
0000-0002-9838-8191]Marcel Neeleman \move@AU\move@AF\@affiliationMax Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 USA
National Radio Astronomy Observatory, Pete V. Domenici Array Science Center, P.O. Box O, Socorro, NM 87801, USA \move@AU\move@AF\@affiliationMax Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
Kavli Institute of Astronomy and Astrophysics at Peking University, No.5 Yiheyuan Road, Haidian District, Beijing, 100871, China
The most distant quasars at z6 are unique probes of galaxy and structure formation in the first Gigayear of the universe. Their rapidly accreting supermassive black holes make these quasars the most luminous (non–transient) sources known at this cosmic epoch (e.g., Fan et al., 2006; Venemans et al., 2015; Bañados et al., 2016; Wu et al., 2015). The interstellar medium of more than 30 host galaxies of z6 quasars has now been detected through the [C II] 158 m line (hereafter: [C II]) and far–infrared (FIR) continuum emission. This suggests that the host galaxies have intense star formation rates at many 100’s of Myr occuring coevally with the growth of the accreting central black holes (e.g., Bertoldi et al., 2003; Walter et al., 2003, 2009; Wang et al., 2013; Willott et al., 2013, 2015; Decarli et al., 2018; Venemans et al., 2018).
Using ALMA observations, Decarli et al. (2017) found dust–obscured companions to 4 quasars in a [C II] survey of 27 z 6 quasars. These massive star-forming companion galaxies show no signs of AGN activity and were detected serendipitously within 8â-60 kpc at the same redshifts in four systems, suggesting that they are physically associated with overdensities or proto–clusters of galaxies at z6. These companion galaxies rival the FIR and [C II] luminosities of the neighboring quasars, in some cases even surpassing their (already extreme) luminosity. These quasar–galaxy pairs thus provide a unique opportunity to efficiently study the interstellar medium in distinctly different galactic environments in the first Gyr of the Universe: pure starbursts, and those in the presence of a powerful AGN.
To date, the [C II] line has been the workhorse line of ISM studies at z6 (e.g., Carilli & Walter, 2013), but it is only one of several fine structure lines that help to charaterize the physical properties of the interstellar medium (e.g., Díaz-Santos et al., 2017). Of these, the [O III] 88 m line (PP) (hereafter: [O III]), emitted by doubly ionized oxygen (O), is a particularly important diagnostic of the star formation process (e.g., Ferkinhoff et al., 2010, 2011, 2015; Vishwas et al., 2018): Given the O ionization potential of 35.1 eV, O is almost exclusively found in dense HII regions around O–type stars or in AGN environments. Indeed, an enhancement of the optical 5008Å [O III] line is a key diagnostic of AGN activity (e.g., Baldwin et al., 1981). In contrast, neutral carbon has an ionization potential of only 11.3 eV, and [C II] is therefore be emitted both in the neutral and ionized phase. So far, [O III] detections have been reported in a number of z 6 galaxies (e.g., Inoue et al., 2016; Laporte et al., 2017; Carniani et al., 2017; Tamura et al., 2018; Hashimoto et al., 2018a; Hashimoto et al., 2018b; Marrone et al., 2018). However, in this redshift range, only 2 systems have been detected in [O III], [C II], and the underlying dust continuum: the Lyman break galaxy B14–65666 at z=7.15 (Hashimoto et al., 2018b) and the lensed galaxy SPT0311–58 at z=6.90 (Marrone et al., 2018, the latter actually consists of two sources).
In this paper we report [O III] and underlying dust continuum observations with ALMA of the z=6.08 quasar CFHQS J2100–1715 (hereafter J2100–Q) and its dust–enshrouded companion starburst (J2100–SB). J2100–Q was discovered by Willott et al. (2010a), and a black hole mass of 9.42.610 M was reported in Willott et al. (2010b). J2100–SB was discovered by Decarli et al. (2017) in the far–infrared using ALMA. Multi–wavelength follow–up observations of the companion galaxy using Spitzer (3.6m, 4.5m) and HST (WFC3/140W) as well as MUSE and LBT spectroscopy did not detect the companion in the rest–frame UV/optical, with an implied ratio of obscured to un–obscured star formation of 99% (Mazzucchelli et al. in prep.). No metallicity measurement of our targets exists, but they must be significantly enriched with heavy elements as they harbor significant amounts of dust, and have high (100 M yr) star formation rates.
In Section 2 we describe our ALMA band 8 observations, in Sec. 3 we present our results, followed by a discussion and a summary in Secs. 4 and 5. Throughout this paper we use cosmological parameters H = 70 km s Mpc, = 0.3, and = 0.7, in agreement with Planck Collaboration XVI (2014), leading to a scale of 5.67 kpc per arcsec at the redshift of our sources.
We observed the [O III] line ( = 3393.0062 GHz) of J2100–Q and J2100–SB redshifted to =479.2GHz at z=6.08 with ALMA in band 8 and configuration C43–1 on 2018 May 20 (two executions) and 2018–May–23 (one execution) for a total of 35200 s or 4.3 h, of which 2.5 h were spent on source. The redshifts of the sources were known from the earlier detections of [[C II]] line emission (z=6.08060.0011 for the quasar, 6.07960.0008 for the companion). The two sources are very close in redshift (=41 km s, Decarli et al. 2017), meaning that the [O III] emission in the two objects could be covered with the same frequency setup. However, the projected separation on the sky (11) required two pointings to cover the sources in band 8. The software package CASA (McMullin et al., 2007) was used for data reduction and imaging. For optimal sensitivity, we employed natural weighting when imaging the data, leading to a synthesized beam size of 0.730.56 (3.7 kpc at z=6.08), and an rms of 0.45 mJy beam per 50 km/s (80 MHz) channel. The corresponding sensitivity in the continuum (excluding the channels that contain line emission) is 54 Jy beam over an effective bandwidth of 5.625 GHz. For all analysis of the [O III] line data, the continuum has been subtracted in the uv–plane using line–free channels.
We also re–analyze the [C II] and underlying continuum ALMA data published by Decarli et al. (2017) and Venemans et al. (2018). The beam size of these [C II] data (0.730.57) is almost identical to that presented here. Therefore, no additional beam matching was required. In order to optimize the S/N in the measurements, and given the fact that any extent, if present, is marginal, we extract [C II] line and underlying continuum fluxes at the peak position of both sources. We report the [C II] values adopted in this study in Tab. 1.333We note that both Decarli et al. (2017) and Venemans et al. (2018) report [C II] and FIR luminosities that are slightly higher than the values adopted here, as their reported fluxes are based on 2D fitting of the line emission in the image domain.
3.1 Continuum emission
In Figure 1 we show the rest–frame 88 m (observed: 475 GHz) continuum map of J2100–Q and J2100–SB. This map is based on those channels in the data cube that do not contain the [O III] line. The sources are detected at high significance and are unresolved at the resolution of our beam (0.730.57, flux densities are reported in Tab. 1). The coordinates (RA/DEC in J2000.0) of the sources, 21:00:54.70, –17:15:21.9 (J2100–Q), and 21:00:55.45, –17:15:21.7 (J2100–SB) are in agreement with Decarli et al. (2017).
Together with our earlier restframe 158 m continuum measurement around the [C II] line (Decarli et al. 2017, see Tab. 1 for slightly updated numbers) we can use the new restframe 88 m observations to further constrain the properties of the FIR continuum emission. The IR and FIR luminosities are obtained by fitting the 88 m and 158 m data points (assuming they are optically thin) with a modified blackbody, taking the effect of the cosmic microwave background into account (da Cunha et al., 2013). These results are reported in Tab. 1. The resulting L–based star formation rates (Kennicutt & Evans, 2012) are 12514 M yr (J2100–Q) and 28430 M yr (J2100–SB).
3.2 [O Iii] line emission
We show the continuum–subtracted [O III] spectra of both sources in Fig. 2 as extracted from the peak positions of the continuum map. The Gaussian fit results are reported in Tab. 1. For reference, we overplot scaled, continuum–subtracted [C II] spectra from Decarli et al. (2017) in greyscale. Both the redshifts and the linewidths of the [O III] and [C II] detections agree within the uncertainties. The [O III] line width of the quasar is much smaller than that of emission lines seen in the broad line region of the quasar (FWHM: 3600 km s, Willott et al. 2010b). In Fig. 3 we show the integrated [O III] maps of both sources, after continuum subtraction. 2D–Gaussian fitting shows that the [O III] line emission is unresolved at our resolution like the continuum emission. In Tab. 1 we summarize the resulting [O III] luminosities. Using the [O III]–SFR scaling relations of De Looze et al. (2014) for their ‘high–redshift sample’ we derive [O III]–based star formation rates of approximately 100 M yr and 360 M yr for the quasar and the companion, respectively, consistent with the FIR–based estimates (Sec. 3.1).
3.3 Line/FIR luminosity ratio
We now combine the [O III] and rest–frame 88m continuum measurements with the earlier [C II] line measurements to constrain the properties of the interstellar medium in our targets. This is particularly interesting, as our sources are physically quite different. In Fig. 4 (left panel) we plot the [O III]/FIR luminosity ratios vs. the continuum flux density ratios (S/S) of our sources, the latter being a proxy for the temperature of the dust (color–coded by FIR luminosity), and compare them to the recent compilation of local luminous infrared galaxies (LIRGs) by Díaz-Santos et al. (2017). The mid–infrared emission in the majority of the local sample is not dominated by AGN emission (circles vs. stars in Fig. 4).
Our two sources have luminosities at the high end of the LIRG sample studied in Díaz-Santos et al. (2017), but their [O III]/FIR luminosity ratios lie only slightly above the average LIRG value of similar dust temperature, within the scatter found locally. Perhaps unexpectedly, despite the presence of an accreting 10 M, the quasar does not show enhanced [O III]emission compared to either the companion or the local comparison sample.
It is instructive to also compare the [C II]/FIR luminosity ratios for our sources to the local LIRGs (middle panel of Fig. 4). The [C II]/FIR luminosity ratios of our targets (1.610) are located within the distribution observed for the local LIRGS at the same temperature (around 110). The right panel in Fig. 4 shows the [O III]/[C II] luminosity ratios as a function of dust temperature. Again, our two targets fall within the distribution seen in local LIRGS, and the [O III]/[C II] ratio is actually lower in the quasar compared to the companion.
4.1 Origin of the [C II] and [O III] emission
The line–widths and positions of the [O III] and [C II] emission are identical within the uncertainties, which indicates that they trace the same gravitational potential on galactic scales (cf. Carniani et al. 2017). Díaz-Santos et al. (2017) show that the ratio of the [C II] emission emerging from PDRs, to that from ionized gas, is correlated with S/S (see their Fig. 3). We determine the S/S ratio from our S/S measurements assuming the same blackbody parameters as used in Sec. 3.1 (see also Tab. 1), finding S/S = 1.7 for the quasar, and 1.8 for the companion, respectively.
If the correlations of Díaz-Santos et al. (2017) hold, then more than 80% of the [C II] emission in both the quasar and the companion are from the neutral phase of the interstellar medium (L 0.2 L). This is in contrast to the [O III] emission, which only traces the ionized phase, because the O ionization potential is 35.1 eV (the next ionization level is at 54.9 eV, e.g., Kramida et al. 2018). Following Ferkinhoff et al. (2010, 2011) and Vacca et al. (1996), a blackbody with a temperature 35,000 K (45,000 K) produces significant emission of photons with energies eV (55 eV). These temperatures correspond to stars of type O7V and O3V, respectively.
The presence of a FUV–bright central source (i.e., the accreting black hole in J2100–Q) will will lead to optical 5008 Å [O III] line emission in the broad line region (e.g., Vanden Berk et al. 2001) and the narrow line region (e.g., Baldwin et al. 1981) that will become accessible to the James Webb Space Telescope at z = 6.08. The FIR [O III] line on the other hand is likely to arise from stellar HII regions given that it has the same systemic velocity and line width as the [C II] line, and is much narrower than the optical lines emerging from the BLR. Additionally, the [O III]/FIR and [O III]/[C II] ratios are similar to those seen for systems where the emission is known to arise from HII regions.
4.2 Minimum ionized gas mass
We can infer the mass of doubly–ionized oxygen in our targets from the observed [O III] luminosity (see, e.g., Ferkinhoff et al., 2011) to be (O)/M= /L = M for the quasar, and M for the companion. Assuming that most the oxygen is in doubly–ionized form, and that the oxygen abundance is (O/H)= (Savage & Sembach, 1996), this implies a minimum mass of ionized gas of (H)= M for the J2100–Q, and (H)= M for J2100–SB. Other high-z systems have shown O/O fractions of 10% (Ferkinhoff et al., 2011), so that the actual ionized gas mass could be a factor of 10 higher than our minimum value here. Using the dust mass estimates of 3.210 M and 5.510 M derived for J2100–Q and J2100–SB, respectively (Decarli et al., 2017), and assuming a standard gas–to–dust ratio of 100 (Draine et al., 2007; Sandstrom et al., 2013), the molecular gas content (M(H)) in both sources is 10 M. A ratio of /M(H) 10 is similar to what is found in nearby galaxies (Brauher et al. 2008, see also the compilation in Ferkinhoff et al. 2011, their figure 3). The total ionized gas mass is thus only a small fraction of the total mass of the dense interstellar medium (and thus baryonic mass) in the two objects.
4.3 Minimum number of ionizing photons and stars
Following Ferkinhoff et al. (2010, 2011, 2015) and Vishwas et al. (2018), we calculate minimum O–ionizing photon rates for J2100–Q and the J2100–SB of s and s, respectively. It is interesting to compare this number with back-of-the-envelope expectations based on the observed SFR. The companion galaxy has a SFR300 M yr. The stars that provide significant radiation energetic enough to doubly–ionize oxygen live only for 5 Myr. Assuming that the SFR is constant, we infer the total stellar mass created in the last 5 Myr, and distribute it over different stellar mass bins based on a Chabrier (2003) Initial Mass Function (IMF). The stellar main sequence provides us with an estimate of the luminosity and effective temperature of stars (assuming a simple blackbody) based on their mass, thus we can compute for each stellar mass bin how many photons are produced with energy 35.1 eV, needed to ionize OII. By integrating over the stellar mass distribution, scaled to match the total mass of stars produced in 5 Myr, we obtain s (this number changes to s for a Kroupa 2001 IMF). This number is slightly higher than the one inferred from our [O III] observations (we get similar results for J2100–Q). This suggests that star formation with a standard IMF is sufficient to account for the observed [O III] luminosity, with no need to invoke a top-heavy IMF or a contribution by a ‘buried’ quasar.
We present ALMA [O III] observations of the z=6.08 quasar J2100–Q and its dust–enshrouded starburst companion J2100–SB. This system is unique, as it offers the possibility to study the physical properties of the interstellar medium of two starbursts that are physically associated with each other within 1 Gyr of the Big Bang. They display distinct galactic environments: one source shows the presence of a powerful AGN, while the other does not. Interestingly, we find that the [O III]/FIR luminosity ratio of the starburst companion J2100–SB is higher than that of the quasar J2100–Q, even though the accreting supermassive black hole in the latter provides additional photons with energies high enough to produce O (which typically leads to strong optical 5008 Å[O III] emission from the broad and narrow line regions, measureable with JWST in the future). The fact that we do not see enhanced FIR [O III] emission in J2100–Q indicates that its [O III] emission is dominated by star formation, and not AGN activity. This is supported by the finding that the linewidths and positions of the [O III] and [C II] line are the same in both sources and consistent with previous findings that the interstellar medium properties in distant quasar host galaxies are predominantly powered by star formation (e.g., Leipski et al., 2014; Barnett et al., 2015; Venemans et al., 2017). One caveat is that J2100–SB could in principle also host an obscured accreting supermassive black hole. In that case, the [O III] emission in the companion may not be solely due to star formation. However, the fact that the [O III]/FIR and [O III]/[C II] luminosity ratios of J2100–SB lie within the range of local LIRGs (including those that are not dominated by an AGN) argues against such a scenario. This latter finding, together with our analysis of the number of photons needed for the creation of O, implies that no extreme (top–heavy) initial stellar mass functions are needed to explain the [O III] luminosity in our sources.
We thank the referee for helpful comments that improved our manuscript and Gordon Stacey for useful suggestions. FW, BV, MN, and MaN acknowledge support from the ERC Advanced Grant 740246 (Cosmic_Gas). DR acknowledges support from the National Science Foundation under grant number AST–1614213. RW acknowledge supports from the National Science Foundation of China (11721303). This paper makes use of the following ALMA data: ADS/JAO.ALMA#2017.1.00118.S,
ADS/JAO.ALMA#2015.1.01115.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.
Facilities: ALMA Software: CASA (McMullin et al., 2007)
- Bañados et al. (2016) Bañados, E., Venemans, B. P., Decarli, R., et al. 2016, ApJS, 227, 11
- Baldwin et al. (1981) Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5
- Barnett et al. (2015) Barnett, R., Warren, S. J., Banerji, M., et al. 2015, A&A, 575, A31
- Bertoldi et al. (2003) Bertoldi, F., Cox, P., Neri, R., et al. 2003, A&A, 409, L47
- Brauher et al. (2008) Brauher, J. R., Dale, D. A., & Helou, G. 2008, The Astrophysical Journal Supplement Series, 178, 280
- Carilli & Walter (2013) Carilli, C. L., & Walter, F. 2013, Annual Review of Astronomy and Astrophysics, 51, 105
- Carniani et al. (2017) Carniani, S., Maiolino, R., Pallottini, A., et al. 2017, A&A, 605, A42
- Chabrier (2003) Chabrier, G. 2003, Publications of the Astronomical Society of the Pacific, 115, 763
- da Cunha et al. (2013) da Cunha, E., Groves, B., Walter, F., et al. 2013, ApJ, 766, 13
- De Looze et al. (2014) De Looze, I., Cormier, D., Lebouteiller, V., et al. 2014, A&A, 568, A62
- Decarli et al. (2017) Decarli, R., Walter, F., Venemans, B. P., et al. 2017, Nature, 545, 457
- Decarli et al. (2018) —. 2018, ApJ, 854, 97
- Díaz-Santos et al. (2017) Díaz-Santos, T., Armus, L., Charmandaris, V., et al. 2017, ApJ, 846, 32
- Draine et al. (2007) Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, 663, 866
- Fan et al. (2006) Fan, X., Strauss, M. A., Becker, R. H., et al. 2006, AJ, 132, 117
- Ferkinhoff et al. (2015) Ferkinhoff, C., Brisbin, D., Nikola, T., et al. 2015, ApJ, 806, 260
- Ferkinhoff et al. (2010) Ferkinhoff, C., Hailey-Dunsheath, S., Nikola, T., et al. 2010, ApJ, 714, L147
- Ferkinhoff et al. (2011) Ferkinhoff, C., Brisbin, D., Nikola, T., et al. 2011, ApJ, 740, L29
- Hashimoto et al. (2018a) Hashimoto, T., Laporte, N., Mawatari, K., et al. 2018a, Nature, 557, 392
- Hashimoto et al. (2018b) Hashimoto, T., Inoue, A. K., Mawatari, K., et al. 2018b, ArXiv e-prints, arXiv:1806.00486
- Inoue et al. (2016) Inoue, A. K., Tamura, Y., Matsuo, H., et al. 2016, Science, 352, 1559
- Kennicutt & Evans (2012) Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531
- Kramida et al. (2018) Kramida, A., Ralchenko, Y., Nave, G., & Reader, J. 2018, in APS Division of Atomic, Molecular and Optical Physics Meeting Abstracts, Vol. 2018, M01.004
- Kroupa (2001) Kroupa, P. 2001, MNRAS, 322, 231
- Laporte et al. (2017) Laporte, N., Nakajima, K., Ellis, R. S., et al. 2017, ApJ, 851, 40
- Leipski et al. (2014) Leipski, C., Meisenheimer, K., Walter, F., et al. 2014, ApJ, 785, 154
- Marrone et al. (2018) Marrone, D. P., Spilker, J. S., Hayward, C. C., et al. 2018, Nature, 553, 51
- McMullin et al. (2007) McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell, Vol. 376, 127
- Sandstrom et al. (2013) Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5
- Savage & Sembach (1996) Savage, B. D., & Sembach, K. R. 1996, ApJ, 470, 893
- Tamura et al. (2018) Tamura, Y., Mawatari, K., Hashimoto, T., et al. 2018, ArXiv e-prints, arXiv:1806.04132
- Vacca et al. (1996) Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914
- Vanden Berk et al. (2001) Vanden Berk, D. E., Richards, G. T., Bauer, A., et al. 2001, AJ, 122, 549
- Venemans et al. (2015) Venemans, B. P., Verdoes Kleijn, G. A., Mwebaze, J., et al. 2015, MNRAS, 453, 2259
- Venemans et al. (2017) Venemans, B. P., Walter, F., Decarli, R., et al. 2017, ApJ, 845, 154
- Venemans et al. (2018) Venemans, B. P., Decarli, R., Walter, F., et al. 2018, ApJ, 866, 159
- Vishwas et al. (2018) Vishwas, A., Ferkinhoff, C., Nikola, T., et al. 2018, ApJ, 856, 174
- Walter et al. (2009) Walter, F., Riechers, D., Cox, P., et al. 2009, Nature, 457, 699
- Walter et al. (2003) Walter, F., Bertoldi, F., Carilli, C., et al. 2003, Nature, 424, 406
- Wang et al. (2013) Wang, R., Wagg, J., Carilli, C. L., et al. 2013, in Molecular Gas, Dust, and Star Formation in Galaxies, ed. T. Wong & J. Ott, Vol. 292, 184–187
- Willott et al. (2015) Willott, C. J., Bergeron, J., & Omont, A. 2015, ApJ, 801, 123
- Willott et al. (2013) Willott, C. J., Omont, A., & Bergeron, J. 2013, ApJ, 770, 13
- Willott et al. (2010a) Willott, C. J., Delorme, P., Reylé, C., et al. 2010a, AJ, 139, 906
- Willott et al. (2010b) Willott, C. J., Albert, L., Arzoumanian, D., et al. 2010b, AJ, 140, 546
- Wu et al. (2015) Wu, X.-B., Wang, F., Fan, X., et al. 2015, Nature, 518, 512