Unveiling Sizes of Compact AGN Hosts with ALMA

, , , , , , , , , and

Published 2020 January 7 © 2020. The American Astronomical Society. All rights reserved.
, , Citation Yu-Yen Chang et al 2020 ApJ 888 44 DOI 10.3847/1538-4357/ab595b

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/888/1/44

Abstract

We present rest-frame far-infrared (FIR) and optical size measurements of active galactic nucleus (AGN) hosts and star-forming galaxies (SFGs) in the COSMOS field, enabled by high-resolution Atacama Large Millimeter/submillimeter Array (ALMA)/1 mm (0farcs1–0farcs4) and Hubble Space Telescope (HST)/F814W imaging (∼0farcs1). Our sample includes 27 galaxies at z < 2.5, classified as infrared-selected AGN (three sources), X-ray selected AGN (four sources), and non-AGN SFGs (20 sources), for which high-resolution Band 6/7 ALMA images are available at 1 mm from our own observing program as well as archival observations. The sizes and star formation rate surface densities measured from both ALMA/1 mm and HST/F814W images show that obscured AGN host galaxies are more compact than non-AGN SFGs at similar redshift and stellar mass. This result suggests that the obscured accretion phase may be related to galaxies experiencing a compaction of their gaseous component, which could be associated with enhanced central star formation before a subsequent quenching driving the formation of compact passive galaxies. Moreover, most of the detected and stacked rest-frame FIR sizes of AGNs in our sample are similar or more compact than their rest-frame optical sizes, which is consistent with recent results of ALMA-detected sources. This might be explained by the fact that the dusty starbursts take place in the compact regions, and suggests that the star formation mechanisms in the compact regions of AGN hosts are similar to those observed in SFGs observed with ALMA.

Export citation and abstract BibTeX RIS

1. Introduction

The processes driving the coevolution of galaxies and their supermassive black holes remain a largely debated issue in extragalactic astrophysics. Specifically, the connection between the formation of stars in galaxies and the fueling of their central black holes is still not fully understood. Previous results revealed that X-ray and optically selected active galactic nuclei (AGN) reside in galaxies harboring sustained activity of star formation (e.g., Harrison et al. 2012; Mullaney et al. 2012; Juneau et al. 2013; Rosario et al. 2013). Furthermore, although high luminosity AGNs seem to be connected to violent events (i.e., quasars) and may reside in galaxies more often involved in major mergers (e.g., Treister et al. 2009; Kartaltepe et al. 2012), the majority of moderate-luminosity X-ray selected AGNs appear to live mostly in disk-dominated isolated systems (Mainieri et al. 2011; Silverman et al. 2011; Fan et al. 2014; Villforth et al. 2014). Also, major mergers do not appear to be the dominating triggering mechanism for luminous unobscured AGN (Mechtley et al. 2016; Villforth et al. 2017). Their morphology should thus be very similar to that of sources populating the so-called main sequence (stellar mass versus star formation rate (SFR); e.g., Brinchmann et al. 2004; Daddi et al. 2007; Elbaz et al. 2007; Noeske et al. 2007; Ilbert et al. 2015; Schreiber et al. 2015; Chang et al. 2015) of star-forming galaxies (SFGs; e.g., Gabor et al. 2009; Cisternas et al. 2011; Kocevski et al. 2012), implying that the bulk of supermassive black hole accretion is likely driven by internal processes and not by major mergers. Nonetheless, it is still an open question whether the internal structure of AGN hosts is comparable to or different from that of normal SFGs.

Most current studies have focused on samples of SFGs. Recent observations from the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to measure far-infrared (FIR) sizes with high angular resolution and sensitivity (e.g., Ikarashi et al. 2015; Simpson et al. 2015a, 2015b; Harrison et al. 2016; Hodge et al. 2016). For size comparison, it has been discussed that rest-frame FIR sizes of massive SFGs are smaller than their rest-frame optical sizes (e.g., Chen et al. 2015; Barro et al. 2016; Fujimoto et al. 2017, 2018; Tadaki et al. 2017; Elbaz et al. 2018; Lang et al. 2019). A possible explanation of the central compact component is a compact dusty bulge which is related to bulge formation in the center (Fujimoto et al. 2017). Moreover, it has been discussed that massive galaxies turn out to exhibit compact star formation (Barro et al. 2016; Elbaz et al. 2018), and transform into compact quiescent galaxies at z ∼ 2 (e.g., van der Wel et al. 2014). There have been several attempts to estimate FIR sizes for faint submillimeter galaxies (Rujopakarn et al. 2016; Fujimoto et al. 2017; González-López et al. 2017), but the uncertainties are still constrained by the small number statistics. Moreover, it has been discussed that both high sensitivities and high angular resolutions are required to measure FIR sizes for less massive objects (Fujimoto et al. 2018). On the other hand, it is still unclear whether the existence of AGNs affects the structures of their host galaxies, and whether AGN host galaxies show different rest-frame FIR sizes compared to other SFGs.

Chang et al. (2017b) investigated AGN and non-AGN host morphologies through Hubble Space Telescope (HST)/F814W images, and found that obscured IR-selected AGN hosts are more compact than normal star-forming sources at a given stellar mass at z ∼ 1. Following that work, in this paper we resolve the dust emission of AGN host galaxies using high-resolution ALMA imaging. We show that their compact sizes are consistent with what we found using HST/F814W-band imaging, and the FIR (ALMA/ 1 mm) sizes are yet even smaller than the optical sizes. We compare all available FIR sizes of AGN and non-AGN hosts in the Cosmic Evolution Survey (COSMOS; Scoville et al. 2007) field. For comparison, we also measure the FIR sizes of a control sample of non-AGN SFGs in the COSMOS field for which we find complementary high-resolution ALMA imaging in the archive. The structure of this paper is as follows. We describe the data and sample selections in Section 2. We analyze the physical and structural properties in Section 3. We discuss the results in Section 4 and summarize in Section 5. Throughout the paper, we use AB magnitudes, adopt the cosmological parameters (ΩM, ΩΛ, h) = (0.30, 0.70, 0.70), and assume the stellar initial mass function of Chabrier (2003).

2. Samples and Observations

2.1. Our ALMA Observations

Our ALMA Band 7 (275–373 GHz) observations (#2016.1.01355.S) contain 10 obscured AGN host galaxies selected in the infrared (based on identifying mid-IR power-law signatures which are characteristic of AGNs as opposed to star-forming dominated sources; more details in Chang et al. 2017a) at 0.5 < z < 1.0 in the COSMOS field. Our targets were selected along the main sequence within 0.3 dex of star formation, as well as a star-forming sample in the rest-frame UVJ plot (Williams et al. 2009), with stellar masses $\mathrm{log}\,{M}_{* }\gt 10.5$ in a 24 μm selected sample in the COSMOS field. Reliable spectroscopic redshifts are available. In order to ensure a robust estimate of the dust-obscured star formation, we restricted our sample to Herschel-detected sources. We confirmed that these objects show compact F814W-band radial profiles as described in Chang et al. (2017b).

We proposed the standard set-up for single continuum centered at 343.5 GHz of Band 7 over a 7.5 GHz bandwidth. The proposed angular resolution is 0farcs1, which is close to the resolution of the F814W-band imaging. The array configuration was C40-6 with 40 antennas, and the proposed rms value is 0.1 mJy/beam per 7.5 GHz. The data were reduced with the Common Astronomy Software Applications package (CASA; McMullin et al. 2007), and the steps of data reduction and calibration were performed with the scripts provided by the ALMA observatory. The entire map was produced by the CLEAN algorithm with the TCLEAN task. The final images are characterized by a synthesized beam size of 0farcs12 × 0farcs1012 as shown in Table 3. One of the objects (Source name: 972308) is detected, and is shown in Figure 1 along with the remaining detections from the ALMA archive. The rms value is 0.096 mJy/beam per 8 GHz bandwidth (∼0.099 mJy/beam per 7.5 GHz). For the remaining nine objects, we stacked their imaging and computed the median values in each pixel to create a stacked imaging as shown in Figure 2. The peak flux is 0.09 mJy ± 0.03 mJy (3.39σ; rms value is calculated from the stacked image).13 We also stacked these objects by STACKER (Lindroos et al. 2015) in the uv plane. The peak flux is 0.08 mJy ± 0.03 mJy (2.24σ; rms value is calculated from the output image of STACKER), which is consistent with the stacking results from the image plane. For point sources, nine objects with an rms of 0.1 mJy/beam can produce at best a stacked image with an rms of 0.03 mJy/beam, and given that the stacked signal is 0.09 the signal-to-noise ratio (S/N) would be around 3. Our stacking result from the image plane shows slightly large signal-to-noise ratios, which could be explained by the extended profile and the beam size. It is not straightforward to convert a stacked image in the uv plane from angular scale to physical scale, so we adopted the FIR stacked image in the image plane in this paper.

Figure 1.

Figure 1. HST/F814W-band images for the 27 detected sources. The green contours show the ALMA band 6/7 intensities from the 3σ to 24σ with a 3σ step. The yellow circles represent the beam size. The white circles are the 2'' radii of the optical objects. The numbers are the COSMOS2015 IDs. The infrared-selected AGNs are labeled as IRAGN, and the X-ray selected AGNs are labeled as XAGN. Our detected source (972308) is in the blue frame.

Standard image High-resolution image
Figure 2.

Figure 2. Stacked image of the nine undetected sources. The background is the HST/F814W image, and the green contour shows the ALMA/1 mm intensities of the 3σ detection. The yellow circle represents the beam size. The box size is 5 kpc × 5 kpc.

Standard image High-resolution image

2.2. Complementary Sample from the Archive

We searched all available ALMA public data by 2019 May to acquire a complementary sample. We found 27 ALMA-detected sources as shown in Figure 1, including AGNs and non-AGN galaxies in Band 6/7 with angular resolution smaller than 0farcs4 and peak flux values larger than 3σ at z < 2.5, and the one in our program. Table 3 shows the properties of these archival data (#2013.1.00034.S, #2015.1.00026.S, #2016.1.00478.S, #2016.1.00735.S, #2016.1.00778.S, #2016.1.01355.S, #2016.1.01426.S, and #2016.1.01604.S.). The angular resolution is in the range of 0farcs1–0farcs4. We simulated ALMA images by assuming Gaussian profiles, observed the images by SIMOBSERVE, and analyzed the observations by SIMANALYZE with similar inputs as our program (single continuum centered at 343.5 GHz of Band 7 over a 7.5 GHz bandwidth) in CASA. Then we used the UVMODELFIT task to fit single component Gaussian models to the UV data, and calculated the fitted sizes. We assumed sizes in a range of 0farcs10–0farcs40 (0.8 kpc–3.2 kpc at z = 1; 0.9 kpc–3.3 kpc at z = 2), as well as fluxes in a range of 0.25 mJy–2.00 mJy (S/Ns in a range of 4–50 with 300 s total time). For high-resolution observations (0farcs1–0farcs4) close to our program, the size measurements are increasing from underestimation to overestimation. The underestimation can be caused by over-resolving low S/N objects, and the overestimation can be caused by insufficient resolutions. Nevertheless, the measured size differences are within 20% of the simulated galaxies if the resolution is around 0farcs1–0farcs4. For spatial resolution larger than 0farcs5, the size measurements from UVMODELFIT are getting larger and overestimated by more than 30% size differences, and the profiles sometimes fail to be fitted, especially for small sources or low signal-to-noise ratio objects. We also did a simulation by IMSMOOTH function in CASA to create images with 0farcs4 angular resolution for images with 0farcs1 angular resolution. Most sources (>90%) are well-fitted and their differences of sizes are within 10%, so we confirmed that these data can be complementary to our program.14 We matched these sources with the COSMOS2015 catalog (Laigle et al. 2016). The redshifts we adopted are first taken from the Chandra COSMOS Legacy Survey (14 sources; Salvato et al. 2011; Civano et al. 2016; Marchesi et al. 2016), then from the COSMOS2015 photometric redshifts (13 sources; Laigle et al. 2016). Among the 27 sources as shown in Table 3, 13 of them have reliable spectroscopic redshifts, 23 of them are detected by Herschel (S/N > 3 in any Herschel bands), 3 of them are X-ray selected AGNs, and 3 of them are infrared-selected AGNs according to Chang et al. (2017a).

2.3. Spectral Energy Distribution (SED) Fitting Results

We fitted the SEDs of these 27 sources from the UV to infrared with an updated version of the MAGPHYS fitting technique (da Cunha et al. 2008, 2015). The version we used accounts for a possible AGN contribution assuming empirical AGN components from a set of empirical templates as described in Chang et al. (2017a), Mullaney et al. (2011, type-2), Richards et al. (2006, QSO), Prieto et al. (2010, QSO), and Polletta et al. (2007, Seyfert 1). We included ALMA detections or upper limits in addition to the photometry in Chang et al. (2017a). The stellar masses of these 27 sources are all above $\mathrm{log}\,{M}_{* }=10.5$, which is complementary to our selected sample. It has been known that the energy balance for UV and FIR in SED can break down because the dominant optical and IR emission of these galaxies do not come from the same physical regions (see Simpson et al. 2017, for a recent example). Therefore, we calculated the decomposed infrared and UV luminosity from the star formation component after subtracting the AGN contribution by our SEDs. Then we calculated the SFRs by the sum of these decomposed infrared and UV SFRs according to Kennicutt (1998) and Kennicutt & Evans (2012). We found that the SFRs calculated by this method are about 0.17 ± 0.02 dex higher than the original SFRs estimated by SED fitting. However, estimating sizes and estimating SFR are two different things that do not impact each other, because the SFR differences for IRAGN (0.20 ± 0.06 dex), XAGN (0.16 ± 0.04 dex), and the SF (0.19 ± 0.03 dex) galaxies are very close. We also visually inspected the SED fitting results, and found that the fit quality was adequate (i.e., no major failures). As a result, we adopt the estimation above in this paper.

2.4. Size Measurements

The effective radius in F814W-band is about 3 kpc (∼0farcs4 at z = 0.75) and the FWHM of the HST/F814W imaging is 0farcs09, so high-resolution ALMA observations (<0farcs4) are required to resolve the dust distribution at a complementary angular resolution of the HST imaging.

First, we measured the rest-frame FIR sizes with the IMFIT task in the image plane, which fitted two-dimensional elliptical Gaussian components, in CASA (version 5.1.0). The effective radius by IMFIT is calculated by ${R}_{e}=\sqrt{{R}_{\mathrm{Maj}}\times {R}_{\mathrm{Min}}}/2$, where RMaj and RMin are the deconvolved FWHM major and minor axes. Then we used the UVMODELFIT task in the uv plane, which fits a single component source model to the UV data. We used the flux density and the size measurements obtained by the IMFIT task as the initial values, and fitted Gaussian models with five iterations by UVMODELFIT task. The effective radius by UVMODELFIT is also calculated by the fitted major and minor axes. We compared the measurements by IMFIT for the image plane and UVMODELFIT for the uv plane. The difference of the sizes are with a 2% systematic difference (UVMODELFIT gives lower values). To avoid the question about the fidelity of the ALMA imaging, we adopted the size measurements from UVMODELFIT for single objects in this paper. We also measured the rest-frame FIR sizes with GALFIT (Peng et al. 2010), assuming galaxy models with single Sérsic profiles. The sizes measured from GALFIT are slightly larger (systematic difference ∼10%) than the sizes measured from IMFIT. These comparisons show that our size measurements are consistent with each other.

There are 23 out of 27 sources with available HST/F814W images. We measured the F814W-band light distribution of our COSMOS galaxy sample with GALFIT, assuming galaxy models with the single Sérsic profiles. Among the 27 sources with reliable redshifts, 2 of them can be found in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS, Grogin et al. 2011; Koekemoer et al. 2011) region. Therefore, we used the single Sérsic profile measurement of available F125W and F160W imaging from van der Wel et al. (2012) in the CANDELS/COSMOS (Grogin et al. 2011; Koekemoer et al. 2011) field. The effective radius from single Sérsic profile represents the rest-frame optical size.

3. Analysis

We separated the sample between non-AGN SFGs, X-ray selected AGNs (XAGNs: LX(2–10 keV) > 1042 ergs/s), and infrared-selected AGNs (IRAGNs: using the mid-infrared colors from the IRAC 3.6, 5.8, 4.5, and 8.0 μm; obscured IRAGNS: ${L}_{\mathrm{IR},\mathrm{AGN}}/{L}_{{\rm{X}},\mathrm{AGN}}\gt 20$), that are identified by previous work (Civano et al. 2016; Marchesi et al. 2016; Chang et al. 2017a, 2017b). Figure 3 shows the SFR versus stellar mass diagrams in the two redshift bins. Most of the sources in the sample from the ALMA archive are massive and have high SFRs, and they are in general on and above the star-forming sequence. Our star-forming sequences are estimated by the SED fitting results for $\mathrm{NUV}\mbox{--}r$ versus rJ selected galaxies. We also show the main-sequence relations from Speagle et al. (2014) and Whitaker et al. (2014). Our SFR estimations are generally lower than those in the literature by ∼0.1–0.3 dex at z ∼ 2. The main reason might be selection of the star-forming sample. With the same sample selection, our main sequences are very close (<0.1 dex) to the main sequences calculated by the stellar masses and SFRs in the COSMOS2015 catalogs. To avoid bias from different selection criteria, we use the $\mathrm{NUV}-r$ versus r − J selection, but show the literature in Figure 3 for comparison. The SFR differences between the various works are more significant at the high stellar mass end, so it is important to compare AGN hosts with non-AGN galaxies using SFRs derived with the same method. Comparing to the archival data at z ∼ 1, our AGN sample have low SFRs because they were selected to be within ±0.3 dex of the main sequence. This may be the reason why our sample was mostly undetected with ALMA. The only one IRAGN at z ∼ 1 has lower stellar mass (∼0.14 dex) and higher SFR (∼0.20 dex) than non-AGN galaxies. There is only one IRAGNs at z ∼ 2, with lower stellar masses (∼0.04 dex) and lower SFRs (∼0.35 dex) than the median of the non-AGN galaxies.

Figure 3.

Figure 3. SFR vs. stellar mass plots at z ∼ 1 and z ∼ 2. The black solid and dashed lines are the star-forming sequence defined by our SED fitting results. The detected obscured IRAGN (pink square), undetected obscured IRAGN (purple square), and stacked values of the undetected obscured IRAGNs (star) in our ALMA sample are shown. The blue dots are star-forming galaxies (SFGs), green dots are X-ray selected AGN hosts (XAGN), and red dots are infrared-selected AGN hosts (IRAGN) from the ALMA archive. Most of our obscured IRAGNs have lower stellar masses and SFRs. The main sequence derived by Speagle et al. (2014) (yellow dotted lines) and Whitaker et al. (2014; yellow dash lines) are also shown.

Standard image High-resolution image

Figure 4 shows the FIR size as a function of redshift (left), stellar mass (middle), and the distance to the star-forming sequence (right) as shown in Figure 3. The sizes of our detected (0.43 ± 0.28 kpc) and stacked (0.90 ± 0.21 kpc) obscured AGN hosts are very compact compared to those from the archive data (mean = 1.76 ± 0.25 kpc and median = 1.67 ± 0.30 kpc at z ∼ 1 and $\mathrm{log}\,{M}_{* }\sim 11$). The uncertainties here and below are estimated from a bootstrapping analysis. As shown in Table 1, for $\mathrm{log}\,{M}_{* }\sim 11$ sources at z ∼ 1, the median (mean) FIR sizes of the SFGs, XAGNs, IRAGNs, are 1.74 ± 0.50 (1.62 ± 0.50) 2.46 (one object),15 0.77 ± 0.50 (0.76 ± 0.48), respectively. As shown in Table 1, for $\mathrm{log}\,{M}_{* }\sim 11$ sources at z ∼ 2, the median (mean) FIR sizes of the SFGs, XAGNs, and IRAGNs, are 1.25 ± 0.73 (1.30 ± 0.39), 0.69 ± 0.25 (0.70 ± 0.24), 0.70 (one object), respectively. We find systematically smaller galaxy sizes at z ∼ 2 relative to the corresponding samples at z ∼ 1. The measurements of HST/F814W images are rest-frame optical sizes for z ∼ 1 sources but rest-frame UV sizes for z ∼ 2 sources (e.g., van der Wel et al. 2014), and color gradient which shows redder centers should also be taken into account (e.g., Wuyts et al. 2012). Most importantly, the number of objects in each sample is limited, especially for the AGN sample. Our IRAGN sample was selected to be within ±0.3 dex of the main sequence, so we also consider $\mathrm{log}\,{M}_{* }\sim 11$ SFGs on the main sequence. For $\mathrm{log}\,{M}_{* }\sim 11$ SFGs on the main sequence (the distance to the star-forming sequence <0.5 dex; a less strict range because of the limited sample size) at z ∼ 1, the median (mean) sizes are 1.97 ± 0.44 kpc (1.83 ± 0.49 kpc). For $\mathrm{log}\,{M}_{* }\sim 11$ SFGs on the main sequence (the distance to the star-forming sequence <0.5 dex) at z ∼ 2, the median (mean) sizes are 1.26 ± 1.43 kpc (1.27 ± 1.70 kpc). As a result, we find that IRAGN hosts are in general smaller than SFGs on the main sequence at a given stellar mass and redshift bin. However, it is difficult to find a reasonable parent sample for our sources on the main sequence because all SFGs are above the main sequence and mainly in the starburst regime. Therefore, future ALMA observations with better matched control samples are necessary to expand the parameter spaces.

Figure 4.

Figure 4. Rest-frame FIR size (effective radius) vs. redshift, stellar mass, and distance to the star-forming sequence according to our SED fitting results. The symbols are the same as those described in Figure 3.

Standard image High-resolution image

Table 1.  The Median (mean) FIR Sizes for $\mathrm{log}\,{M}_{* }\sim 11$ Sources in Kiloparsecs

Type z ∼ 1 z ∼ 2
SFG 1.74 ± 0.50 (1.62 ± 0.50) 1.25 ± 0.73 (1.30 ± 0.39)
XAGN 2.46 (one object) 0.69 ± 0.25 (0.70 ± 0.24)
IRAGN 0.77 ± 0.50 (0.76 ± 0.48) 0.70 (one object)

Download table as:  ASCIITypeset image

Figure 5 shows the histograms of HST/F814W, and ALMA/1 mm sizes. As shown in Table 2, the mean (median) optical sizes of the SFGs, XAGNs, and IRAGNs, are 1.95 ± 0.74 (2.16 ± 0.74), 1.58 ± 0.50 (1.56 ± 0.50), and 0.42 ± 0.42 (0.41 ± 0.42), respectively. And the mean (median) FIR sizes of the SFGs, XAGNs, and IRAGNs, are 1.58 ± 0.55 (1.48 ± 0.39), 1.07 ± 0.52 (0.89 ± 0.56), 0.75 ± 0.39 (0.68 ± 0.45), respectively. Considering the uncertainties of the sizes and the different measurement techniques as described in Section 2, these numbers imply that the rest-frame FIR sizes of all galaxies seem to be similar or more compact than their rest-frame optical sizes. In particular, the FIR size (0.90 ± 0.21 kpc) of our stacked sources is smaller than the optical one (1.58 ± 0.01). Though there are few IRAGNs, most of them show both very compact optical and FIR sizes. The rest-frame optical and FIR sizes of XAGNs are also relatively small or similar as shown in Figure 5.

Figure 5.

Figure 5. Upper left: histograms of ALMA/1 mm size. Lower right: histogram of HST/F814W size. The magenta lines show the effective radii of the stacked images from the nine undetected obscured AGN sources. Here the detected obscured AGN is included in the IRAGN sample. Lower left: HST/F814W size to ALMA/1 mm size plot. The blue color represents the density of SFGs in each size bin. This subplot only includes objects that have measurements of both HST/F814W and ALMA/1 mm sizes.

Standard image High-resolution image

Table 2.  The Median (mean) Optical (HST/F814W) and FIR (ALMA/1 mm) Sizes in Kiloparsecs

Type HST/F814W Size ALMA/1 mm Size
SFG 1.95 ± 0.74 (2.16 ± 0.74) 1.58 ± 0.55 (1.48 ± 0.39)
XAGN 1.58 ± 0.50 (1.56 ± 0.50) 1.07 ± 0.52 (0.89 ± 0.56)
IRAGN 0.42 ± 0.42 (0.41 ± 0.42) 0.75 ± 0.39 (0.68 ± 0.45)

Download table as:  ASCIITypeset image

There are only two objects in the CANDELS field as shown in Table 3, and one XAGN is at z ∼ 2. The HST/F160W measurements are better tracers for optical sizes than the HST/F814W ones at z ∼ 2. The current data shows that the two sizes in HST/F160W are comparable but larger than most sizes in HST/F814W. HST/F160W traces rest-frame the optical size and HST/F814W traces the UV size for the same object, so it is difficult to make a direct comparison. We need size measurements using observations at the same (rest-frame) wavelength for larger numbers of objects to draw a statistical conclusion at z ∼ 2.

Table 3.  Properties of 27 Detected (Upper) and 9 Undetected (Lower) ALMA Band 6/7 High-resolution (0farcs1–0farcs4) Sources

Project Code Source Name Band Beam Size ${\mathrm{SN}}_{\mathrm{peak}}$ ${{\rm{S}}}_{\mathrm{peak}}$ ${{\rm{S}}}_{\mathrm{int}}$ ${\mathrm{ID}}_{\mathrm{COSMOS}}$ R.A. Decl. z ${{\rm{M}}}_{* }$ SFR ${{\rm{R}}}_{\mathrm{ACS}({\rm{g}})}$ ${{\rm{R}}}_{\mathrm{WFC}3({\rm{g}})}$ ${{\rm{R}}}_{\mathrm{ALMA}({\rm{g}})}$ ${{\rm{R}}}_{\mathrm{ACS}({\rm{c}})}$ ${{\rm{R}}}_{\mathrm{ALMA}({\rm{c}})}$ ${{\rm{R}}}_{\mathrm{ALMA}({\rm{d}})}$ ${{\rm{R}}}_{\mathrm{ALMA}(\mathrm{UV})}$ AGN
      $(^{\prime\prime} \times ^{\prime\prime} )$   (mJy/beam) (mJy) (deg) (deg)   (M) (M yr−1) (kpc) (kpc) (kpc) (kpc) (kpc) (kpc) (kpc)  
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21)
2013.1.00034.S lowz_cell8_6496 7 0.33 × 0.24 3.54 0.39 ± 0.11 2.17 ± 0.72 239908 149.97412 1.64352 1.04 11.2 ± 0.1 2.2 ± 0.0 1.4 ± 0.2 0.6 ± 0.0 2.7 ± 0.5 2.4 ± 0.8 2.5 ± 0.8 XAGN
2016.1.01604.S 317108 7 0.39 × 0.39 14.61 0.71 ± 0.05 1.64 ± 0.15 317108 150.15256 1.76921 1.56 11.2 ± 0.1 2.3 ± 0.1 1.9 ± 0.1 3.6 ± 0.1 2.5 ± 0.1 1.8 ± 0.2 1.4 ± 0.2
2016.1.01604.S 360903 7 0.34 × 0.29 5.67 0.30 ± 0.05 0.85 ± 0.20 360903 149.66491 1.83649 1.35 11.5 ± 0.1 1.9 ± 0.1 1.2 ± 0.6 3.6 ± 0.0 2.2 ± 0.3 1.7 ± 0.4 2.7 ± 0.4
2016.1.01604.S 405986 7 0.33 × 0.29 8.04 0.36 ± 0.04 2.03 ± 0.29 405986 149.79795 1.90880 1.35 11.3 ± 0.1 2.0 ± 0.1 0.6 ± 0.5 3.0 ± 0.1 2.4 ± 0.0 3.2 ± 0.3 2.8 ± 0.3 3.2 ± 0.3
2016.1.01604.S 495291 7 0.39 × 0.39 26.43 1.07 ± 0.04 1.75 ± 0.10 495291 150.62793 2.05289 1.30 11.1 ± 0.1 2.2 ± 0.1 3.4 ± 0.7 3.1 ± 0.1 2.1 ± 0.1 1.3 ± 0.1 1.1 ± 0.1
2016.1.01604.S 500431 7 0.39 × 0.39 13.05 0.53 ± 0.04 1.43 ± 0.15 500431 150.13314 2.06073 1.39 11.0 ± 0.1 1.9 ± 0.1 2.5 ± 0.9 3.5 ± 0.1 2.7 ± 0.2 1.8 ± 0.3 1.8 ± 0.3
2013.1.00034.S lowz_cell11_120263 7 0.40 × 0.24 9.55 1.73 ± 0.18 2.92 ± 0.45 504172 149.83252 2.06602 1.14 11.4 ± 0.0 2.4 ± 0.0 2.2 ± 0.6 2.8 ± 0.0 1.6 ± 0.1 0.8 ± 0.3 0.7 ± 0.3
2016.1.00735.S CXID-36 7 0.11 × 0.10 15.19 0.66 ± 0.04 2.69 ± 0.22 553928 150.15841 2.13958 1.83 11.8 ± 0.7 3.0 ± 0.0 0.0 ± 8.3 0.4 ± 0.0 0.9 ± 0.0 0.8 ± 0.1 0.9 ± 0.1 XAGN
2016.1.01426.S PACS_819 6 0.13 × 0.09 13.47 0.14 ± 0.01 0.56 ± 0.05 628146 149.98145 2.25320 1.40 10.8 ± 0.0 2.6 ± 0.0 1.0 ± 0.4 7.0 ± 0.1 0.9 ± 0.0 0.8 ± 0.1 0.4 ± 0.1
2016.1.01604.S 659927 7 0.34 × 0.29 8.05 0.41 ± 0.05 1.17 ± 0.19 659927 150.34101 2.30099 1.44 11.4 ± 0.1 2.0 ± 0.1 2.6 ± 0.6 3.8 ± 0.0 2.3 ± 0.2 1.6 ± 0.3 2.0 ± 0.3
2016.1.01604.S 766157 7 0.33 × 0.29 13.39 0.84 ± 0.06 1.51 ± 0.16 766157 150.00352 2.46146 1.40 11.3 ± 0.1 2.2 ± 0.0 1.1 ± 101.3 2.5 ± 0.1 1.8 ± 0.1 1.2 ± 0.2 1.2 ± 0.2
2013.1.00034.S midz_cell8_254150 7 0.33 × 0.30 15.95 2.37 ± 0.15 3.09 ± 0.31 794848 150.09342 2.50734 2.20 10.8 ± 0.1 0.8 ± 0.2 1.9 ± 0.7 3.1 ± 0.1 2.0 ± 0.1 1.5 ± 0.1 0.7 ± 0.3 0.6 ± 0.3 XAGN
2016.1.00778.S J100054.83 + 023126.2 6 0.45 × 0.32 9.68 0.33 ± 0.03 0.48 ± 0.08 806579 150.23118 2.52466 1.95 10.9 ± 0.1 1.6 ± 0.2 22.3 ± 22.4 1.9 ± 0.1 0.6 ± 0.7 0.6 ± 3.3
2013.1.00034.S lowz_cell3_258881 7 0.52 × 0.36 7.44 0.97 ± 0.13 1.64 ± 0.33 811432 150.04247 2.52656 1.23 11.2 ± 0.1 2.1 ± 0.0 1.8 ± 1.1 2.3 ± 1.0 4.5 ± 0.1 2.4 ± 0.2 1.2 ± 0.7 1.5 ± 0.5
2013.1.00034.S midz_cell11_262500 7 0.37 × 0.34 26.26 3.72 ± 0.14 6.70 ± 0.37 815012 150.60329 2.53654 2.19 11.3 ± 0.1 3.2 ± 0.0 1.1 ± 2.0 0.8 ± 0.0 2.0 ± 0.1 1.3 ± 0.1 1.2 ± 0.1
2016.1.01604.S 816106 7 0.39 × 0.39 19.19 0.86 ± 0.05 1.81 ± 0.13 816106 150.08329 2.53609 1.36 11.3 ± 0.1 2.3 ± 0.1 2.3 ± 0.7 3.5 ± 0.0 3.2 ± 0.0 2.4 ± 0.1 1.7 ± 0.1 1.5 ± 0.1
2013.1.00034.S lowz_cell11_268592 7 0.53 × 0.36 3.45 0.36 ± 0.10 2.24 ± 0.75 831023 150.41785 2.55852 1.21 11.0 ± 0.1 2.4 ± 0.0 2.5 ± 0.3 2.7 ± 0.0 4.5 ± 0.9 4.1 ± 1.2 3.4 ± 1.4
2015.1.00026.S SHIZELS-14 6 0.25 × 0.21 12.96 0.24 ± 0.02 2.33 ± 0.20 831167 150.21496 2.55952 2.08 11.2 ± 0.1 2.7 ± 0.0 2.0 ± 0.8 2.3 ± 0.0 3.0 ± 0.2 2.7 ± 0.2 1.6 ± 0.2
2016.1.00478.S AzTECC2a 7 0.20 × 0.17 7.72 1.26 ± 0.16 1.78 ± 0.36 842595 149.99797 2.57823 2.42 10.8 ± 0.0 3.0 ± 0.0 0.7 ± 175.9 3.6 ± 0.1 0.9 ± 0.1 0.5 ± 0.2 0.7 ± 0.2 XAGN
2013.1.00034.S lowz_cell11_286398 7 0.54 × 0.37 6.84 0.78 ± 0.11 4.03 ± 0.70 872762 150.12960 2.62144 1.41 11.4 ± 0.1 2.3 ± 0.1 8.2 ± 0.2 4.3 ± 0.4 3.6 ± 0.6 2.4 ± 0.6
2016.1.01604.S 888769 7 0.39 × 0.39 23.00 1.02 ± 0.04 1.62 ± 0.11 888769 150.10745 2.64591 1.54 10.9 ± 0.2 2.2 ± 0.1 2.0 ± 1.1 1.5 ± 0.0 2.1 ± 0.1 1.3 ± 0.1 1.2 ± 0.1
2016.1.00478.S AzTECC22a 7 0.19 × 0.18 19.99 3.56 ± 0.18 6.41 ± 0.47 902320 150.03729 2.66960 1.85 11.0 ± 0.1 3.0 ± 0.0 2.3 ± 0.6 1.9 ± 0.0 1.0 ± 0.0 0.7 ± 0.1 0.7 ± 0.1
2016.1.01604.S 917430 7 0.34 × 0.29 13.33 0.70 ± 0.05 0.86 ± 0.11 917430 149.91774 2.69211 1.60 11.1 ± 0.1 2.3 ± 0.1 0.8 ± 0.3 0.6 ± 0.0 1.5 ± 0.1 0.6 ± 0.2 0.7 ± 0.2 IRAGN
2013.1.00034.S lowz_cell10_324567 7 0.56 × 0.38 5.53 0.73 ± 0.13 1.45 ± 0.38 960580 149.91909 2.76053 1.13 10.7 ± 0.1 2.1 ± 0.0 2.3 ± 0.3 3.0 ± 0.0 2.7 ± 0.4 1.9 ± 0.7 1.8 ± 0.8
2013.1.00034.S lowz_cell10_328266 7 0.57 × 0.38 5.81 0.80 ± 0.14 1.61 ± 0.40 969105 150.53066 2.77548 1.35 11.2 ± 0.0 2.4 ± 0.0 0.1 ± 0.3 0.6 ± 0.0 2.8 ± 0.3 2.0 ± 0.6 1.2 ± 0.6 IRAGN
2016.1.01355.S 972308 7 0.12 × 0.10 5.79 0.52 ± 0.09 1.37 ± 0.31 972308 150.40459 2.78053 0.60 10.6 ± 0.1 1.4 ± 0.0 0.0 ± 0.0 0.4 ± 0.0 0.6 ± 0.1 0.5 ± 0.1 0.5 ± 0.1 IRAGN
2013.1.00034.S midz_cell6_332275 7 0.37 × 0.33 19.46 3.00 ± 0.15 7.19 ± 0.50 980250 150.01611 2.79235 1.76 11.3 ± 0.1 2.6 ± 0.0 0.7 ± 0.0 2.3 ± 0.1 1.7 ± 0.1 1.6 ± 0.1
2016.1.01355.S 263569 7 0.12 × 0.10 263569 149.95048 1.68282 0.78 10.8 ± 0.1 1.3 ± 0.1 2.4 ± 0.4 IRAGN
2016.1.01355.S 334665 7 0.12 × 0.10 334665 150.43685 1.79439 0.56 10.9 ± 0.1 1.2 ± 0.1 1.8 ± 0.1 IRAGN
2016.1.01355.S 335028 7 0.12 × 0.10 335028 149.50325 1.79492 0.65 11.3 ± 0.1 1.8 ± 0.1 1.5 ± 0.1 IRAGN
2016.1.01355.S 410490 7 0.12 × 0.10 410490 150.02499 1.91477 0.97 11.2 ± 0.1 1.9 ± 0.0 IRAGN
2016.1.01355.S 414248 7 0.12 × 0.10 414248 149.84769 1.92058 0.71 10.7 ± 0.1 1.2 ± 0.2 0.0 ± 0.0 IRAGN
2016.1.01355.S 417645 7 0.12 × 0.10 417645 150.02586 1.92639 0.69 10.7 ± 0.1 1.4 ± 0.1 0.7 ± 0.1 IRAGN
2016.1.01355.S 456808 7 0.12 × 0.10 456808 149.57039 1.99058 1.00 11.2 ± 0.1 1.1 ± 0.2 1.1 ± 0.2 IRAGN
2016.1.01355.S 711692 7 0.12 × 0.10 711692 149.99343 2.37723 0.93 10.7 ± 0.1 1.6 ± 0.1 0.9 ± 0.3 IRAGN
2016.1.01355.S 755017 7 0.12 × 0.10 755017 150.76449 2.44222 0.70 10.9 ± 0.1 1.8 ± 0.1 3.4 ± 0.2 IRAGN

Note. (1) ALMA project code. (2) ALMA source name. (3) ALMA bands. (4) ALMA beam size. (5) ALMA S/N. (6) ALMA peak flux from IMFIT. (7) Integrated total flux from IMFIT. (8) COSMOS ID. (9) COSMOS R.A. (10) COSMOS decl. (11) Best redshift in COSMOS. (12) Stellar mass estimation from MAGPHYS. The uncertainties of stellar masses are from the probability density functions, which depend on the model templates. The real uncertainty can be larger than our estimation. (13) SFR estimation from MAGPHYS. (14) Sizes measured from HST/F814W images with GALFIT. (15) Sizes measured from HST/F160W images with GALFIT (16) Sizes measured from HST/F160W images with GALFIT. (17) Sizes measured from HST/F814W images with IMFIT (convolved). (18) Sizes measured from ALMA/ 1mm images with IMFIT (convolved). (19) Sizes measured from ALMA/ 1mm images with IMFIT (deconvolved). (20) Sizes measured from ALMA/ 1mm images with UVMODELFIT. (21) IRAGN, XAGN, or non-AGN galaxies.

Download table as:  ASCIITypeset image

Figure 6 shows the ALMA/1 mm sizes to the infrared luminosity relation. Our galaxies fall more or less on top of the Fujimoto et al. (2017) line, but it is difficult to show a clear relation because of the scatter and sample size. Here we decompose the AGN and star formation components according to the SEDs by MAGPHYS+AGN as described in Chang et al. (2017a). Therefore, we are able to derive the infrared luminosity from the star formation and AGN components, respectively. Because FIR luminosity traces dusts, this decomposition can separate dust emission of AGNs from their host galaxies. Figure 6 does not show an obvious correlation between the FIR size and the IR luminosity from the AGN component. According to our SED decomposition, the AGN contribution in the FIR range is negligible (<1% of the total flux), so the ALMA/ 1mm sizes are dominated by star-forming populations. Therefore, the middle panel in Figure 6 is approximate to the FIR size of hosts versus the infrared luminosity from the star-forming part. It is possible that there are correlations between the FIR sizes of the AGNs and the infrared luminosity from the AGN part. However, deeper and higher resolution data are required to resolve the AGN component at 1mm with ALMA.

Figure 6.

Figure 6. FIR size vs. total infrared luminosity, infrared luminosity from star formation, and AGN contributions. The black lines are from Fujimoto et al. (2017). The symbols are the same as those described in Figure 3.

Standard image High-resolution image

Figure 7 shows the specific SFR (sSFR = SFR/M*) to SFR surface density (${{\rm{\Sigma }}}_{\mathrm{SFR}}=\mathrm{SFR}/(\pi {R}_{e}^{2})$, where Re is the rest-frame FIR size) plots. Because there are more massive and star-forming sources at high redshifts, it is clear to see that high redshift sources having higher SFR will naturally imply higher SFR surface density. It was also mentioned earlier that the sizes were systematically smaller at higher redshift, which would also increase the surface densities even at a fixed SFR or a fixed sSFR. At z ∼ 1, the SFR surface density of our obscured IRAGN hosts are high compared to archived ALMA observations at a given sSFR. At z ∼ 2, the SFR surface density of both XAGN and IRAGN hosts seem to be similar or slightly higher than that of SFGs, but the uncertainties of the SFR surface density are large. In Figure 7, SFR surface densities of our IRAGN hosts are much higher than SFGs on the main sequence at both z ∼ 1. This suggests that the compactness of IRAGN hosts will be significant if we have a large control sample of SFGs. Therefore, more data with high resolution are needed to reach a statistically meaningful conclusion.

Figure 7.

Figure 7. Specific star formation rate to SFR surface density plot in the two redshift bins. Our obscured IRAGN hosts are very compact compared to archival ALMA observations, especially for those on the star-forming sequence. The symbols are the same as those described in Figure 3.

Standard image High-resolution image

4. Discussion

4.1. Compact Rest-frame FIR and Optical sizes

We used high-resolution ALMA images to measure the rest-frame FIR sizes, and compared the dust maps with the HST optical images for non-AGN SFGs, X-ray selected AGNs, and IR-selected AGNs. We found that the rest-frame FIR sizes of IRAGN hosts are as compact (<1 kpc) as their rest-frame optical sizes, which suggests that they are indeed compact SFGs which have nuclear activities (also see Lutz et al. 2018) due to compaction (e.g., Dubois et al. 2015; Zolotov et al. 2015; Tacchella et al. 2016). It is possible that the detected sources may be more concentrated galaxies than undetected ones because we may not detect galaxies that are too extended in the high-resolution ALMA observations given the sensitivity limit. However, we apply the analyses equally to AGN and non-AGN galaxies, so the conclusion would not be changed. A possible explanation might be that obscured AGNs are triggered in compact galaxies. Such compact profiles could also be explained if obscured AGNs are caught at a special moment when their host galaxies still harbor intense star formation but simultaneously experience a phase of compaction, as what was recently suggested from simulations of gas-rich galaxies undergoing violent disk instabilities. The possible migration of their giant star-forming clumps toward their central region could lead to enhanced AGN fueling (Bournaud et al. 2012; Dubois et al. 2015; Tacchella et al. 2016) as well as to the quick formation of a dense stellar component (Dekel et al. 2009; Ceverino et al. 2010). According to Tacchella et al. (2016), post-compaction blue nuggets are expected to be found above the main-sequence ridge. Most of the archival sources seem to be biased on the higher side of the main sequence and above, so there is an important selection effect here when comparing the samples. Our sources were selected to lie on the MS relation within ±0.3 dex (Figure 3), so by selection, they do not have the same bias toward higher SFR. That might also explain why they do not lie in the expected blue nugget location. The lack of a comparison sample at the same stellar mass, SFR, and redshift bin might also be the reason why we did not see much higher SFR surface density of our obscured AGNs in Figure 7. Our results may imply that the dusty star formation of AGN hosts happens in compact regions, and forms the stellar component of the host galaxy. However, more FIR data with a better matched control sample are needed to confirm the simulation predictions.

4.2. FIR versus Optical Sizes in AGN Hosts

It has been found that the rest-frame FIR sizes are in general smaller than the rest-frame optical sizes (e.g., Barro et al. 2016; Fujimoto et al. 2017, 2018; Tadaki et al. 2017; Elbaz et al. 2018), which can be explained by compact dusty starbursts. Though the size estimation can be sensitive to the measurement methods, model profiles, and initial values (e.g., Elbaz et al. 2018; Fujimoto et al. 2018), the main results of our analysis do not change even if we use different ways to measure the sizes. We cross-matched our sample with Fujimoto et al. (2018), who showed consistent size measurement with previous ALMA results (Ikarashi et al. 2015; Simpson et al. 2015b; Barro et al. 2016; Hodge et al. 2016; Rujopakarn et al. 2016; González-López et al. 2017; Tadaki et al. 2017). We found 20 matches out of our 27 sources. Among them, there is no significant difference (systematic difference ∼0.05 ± 0.10 kpc; our value is smaller but close to Fujimoto et al. 2018) in between our size estimation with UVMODELFIT and their FWHM estimated with the UVMULTIFIT task of CASA. Moreover, Elbaz et al. (2018) show that the semimajor axis from uv (UVMODELFIT in CASA) and image (GALFIT) planes agree with only a 5% systematic difference. Therefore, our size measurements should be robust against those of other size measurement methods. A more sophisticated analysis might be performed with GALFIT (e.g., bulge-disk decomposition for the PHIBSS2 galaxies; Tacconi et al. 2018; Freundlich et al. 2019), but more HST/F160W images or higher resolution ALMA images (∼0farcs1) of the complementary sample at z < 1 are required. Our results show that the rest-frame FIR sizes of both AGN and SFG seem to be similar or more compact than their rest-frame optical sizes. This may suggest that the star formation mechanisms in the compact regions of AGN hosts are similar to those of other ALMA-detected sources. However, more detected AGN hosts with ALMA observations are needed to confirm this result.

4.3. SFR Surface Densities of Obscured IRAGNs

Our galaxies fall more or less on top of the Fujimoto et al. (2017) line (also see Harrison et al. (2016) for an X-ray selected AGN sample), but it is difficult to show a clear relation because of the scatter and sample size. Our measurements are also above the maximum SFR surface density ($\sim 1000\,{M}_{\odot }\,{\mathrm{yr}}^{-1}\,{\mathrm{kpc}}^{-2}$ as shown in Figure 7) which was estimated by Simpson et al. (2015b). The decomposed IRAGN luminosity may imply that AGN activity has little influence on the FIR sizes. However, more high-resolution imaging from future ALMA observations are required to show a statistical result. On the other hand, the median SFR surface density of the obscured AGN hosts is more than two times higher than those of other ALMA-detected sources at $\mathrm{log}\,{M}_{* }/{M}_{\odot }\sim 11$ and z ∼ 1. This implies that obscured AGN host galaxies may be in a transition phase of compact SFGs, and may evolve to compact quiescent galaxies at later time (Dubois et al. 2015; Zolotov et al. 2015). The obscured black hole accretion phase may be related to galaxies experiencing a compaction of their gaseous component, which could play a key role in the final quenching mechanisms driving the formation of passive red sources. Nevertheless, there are still very few complementary samples with matched stellar mass, SFR, and redshift. Most of the archival ALMA observations focus on massive galaxies above the star-forming sequence. Therefore, future high-resolution ALMA observations would be crucial to investigating further about these compact AGN host galaxies.

