We make use of the ASAGAO, deep 1.2 mm continuum observations of a 26 arcmin$^2$ region in the GOODS-South field obtained with ALMA, to probe dust-enshrouded star formation in $K$-band selected (i.e., stellar mass selected) galaxies, which are drawn from the ZFOURGE catalog. Based on the ASAGAO combined map, which was created by combining ASAGAO and ALMA archival data in the GOODS-South field, we find that 24 ZFOURGE sources have 1.2 mm counterparts with a signal-to-noise ratio $>$ 4.5 (1$\sigma\simeq$ 30 - 70 $\mu$Jy beam$^{-1}$ at 1.2 mm). Their median redshift is estimated to be $z_\mathrm{median}=$ 2.38 $\pm$ 0.14. They generally follow the tight relationship of the stellar mass versus star formation rate (i.e., the main sequence of star-forming galaxies). ALMA-detected ZFOURGE sources exhibit systematically larger infrared (IR) excess (IRX $\equiv L_\mathrm{IR}/L_\mathrm{UV}$) compared to ZFOURGE galaxies without ALMA detections even though they have similar redshifts, stellar masses, and star formation rates. This implies the consensus stellar-mass versus IRX relation, which is known to be tight among rest-frame-UV-selected galaxies, can not fully predict the ALMA detectability of stellar-mass-selected galaxies. We find that ALMA-detected ZFOURGE sources are the main contributors to the cosmic IR star formation rate density at $z$ = 2 - 3.

We are developing an ultra-wideband spectroscopic instrument, DESHIMA (DEep Spectroscopic HIgh-redshift MApper), based on the technologies of an on-chip filter-bank and Microwave Kinetic Inductance Detector (MKID) to investigate dusty star-burst galaxies in the distant universe at millimeter and submillimeter wavelength. An on-site experiment of DESHIMA was performed using the ASTE 10-m telescope. We established a responsivity model that converts frequency responses of the MKIDs to line-of-sight brightness temperature. We estimated two parameters of the responsivity model using a set of skydip data taken under various precipitable water vapor (PWV, 0.4-3.0 mm) conditions for each MKID. The line-of-sight brightness temperature of sky is estimated using an atmospheric transmission model and the PWVs. As a result, we obtain an average temperature calibration uncertainty of $1\sigma=4$%, which is smaller than other photometric biases. In addition, the average forward efficiency of 0.88 in our responsivity model is consistent with the value expected from the geometrical support structure of the telescope. We also estimate line-of-sight PWVs of each skydip observation using the frequency response of MKIDs, and confirm the consistency with PWVs reported by the Atacama Large Millimeter/submillimeter Array.

Ultra-wideband 3D imaging spectrometry in the millimeter-submillimeter (mm-submm) band is an essential tool for uncovering the dust-enshrouded portion of the cosmic history of star formation and galaxy evolution. However, it is challenging to scale up conventional coherent heterodyne receivers or free-space diffraction techniques to sufficient bandwidths ($\geq$1 octave) and numbers of spatial pixels (>$10^2$). Here we present the design and first astronomical spectra of an intrinsically scalable, integrated superconducting spectrometer, which covers 332-377 GHz with a spectral resolution of $F/\Delta F \sim 380$. It combines the multiplexing advantage of microwave kinetic inductance detectors (MKIDs) with planar superconducting filters for dispersing the signal in a single, small superconducting integrated circuit. We demonstrate the two key applications for an instrument of this type: as an efficient redshift machine, and as a fast multi-line spectral mapper of extended areas. The line detection sensitivity is in excellent agreement with the instrument design and laboratory performance, reaching the atmospheric foreground photon noise limit on sky. The design can be scaled to bandwidths in excess of an octave, spectral resolution up to a few thousand and frequencies up to $\sim$1.1 THz. The miniature chip footprint of a few $\mathrm{cm^2}$ allows for compact multi-pixel spectral imagers, which would enable spectroscopic direct imaging and large volume spectroscopic surveys that are several orders of magnitude faster than what is currently possible.

Our current knowledge of cosmic star-formation history during the first two billion years (corresponding to redshift z >3) is mainly based on galaxies identified in rest-frame ultraviolet light. However, this population of galaxies is known to under-represent the most massive galaxies, which have rich dust content and/or old stellar populations. This raises the questions of the true abundance of massive galaxies and the star-formation-rate density in the early universe. Although several massive galaxies that are invisible in the ultraviolet have recently been confirmed at early epochs, most of them are extreme starbursts with star-formation rates exceeding 1000 solar masses per year, suggesting that they are unlikely to represent the bulk population of massive galaxies. Here we report submillimeter (wavelength 870um) detections of 39 massive star-forming galaxies at z > 3, which are unseen in the spectral region from the deepest ultraviolet to the near-infrared. With a space density of about $2 \times 10^{-5}$ per cubic megaparsec (two orders of magnitudes higher than extreme starbursts) and star-formation rates of 200 solar masses per year, these galaxies represent the bulk population of massive galaxies that have been missed from previous surveys. They contribute a total star-formation-rate density ten times larger than that of equivalently massive ultraviolet-bright galaxies at z >3. Residing in the most massive dark matter halos at their redshifts, they are probably the progenitors of the largest present-day galaxies in massive groups and clusters. Such a high abundance of massive and dusty galaxies in the early universe challenges our understanding of massive-galaxy formation.

