Gas Rich Starbursts in Luminous Reddened Quasars

The Discovery of Gas Rich, Dusty Starbursts in Luminous Reddened Quasars at with ALMA

Manda Banerji, C. L. Carilli, G. Jones, J. Wagg, R. G. McMahon, P. C. Hewett, S. Alaghband-Zadeh, C. Feruglio
Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
Kavli Institute of Cosmology Cambridge, Madingley Road, Cambridge CB3 0HA, UK
National Radio Astronomy Observatory, 1003 Lopezville Road, Socorro, NM 87801, USA
Battcock Centre for Experimental Astrophysics, Cavendish Laboratory, Cambridge CB3 0HE, UK
Physics Department, New Mexico Institute of Mining and Technology, Socorro, NM 87801, USA
SKA Organization, Lower Withington Macclesfield, Cheshire SK11 9DL, UK
Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy
E-mail: mbanerji@ast.cam.ac.uk
Abstract

We present ALMA observations of cold dust and molecular gas in four high-luminosity, heavily reddened (A mag) Type 1 quasars at with virial MM, to test whether dusty, massive quasars represent the evolutionary link between submillimetre bright galaxies (SMGs) and unobscured quasars. All four quasars are detected in both the dust continuum and in the CO(3-2) line. The mean dust mass is 610M assuming a typical high redshift quasar spectral energy distribution (T=41K, =1.95 or T=47K, =1.6). The implied star formation rates are very high - 1000 M yr in all cases. Gas masses estimated from the CO line luminosities cover 1-5()M and the gas depletion timescales are very short - Myr. A range of gas-to-dust ratios is observed in the sample. We resolve the molecular gas in one quasar - ULASJ23150143 () - which shows a strong velocity gradient over 20 kpc. The velocity field is consistent with a rotationally supported gas disk but other scenarios, e.g. mergers, cannot be ruled out at the current resolution of these data. In another quasar - ULASJ1234+0907 () - we detected molecular line emission from two millimetre bright galaxies within 200 kpc of the quasar, suggesting that this quasar resides in a significant over-density. The high detection rate of both cold dust and molecular gas in these sources, suggests that reddened quasars could correspond to an early phase in massive galaxy formation associated with large gas reservoirs and significant star formation.

keywords:
galaxies:evolution – galaxies:formation – galaxies:high-redshift – galaxies:starburst

1 Introduction

Since the discovery of the correlation between the stellar bulge mass in galaxies and the mass of their supermassive black holes (Magorrian:98; Kormendy:13), it is now widely acknowledged that these central supermassive black holes play an important role in governing the formation and evolution of their host galaxies. In our current picture of galaxy formation, galaxies and their supermassive black-holes co-evolve. Molecular gas and the dusty interstellar medium (ISM) in galaxies, acts as the fuel for both star formation and black hole accretion. Therefore, studying the distribution of dust and gas in star forming galaxies with actively accreting supermassive black holes can help shed light on the exact physical processes that are driving co evolution.

Both star formation and black hole accretion activity in the Universe peak at (e.g. Madau:98; Richards:06) and this represents the key epoch in the Universe’s history at which to study the feeding and feedback processes in galaxies and quasars that are the main drivers of co-evolution. Over the last two decades, surveys starting with the SCUBA bolometer on the James Clerk Maxwell Telescope (JCMT) and, more recently, using the Herschel space satellite, have led to the discovery of large populations of far infrared (FIR) and submillimetre bright starburst galaxies (SMGs) at the main epoch of galaxy formation (e.g. Smail:97; Blain:99; Casey:14). Molecular gas detections are now being assembled for many tens to hundreds of such galaxies (e.g. Greve:05; Tacconi:06; Bothwell:13; Genzel:15). Several studies of molecular gas in galaxies have also taken advantage of gravitational lensing to amplify the molecular line signal (e.g. Coppin:07; Danielson:11; Aravena:16). While clearly an advantage for probing intrinsically fainter galaxies, differential magnification between the different dust and gas components can make interpretations regarding their actual spatial distribution complicated in such lensed systems.

Quasars, on account of being among the brightest extragalactic sources in our Universe, were some of the first sources to be observed in molecular gas (e.g. Solomon:05; Barvainis:97). Using the compilations of extragalactic sources with CO detections in Carilli:13 and Heywood:13, we find there are a total of 89 quasars that have been detected in CO all the way from the local Universe out to the highest redshifts of . Many of the recent efforts have focussed on quasars (e.g. Carilli:02; Walter:04; Wang:10) and, as a consequence, the numbers of quasars with CO detections at the peak epoch of galaxy formation (1.5) is considerably smaller - 20 quasars. Many of the quasars that have been observed in molecular gas at these epochs are gravitationally lensed systems (e.g. Riechers:11b; Sharon:16). Excluding these lensed systems and focussing only on unlensed quasars where spatially mapping the gas and dust distributions is considerably easier, there are only 13 quasars at with CO detections (Coppin:08; Simpson:12; Schumacher:12; Willott:07). Detecting molecular gas in even a small number of unlensed quasars at these redshifts, to synchronously study the fuelling of star formation and black hole accretion at the epoch when both are at their peak, is therefore valuable.

Most searches for molecular gas in quasars have started from samples of ultraviolet (UV) luminous quasars (e.g. Engels:98; Hewett:95; Schneider:10) that were already known to be FIR/millimetre bright (Omont:03; Priddey:03; Stevens:05). The studies by Priddey:03 and Omont:03 find that only 9 out of 53 UV luminous quasars at (i.e. 17 per cent) are detected at 850m down to 3 flux density limits of 7-9 mJy and 9 out of 26 UV luminous quasars at (i.e. 35 per cent) are detected at 1.2 mm down to 3 flux density limits of 1.8-4 mJy. In the galaxy formation scenario first advocated by Sanders:88 and now commonly adopted in galaxy formation models (e.g. Hopkins:08; Narayanan:10), the most highly star forming galaxies at high redshift e.g. the SMGs, will eventually evolve into UV luminous quasars. The relatively small fraction of FIR/mm-bright quasars argues for a quick transition between the two populations and recent studies suggest that the transition phase could be as short lived as 1 Myr (Simpson:12). In such a picture of galaxy formation, transition populations of hybrid SMG-quasars should exist and a larger fraction of these transition quasars should be FIR bright, gas rich systems compared to the UV luminous quasars. Molecular gas observations of SMGs and UV luminous quasars have also revealed some differences between the two populations. While SMGs appear to have substantial amounts of low excitation gas (Riechers:11a; Bothwell:13), this does not seem to be the case in quasars (Riechers:11b) which could therefore represent a later evolutionary stage when the extended gas reservoirs in the galaxy have been significantly depleted. Once again, if transition objects between these two populations do exist, their gas fractions might be expected to be intermediate between the SMGs and UV luminous quasars.

Searches for these so called “transition” galaxies have often focussed on obscured AGN populations with different selection methods picking out AGN and quasars with a wide variety of dust extinctions and luminosities (e.g. Urrutia:08; Brusa:10; Banerji:12; Banerji:13; Banerji:15; Glikman:12; Eisenhardt:12; Tsai:15). While in principle, any or all of these populations could represent the missing evolutionary link between star forming galaxies and optical quasars, in practice further multiwavelength observations are necessary to establish that these obscured quasars are indeed distinct in terms of their physical properties from matched control samples of unobscured quasars. Our searches for obscured quasars have been focussed at the highest luminosities and at the peak epoch of galaxy formation at with colour cuts deliberately chosen to isolate high luminosity quasars with the same levels of dust extinction as seen in SMGs at similar redshifts. To date, we have spectroscopically confirmed a new sample of almost 60 heavily reddened, luminous broad line quasars with a median A 2.5 mags (c.f. A for SMGs; Takata:06). We have already found evidence that a higher fraction of our reddened quasars are actively star forming relative to optical quasars (Banerji:14; Alaghband-Zadeh:16) consistent with these being transition objects seen as they are both rapidly forming stars and rapidly growing their black holes. Our reddened quasars are also among the most luminous quasars known at these epochs with recent results suggesting that their space density actually exceeds that of unobscured, UV luminous quasars at the highest luminosities (Banerji:15). The bolometric luminosities are 10 erg s and black hole masses are 10-10 M inferred from the broad H emission lines in the NIR spectra, consistent with the most luminous quasars found in optical surveys such as the Sloan Digital Sky Survey (SDSS).

In order to determine whether the reddened broad line quasars are indeed a transition population between SMGs and optical quasars, direct comparison of their molecular gas and ISM properties to both SMGs and optical quasars is now necessary. Here we discuss the first ALMA observations of four heavily reddened quasars from Banerji:12 and Banerji:15 (B12 and B15 hereafter). Our new data is supplemented with two more detections of the CO(3-2) line in two quasars from our B12 parent sample (Feruglio:14; Brusa:15) to give a sample of six heavily reddened quasars with molecular gas measurements. In Section 2 we present the ALMA observations of the four new quasars. Section 3 discusses both the dust continuum and CO(3-2) line properties of the quasars and presents physical properties of the quasar host galaxies including star formation rates, dust and molecular gas masses. In Section 4 we attempt to put our results in context with previous observations of molecular gas in both high redshift SMGs and optical quasars to determine whether the dust and gas properties of the reddened quasars are in fact intermediate between the two. Throughout this paper we assume a flat CDM cosmology with h=0.7, , .

2 Data

2.1 Sample Selection & ALMA Observations:

In this work we make use of ALMA data obtained as part of the Cycle 3 project 2015.1.01247.S (PI:Banerji). The aim is to detect, for the first time, the molecular gas reservoirs in four heavily reddened quasars from B12 and B15 via the CO(3-2) emission line. We refer readers to B12, B15 and Banerji:13 for details of the parent sample of reddened (A) broad emission line quasars from which our targets are drawn. Three of the four quasars selected for ALMA observations - ULASJ0123+1525, ULASJ1234+0907 and ULASJ2315+0143 - correspond to three of the reddest quasars in our sample with inferred dust extinctions towards the quasar continuum of A=4.0, 6.0 and 3.4 mags respectively. Assuming the extinction is connected to star formation in the quasar host galaxy as discussed extensively in B12 and B15, these broad line quasars would be the best candidates for highly star forming gas rich hosts. The fourth quasar in this paper - VHSJ2101-5943 - represents a quasar with average extinction (A=2.5 mags) in our sample. As such, it can be taken to be representative of the gas and dust properties of the larger sample of heavily reddened broad line quasars. We searched for radio counterparts to these quasars in the VLA-FIRST, NVSS and SUMSS radio catalogues and did not find any matches within 10. The H derived redshifts, AGN luminosities (at rest-frame 6m) and black hole masses for these quasars can be seen in Table 1.

We emphasise that apart from ULASJ1234+0907, which was detected in the Herschel PACS and SPIRE bands (Banerji:14), none of the other three targets were known a-priori to be FIR or millimetre bright. This is in contrast to the targeting strategies for detecting molecular gas in UV luminous, unobscured quasars where only UV luminous quasars already known to be FIR bright have typically been observed in molecular gas. As discussed in B12 and Banerji:14, it is unlikely that our reddened quasars are lensed. Optical and near infrared band images of our quasars demonstrate that the reddened quasars are extremely faint or completely invisible at these wavelengths and, given the depths of these imaging data, a lensing galaxy should have been visible out to .

ULASJ01231525 ULASJ12340907 VHSJ21015943 ULASJ23150143
R.A. (J2000) 01:23:12.52 12:34:27.52 21:01:19.46 23:15:56.23
Dec (J2000) 15:25:22.7 09:07:54.2 59:43:44.8 01:43:50.4
K 18.59 18.05 16.68 18.38
z 2.629 2.503 2.313 2.560
A / mag 4.0 6.0 2.5 3.4
log(L / erg s) 47.3 46.8 46.9 47.0
log(M / M) 9.7 10.4 10.5 10.1

A=R=3.1

Table 1: Summary of the properties - band magnitudes, redshifts, line-of-sight extinctions, AGN luminosities and black hole masses - for the four reddened quasars observed with ALMA. The properties have been derived from NIR photometry and spectroscopy and presented in B12 and B15.

Observations were conducted between 2015 December and 2016 January and, as a detection experiment, made use of the most compact configuration available in Cycle 3, which corresponds to the lowest angular resolution (2-3 arcsec). The number of useable antennae increased over the duration of the observations resulting in increased sensitivity for observations conducted at a later date. The correlator was configured to four dual polarization bands of 2 GHz (1.875 GHz effective) bandwidth each, providing a channel width of 15.6 MHz. The basebands were set up such that the CO(3-2) emission line is positioned in the first band with the other three bands providing a measurement of the dust continuum at observed wavelengths of 3 mm. The central frequencies of the four bands for each quasar, together with other details of the observations, can be found in Table 2.

ULASJ01231525 ULASJ12340907 VHSJ21015943 ULASJ23150143
Date Completed 2016-01-10 2015-12-26 2016-01-21 2016-01-07
Exposure Time / sec 3629 1058 4052 3145
Number of Antennae 46 34 48 47
Beam Size (line) / arcsec 3.02.5 3.22.6 3.62.5 2.92.2
Observed Frequencies / GHz 95.286, 94.286 98.714, 97.714 104.375, 103.267 97.134, 96.134
107.286, 106.286 86.714, 85.714 92.375, 91.375 109.134, 108.134
Channel RMS / mJy / beam 0.13 0.21 0.13 0.11
Continuum RMS / Jy / beam 16 20 11 15
Continuum Flux Density / Jy / beam 9816 6721 4112 25917
z 2.62970.0012 2.50260.0012 2.31130.0034 2.56140.0028
CO(3-2) Line Intensity / Jy km s 1.400.08 0.970.18 0.460.03 0.910.05

Median RMS over 15.6 MHz channels

Table 2: Details of the ALMA observations for the four reddened quasars presented in this work.

2.2 Data Reduction

The ALMA data were calibrated using the ALMA pipeline in the Common Astronomy Software Applications package, casa (v4.5.1-4.5.2) by executing the appropriate ALMA calibration scripts corresponding to the release data of the observations. Time dependent amplitude and phase variations were calibrated using nearby quasars and radio galaxies. Flux calibrations made use of observations of Uranus and Neptune. The typical calibration uncertainties are of the order of 10-20 per cent. The beam sizes as well as the root mean square (RMS) sensitivity of both the CO(3-2) (median over 15.6 MHz channels) and the dust continuum observations can be seen in Table 2. Dust continuum images were produced for each quasar from the calibrated visibilities, by combining the line free channels from the three continuum spectral windows in multi-frequency synthesis mode using the CASA task clean and a natural weighting scheme to maximise the sensitivity. These same channels were also used to produce a UV plane model of the continuum emission, which was then subtracted from the first spectral window containing the line using the CASA task uvcontsub. The continuum subtracted line visibilities were then imaged using clean, once again employing a natural weighting scheme in order to produce the final line cubes. Finally, both the continuum images and line cubes were corrected for the primary beam response.

The spectral profile in the data cubes is a convolution of the intrinsic spectrum with a spectral resolution function. By default, the ALMA pipeline calibration performs Hanning smoothing of the data cubes resulting in a spectral resolution of 2 the channel spacing. We have corrected for this effective reduction in resolution in all the spectral fits presented in this paper, although the effect is small. We have detected the dust continuum as well as the CO(3-2) line in all four quasars (Figs. 2 and 3) and the dust continuum flux densities, CO line intensities and CO derived redshifts are presented in Table 2.

3 Results & Analysis