5. Summary

In this paper, we have investigated rest-frame FIR and optical sizes of non-AGN SFGs, infrared-selected AGNs, and X-ray selected AGNs at z < 2.5 by using high-resolution ALMA (0farcs10–0farcs85) and HST imaging (∼0farcs1) in the COSMOS field. Our main findings are as follows:

  • 1.  
    The sizes measured from both ALMA/1 mm and HST/F814W images show that obscured AGN host galaxies at z ∼ 1 are more compact than SFGs at a given redshift and stellar mass.
  • 2.  
    Most detected rest-frame FIR sizes of the sources are similar or more compact than their rest-frame optical sizes, which is consistent with recent results for ALMA-detected sources. This can be explained by the fact that the dusty starbursts take place in the compact regions, and suggests that the star formation mechanisms in the compact regions of AGN hosts are similar to those taking place in SFGs with no evidence of an AGN at their center.
  • 3.  
    The decomposed IRAGN luminosity may imply that AGN activity has little influence on the FIR sizes. However, more high-resolution imaging from future ALMA observations are needed to confirm this result.
  • 4.  
    The median SFR surface density of the obscured AGN hosts is higher than those of other ALMA-detected sources at $\mathrm{log}\,{M}_{* }/{M}_{\odot }\sim 11$ and z ∼ 1. This implies that obscured AGN host galaxies may be in a transition phase of compact SFGs, and may evolve into compact quiescent galaxies at later time.
  • 5.  
    Future high-resolution ALMA observations would be crucial to investigating these compact AGN host galaxies further.

We thank the referee, as well as J. Freundlich for helpful comments and discussions. Y.Y.C. acknowledges financial support from the Ministry of Science and Technology of Taiwan grant (105-2112-M-001-029-MY3 and 108-2112-M-001-014-), and the Agence Nationale de la Recherche (contract #ANR-12-JS05-0008-01). E.d.C. gratefully acknowledges the Australian Research Council as the recipient of a Future Fellowship (project FT150100079). This paper makes use of the following ALMA data: ADS/JAO.ALMA #2013.1.00034.S, #2015.1.00026.S, #2016.1.00478.S, #2016.1.00735.S, #2016.1.00778.S, #2016.1.01355.S, #2016.1.01426.S, and #2016.1.01604.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), MOST 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. We gratefully acknowledge the contributions of the entire COSMOS collaboration.

Footnotes

  • 12 

    The synthesized beam sizes are slightly different if we use different weighting. A synthesized beam size of 0farcs12 × 0farcs10 corresponds to Briggs weighting with 0.5 robustness parameter in CASA.

  • 13 

    In order to have a high signal-to-noise ratio, we used natural weighting with TCLEAN in CASA. As a result, the synthesized beam size becomes 0farcs14 × 0farcs12.

  • 14 

    For a choice of 1farcs0 (0farcs8) simulated images, the differences of sizes are 40% (30%). Moreover, >90% of sources in the range of 0farcs8–1farcs0 fail to be well-fitted in the image plane, and >50% of sources in the range of 0farcs4–0farcs8 fail to be well-fitted in the uv plane. Considering the simulation and uncertainties of the fitting results, we decided to have a range of 0farcs1–0farcs4 as our complementary sample.

  • 15 

    For only one object in the stellar mass and redshift bin, the uncertainties are not available.

Please wait… references are loading.
10.3847/1538-4357/ab595b