State-of-the-art rest-frame UV and FIR photometric and spectroscopic observations are now pushing the redshift frontiers of galaxy formation studies up to $z\sim9-11$ and beyond. Recent HST observations unveiled the presence of a star-forming galaxy exhibiting the Lyman break at $\lambda_{\rm obs}=1.47\pm0.01$ $\mu$m, i.e., a $z=11.09^{+0.08}_{-0.12}$ galaxy with a stellar mass of $\sim10^9 M_\odot$, demonstrating that galaxy build-up was well underway early in the epoch of reionization (EoR) at $z>10$. Targeted spectroscopy of a lensed Lyman break galaxy uncovers the earliest metals known to date up to $z=9.1096\pm0.0006$ by detecting the bright [OIII] 88~$\mu$m nebular line, indicating the onset of star formation 250 million years after the Big Bang, i.e., corresponding to a redshift of $z\sim15$. These latest findings lead us to a number of key questions: How and when metal enrichment happened in the EoR? What was the nature of the earliest-epoch star-forming galaxies at $z=10-15$? What was the spatial distribution of such galaxies, and what was the relation to the putative large-scale ionization bubbles during the EoR? What were the dark-halo masses of such earliest-epoch star-forming galaxies? To address all these questions, we need to uncover a statistically large number of $z=10-15$ galaxies in the pre-reionization era. Here we argue two possible pathways: (1) a wide-area, sensitive blind spectroscopic survey of [OIII] 88 $\mu$m line-emitting galaxies at submillimeter wavelengths, and (2) an ultra-wide-area, high-cadence photometric survey of transient sources at radio-to-(sub)millimeter wavelengths, together with the immediate follow-up spectroscopy with an ultra-wide-band spectrograph, to catch the pop-III $\gamma$-ray bursts.

We report detections of two 1.2 mm continuum sources ($S_\mathrm{1.2mm}$ ~ 0.6 mJy) without any counterparts in the deep $H$- and/or $K$-band image (i.e., $K$-band magnitude $\gtrsim$ 26 mag). These near-infrared-dark faint millimeter sources are uncovered by ASAGAO, a deep and wide-field ($\simeq$ 26 arcmin$^2$) Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm survey. One has a red IRAC (3.6 and 4.5 $\mu$m) counterpart, and the other has been independently detected at 850 and 870 $\mu$m using SCUBA2 and ALMA Band 7, respectively. Their optical to radio spectral energy distributions indicate that they can lie at $z \gtrsim$ 3-5 and can be in the early phase of massive galaxy formation. Their contribution to the cosmic star formation rate density is estimated to be ~ 1 $\times$ 10$^{-3}$ $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ if they lie somewhere in the redshift range of $z$ ~ 3-5. This value can be consistent with, or greater than that of bright submillimeter galaxies ($S_\mathrm{870\mu m}>$ 4.2 mJy) at $z$ ~ 3-5. We also uncover 3 more candidates near-infrared-dark faint ALMA sources without any counterparts ($S_\mathrm{1.2mm}$ ~ 0.45-0.86 mJy). These results show that an unbiased ALMA survey can reveal the dust-obscured star formation activities, which were missed in previous deep optical/near-infrared surveys.

We present the survey design, data reduction, construction of images, and source catalog of the Atacama Large Millimeter/submillimeter Array (ALMA) twenty-six arcmin^2 survey of GOODS-S at one-millimeter (ASAGAO). ASAGAO is a deep (1sigma ~ 61 uJy/beam for a 250 klambda-tapered map with a synthesized beam size of 0.51" x 0.45") and wide area (26 arcmin^2) survey on a contiguous field at 1.2 mm. By combining with ALMA archival data in the GOODS-South field, we obtained a deeper map in the same region (1sigma ~ 30 uJy/beam for a deep region with a 250 klambda-taper, and a synthesized beam size of 0.59" x 0.53"), providing the largest sample of sources (25 sources at >=5.0sigma, 45 sources at >=4.5sigma) among ALMA blank-field surveys to date. The number counts shows that 52(+11 -8)% of the extragalactic background light at 1.2 mm is resolved into discrete sources at S1.2m > 135 uJy. We create infrared (IR) luminosity functions (LFs) in the redshift range of z = 1-3 from the ASAGAO sources with KS-band counterparts, and constrain the faintest luminosity of the LF at 2.0 < z < 3.0. The LFs are consistent with previous results based on other ALMA and SCUBA-2 observations, which suggest a positive luminosity evolution and negative density evolution with increasing redshift. We find that obscured star-formation of sources with IR luminosities of log(L(IR)/Lsun)} ~> 11.8 account for ~~60%-90% of the z ~ 2 cosmic star-formation rate density.

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges.