{'fraction_non_alphanumeric': 0.04556170097016955, 'fraction_numerical': 0.05417535587995285, 'mean_word_length': 4.388128969223254, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Karl G. Jansky Very Large Array observations of cold dust and molecular gas in starbursting quasar host galaxies at z~4.5

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%.


Introduction
The high-redshift formation of some of the most massive present-day galaxies is often accompanied by episodes of extreme far-infared (FIR) luminosity, as high-lighted by the phenomena of submm/mm bright quasar host galaxies and starbursting submm galaxies (SMGs). The interpretation of dust-heated star-formation to explain the rest-frame FIR continuum properties of these populations is supported by the detection of redshifted [CII] line emission (e.g. Maiolino et al. 2005;Stacey et al. 2010;Venemans et al. 2012;Wang et al. 2013), as this is typically the strongest cooling line in the FIR spectrum of nearby galaxies (e.g. Stacey et al. 1991). Fueling the star-formation and active galactic nucleus (AGN) activity requires significant molecular gas reservoirs, most efficiently studied through observations of redshifted CO line emission (see  for a recent review). Together, these line and continuum observations can probe the physical conditions and kinematic properties of the interstellar medium in galaxies over much of the history of the Universe.
Recent studies of gas and dust in high-redshift galaxies have used photoionization models along with measured values of the relative intensity between the FIR continuum, low-J CO and [CII] line emission to constrain the ionization rate and density of the interstellar medium (e.g. Stacey et al. 2010;Ivison et al. 2010). With new facilities like the Atacama Large Millimeter/submillimeter Array (ALMA) and the upgraded VLA it is now possible to perform similar analyses at kpc scale resolution with some of the most distant starburst galaxies and AGN, while also constraining their kinematic properties (e.g. Riechers et al. 2011;Ivison et al. 2011;Wang et al. 2013;. Furthermore, with the combined continuum sensitivity of these interferometers one can constrain the emissivity of the thermal dust continuum emission and obtain independent estimates of the star-formation rates (SFRs) through detection of thermal free-free emission.
The FIR luminous quasar host galaxies, BRI1335-0417 and BR1202-0725, are two wellstudied cases of massive, unlensed galaxy formation roughly 1.5 Gyr after the big bang. BRI1335-0417 is believed to be a "wet" merger at z = 4.4 with an estimated FIR luminosity of L F IR ∼ 3 × 10 13 L ⊙ fueled by ∼ 9 × 10 10 M ⊙ of cold molecular gas traced by CO line emission (Storrie-Lombardi et al. 1996;Ohta et al. 1996;Guilloteau et al. 1997;Benford et al. 1999;Carilli et al. 1999Carilli et al. , 2002Riechers et al. 2008). Similarly, BR1202-0725 at z = 4.7 is a quasar host galaxy discovered in the APM survey (Irwin et al. 1991) and has been shown to be comprised of two FIR luminous radio sources separated by ∼3.8 ′′ (Omont et al. 1996;Yun et al. 2000). Both the quasar host galaxy BR1202-0725 South and the optically faint SMG companion to the North have FIR luminosities in excess of 10 13 L ⊙ and a combined cold molecular gas mass of ∼10 11 M ⊙ (Carilli et al. 2002;Iono et al. 2006;Riechers et al. 2006). Within the same region, Hu et al. (1996) identified two fainter Lyα companions, both of which are thought to exhibit faint 157.7µm [CII] line emission . Submm interferometric observations of [CII] and dust continuum emission with ALMA show that the Lyman-α companion to the SW of the quasar is also luminous in the FIR (Wagg et al. 2012;R. Williams et al. in prep.). It is likely that BR1202-0725 is at an early stage of a merger between multiple gas-rich systems (Salomé et al. 2012).

Observations and data reduction
Observations were carried out with 27 elements of the VLA on four dates in January and February of 2013, while the array was in the most compact D-configuration. The rest frequency of the CO J =2-1 line is 230.538 GHz, so that at the redshifts of our targets this line can be observed at ∼40.5 GHz for the two most luminous galaxies in the BR1202-0725 system, and at 42.6 GHz in BRI1335-0417. Each field was observed for a total of 8 hours including overheads associated with complex gain, bandpass and flux density calibration. In the case of each target field, 3C286 was used to calibrate the flux density and bandpass shape, while J1229+0203 and J1354-0206 were used to calibrate the complex gains for BR1202-0725 and BRI1335-0417, respectively. Using the recently commissioned 8 GHz bandwidth correlator mode with the Q-band receivers our observations covered 40 to 48 GHz with a spectral resolution of 2 MHz. The frequency coverage in this mode is not contiguous, and small 8 MHz gaps between sub-bands mean that two slightly offset tunings were used to fully sample the entire range.
The data reduction was performed using the CASA software package 1 . Standard calibration steps are carried out using a new set of CASA pipeline reduction procedures developed by NRAO staff. Including only the spectral windows without line emission, we then use natural weighting and multi-frequency synthesis imaging to generate the continuum images, achieving sensitivities of 5.9µJy and 6.5µJy beam −1 for BR1202-0725 and BR1335-0417, respectively. For the spectral line cubes, we also used natural weighting in the imaging and the synthesized beamsizes are 2.55 ′′ × 1.95 ′′ (PA = 2.3 • ) and 2.57 ′′ × 1.75 ′′ (PA = 7.4 • ) in the two images. The spectral line maps are resampled to 8 MHz resolution, corresponding to velocity channel widths of 59 and 56 km s −1 at the redshifted frequencies of the CO J =2-1 line in BR1202-0725 and BRI1335-0417, while the rms per channel is 160 and 170µJy, respectively. Figure 1 shows the 44 GHz continuum images of the two target fields, where we also overlay contours of the 340 GHz continuum observed in BR1202-0725 with ALMA. The two most luminous submm sources in the BR1202-0725 field are detected at 44 GHz, with peak flux densities of 53 ± 6µJy and 24 ± 6µJy, for the Northern and Southern components. The integrated flux densities are similar to the peak values for both components so they are unresolved. BRI1335-0417 is also detected with a peak flux density of 40 ± 7µJy, which is in a The CO J=2-1 redshifts and line widths (∆V F W HM ) are determined from the best-fit parameters of a Gaussian fit to the spectra. agreement with the prediction by Carilli et al. (1999) who assume a thermal dust emissivity index, β = 1.5. Future higher spatial resolution continuum observations will determine the size of this emission region.

Results
From the spectral line data cubes we extract spectra of the CO J=2-1 line emission at the positions of our three primary targets ( Figure 2). For comparison, we show the [CII] line profiles for the BR1202-0725 SMG and quasar host galaxy measured from recent ALMA commissioning observations (Wagg et al. 2012;Lentati et al. 2013). Table 1 gives the redshift, line widths and integrated intensities of the CO J=2-1 line emission observed here with the VLA, which have been calculated assuming the data can be decribed by a single Gaussian line profile. However, based on the [CII] line observations and the CO J=5-4 and CO J=7-6 line emission observed with the Plateau de Bure Interferometer (PdBI) by Salomé et al. (2012), the line profiles of BR1202-0725 North appear to be best described by two Gaussian components. The CO J=2-1 redshift and line width for the quasar BR1202-0725 South are in excellent agreement with previous observations of [CII] and CO line emission (Omont et al. 1996;Carilli et al. 2002;Wagg et al. 2012), however these data hint at the presence of broad wings in both sources, and more sensitive observations are needed to confirm this.
In the case of BRI1335-0417, the CO J=2-1 line profile measured here is slightly narrower than that of the previous [CII] and CO J=5-4 lines detected (∼ 430 km s −1 ; Guilloteau et al. 1997;Wagg et al. 2010), and also blueshifted by ∼50 km s −1 . The integrated intensity of the CO J=2-1 line emission measured here (0.38±0.03 Jy km s −1 ) is in good agreement with the previous observations reported by Riechers et al. (2008), who measure 0.43±0.02 Jy km s −1 using the old VLA correlator. This implies that the narrower bandwidth of past observations was wide enough to encompass the velocity width of the line, suggesting that the differences in the linewidths can be attributed to the low signal-to-noise of the [CII] and CO J=5-4 spectra.

Continuum emission
We use archival Herschel SPIRE submm-wavelength and lower frequency, 1.4 and 5 GHz VLA observations of our targets (project AC878) to better sample their (rest-frame) FIR through cm-wavelength spectral energy distributions. For BR1202-0725 the FIR spectral energy distribution has been observed with ALMA, PdBI and Herschel SPIRE. From the level 2 pipeline processed SPIRE maps we extract flux densities at the pixels corresponding to the 44 GHz peak position of BRI1335-0417, measuring, S 250µm = 33 ± 6 mJy, S 350µm = 41 ± 5 mJy, and S 500µm = 44 ± 6 mJy, while in the case of BR1202-0725 the total 350 µm flux density of both components measured in the SPIRE beam is S 350µm = 78 ± 6 mJy. Figure 3 shows the observed spectral energy distributions for BR1202-0725 and BRI1335-0417.
In order to model the spectral energy distribution we assume the data can be described by a power-law synchrotron component dominating at frequencies below ∼30 GHz, a thermal free-free component with a power-law index of α f f = −0.1, and a single temperature greybody component which is described by a dust temperature, T d and emissivity index, β. Following Yun & Carilli (2002), we use the linear relationships between the total SFR and the three components describing the radio-to-FIR spectral energy distribution (Condon 1992) to obtain simultaneous constraints on all parameters in the model. The parameters we fit are therefore synchrotron spectral index, ν −α , a scaling factor for the non-thermal synchrotron emission (f nth ; see equation (13) of Yun & Carilli 2002), and finally the SFR, T d and β. Our fitting method is based on the principles of Bayesian inference, which provides a consistent approach to the estimation of a set of parameters Θ in a model or hypothesis H given the data, D. Bayes' theorem states that: In parameter estimation, the normalizing evidence factor is usually ignored, since it is independent of the parameters Θ. Inferences are therefore obtained by taking samples from the (unnormalised) posterior using, for example, standard Markov chain Monte Carlo (MCMC) sampling methods. An alternative to MCMC is the nested sampling approach (Skilling 2004), a Monte-Carlo method targeted at the efficient calculation of the evidence, that also produces posterior inferences as a by-product. In Feroz et al. (2008Feroz et al. ( , 2009) this nested sampling framework was developed further with the introduction of the MultiNest algorithm, which provides an efficient means of sampling from posteriors that may contain multiple modes and/or large (curving) degeneracies, and also calculates the evidence. We make use of the MultiNest algorithm to obtain our estimates of the posterior probability distributions for the spectral energy distribution parameters. Naturally, the model parameters are correlated. For example, when fitting thermal dust continuum emission there are degeneracies between the emissivity index, β, and T d (e.g. Priddey & McMahon 2001;Blain et al. 2003). These correlations can be seen in the  Far-infared spectral energy distributions of BR1202-0725 North and South along with the best-fitting greybody, synchrotron and thermal free-free models described in the text and parameters given in Table 2. Observed data are from Carilli et al. (2002), this work, Salomé et al. (2012), and Wagg et al. (2012). The solid curves show the best-fit models while the dotted lines show the synchrotron and thermal dust components. The thermal freefree component of the model is not shown but does make a small contribution to the solid line. bottom: The observed spectral energy distribution of BRI1335-0417 plotted with the best-fitting model (solid line) composed of the thermal dust continuum and synchrotron emission curves (dotted lines), and only a small contribution from the thermal free-free emission. Observed data are from Carilli et al. (1999), this work, and Guilloteau et al. (1997). probability density plots shown in Figure 4 and the mean of the posterior parameter values are given in Table 2. Although we impose the prior that the total 350 µm emission in BR1202-0725 North and South should be equal to the flux density measured by SPIRE, the relative temperatures of these two components is unknown. In an attempt to overcome this, we include priors on the dust temperatures based on recent submm-wavelength observations of high-redshift starburst galaxies (T d = 40 K; e.g. Magnelli et al. 2012) and quasar host galaxies (T d = 50 K; Beelen et al. 2006;Wang et al. 2008). However, we find that the resulting fits do a poor job of reproducing the measured 44 GHz flux densities, and so we adopt the values derived in the free floating parameter analysis.
Evidence for thermal free-free emission has recently been observed in gravitationally lensed, FIR luminous starburst galaxies (Thomson et al. 2012;Aravena et al. 2013). The models provide an estimate of the contribution from such free-free emission to the observed 44 GHz continuum emission in our targets. Based on the fitting analysis, the main contribution to the 44 GHz flux density in each source is thermal dust continuum emission. In the case of BRI1335-0417, that contribution is 54%, with 30% coming from synchrotron emission and only 16% from thermal free-free emission. In the case of the SMG BR1202-0725 North, the thermal dust emission contribution from our model fit is 75%, while 17% and 8% come from synchrotron and free-free emission, respectively. The model fit to the spectral energy distribution of the quasar BR1202-0725 South has a 67% contribution to the 44 GHz flux density from thermal dust emission, while only 12% comes from synchrotron emission and 21% is found to come from free-free emission. Although it is not considered in our models, another possibility is that the 44 GHz emission arises from a flat-spectrum AGN embedded in the starburst host galaxies, however higher resolution imaging is needed before we can determine the likelihood of this.

CO J=2-1 Line Emission
Following the definition of line luminosity given by Solomon et al.( 1992), we calculate the line luminosities to be L ′ CO = (8.7 ± 0.8) × 10 10 K km s −1 pc 2 and L ′ CO = (6.0 ± 0.5) × 10 10 K km s −1 pc 2 for BR1202-0725 North and South, respectively, and L ′ CO = (7.3 ± 0.6) × 10 10 K km s −1 for BRI1335-0417. For BR1202-0725, these luminosities are in excellent agreement with previous lower spectral resolution observations with the VLA (Carilli et al. 2002), and with the high-J CO lines measured by Omont et al. (1996) and Salomé et al. (2013), and supports the claim that the line emission is thermalized. This is also the case for BRI1335-0417, for which the CO J=5-4 line luminosity, L ′ CO = (8.5±0.9)×10 10 K km s −1 pc 2 measured by Guilloteau et al. (1997) is similar to what we measure in CO J=2-1 . ALMA observations reveal [CII] line emission in the two LAEs at z ∼ 4.7 associated with the BR1202-0725 system . We search the Q-band spectral line data cubes at the positions of Lyα-1 and Lyα-2 and the expected frequencies of redshifted CO J=2-1 , but this emission is undetected. Assuming Gaussian line profiles we set 3-σ upper limits to the line luminosities, L ′ CO < 4.8 × 10 9 K km s −1 pc 2 and L ′ CO < 1.2 × 10 10 K km s −1 pc 2 for Lyα-1 and Lyα-2 assuming FWHM linewidths of 56 and 338 km s −1 .

Discussion
Although the uncertainties are large, the star-formation rates estimated for the three targets are similar to previous estimates derived from the FIR luminosities (e.g. Wagg et al. 2010;Salomé et al. 2012;Carniani et al. 2013). The FIR emission model fit to the SMG BR1202-0725 North favours an emissivity index of β ∼ 1.5, while the quasars BR1202-0725 South and BRI1335-0417 appear to have values closer to β ∼ 2.0. Both of these values are similar to what has been measured in 60µm selected galaxies in the nearby Universe (e.g. Dunne et al. 2000). More recent studies of the full FIR through submm-wavelength emission in normal star-forming galaxies from the KINGFISH sample show that nearby galaxies can exhibit a broad range in emissivity indices, β ∼0.8-2.5 (Galametz et al. 2012). Our measured values are therefore not unusual. If free-free emission does not contribute significantly to the 44 GHz emission observed in BR1202-0725 North, then these data suggest that the cold dust mass is a factor of ∼2.5× higher than that in the quasar BR1202-0725 South.
In all three sources, the contribution from thermal free-free emission to the rest-frame 250 GHz continuum is < ∼ 20%. Free-free emission typically contributes less than 10% of the 230 GHz continuum emission observed in normal, nearby galaxies (e.g. Albrecht et al. 2007), and so our results are consistent with what is observed in these more quiescent starforming galaxies. Strong free-free emission has been detected in lensed SMGs (Aravena et al. 2013), and the SFRs estimated are in good agreement with those estimated from the FIR luminosities.
Our measured CO J=2-1 line luminosities are in good agreement with previous observations of both the low and high-J CO emission observed in these objects (Guilloteau et al. 1997;Carilli et al. 2002;Riechers et al. 2006Riechers et al. , 2008. The luminosities imply cold molecular gas masses ∼10 11 M ⊙ assuming the recently measured value α CO = 0.6 M ⊙ [K km s −1 pc 2 ] −1 for the CO-to-H 2 conversion factor in ULIRGs (Papadopoulos et al. 2012), similar to the original value of Downes & Solomon (1998). However, we note that this factor can be uncertain by ∼2. A recent analysis of the [CII] line emission observed in BR1202-0725 infers dynamical masses for the quasar and SMG based on the kinematic properties of the gas (Carniani et al. 2013). They calculate dynamical masses of M dyn ∼ 4 × 10 10 M ⊙ and ∼ 6 × 10 10 M ⊙ , for the quasar and SMG, implying very high molecular gas mass fractions of ∼ 80 − 90%. Such high gas fractions are inferred for some local ultraluminous infared galaxies (e.g. Papadopoulos et al. 2012). Planned sub-arcsecond imaging of this cold molecular gas in these targets would be complemented by higher resolution imaging of the (rest-frame) FIR continuum and [CII] line emission with ALMA in order to infer the ionizing radiation field and gas density on kpc scales.

Summary
We present VLA observations of redshifted CO J=2-1 line and thermal dust continuum emission in high-redshift quasar host galaxies, observing BR1202-0725 and BRI1335-0417. The luminosities measured in the CO J=2-1 line detected in all three targets are similar to previous observations of CO J=5-4 line emission, implying that the gas kinetic temperature is close to the excitation temperature. The observations reveal strong (rest-frame) 250 GHz continuum emission associated with all three FIR luminous galaxies, with stronger emission observed in the two starburst galaxies than the quasar BR1202-0725 South. Although synchrotron and thermal free-free likely contribute to some of this emission, we intrepret the remainder as due to thermal dust emission. Future observations at frequencies between 10 and 40 GHz will allow us to better constrain the relative importance of the three processes.
We thank Maud Galametz, Rob Kennicutt, Thomas Greve and Padelis Papadopoulos for helpful discussions, and the anonymous referee for their comments on the submitted manuscript. This work was co-funded under the Marie Curie Actions of the European Commission (FP7-COFUND). We thank all those involved in the VLA project for making these observations possible (project code 13A-012). This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00006.SV. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), 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.
This preprint was prepared with the AAS L A T E X macros v5.2.   Table 2.