Search

We present confusion-limited SCUBA-2 450-$\mu$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $\mu$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $\sigma_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $\sigma_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The corresponding confusion noise levels are estimated to be 0.65 and 0.36 mJy beam$^{-1}$. Above the 4 (3.5) $\sigma$ threshold, we detected 360 (479) sources at 450 $\mu$m and 237 (314) sources at 850 $\mu$m. We derive the deepest blank-field number counts at 450 $\mu$m, covering the flux-density range of 2 to 43 mJy. These are in agreement with other SCUBA-2 blank-field and lensing-cluster observations, but are lower than various model counts. We compare the counts with those in other fields and find that the field-to-field variance observed at 450 $\mu$m at the $R=6^\prime$ scale is consistent with Poisson noise, so there is no evidence of strong 2-D clustering at this scale. Additionally, we derive the integrated surface brightness at 450 $\mu$m down to 2.1 mJy to be $57.3^{+1.0}_{-6.2}$~Jy deg$^{-2}$, contributing to (41$\pm$4)\% of the 450-$\mu$m extragalactic background light (EBL) measured by COBE and Planck. Our results suggest that the 450-$\mu$m EBL may be fully resolved at $0.08^{+0.09}_{-0.08}$~mJy, which extremely deep lensing-cluster observations and next-generation submillimeter instruments with large aperture sizes may be able to achieve., Comment: 29 pages, 14 figures, accepted for publication in ApJ

We present ALMA 870um and JCMT SCUBA2 850um dust continuum observations of a sample of optically dark and strongly lensed galaxies in the cluster fields. The ALMA and SCUBA2 observations reach a median rms of about 0.11 mJy and 0.44 mJy, respectively, with the latter close to the confusion limit of the data at 850um. This represents one of the most sensitive searches for dust emission in optically dark galaxies. We detect the dust emission in 12 out of 15 galaxies at >3.8 sigma, corresponding to a detection rate of 80 per cent. Thanks to the gravitational lensing, our observations reach a deeper limiting flux than previous surveys in blank fields by a factor of 3. We estimate delensed infrared luminosities in the range log(LIR)=11.5-12.7 Lsun, which correspond to dust-obscured star formation rates (SFRs) of 30 to 520 Msun per year. Stellar population fits to the optical-to-NIR photometric data yield a median redshift z=4.26 and de-lensed stellar mass log(Mstar)=10.78 Msun. They contribute a lensing-corrected star-formation rate density at least an order of magnitude higher than that of equivalently massive UV-selected galaxies at z>3. The results suggest that there is a missing population of massive star-forming galaxies in the early Universe, which may dominate the SFR density at the massive end. Five optically dark galaxies are located within r<50 arcsec in one cluster field, representing a potential overdensity structure that has a physical origin at a confidence level >99.974% from Poisson statistics. Follow-up spectroscopic observations with ALMA and JWST are crucial to confirm whether it is associated with a protocluster at similar redshifts., Comment: Accepted for publication in The Astrophysical Journal, 19 pages, 13 figures, 1 table

We use machine learning techniques to investigate their performance in classifying active galactic nuclei (AGNs), including X-ray selected AGNs (XAGNs), infrared selected AGNs (IRAGNs), and radio selected AGNs (RAGNs). Using known physical parameters in the Cosmic Evolution Survey (COSMOS) field, we are able to well-established training samples in the region of Hyper Suprime-Cam (HSC) survey. We compare several Python packages (e.g., scikit-learn, Keras, and XGBoost), and use XGBoost to identify AGNs and show the performance (e.g., accuracy, precision, recall, F1 score, and AUROC). Our results indicate that the performance is high for bright XAGN and IRAGN host galaxies. The combination of the HSC (optical) information with the Wide-field Infrared Survey Explorer (WISE) band-1 and WISE band-2 (near-infrared) information perform well to identify AGN hosts. For both type-1 (broad-line) XAGNs and type-1 (unobscured) IRAGNs, the performance is very good by using optical to infrared information. These results can apply to the five-band data from the wide regions of the HSC survey, and future all-sky surveys., Comment: accepted for publication in ApJ

We examine the robustness of the color-color selection of quiescent galaxies (QGs) against contamination of dusty star-forming galaxies using the latest submillimeter data. We selected 18,304 QG candidates out to $z\sim$ 3 using the commonly adopted $NUV-r-J$ selection based on the high-quality multi-wavelength COSMOS2015 catalog. Using extremely deep 450 and 850 $\mu$m catalogs from the latest JCMT SCUBA-2 Large Programs, S2COSMOS, and STUDIES, as well as ALMA submillimeter, VLA 3 GHz, and $Spitzer$ MIPS 24 $\mu$m catalogs, we identified luminous dusty star-forming galaxies among the QG candidates. We also conducted stacking analyses in the SCUBA-2 450 and 850 $\mu$m images to look for less-luminous dusty galaxies among the QG candidates. By cross-matching to the 24 $\mu$m and 3 GHz data, we were able to identify a sub-group of "IR-radio-bright" QGs who possess a strong 450 and 850 $\mu$m stacking signal. The potential contamination of these luminous and less-luminous dusty galaxies accounts for approximately 10% of the color-selected QG candidates. In addition, there exists a spatial correlation between the luminous star-forming galaxies and the QGs at a $\lesssim60$ kpc scale. Finally, we found a high QG fraction among radio AGNs at $z<$ 1.5. Our data show a strong correlation between QGs and radio AGNs, which may suggest a connection between the quenching process and the radio-mode AGN feedback., Comment: This paper is accepted for publication on ApJ