3.1 Dust Continuum: Dust Masses, FIR Luminosities & Star Formation Rates

The dust continuum detections at observed frame wavelengths of 2.9-3.5 mm (rest frame 800-900m), suggest that there is a significant amount of cold dust in all four quasar host galaxies. The dust continuum emission is unresolved in all cases. The ALMA continuum detections probe the Rayleigh Jeans tail of the dust spectral energy distribution (SED) for these quasars, where the emission from the dust is expected to be optically thin. Under these assumptions dust masses can be calculated from the single photometric ALMA data point, after assuming a typical dust temperature, T and dust emissivity index, as follows:

(1)

where the mass absorption coefficient of the dust:

(2)

is the luminosity distance and is the Planck function. We assume =0.045 m kg at =250 GHz (Greve:12). The uncertainty in the dust mass therefore depends primarily on the uncertainty in the dust temperature. In order to estimate the dust temperature, T, and emissivity index, , that is most appropriate for our quasars, we begin by using the fact that for one of our quasars, ULASJ1234+0907, we already have a complete sampling of the FIR spectral energy distribution (SED) from Herschel observations (Banerji:14). The Herschel data traces rest-frame wavelengths of 70140m in this quasar where the optically thin assumption tends to break down. We therefore fit a modified single temperature greybody to the Herschel observations combined with the ALMA 3mm continuum detection assuming the emission is optically thick. The best-fit SED has parameters: T=314K and =2.50.3.

Previous investigations of the dust SEDs of high redshift quasars have found T=47K and =1.6 (Beelen:06; Wang:08) or T=41K and =1.95 (Priddey:01). In Fig. 1, we show all three assumptions for single temperature greybodies overlaid on the photometry for ULASJ12340907. As we shall see in Section 3.4.1, the Herschel fluxes for ULASJ12340907 could be boosted by the presence of other millimetre bright galaxies within the Herschel beam and the dust temperature and emissivity index derived from fitting the Herschel + ALMA data may therefore not be accurate. Throughout this paper we will attempt to relate our measurements to the measured dust properties of other high redshift quasars. To be consistent with these previous works, from hereon we therefore choose to adopt the average of two commonly used assumptions for the dust SED of high redshift quasars: (i) T=47K, =1.6 and (ii) T=41K, =1.95. The dust masses in Table 3 therefore correspond to dust temperatures of 4147 K.

Under the same assumptions for the dust SED, we can also calculate the FIR luminosities and star formation rates by simply scaling the optically thick single temperature greybodies to match the dust continuum detections at 3mm. The FIR luminosity here is defined between 40 and 300m and we use the Kennicutt:12 relation to convert this FIR luminosity to a star formation rate:

(3)

Once again, these values can be found in Table 3 and correspond to the average of the T=47K, =1.6 and T=41K, =1.95 greybodies. These assumptions should be kept in mind when interpreting the physical properties of the quasar host galaxies.

ULASJ01231525 ULASJ12340907 VHSJ21015943 ULASJ23150143 ULASJ15390557 ULASJ10020137
log(L/L) 13.220.07 13.220.12 12.910.15 13.610.03 13.36 12.47
SFR / M yr 2500500 2500900 1200400 6100400 3400 270
log(M/M) 8.67 8.68 8.40 9.08 8.89
M / 10 M 4.30.9 4.41.8 2.30.9 110.8 7.0
L / 10 K km s pc 6.50.3 4.20.8 1.70.1 4.10.2 6.41.1 2.40.5
M / 10 M 5.20.3 3.40.7 1.40.1 3.30.2 5.10.9 1.90.4
Gas-to-dust Ratio 110 70 55 28 24
/ M (K km s pc) 0.7 1.0 1.3 2.6 2.9
FWHM / km s 52040 700120 19020 35050 840 550200
M sin(i) / 10 M 0.4-3.1 0.8-5.7 0.07-0.5 0.3-1.8 0.6 4.5

This Work; Feruglio:14; Brusa:15

(a) Assuming the average of T=47K, =1.6 and T=41K, =1.95 for the four quasars in this work as well as ULASJ1539+0557; For ULASJ1002+0137 the FIR luminosity and SFR are taken directly from Brusa:15.

(b) Assuming gas-to-dust ratio of 91 (Sandstrom:13);

(c) Assuming M (K km s pc);

(d) Assuming M=M;

(e) Assuming M=M;

(f) Corrected for instrumental resolution due to Hanning smoothing of spectra;

(g) Full range in values quoted using the different estimators detailed in Section 3.3 for the four quasars in this work. For ULASJ1539+0557 and ULASJ1002+0137, the dynamical masses are taken directly from Feruglio:14 and Brusa:15.

Table 3: FIR luminosities, star formation rates, gas and dust masses for the four reddened quasars observed with ALMA and presented in this work, as well as the two reddened quasars from the same parent sample studied in Feruglio:14 and Brusa:15. In the latter cases, the CO(1-0) luminosities and gas masses have been re-derived assuming the same CO excitation and CO-to-H conversion factor as has been assumed for the quasars in this work.
Figure 1: Herschel SPIRE photometry at observed wavelengths of 250, 350 and 500m together with the ALMA 3mm continuum photometry for ULASJ1234+0907 (). The best fit power-law + single temperature greybody fit to these data points is also shown (see Banerji:14). We also plot three different single temperature, optically thick greybodies with T=31K, =2.5 (best fit), T=41K, (Priddey:01) and T=47K, (Beelen:06; Wang:08). In Section 3.4.1 we show that the Herschel fluxes are likely contaminated by other galaxies within the Herschel beam. The latter two spectral energy distributions have therefore been used in the paper to calculate FIR luminosities, star formation rates and dust masses for all the reddened quasars in our sample.

3.2 CO(3-2) Line Luminosities and Molecular Gas Masses

All four quasars are also detected in the CO(3-2) line. Most of these quasars were not previously known to be millimetre bright and were identified for follow-up observations based purely on their near infrared colours and implied dust extinctions. This strongly suggests that the near-infrared selection is successfully identifying broad line quasars in gas rich, highly star forming host galaxies at .

Line maps extracted over the full-width-half-maximum (FWHM) velocity interval centred on the line can be seen in Fig. 2 together with the corresponding CO spectra in Fig. 3. We checked by fitting a 2-dimensional (2-D) Gaussian to the line maps that the centroid of the CO line emission was always within 0.5-1 arcsec of the centroid of the quasar emission as measured from the near infrared images in B12 and B15. The largest offset is seen in ULASJ2315+0143 where the peak CO emission lies 1 arcsec from the quasar emission. ULASJ23150143 is spatially resolved in our data and is discussed further below in Section 3.3 and 3.4.2. All CO line properties quoted for ULASJ2315+0143 correspond to those derived from the spatially integrated spectrum of this source.

Figure 2: CO(3-2) (left) and dust continuum (right) maps for all four quasars in our sample. Contours start at 1.5 and are shown in intervals of 1.5. Solid contours denote regions of positive flux while dashed contours denote regions of negative flux. The star marks the position of the quasar from the near infrared observations. The beam size is indicated in the bottom left corner of the map.
ULASJ01231525 (z=2.630) ULASJ12340907 (z=2.503)
VHSJ21015943 (z=2.311) ULASJ23150143 (z=2.561)
Figure 3: CO(3-2) spectra for all four quasars in our sample. The error bars denote the 1 uncertainties on the measured flux. The spectra are Hanning smoothed straight out of the ALMA pipeline. Dashed lines show the best fit Gaussian profile fit to these data. In the case of ULASJ01231525 two Gaussians have been used to fit the data. The velocities are relative to the H redshifts from B12 and B15 and the vertical dashed line marks zero velocity.

Our observations of the dust continuum as well as the CO(3-2) line reveal a range in the CO(3-2)-to-FIR luminosities in the sample. While ULASJ0123+1525 has L/L=2.510, in ULASJ2315+0143 this ratio is 2.210 - i.e. almost an order of magnitude larger. Seen another way, the CO equivalent widths in the sample also cover a large range from 3.510 km s in ULASJ2315+0143 to 1.410 km s in ULASJ0123+1525. For ULIRG type systems, the CO(3-2) line typically comprises 25 per cent of the 850m flux within a 310 km s bandwidth, which corresponds to an equivalent width of the CO(3-2) line of 110 km s (Seaquist:04). The range in CO-to-FIR luminosities and equivalent widths seen in our reddened quasars may indicate a range in ISM properties in the sample and further observations of other molecular lines and CO transitions would help constrain the excitation conditions of the gas in these systems.

The dust and gas detections also allow us to calculate molecular gas masses in two different ways. Recently there have been several studies in the literature advocating the use of dust masses as a proxy for the gas mass assuming a metallicity dependent gas-to-dust ratio (e.g. Scoville:14). By assuming a gas-to-dust ratio of 91 from observations of nearby galaxies (Sandstrom:13), we can convert the dust masses in Table 3 to gas masses and, using the CO line luminosity, we can therefore estimate the CO-to-H conversion factor, . These values can be seen in Table 3 and we find that a range of values are derived for the sample all the way from 0.7 M (K km s pc) appropriate for nuclear starbursts and quasars to 3 (K km s pc) in ULASJ2315+0143, appropriate for disk like galaxies like our own Milky Way as well as high redshift star forming disks. It is interesting to note that ULASJ2315+0143 is also the source with the most extended gas emission. However, we caution that these values are highly dependent on the assumptions made regarding the dust SED of the quasars in Section 3.1.

Gas masses can also be estimated directly from the CO line luminosities, which are also given in Table 3. We assume an excitation ratio between the CO(3-2) and CO(1-0) line of =0.8, which is intermediate between the typical values for SMGs and optical quasars from Carilli:13 and consistent with our interpretation that reddened quasars are transitioning from an SMG-like to an optical quasar-like phase. The corresponding CO(1-0) line luminosities can be seen in Table 3 and these are converted to gas masses assuming (K km s pc) appropriate for nuclear starbursts and quasars. The resultant gas-to-dust ratios presented in Table 3 range from values of 100, similar to the values in nearby galaxies, to 30-50, consistent with other high redshift star forming galaxies (Aravena:16), but are once again dependant on the exact choice of dust SED parameters.

3.3 CO(3-2) Kinematics and Dynamical Masses

Our observations were carried out in low resolution mode corresponding to the most compact ALMA Cycle 3 configuration. However, the gas emission nevertheless appears to be extended in some of our quasars. We tested for evidence for spatially resolved emission using the spectroastrometry method detailed in e.g. Carniani:13. Specifically, we constructed separate line maps from the red and blue side of the line and looked for evidence for a spatial offset between the two maps. These line maps can be seen in Fig. 4 for our most extended source, ULASJ2315+0143. We measure a (1.200.13) arcsec offset between the two maps in Fig. 4 corresponding to physical scales of (9.61.0) kpc. The dust continuum emission from ULASJ2315+0143 is more compact and spatially unresolved in our data, suggesting that the gas reservoirs extend to larger scales than the cold dust in this galaxy. Using the same method, we find a (0.780.16) arcsec (6.41.3 kpc) offset between the red and blue emission in VHSJ21015943 and a (0.580.2) arcsec (4.71.6 kpc) offset in ULASJ1234+0907. In ULASJ0123+1525 the red and blue wing emission map centroids are separated by (0.040.09) arcsec (0.30.7 kpc). In all cases (apart from ULASJ2315+0143 which is discussed in Section 3.4.2), the line emission is therefore spatially unresolved or only marginally resolved.

Figure 4: CO line maps for ULASJ2315+0143 extracted using the blue and red sides of the line separately. Given the high signal-to-noise (S/N) of these data, the typical centroiding errors for the red and blue wings are 0.07 and are therefore too small to be visible on this plot. A clear spatial offset is seen between the two maps with the dots marking the centres of the two maps as calculated using a 2-D (in RA, Dec) Gaussian fitting. This suggests that the source is spatially resolved.

We use either a single or double Gaussian to fit the line profiles in Fig. 3 to characterise the line shape and also compute non-parametric first and second intensity weighted moments to derive line centroids and line widths. The CO redshifts quoted in Table 2 are derived from these line centroids. The CO-derived redshifts are all within 250 km s of the H derived redshifts in B12 and B15. Given the broad H line profiles in these quasars, the typical uncertainties in the H redshift estimates can be several hundred km s.

In all cases, apart from ULASJ0123+1525, a single Gaussian gives an acceptable fit to the data with a reduced value of 1. In the case of ULASJ0123+1525, a double Gaussian leads to a considerably better fit to the data. The parametric and non-parametric values for the line centroid and line widths are consistent within the error bars and in Table 3 we quote the non-parametric widths of the CO emission line.

Dynamical masses can be estimated from the widths of the CO lines in Table 3 assuming a source size. We have already seen that in ULASJ2315+0143 the gas emission is spatially resolved with an implied source radius of 4.8 kpc from the spectroastrometry method. We can use the relation in Gnerucci:11 to derive the corresponding dynamical mass:

(4)

This assumes that the gas emission comes from a rotating disk with the red and blue images shown in Fig. 4 tracing different sides of the disk due to rotation. In such an estimator, does not correspond to a conventional size but is rather the average distance between the redshifted and blueshifted gas weighted by the surface brightness of the CO line. We leave the dependence on inclination explicit in the derivations of dynamical masses. This gives M=1.5M for ULASJ2315+0143. Assuming a typical R=4 kpc for the other three quasars, the dynamical masses are (1.3, 2.3 and 0.2)M for ULASJ0123+1525, ULASJ1234+0907 and VHSJ2101-5943 respectively. For completeness, we also estimate the dynamical masses using two other estimators: 1) the virial estimator: M=2.8210 and the rotating disk estimator (Neri:03): Msin(i)=410 . The values quoted in Table 3 encompass the full range in dynamical masses from these three estimators.

3.4 Notes on Individual Objects

3.4.1 ULASJ1234+0907: A Significant Overdensity of Millimetre Bright Galaxies

Given that our high luminosity, massive quasars are likely to be the progenitors of the most massive galaxies seen in the Universe today, they may be expected to reside in massive halos. Significant overdensities of galaxies around the quasar may therefore be expected even at . The ALMA field-of-view at these frequencies (FWHM of primary beam 1 arcmin) allows us to search for other millimetre-bright galaxies in the vicinity of these quasars. We have detected two more galaxies with molecular line emission - one located 21 (170 kpc at ) to the north of the quasar and the other 11.5 (90 kpc at ) to the south of the quasar. We name these galaxies G1234-N and G1234-S respectively. The northern source G1234-N is also detected in the dust continuum. The CO line profiles and line maps can be seen in Fig. 5 and we list the positions of the two galaxies (as measured from the centroid of the CO emission) together with continuum flux densities and CO line intensities in Table 4. As the galaxies are not located at the phase center, we determine the scaling that needs to be applied to the flux densities in order to account for the primary beam correction. This scaling is 1.02 for G1234-N and 1.01 for G1234-S and the flux densities in Table 4 have been adjusted accordingly. In light of the new detections of these two galaxies, we also note that the Herschel SPIRE fluxes for ULASJ1234+0907 in Banerji:14 are almost certainly overestimated as all three galaxies fall within the Herschel beam at these wavelengths. Reducing the Herschel SPIRE fluxes for ULASJ1234+0907 to 40 per cent of their measured values (to account for the flux density ratio measured between the quasar and the other millimetre bright galaxies with ALMA), means the FIR luminosity of this source decreases by 0.5 dex. Even accounting for this reduction in the FIR luminosity, ULASJ1234+0907 still has a measured SFR of 1000 M yr. Note, the dust mass for ULASJ12340907 in Table 3 is calculated from the ALMA continuum detection only and is therefore not affected by the Herschel blending.

Figure 5: Top: CO line maps constructed by summing the channels containing the CO line in G1234-N (left) and G1234-S (right). The offsets in RA and Dec are relative to the position of the quasar (the red star), ULASJ12340907 which marks the phase center. G1234-N lies north-east of the quasar while G1234-S lies south of the quasar. The G1234-S channel map overlaps in frequency with the CO line in ULASJ1234+0907 and both sources are therefore visible in the map. Bottom: CO(3-2) spectrum for G1234-N (left) and G1234-S (right) with the dashed line showing the best fit Gaussian fit to the line profile. The error bars denote the 1 uncertainties on the flux measurements.

. Object RA Dec z S / Jy I / Jy km s G1234-N 12:34:26.36 09:08:06.1 2.514 975 0.790.14 G1234-S 12:34:27.31 09:07:43.1 2.497 40 (2) 0.390.12

Table 4: Properties of the two millimetre-bright galaxies detected in the field-of-view of the quasar ULASJ1234+0907 at

The CO line profile in G1234-N is well fit by a single Gaussian with FWHM=530100 km s, while the line in G1234-S is fit by a Gaussian with FWHM=540130 km s. Based on these CO and dust continuum detections, we can estimate the dust masses and gas masses for the two companion galaxies. For the gas masses, we assume an excitation ratio of , representative of SMGs, and a CO-to-H conversion factor of =0.8 M (K km s pc). The gas masses for G1234-N and G1234-S are (3.30.6)M and (1.90.5)M respectively. Based on the dust continuum detection in G1234-N and assuming a dust temperature of 30K and a dust emissivity index, (e.g. Casey:14), the dust mass is 1.4M and the SFR is 600 M yr.

The two galaxies G1234-N and G1234-S while clearly not directly interacting with the quasar, are likely to be part of the same overdensity. Several protoclusters have been found at with overdensities of galaxies extending tens of Mpc across e.g. (Venemans:07; Carilli:11; Casey:15; Capak:11).

3.4.2 ULASJ2315+0143: Spatially Resolved Studies of the Molecular Gas

As seen in Section 3.3 and Fig. 4, in ULASJ2315+0143 we find evidence that the gas emission is spatially resolved. The deconvolved source size derived from the collapsed channel map summed across the full velocity extent of the line, is (2.841.36) (0.300.53) arcsec which corresponds to 23 kpc along the major axis at the quasar redshift. At the current resolution of our data, this size measurement is still uncertain. However in Section 4.4 we show that the large inferred size is robust to how the visibilities are weighted to create the final images. We compare this size estimate to the prototype SMG GN20, which has one of the most extended gas reservoirs seen in high redshift galaxies - 144 kpc in diameter as mapped using the CO(2-1) emission (Hodge:12). Thus, the warm, excited gas in ULASJ23150143 appears to be even more extended compared to what is seen in lower excitation CO lines in the most extreme SMGs. However, the constraint on the source size in ULASJ2315+0143, is very similar to the size of the beam, and further higher resolution observations are needed to determine whether this extended gas reservoir is associated with a single source.

We construct zeroth- and first-moment maps of the line emission by fitting a single Gaussian profile to the line emission in every pixel, retaining pixels in the central region (i.e. avoiding the edges of the field-of-view) and where the S/N of the line is at least 3. These maps can be seen in Fig. 6 and a velocity gradient is evident across the source.

Figure 6: Zeroth (top) and first moment (middle) maps of the CO emission in ULASJ2315+0143 derived by fitting a single Gaussian to the line emission in each pixel. The systemic velocity has been subtracted from the velocity field shown in the middle panel, which therefore traces the rotational velocity only. The bottom panel shows the rotation curve derived by fitting a tilted ring gas disk dynamical model to the velocity field shown in the middle panel.

We fitted a tilted ring disk model to the velocity field shown in Fig. 6. In such an analysis, the gas emission is modelled as series of concentric rings with inclination, , position angle, P.A., systemic velocity , centre and rotational velocity, . The line of sight velocity at any position can then be described on a ring with radius, as:

(5)

where the angle, is related to P.A. in the plane of the sky by

(6)
(7)

We used the rotcur task within gipsy in iterative mode for the fitting. The expansion velocity was set to zero and the dynamical centre of all the rings was set to be identical. The rings have widths ranging from 0.21-3. The source has a maximum radial extent of 2.9 or 23 kpc. We fixed the P.A. to be that measured from the velocity field by looking at the axis connecting the maximum and minimum velocities. This P.A. is measured to be 352.3. The rotcur task was run 15 times with one to fifteen rings covering the entire radius with only the systemic velocity and then the inclination angle, , as free parameters. The systemic velocity and were then determined from the average of all successful rings where a successful ring is defined to be one with errors on of 20. We have checked that using an error weighted average instead makes very little difference to the model fit parameters. Finally, we ran the task with as the only free parameter and used the results from the most successful rings (defined as those with uncertainties on ) in a single run to produce a model velocity field and a rotation curve for the source. The best fit rotation curve can be seen in the bottom panel of Fig. 6. The best fit tilted ring model has =578.

Given the resolution of our data, the model fails to constrain the internal dynamics of this source at kpc and the rotation curve shown in Fig. 6 is noisy and should be considered as illustrative. However, given that these data were taken at the lowest spatial resolution available in ALMA Band 3, the fact that dynamical modelling is even possible with this data highlights the enormous promise of gaining detailed insight into the dynamics of such high redshift quasar host galaxies with higher resolution ALMA observations. We defer a more detailed discussion of the dynamical properties of this source to future papers.

We can estimate the dynamical mass from the rotation curve shown in Fig. 6. The maximum rotational velocity, =1256 km s at r=21.5 kpc. Using , we get a dynamical mass of (7.820.02)10M which is a factor of 2 lower than that derived using the spectroastrometry method in Section 3.3. Correcting for the full range in inclination angles allowed by the fit, M=(1.01.4)10M.

3.4.3 Ulasj01231525: Distinct Kinematic Components

ULASJ0123+1525 shows two distinct peaks in its CO line profile, which could indicate the presence of kinematically distinct components in the gas emission from this galaxy. A double Gaussian fit yields a significantly better reduced value compared to a single Gaussian (1.2 versus 2.6 for the single Gaussian fit). The two Gaussians have FWHM=24080 km s and 57070 km s respectively and are separated in velocity by 26060 km s. At the spatial resolution of our data, there is no evidence for a spatial offset between these velocity components. Similar double peaked line profiles have been observed in SMGs (e.g. Greve:05; Tacconi:06; Bothwell:13) and could arise from a galaxy merger, a rotating gas disk with a large inclination angle relative to the sky plane or large scale inflows and outflows affecting the molecular gas. Further high resolution observations of the molecular gas would help discriminate between these scenarios.

4 Discussion

4.1 Comparison to Unobscured “Blue” Quasars and X-Ray AGN

In order to understand whether our reddened quasars are indeed a distinct population in terms of their molecular gas and dust continuum properties, we now draw direct comparisons to observations of the molecular gas and dust continuum in samples of blue, unobscured quasars with comparable luminosities and redshifts. To date, the largest surveys for cold dust in high luminosity quasars at , are still those conducted using the SCUBA bolometer (at 850m) on the James Clerk Maxwell Telescope and using the MAMBO bolometer (at 1200m) on the IRAM 30-m telescope (e.g. McMahon:94; McMahon:99; Isaak:02; Priddey:03; Omont:03). The SCUBA surveys have typically probed down to 3 flux limits of 7 mJy at 850m while the MAMBO survey reaches sensitivities of 2-4 mJy (3).

