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 report Atacama Large Millimetre/submillimeter Array (ALMA) Cycle-5 follow-up observations of two candidate [CII] emitters at z ~ 6 in the ALMA deep field in SSA22 (ADF22). The candidates were detected blindly in a Cycle-2 ALMA survey covering ~ 5 square arcmins, with a single tuning, along with two CO lines associated with galaxies at lower redshifts. Various tests suggested at least one of the two > 6-sigma [CII] candidates should be robust (Hayatsu et al. 2017). Nevertheless, our new, deeper observations recover neither candidate, demonstrating a higher contamination rate than expected. The cause of the spurious detections is under investigation but at present it remains unclarified.

The matched filter (MF) represents one of the main tools to detect signals from known sources embedded in the noise. In the Gaussian case the noise is assumed to be the realization of a Gaussian random field (GRF). The most important property of the MF, the maximization of the probability of detection subject to a constant probability of false detection or false alarm (PFA), makes it one of the most popular techniques. However, the MF technique relies upon the a priori knowledge of the number and the position of the searched signals in the GRF which usually are not available. A typical way out is to assume that the position of a signal coincides with one of the peaks in the matched filtered data. A detection is claimed when the probability that a given peak is due only to the noise (i.e. the PFA) is smaller than a prefixed threshold. In this case the probability density function (PDF) of the amplitudes has to be used for the computation of the PFA, which is different from the Gaussian. Moreover, the probability that a detection is false depends on the number of peaks present in the filtered GRF, the greater the number of peaks in a GRF, the higher the probability of peaks due to the noise that exceed the detection threshold. If not taken into account, the PFA can be severely underestimated. Many solutions proposed to this problem are non-parametric hence not able to exploit all the available information. This limitation has been overcome by means of two efficient parametric approaches, one based on the PDF of the peak amplitudes of a smooth and isotropic GRF whereas the other uses the Gumbel distribution (the asymptotic PDF of the corresponding extreme). Simulations and ALMA maps show that, although the two methods produce almost identical results, the first is more flexible and allows us to check the reliability of the detection procedure.

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.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

To search for dust-obscured star-formation activity in the early Universe, it is essential to obtain a deep and wide submillimeter/millimeter map. The advent of the Atacama Large Millimeter/submillimeter Array (ALMA) has enabled us to obtain such maps at sufficiently high spatial resolution to be free from source confusion. We present a new 1.1mm map obtained by ALMA in the SSA22 field. SSA22 contains a remarkable proto-cluster at z=3.09 and is therefore an ideal region to investigate the role of large-scale cosmic web on dust-obscured star formation. The typical 1sigma depth of our map is 73 uJy/beam at a 0.5 resolution; combined with earlier, archived observations, we map an area of 20 arcmin^2 (71 comoving Mpc^2 at z=3.09). Within the combined survey area we have detected 35 sources at a signal-to-noise ratio (S/N) >5, with flux densities, S1.1mm=0.43--5.6 mJy, equivalent to star-formation rates of >=100--1000 Msun/yr at z=3.09, for a Chabrier initial mass function; of these, 17 are new detections. The cumulative number counts show a factor 3--5x excess compared to blank fields. The excess suggests enhanced dust-enshrouded star-formation activity in the proto-cluster on a 10 comoving Mpc scale, indicating accelerated galaxy evolution in this overdense region.

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.

We present the results of a blind millimeter line emitter search using ALMA Band 6 data with a single frequency tuning toward four gravitational lensing clusters (RXJ1347.5-1145, Abell S0592, MACS J0416.1-2403, and Abell 2744). We construct three-dimensional signal-to-noise ratio (S/N) cubes with 60 MHz and 100 MHz binning, and search for millimeter line emitters. We do not detect any line emitters with a peak S/N > 5, although we do find a line emitter candidate with a peak S/N ~ 4.5. These results provide upper limits to the CO(3-2), CO(4-3), CO(5-4), and [CII] luminosity functions at z ~ 0.3, 0.7, 1.2, and 6, respectively. Because of the magnification effect of gravitational lensing clusters, the new data provide the first constraints on the CO and [CII] luminosity functions at unprecedentedly low luminosity levels, i.e., down to $<~ 10^{-3}-10^{-1} $Mpc$^{-3}$ dex$^{-1}$ at $L'_{CO} ~ 10^8-10^{10}$ K km $s^{-1}$ pc$^2$ and $<~ 10^{-3}-10^{-2}$ Mpc$^{-3}$ dex$^{-1}$ at $L_[CII] ~ 10^8-10^{10} L_\odot$, respectively. Although the constraints to date are not stringent yet, we find that the evolution of the CO and [CII] luminosity functions are broadly consistent with the predictions of semi-analytical models. This study demonstrates that the wide observations with a single frequency tuning toward gravitational lensing clusters are promising for constraining the CO and [CII] luminosity functions.

We report the identification of four millimeter line emitting galaxies with the Atacama Large Milli/submillimeter Array (ALMA) in SSA22 Field (ADF22). We analyze the ALMA 1.1 mm survey data, with an effective survey area of 5 arcmin$^2$, a frequency range of 253.1--256.8 and 269.1--272.8 GHz, angular resolution of 0".7 and RMS noise of 0.8 mJy beam$^{-1}$ at 36 km s$^{-1}$ velocity resolution. We detect four line emitter candidates with significance levels above $6 \sigma$. We identify one of the four sources as a CO(9-8) emitter at $z = 3.1$ in a member of the proto-cluster known in this field. Another line emitter with an optical counterpart is likely a CO(4-3) emitter at $z = 0.7$. The other two sources without any millimeter continuum or optical/near-infrared counterpart are likely to be [CII] emitter candidates at $z = 6.0$ and $6.5$. The equivalent widths of the [CII] candidates are consistent with those of confirmed high-redshift [CII] emitters and candidates, and are a factor of 10 times larger than that of the CO(9-8) emitter detected in this search. The [CII] luminosity of the candidates are $4-7 \times 10^8~\rm L_\odot$. The star formation rates (SFRs) of these sources are estimated to be $10-20~\rm M_{\odot}~yr^{-1}$ if we adopt an empirical [CII] luminosity - SFR relation. One of them has a relatively low-S/N ratio, but shows features characteristic of emission lines. Assuming that at least one of the two candidates is a [CII] emitter, we derive a lower limit of [CII]-based star formation rate density (SFRD) at $z~\sim~6$. The resulting value of $> 10^{-2}$ $\rm M_\odot yr^{-1} Mpc^{-3}$ is consistent with the dust-uncorrected UV-based SFRD. Future millimeter/submillimeter surveys can be used to detect a number of high redshift line emitters, with which to study the star formation history in the early Universe.