Wide-field sub-millimetre surveys have driven many major advances in galaxy evolution in the past decade, but without extensive follow-up observations the coarse angular resolution of these surveys limits the science exploitation. This has driven the development of various analytical deconvolution methods. In the last half a decade Generative Adversarial Networks have been used to attempt deconvolutions on optical data. Here we present an autoencoder with a novel loss function to overcome this problem in the sub-millimeter wavelength range. This approach is successfully demonstrated on Herschel SPIRE 500$\mu$m COSMOS data, with the super-resolving target being the JCMT SCUBA-2 450$\mu$m observations of the same field. We reproduce the JCMT SCUBA-2 images with high fidelity using this autoencoder. This is quantified through the point source fluxes and positions, the completeness and the purity., Comment: Published by MNRAS in Volume 507 issue 1 October 2021, 12 pages, 7 figures. https://doi.org/10.1093/mnras/stab2195

We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies., Comment: ApJ accepted

We present far-infrared (FIR) properties of an extremely luminous infrared galaxy (ELIRG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1 (WISE1013+6112). This ELIRG is selected as an IR-bright dust-obscured galaxy (DOG) based on the photometry from the Sloan digital sky survey (SDSS) and wide-field infrared survey explorer (WISE). In order to derive its accurate IR luminosity, we perform follow-up observations at 89 and 154 $\mu$m using the high-resolution airborne wideband camera-plus (HAWC+) on board the 2.7-m stratospheric observatory for infrared astronomy (SOFIA) telescope. We conduct spectral energy distribution (SED) fitting with CIGALE using 15 photometric data (0.4-1300 $\mu$m). We successfully pin down FIR SED of WISE1013+6112 and its IR luminosity is estimated to be $L_{\rm IR}$ = (1.62 $\pm$ 0.08) $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the universe. We determine the dust temperature of WISE1013+6112 is $T_{\rm dust}$ = 89 $\pm$ 3 K, which is significantly higher than that of other populations such as SMGs and FIR-selected galaxies at similar IR luminosities. The resultant dust mass is $M_{\rm dust}$ = (2.2 $\pm$ 0.1) $\times 10^{8}$ $M_{\odot}$. This indicates that WISE1013+6112 has a significant active galactic nucleus (AGN) and star-forming activity behind a large amount of dust., Comment: 9 pages, 4 figures, and 2 tables, accepted for publication in ApJ

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$., Comment: Published in the ApJ

Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$\beta$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution., Comment: 13 pages, 10 figures, published in MNRAS

We present structural parameters and morphological properties of faint 450-um selected submillimeter galaxies (SMGs) from the JCMT Large Program, STUDIES, in the COSMOS-CANDELS region. Their properties are compared to an 850um selected and a matched star-forming samples. We investigate stellar structures of 169 faint 450-um sources (S450=2.8-29.6mJy; S/N>4) at z<3 using HST near-infrared observations. Based on our spectral energy distribution fitting, half of such faint SMGs (LIR=10^11.65+-0.98Lsun) lie above the star-formation rate (SFR)/stellar mass plane. The size-mass relation shows that these SMGs are generally similar to less-luminous star-forming galaxies selected by NUV-r vs. r-J colors. Because of the intrinsic luminosity of the sample, their rest-frame optical emission is less extended than the 850um sources (S850>2mJy), and more extended than the star-forming galaxies in the same redshift range. For the stellar mass and SFR matched sample at z~=1 and z~=2, the size differences are marginal between faint SMGs and the matched galaxies. Moreover, faint SMGs have similar Sersic indices and projected axis ratios as star-forming galaxies with the same stellar mass and SFR. Both SMGs and the matched galaxies show high fractions (~70%) of disturbed features at z~=2, and the fractions depend on the SFRs. These suggest that their star formation activity is related to galaxy merging, and the stellar structures of SMGs are similar to those of star-forming galaxies. We show that the depths of submillimeter surveys are approaching the lower luminosity end of star-forming galaxies, allowing us to detect galaxies on the main sequence., Comment: accepted for publication in ApJ

We present the discovery of an extremely-luminous dust-obscured galaxy (DOG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1. This DOG is selected as a candidate of extremely-luminous infrared (IR) galaxies based on the photometry from the Sloan Digital Sky Survey and Wide-field Infrared Survey Explorer. In order to derive its accurate IR luminosity, we perform follow-up observations at 450 and 850 $\mu$m using the Submillimetre Common User Bolometer Array 2 on the James Clerk Maxwell Telescope, and at 870 and 1300 $\mu$m using the Submillimeter Array, which enable us to pin down its IR Spectral Energy Distribution (SED). We perform SED fitting using 14 photometric data (0.4 - 1300 $\mu$m) and estimate its IR luminosity, $L_{\rm IR}$ (8-1000 $\mu$m), to be $2.2^{+1.5}_{-1.0}$ $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the Universe. The energy contribution from an active galactic nucleus (AGN) to the IR luminosity is $94^{+6}_{-20}$%, which indicates it is an AGN-dominated DOG. On the other hand, its stellar mass ($M_*$) and star formation rate (SFR) are $\log \,(M_\ast/M_{\odot})$ = $11.2^{+0.6}_{-0.2}$ and $\log \,({\rm SFR}/M_{\odot}\,{\rm yr}^{-1}$) = $3.1^{+0.2}_{-0.1}$, respectively, which means that this DOG can be considered as a starburst galaxy in $M_*$--SFR plane. This extremely-luminous DOG shows significant AGN and star forming activity that provides us an important laboratory to probe the maximum phase of the co-evolution of galaxies and supermassive black holes., Comment: 10 pages, 6 figures, and 1 table, accepted for publication in ApJ

We present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available., Comment: Accepted for publication in ApJS

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies., Comment: ApJ accepted. revised version uploaded