In the SCUBA studies, only 20 per cent of the optically unobscured, blue quasars are “submillimetre loud” with 850m flux densities of 7 mJy. We can estimate the 850m fluxes for the reddened quasars using the assumptions regarding the dust SED set out in Section 3.1. These 850m flux densities are compared to other samples from the literature in Fig. 7 and the errorbars reflect the difference in 850m flux between assuming T=47K and =1.6 (Beelen:06; Wang:08) and T=41K and =1.95 (Priddey:01). Three of the four reddened quasars lie near or below the 2 limit for the SCUBA surveys. From this initial sample of reddened quasars, and with the assumptions regarding their dust SED set out above, 25 per cent of the reddened quasars appear to be “submillimetre loud”, similar to the submillimetre loud fraction of blue quasars. We can also compare the submillimetre and millimetre flux densities for our reddened quasars with the stacked results from the SCUBA surveys, which suggest that the average 850m flux density of a flux limited sample of blue quasars is 1.90.4 mJy at (Priddey:03) and 2.00.6 mJy at (Isaak:02). Based on our ALMA detections, our reddened quasars would all have 850m fluxes that are higher than these average properties for the blue quasars as can be seen in Fig. 7. The larger 850m fluxes could suggest that luminous reddened quasars have higher dust luminosities and are therefore, on average, more likely to be hosted in starburst galaxies relative to their unobscured blue counterparts. However, larger samples of reddened quasars with cold dust continuum measurements and better constraints on their dust SED are clearly required before any robust conclusions can be drawn.

We also compare the submillimetre fluxes of our reddened quasars to those of X-ray selected AGN. X-ray surveys of high redshift AGN have generally covered smaller areas on the sky compared to optical surveys and therefore probe quasars and AGN at lower luminosities. Considering only the most X-ray luminous, high redshift quasars in the COSMOS field with SCUBA-2 850m observations from the SCUBA-2 Cosmology Legacy Survey, the stacked 850m fluxes are 0.70.2 mJy for Type 1 AGN and 1.30.3 mJy for Type 2 AGN (Banerji:15b), once again lower than those measured for the reddened quasars in this work. There are however a number of X-ray AGN that are individually detected at 850m in Banerji:15b with 850m flux densities of 4-10 mJy, comparable with the inferred submillimetre fluxes of our reddened quasars. Previous SCUBA surveys of X-ray absorbed and X-ray unabsorbed quasars have suggested a link between cold dust emission and X-ray obscuration (Page:01; Stevens:05) and we note that for those reddened quasars where we have X-ray spectra (e.g. Banerji:14), moderate column densities of 10 cm are present.

Surveys for molecular gas in blue quasars have focussed on the subset of the population that are (sub)millimetre loud, in contrast to our “blind” search for CO in the reddened quasar population. So far six of our near infrared selected, heavily reddened quasars from the B12 and B15 samples have been observed in CO - four in this paper and one each in Feruglio:14 and Brusa:15. Despite not being selected a-priori to be FIR/mm bright, all six quasars have been detected in CO. As stated earlier, most optical quasars that have been observed in CO were known a-priori to be FIR luminous therefore increasing the probability of a CO detection. Coppin:08 have looked for CO in a sample of 10 quasars of comparable redshift to our reddened quasar sample. Three quasars in their sample were detected in blank field submm surveys over much smaller areas of sky than the wide-field near infrared surveys that have been used to select our reddened quasars (B12,B15), and are therefore lower luminosity. The remaining seven quasars in the Coppin:08 sample, are brighter than the -band flux limit imposed on our reddened quasar searches in B12 and B15 and can therefore be directly compared to our sample. From these seven FIR bright optical quasars, five have been detected in CO with the limits on the CO line luminosity for the remaining two placing them below the CO line luminosities of all six of our reddened quasars. Presumably, starting from a complete sample of unobscured quasars (i.e. not just those known a-priori to be millimetre bright) would result in a smaller fraction being detected in CO. The black hole masses of our reddened quasars and the UV-luminous Coppin:08 quasars are comparable given the typical uncertainties in these black hole mass estimates. Hence, the CO-detection fraction in our sample suggests that a larger fraction of luminous reddened quasars have significant molecular gas reservoirs compared to UV luminous blue quasars of comparable luminosity and mass.

In B15 we demonstrated that our reddened quasars outnumber blue quasars at the highest luminosities whereas the reddened quasar population is subdominant as we approach more typical luminosities around L. One interpretation of these trends put forward in B15 is that the most luminous, massive black holes have a longer duty cycle associated with their assembly or growth phase and therefore are seen as dust obscured for a longer period of time. As all our ALMA targets are drawn from the high luminosity, massive end of the reddened quasar population, the detection of significant amounts of cold dust and molecular gas in these systems fits in with this overall picture. In such an evolutionary sequence, (sub)millimetre-loud blue quasars are presumably a later, more evolved phase when significant gas reservoirs and star formation are still detectable but some of the dust has been cleared out. The (sub)millimetre-faint blue quasars may then correspond to an even more advanced phase when significant star formation has ceased either because the gas supply has been depleted, or as a result of AGN feedback.

4.2 Comparison to Mid-Infrared Luminous AGN

Under the premise that dust reddened quasars are a phase in massive galaxy evolution preceding the unobscured quasar phase, and that this phase is associated with high levels of star formation and high rates of accretion onto the supermassive black hole, we would predict that populations of high luminosity AGN that are even more obscured than our reddened quasars may have similar or higher levels of star formation in their host galaxies. These most highly obscured AGN might therefore be expected to be characterised by high 850m fluxes. The most suitable sample of mid-infrared luminous, heavily obscured AGN with comparable AGN luminosities and redshifts as our reddened quasars are the so called AGN dominated Hot Dust Obscured Galaxies or HotDOGs that have emerged from the WISE All Sky Survey (e.g. Eisenhardt:12). We use as the basis for comparison, the SCUBA-2 850m and ALMA 870m flux densities of the spectroscopically confirmed WISE HotDOGs at from Jones:14; Jones:15 and Lonsdale:15. These studies provide a sample of 42 HotDOGs with secure spectroscopic redshifts overlapping the redshift range of both our reddened quasar sample and the Priddey:03 optical quasars. Once again we extrapolate our reddened quasar 3mm flux densities using the SED assumptions above to predict their 850m fluxes. A comparison of the 850m flux densities can be seen in Fig. 7.

Figure 7: Redshift versus 850m flux density for blue quasars (Priddey:03), reddened quasars (This Work, Feruglio:14; Brusa:15) and WISE HotDOGs (Lonsdale:15; Jones:14; Jones:15). In all cases downward triangles denote upper limits and the large blue open circle shows the stacked 850m flux of the blue quasars from Priddey:03. The dashed lines mark the range in 2 flux limits for the SCUBA survey (Priddey:03) and the ALMA survey (Lonsdale:15).

Out of the 42 WISE HotDOGs, 19 are detected at high significance at 850/870m - a detection fraction of 45 per cent. If we consider the fraction detected to the SCUBA limits of 7 mJy, the number drops to only 3 - i.e. just 7 per cent of the population. As seen in Fig. 7, in many cases the limits on the submillimetre fluxes obtained for the WISE HotDOGs, place them well below the flux densities of our reddened quasars. Thus, as concluded in Lonsdale:15 and several other recent WISE papers, these HotDOGs on average do not show evidence for the dust luminosities being dominated by cold dust heated by star formation and instead, it is likely that most of the dust heating in these galaxies comes from the AGN. The WISE HotDOGs therefore do not seem to fit our simple evolutionary picture linking obscured and unobscured quasars. As suggested in Eisenhardt:12, these WISE AGN could perhaps represent a class of objects where the peak of AGN activity precedes that of star formation. Alternatively, the higher line-of-sight extinction seen in these sources could simply reflect a difference in viewing angle between these WISE AGN and both the reddened and unobscured broad line quasar populations.

There have been relatively few studies of CO molecular gas in the WISE population so far but Wu:14 have looked for CO(3-2) and CO(4-3) in two of these galaxies. Neither source is detected in their observations but the limits on the gas mass from that study are consistent with the low end of the gas masses we have derived for the reddened quasar population in this paper.

It is clear from Fig. 7 that the number of dust continuum observations for optical and WISE-selected AGN in much larger than the number available for reddened quasars. With a large spectroscopic sample of the reddened quasars now in place, it is imperative that we extend the investigations carried out in this work to larger samples in order to gain a coherent picture of the cold dust and molecular gas emission and star formation in these different quasar populations covering a very wide range in line-of-sight extinction.

4.3 Star Formation Efficiencies and Gas Depletion Timescales

Having detected both cold dust and molecular gas in dust reddened quasars, we can attempt to compare the star formation efficiencies and gas depletion timescales to those measured in high redshift starburst galaxies and quasars. We use as the basis for our comparisons, the compilations of CO detections in high redshift galaxies and quasars published by Carilli:13 and Heywood:13. While these compilations are not necessarily exhaustive and complete, they can be taken to be representative of current observational samples. Throughout this analysis we work with the full sample of six heavily reddened quasars with CO detections - i.e. the four in this paper plus the two from Feruglio:14 and Brusa:15. For the last two quasars we have re-derived the CO line luminosities using the same assumptions regarding the CO excitation as in this work (see Table 3). In Fig. 8 we show the CO line luminosity, L, versus the infrared luminosity for all extragalactic sources from Heywood:13 and the six heavily reddened quasars in Table 3. Also shown is the best-fit relation from Carilli:13:

(8)

The reddened quasars in general have higher infrared luminosities or star formation rates for a given gas mass, than predicted by this relation, which has been derived from the entire distribution of high redshift galaxies in the literature. This is consistent with a high star formation efficiency in these reddened quasar hosts compared to other extragalactic sources The infrared luminosities are however broadly consistent with those of other high redshift quasar host galaxies and the scatter in the observed relation shown in Fig. 8 is large. We caution also that the infrared luminosities in the literature have been calculated using heterogenous methods both in terms of the wavelength range used to compute these luminosities and the assumptions made regarding the form of the dust SED.

To gain further insight into the similarities and differences in the star formation efficiencies between SMGs, quasars and reddened quasars, we now compare both the FIR luminosity distributions and CO luminosity distributions over the same redshift range, using a KS-test to determine the probability that they are drawn from the same parent population. There is no evidence for a difference in the FIR luminosity distributions of the three samples when the same assumptions are made regarding the dust temperature and emissivity index of the FIR SED. Comparing the L luminosities for the three samples however, we find that quasars and SMGs have a 3 per cent probability of being drawn from the same distribution, highlighting that the molecular gas reservoirs may indeed be different in the two samples. Reddened quasars on the other hand have a 32 per cent probability of having CO luminosities drawn from the same luminosity distribution as SMGs and a 54 per cent probability of being drawn from the same luminosity distribution as optical quasars so their gas properties are consistent with both SMGs and optical quasars.

We calculate typical gas depletion timescales using the gas masses and star formation rates in Table 3. Regardless of whether the CO based or dust based gas masses are used, the gas depletion timescales are very short - Myr. The black hole mass accretion rates of 100-500 Myr (B12,B15) imply black hole gas consumption timescales that are at least an order of magnitude longer - 100-1000 Myr. The black holes are already very massive - M - suggesting we are witnessing a phase when the galaxies are running out of gas (to further fuel the black hole) following the main period of black-hole mass assembly. The reddened quasars could thus be plausible progenitors of the largest supermassive black holes seen in the Universe today, which have measured black hole masses of 10M (e.g. NGC 3842 and NGC 4899; McConnell:11).

Figure 8: FIR luminosity versus CO line luminosity for extragalactic sources from the literature including quasars (Heywood:13), compared to our sample of reddened quasars presented in this paper and in Brusa:15 and Feruglio:14. The infrared luminosities for the reddened quasars have been calculated as discussed in Section 3.1. The lines show the best-fit relation from Carilli:13 given in Eq. 8 together with the dispersion in this relation.

4.4 CO(3-2) Sizes

A detailed discussion of the sizes of the CO emitting regions in reddened Type 1 quasars is beyond the scope of this paper and will have to wait for higher resolution observations of the gas distribution in these sources. However, we have seen that the CO emission is already resolved in one of our reddened quasars - ULASJ23150143 (Section 3.4.2) - with an implied source diameter of 20 kpc. There are hints also that the gas emission is marginally resolved in ULASJ1234+0907 and VHSJ21015943. We used casa to fit a 2-D Gaussian (in RA, Dec) to the CO maps constructed by summing all the channels covering the full extent of the line. The deconvolved source sizes along the major axis are (2.530.52) arcsec (204 kpc) and (1.340.29) arcsec (112) kpc for ULASJ1234+0907 and VHSJ21015943 respectively. We note that the images used for the size estimates were generated using natural weighting, which results in strong sidelobes in the image plane. We have therefore checked that the sizes quoted here are consistent within the error bars with the deconvolved sizes estimated from images generated using a Briggs weighting (robust=0.5), where the effects of the sidelobes should be minimised. Finally, for our most extended and highest S/N source, ULASJ23150143, we have also checked that the major axis derived from directly fitting to the visibility data using the casa task uvmodelfit, is once again consistent with the deconvolved sizes from the image plane quoted here.

The large implied sizes for the warm gas reservoirs in these quasars, while still highly uncertain, are nevertheless at odds with previous observations of gas in quasar host galaxies. In the few cases where the gas emission has been resolved in high redshift quasar host galaxies, it is distributed in compact regions of few kpc in size (e.g. Carilli:02; Walter:04; Riechers:11b; Wang:13; Willott:13). The conventional interpretation is that the molecular gas in quasar hosts resides in a nuclear region close to the accreting black hole, and its excitation properties can therefore be explained by a single temperature and density gas component. SMGs on the other hand often show evidence for substantially larger reservoirs of low excitation gas 10 kpc in size (e.g. Ivison:11; Riechers:11a; Hodge:12). Although the excitation properties of the gas in reddened quasars remains to be determined, we tentatively conclude that the CO sizes appear to more closely resemble SMGs rather than high redshift quasar hosts.

4.5 Molecular Gas Fractions

With the CO detections providing us with molecular gas masses and the CO velocity widths providing a crude estimate of the dynamical masses, we can estimate the molecular gas fractions and compare directly to gas fractions in both SMGs and optical quasars. A KS-test comparing the observed FWHM of the CO lines in SMGs, optical quasars and reddened quasars reveals no evidence for a difference in the distributions between the three populations. If the three populations indeed form an evolutionary sequence with the optical quasars representing more evolved systems and SMGs representing the earlier gas rich stages of massive galaxy formation, we would expect the reddened quasars to have gas fractions that are intermediate between SMGs and optical quasars. In Fig. 9 we show the FWHM of the CO line - a proxy for the dynamical mass - versus the CO(1-0) luminosities - a proxy for gas mass. We compare our sample to optical quasars from Coppin:08 and Simpson:12 and to the SMGs from Bothwell:13, selecting quasars and SMGs that are unlensed only. We also show the scaling relation between these two quantities derived for unlensed sources by Harris:12, as well as a simple relation that equates the gas mass to the dynamical mass as follows:

(9)

where the constant, depends on the geometry and inclination of the source, represents the line dispersion FWHM/2.35 for a Gaussian line, is the source size and is the gravitational constant. We consider two separate cases: (i) , kpc, =4.6 M (K km s pc), which is appropriate for disk galaxies and (ii) , kpc, =1 M (K km s pc), representing a virialised, spherical source.

Figure 9: FWHM of the CO line (proxy for dynamical mass) versus CO line luminosity (proxy for gas mass) for quasars and SMGs from the literature (Bothwell:13; Coppin:08; Simpson:12) compared to our reddened quasars (This work, Feruglio:14; Brusa:15). All sources shown here are unlensed and therefore no magnification correction needs to be applied to these data. The dashed line shows the best-fit relation to unlensed galaxies from Harris:12 while the dot-dashed lines show the simple relation given by Eq. 9 for both disk like and spherical geometries.

We find that the reddened quasars occupy a region roughly between the majority of the optical quasars and the SMGs in Fig. 9. However, the observed FWHM is also sensitive to inclination with thin rotating disks viewed face-on having narrower FWHM. Bothwell:13 have argued based on the fact that SMGs show a correlation between CO linewidth and line luminosity, that they are primarily thick disks/turbulent ellipsoids where the FWHM is not as dependent on inclination. However, optical quasars, where broad emission lines can be seen in the quasar spectra, may have an inclination bias relative to SMGs although the opening angles needed for detection of broad lines in quasars likely span a large range. If the molecular gas sits in a thin disk in these quasars, the optical quasars with very narrow FWHM that sit above the various lines in Fig. 9, could be seen more face-on. On the other hand, if the gas morphologies are more turbulent and inhomogenous - e.g. as would be expected for merging systems - the FWHM is unlikely to be as dependent on inclination. Further observations are required to established whether the reddened quasars have any inclination bias relative to optical quasars and SMGs. We note that there is a general trend among our objects for the more highly obscured quasars to have larger observed CO FWHM - e.g. ULASJ1234+0907, the dustiest quasar in our sample (A mag; B12) has the broadest CO line whereas VHSJ21015943, the least dusty quasar observed (A mag; B15), has the narrowest CO line, suggesting that some of the reddening in these quasars could be linked to orientation. However, in Section 3.4.2 we have seen that in at least one of our quasars the gas appears to be distributed in a disk at an inclination angle of 57. This is consistent with the inclination angles inferred for SMGs. If SMGs, reddened quasars and optical quasars have similar source sizes and inclinations on average, for a fixed dynamical mass or CO FWHM, reddened quasars have gas masses that are intermediate between SMGs and optical quasars, and are therefore on their way to exhausting their gas supply and evolving into optical quasars. Further spatially resolved observations of the reddened quasars to get rotational dynamical parameters as well as internal velocity dispersions, is clearly the next step forward in order to get robust gas fractions in these sources.

5 Conclusions

We have presented ALMA Band 3 observations of a sample of four heavily reddened (A mag), high luminosity broad line quasars at selected from wide field near infrared surveys. The aim of these observations is to test whether such reddened quasars could represent the evolutionary link between starburst galaxies and unobscured, optically bright quasars at high redshifts. All four quasars are detected in both the dust continuum (at observed frame wavelengths of 3mm) and in the CO(3-2) line. We infer star formation rates, dust and gas masses from these observations and attempt to place these measurements in context with corresponding values for other high redshift galaxy populations:

  • The dust masses range from 2.510M to 1.210M. The FIR luminosities are 10L and the implied star formation rates in the quasar host galaxies is M yr, suggesting that these reddened quasars reside in prodigiously star forming host galaxies with very high dust luminosities.

  • We estimate gas masses from the measured CO(3-2) line luminosities, assuming an excitation ratio =0.8 and a CO-to-H conversion factor of (K km s pc). These gas masses are in the range 10M and the gas-to-dust ratios are in the range 30-110.

  • Instead, adopting a gas-to-dust ratio of 91, as found in nearby galaxies, dust- and CO-based gas masses agree within a factor of 2 for three of the quasars. For ULASJ2315+0143, the dust-based gas mass is 3 times larger. The CO- and dust-based masses could be made consistent for ULASJ2315+0143 if a CO-to-H conversion factor of (K km s pc) applies. Such a value is intermediate between that found in nearby galaxies/high-redshift star-forming disks and in nuclear starbursts/quasar host galaxies.

  • The CO linewidths span kms, comparable to linewidths in other high-redshift SMGs and quasars. Using several different estimators, dynamical masses of Msin(i) 310M to 310M are found.

  • ULASJ2315+0143 () shows evidence for spatially resolved CO emission over 20 kpc, with a strong velocity gradient across the source. We demonstrate that the velocity field is consistent with the presence of a large, rotating gas disk although other scenarios such as mergers cannot be ruled out at the resolution of the current data. ULASJ23150143 is the quasar with the largest discrepancy between its dust and CO based gas masses, and where there is evidence for (K km s pc).

  • The sizes of the CO(3-2) emitting regions appear to be large ( 10 kpc) hinting at the presence of extended reservoirs of warm gas. Gas emission in high-redshift, unobscured quasars is generally compact (few kpc). If extended emission is confirmed via higher resolution observations, the reddened quasars would instead be more similar to the properties of SMGs.

  • ULASJ1234+0907 () has two other millimetre bright galaxies at the same redshift in the field-of-view, with implied gas masses of a few times 10M. The two galaxies are located 90 kpc and 170 kpc from the quasar. ULASJ1234+0907 may reside in a significant over-density at .

  • A similar fraction of reddened quasars are “submillimetre loud” compared to unobscured UV luminous quasars. However the average dust luminosities of reddened quasars would appear to be higher than for the UV luminous quasar population at similar redshifts. Reddened quasars also appear to have higher cold dust luminosities compared to the recently discovered population of Hot Dust Obscured Galaxies from WISE. So far, 100 per cent of our reddened quasars (6/6) have been observed to have very large molecular gas reservoirs. A direct comparison to the unobscured quasar population is not possible as surveys of molecular gas in UV luminous quasars have focussed on the subset of these that were already known to be millimetre bright, therefore increasing the probability of detecting molecular gas. Nevertheless, several UV luminous, millimetre bright quasars with comparable black hole masses and AGN luminosities to our sample, have limits on their molecular gas masses that place them below the gas masses seen so far in the reddened quasar population, indicating that these unobscured quasars could indeed correspond to a later evolutionary phase. Further observations of molecular gas in larger samples of both reddened and unobscured quasars matched in luminosity and black hole mass would help confirm this trend.

  • Using the CO linewidths in our sample as a proxy for the dynamical mass, and the CO line luminosities as a proxy for gas mass, we find that the reddened quasars have molecular gas fractions that are intermediate between SMGs and optical quasars assuming that all three populations have similar source sizes and inclinations. It therefore seems plausible that reddened quasars are being seen as they are transitioning to optical quasars from starburst galaxies. Their supermassive black holes already appear to be fully assembled and, given the rate of ongoing star formation, they will very quickly exhaust their gas reservoirs and presumably blow out the surrounding gas and dust on their way to evolving into some of the most luminous, unobscured quasars - progenitors of the largest supermassive black holes seen in the Universe today.

Overall, our observations demonstrate the high scientific return possible from ALMA observations of luminous quasars at the main epoch of galaxy formation. Studies of the host galaxy properties of intrinsically luminous, unlensed quasars at are still scarce and further observations with ALMA both in terms of assembling larger samples of such sources and obtaining higher spatial resolution detections, will help shed light on the exact nature of the host galaxies of the largest supermassive black holes, seen as they are being assembled in the high redshift Universe.

Acknowledgements

The authors thank the referee for a useful and detailed report which helped improve the paper. MB would like to thank Matt Bothwell and Stefano Carniani for useful discussions. MB acknowledges funding from the UK Science and Technology Facilities Council (STFC) via an Ernest Rutherford Fellowship. GJ is grateful for support from NRAO through the Grote Reber Doctoral Fellowship Program. RGM and PCH acknowledge funding from STFC via the Institute of Astronomy, Cambridge Consolidated Grant. SA-Z acknowledges support from Peterhouse, Cambridge.

This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.01247.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.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
124616
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description