Source author record

Yuichi Matsuda

Yuichi Matsuda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2023arXiv

The JCMT SCUBA-2 Survey of the James Webb Space Telescope North Ecliptic Pole Time-Domain Field

The James Webb Space Telescope Time-Domain Field (JWST-TDF) is an $\sim$14$'$ diameter field near the North Ecliptic Pole that will be targeted by one of the JWST Guaranteed Time Observations programs. Here, we describe our James Clerk Maxwell Telescope SCUBA-2 850 $μ$m imaging of the JWST-TDF and present the submillimeter source catalog and properties. We also present a catalog of radio sources from Karl J. Jansky Very Large Array 3 GHz observations of the field. These observations were obtained to aid JWST's study of the dust-obscured galaxies that contribute significantly to the cosmic star formation at high redshifts. Our deep 850 $μ$m map covers the JWST TDF at a noise level of $σ_{850}$ = 1.0 mJy beam$^{-1}$, detecting 83/31 sources in the main/supplementary signal-to-noise ratio (S/N $>$ 4 / S/N = 3.5 - 4) sample respectively. The 3 GHz observations cover a 24$'$ diameter field with a 1 $σ$ noise of 1$μ$Jy beam$^{-1}$ at a 0$.\!\!^{\prime\prime}$7 FWHM. We identified eighty-five 3 GHz counterparts to sixty-six 850 $μ$m sources and then matched these with multiwavelength data from the optical to the mid-infrared wave bands. We performed spectral energy distribution fitting for 61 submillimeter galaxies (SMGs) matched with optical/near-infrared data, and found that SMGs at S/N $>$ 4 have a median value of $z_{phot} = $2.22 $\pm$ 0.12, star formation rates of 300 $\pm$ 40 M$_{\odot}\,{\rm yr^{-1}}$ (Chabrier initial mass function), and typical cold dust masses of 5.9 $\pm$ 0.7 $ \times$ 10$^{8} $M$_{\odot}$, in line with bright SMGs from other surveys. The large cold dust masses indicate correspondingly large cool gas masses, which we suggest are a key factor necessary to drive the high star formation rates seen in this population

preprint2022arXiv

An AGN with an ionized gas outflow in a massive quiescent galaxy in a protocluster at $\bf z=3.09$

We report the detection of an ionized gas outflow from an $X$-ray active galactic nucleus (AGN) hosted in a massive quiescent galaxy in a protocluster at $z=3.09$ (J221737.29+001823.4). It is a type-2 QSO with broad ($W_{80}>1000$ km s$^{-1}$) and strong ($\log (L_{\rm [OIII]}$ / erg s$^{-1})\approx43.4$) [O {\footnotesize III}]$λλ$4959,5007 emission lines detected by slit spectroscopy in three-position angles using Multi-Object Infra-Red Camera and Spectrograph (MOIRCS) on the Subaru telescope and the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE) on the Keck-I telescope. In the all slit directions, [O {\footnotesize III}] emission is extended to $\sim15$ physical kpc and indicates a powerful outflow spreading over the host galaxy. The inferred ionized gas mass outflow rate is $\rm 22\pm3~M_{\odot}~yr^{-1}$. Although it is a radio source, according to the line diagnostics using H$β$, [O {\footnotesize II}], and [O {\footnotesize III}], photoionization by the central QSO is likely the dominant ionization mechanism rather than shocks caused by radio jets. On the other hand, the spectral energy distribution of the host galaxy is well characterized as a quiescent galaxy that has shut down star formation by several hundred Myr ago. Our results suggest a scenario that QSOs are powered after the shut-down of the star formation and help to complete the quenching of massive quiescent galaxies at high redshift.

preprint2022arXiv

Detection of nitrogen and oxygen in a galaxy at the end of reionization

We present observations of [NII] 205 $μ$m, [OIII] 88 $μ$m and dust emission in a strongly-lensed, submillimeter galaxy (SMG) at $z=6.0$, G09.83808, with the Atacama Large Millimeter/submillimeter Array (ALMA). Both [NII] and [OIII] line emissions are detected at $>12σ$ in the 0.8$"$-resolution maps. Lens modeling indicates that the spatial distribution of the dust continuum emission is well characterized by a compact disk with an effective radius of 0.64$\pm$0.02 kpc and a high infrared surface brightness of $Σ_\mathrm{IR}=(1.8\pm0.3)\times10^{12}~L_\odot$ kpc$^{-2}$. This result supports that G09.83808 is the progenitors of compact quiescent galaxies at $z\sim4$, where the majority of its stars are expected to be formed through a strong and short burst of star formation. G09.83808 and other lensed SMGs show a decreasing trend of the [NII] line to infrared luminosity ratio with increasing continuum flux density ratio between 63 $μ$m and 158 $μ$m, as seen in local luminous infrared galaxies (LIRGs). The decreasing trend can be reproduced by photoionization models with increasing ionization parameters. Furthermore, by combining the [NII]/[OIII] luminosity ratio with far-infrared continuum flux density ratio in G09.83808, we infer that the gas phase metallicity is already $Z\approx 0.5-0.7~Z_\odot$. G09.83808 is likely one of the earliest galaxies that has been chemically enriched at the end of reionization.

preprint2022arXiv

FOREVER22: galaxy formation in protocluster regions

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.

preprint2022arXiv

Where Outflows Meet Inflows: Gas Kinematics in SSA22 Lyman-$α$ Blob 2 Decoded by Advanced Radiative Transfer Modelling

We present new spectroscopic observations of Lyman-$α$ (Ly$α$) Blob 2 ($z \sim$ 3.1). We observed extended Ly$α$ emission in three distinct regions, where the highest Ly$α$ surface brightness (SB) center is far away from the known continuum sources. We searched through the MOSFIRE slits that cover the high Ly$α$ SB regions, but were unable to detect any significant nebular emission near the highest SB center. We further mapped the flux ratio of the blue peak to the red peak and found it is anti-correlated with Ly$α$ SB with a power-law index of $\sim$ -0.4. We used radiative transfer models with both multiphase, clumpy and shell geometries and successfully reproduced the diverse Ly$α$ morphologies. We found that most spectra suggest outflow-dominated kinematics, while 4/15 spectra imply inflows. A significant correlation exists between parameter pairs, and the multiphase, clumpy model may alleviate previously reported discrepancies. We also modeled Ly$α$ spectra at different positions simultaneously and found that the variation of the inferred clump outflow velocities can be approximately explained by line-of-sight projection effects. Our results support the `central powering + scattering' scenario, i.e. the Ly$α$ photons are generated by a central powering source and then scatter with outflowing, multiphase HI gas while propagating outwards. The infalling of cool gas near the blob outskirts shapes the observed blue-dominated Ly$α$ profiles, but its energy contribution to the total Ly$α$ luminosity is less than 10%, i.e. minor compared to the photo-ionization by star-forming galaxies and/or AGNs.

preprint2020arXiv

3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$α$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$α$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$α$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

preprint2020arXiv

ALMA band 8 observations of DLA2233+131 at z=3.150

We present our ALMA Band 8 observations of a damped Ly$α$ absorption (DLA) system at $z$=3.150 observed in the spectrum of the quasar Q2233+131 at $z$=3.295. The optical counterpart of this DLA has been identified and it shows a double-peaked Ly$α$ emission line. Since one possible origin of DLAs at high redshift is an outflowing gas from star-forming galaxies, DLA2233+131 provides a good laboratory to investigate the nature of high-$z$ DLAs. Motivated by this, we have carried out ALMA band 8 observations to study the [C II] line in this system. However, we do not detect any significant emission line in the observed pass bands. Instead, we have serendipitously found three submm continuum sources in the observed sky area. One appears to be the quasar Q2233+131 itself while the other two sources are newly identified submm galaxies (SMGs), called SMG1 and SMG2 in this paper. They are located at a separation of 400.7 and 800.1 from Q2233+131, respectively. Their 646 $μ$m fluxes are 6.35 mJy and 6.43 mJy, respectively, being higher than that of Q2233+131, 3.62 mJy. Since these two SMGs are not detected in the optical images obtained with the Hubble Space Telescope and the Subaru Telescope, they have a very red spectral energy distribution. It is, therefore,suggested that they are high-redshift galaxies or very dusty galaxies at intermediate redshift although we cannot rule out the possibility that they are optically very faint SMG analogs at low redshift. Follow up observations will be necessary to explore the nature of this interesting region.

preprint2020arXiv

ALMA uncovers the [CII] emission and warm dust continuum in a z = 8.31 Lyman break galaxy

We report on the detection of the [CII] 157.7 $μ$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $μ$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $μ$m dust emission (90 $μ$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $μ$m rest-frame) down to 18 $μ$Jy (3$σ$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $μ$Jy continuum emission at 850 $μ$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($β_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $β_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

preprint2020arXiv

Dual supermassive black holes at close separation revealed by the Hyper Suprime-Cam Subaru Strategic Program

The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

preprint2020arXiv

FIR-luminous [CII] emitters in the ALMA-SCUBA-2 COSMOS survey (AS2COSMOS): The nature of submillimeter galaxies in a 10 comoving Mpc-scale structure at z~4.6

We report the discovery of a 10 comoving Mpc-scale structure traced by massive submillimeter galaxies (SMGs) at z~4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources ($S_{850μ{\rm m}}\geq$ 6.2 mJy) across 1.6 degrees$^2$ in the COSMOS field. We identify four [CII] emitting SMGs and two probable [CII] emitting SMG candidates at z=4.60-4.64 with velocity-integrated signal-to-noise ratio of SNR>8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330-550 km s$^{-1}$ and the inferred host dark-matter halo masses are ~2-8 $\times$ 10$^{12}$M$_{\odot}$. From their estimated halo masses and [CII] luminosity function, we suggest that these galaxies have a high (50-100%) duty cycle and high (~0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute $\simeq$2% to the total star-formation rate density at z=4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ~10$^{15}$M$_{\odot}$ cluster at z=0. This survey reveals synchronized dusty starburst in massive halos at z>4, which could be driven by mergers or fed by smooth gas accretion.

preprint2020arXiv

Large Population of ALMA Galaxies at z>6 with Very High [OIII]88um to [CII]158um Flux Ratios: Evidence of Extremely High Ionization Parameter or PDR Deficit?

We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

preprint2020arXiv

Spatially Resolved Molecular Gas Properties of Host Galaxy of Type I Superluminous Supernova SN 2017egm

We present the results of CO(1-0) observations of the host galaxy of a Type I superluminous supernova (SLSN-I), SN2017egm, one of the closest SLSNe-I at z = 0.03063, by using the Atacama Large Millimeter/submillimeter Array. The molecular gas mass of the host galaxy is $M_{\rm gas} = (4.8 \pm 0.3) \times 10^9$ $M_{\odot}$, placing it on the sequence of normal star-forming galaxies in an $M_{\rm gas}$-star-formation rate (SFR) plane. The molecular hydrogen column density at the location of SN2017egm is higher than that of the Type II SN PTF10bgl, which is also located in the same host galaxy, and those of other Type II and Ia SNe located in different galaxies, suggesting that SLSNe-I have a preference for a dense molecular gas environment. On the other hand, the column density at the location of SN2017egm is comparable to those of Type Ibc SNe. The surface densities of molecular gas and the SFR at the location of SN2017egm are consistent with those of spatially resolved local star-forming galaxies and follow the Schmidt-Kennicutt relation. These facts suggest that SLSNe-I can occur in environments with the same star-formation mechanism as in normal star-forming galaxies.

preprint2019arXiv

Planck far-infrared detection of Hyper Suprime-Cam protoclusters at $\bf z\sim4$: hidden AGN and star formation activity

We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~μ$m wavelength range is successfully detected with $>5σ$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $μ$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical.

preprint2018arXiv

The onset of star formation 250 million years after the Big Bang

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.

preprint2016arXiv

Clumpy and Extended Starbursts in the Brightest Unlensed Submillimeter Galaxies

The central structure in three of the brightest unlensed z=3-4 submillimeter galaxies are investigated through 0.015" - 0.05" (120 -- 360~pc) 860 micron continuum images obtained using the Atacama Large Millimeter/submillimeter Array (ALMA). The distribution in the central kpc in AzTEC1 and AzTEC8 are extremely complex, and they are composed of multiple ~200 pc clumps. AzTEC4 consists of two sources that are separated by ~1.5 kpc, indicating a mid-stage merger. The peak star formation rate densities in the central clumps are ~300 - 3000 Msun/yr/kpc^2, suggesting regions with extreme star formation near the Eddington Limit. By comparing the flux obtained by ALMA and Submillimeter Array (SMA), we find that 68-90% of the emission is extended (> 1 kpc) in AzTEC 4 and 8. For AzTEC1, we identify at least 11 additional compact (~200 pc) clumps in the extended 3 - 4 kpc region. Overall, the data presented here suggest that the luminosity surface densities observed at < 150 pc scales are roughly similar to that observed in local ULIRGs, as in the eastern nucleus of Arp 220. Between 10 to 30% of the 860 micron continuum is concentrated in clumpy structures in the central kpc while the remaining flux is distributed over > 1 kpc regions, some of which could also be clumpy. These sources can be explained by a rapid inflow of gas such as a merger of gas-rich galaxies, surrounded by extended and clumpy starbursts. However, the cold mode accretion model is not ruled out.

preprint2016arXiv

Detection of an oxygen emission line from a high redshift galaxy in the reionization epoch

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

preprint2016arXiv

Investigating the Relation between CO (3-2) and Far Infrared Luminosities for Nearby Merging Galaxies Using ASTE

We present the new single dish CO (3-2) emission data obtained toward 19 early stage and 7 late stage nearby merging galaxies using the Atacama Submillimeter Telescope Experiment (ASTE). Combining with the single dish and interferometric data of galaxies observed in previous studies, we investigate the relation between the CO (3-2) luminosity (L'CO(3-2)) and the far Infrared luminosity (LFIR) in a sample of 29 early stage and 31 late stage merging galaxies, and 28 nearby isolated spiral galaxies. We find that normal isolated spiral galaxies and merging galaxies have different slopes (alpha) in the log L'CO(3-2) - log LFIR plane (alpha ~ 0.79 for spirals and ~ 1.12 for mergers). The large slope (alpha > 1) for merging galaxies can be interpreted as an evidence for increasing Star Formation Efficiency (SFE=LFIR/L'CO(3-2)) as a function of LFIR. Comparing our results with sub-kpc scale local star formation and global star-burst activity in the high-z Universe, we find deviations from the linear relationship in the log L'CO(3-2) - log LFIR plane for the late stage mergers and high-z star forming galaxies. Finally, we find that the average SFE gradually increases from isolated galaxies, merging galaxies, and to high-z submillimeter galaxies / quasi-stellar objects (SMGs/QSOs). By comparing our findings with the results from numerical simulations, we suggest; (1) inefficient star-bursts triggered by disk-wide dense clumps occur in the early stage of interaction and (2) efficient star-bursts triggered by central concentration of gas occur in the final stage. A systematic high spatial resolution survey of diffuse and dense gas tracers is a key to confirm this scenario.

preprint2016arXiv

SXDF-ALMA 2 Arcmin^2 Deep Survey: Resolving and Characterizing the Infrared Extragalactic Background Light Down to 0.5 mJy

We present a multi-wavelength analysis of five submillimeter sources (S_1.1mm = 0.54-2.02 mJy) that were detected during our 1.1-mm-deep continuum survey in the SXDF-UDS-CANDELS field (2 arcmin^2, 1sigma = 0.055 mJy beam^-1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be ~ 4.1^{+5.4}_{-3.0} Jy deg^{-2}, which corresponds to ~ 16^{+22}_{-12}% of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multi-wavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity (L_IR ~ 1e12 L_sun or SFR ~ 100 M_sun yr^{-1}). By fitting the spectral energy distributions (SEDs) at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (z_photo) and stellar masses (M_*): z_photo ~ 1.3-2.5, M_* ~ (3.5-9.5)e10 M_sun. We also derived their star formation rates (SFRs) and specific SFRs (sSFRs) as ~ 30-200 M_sun yr^{-1} and ~ 0.8-2 Gyr^{-1}, respectively. These values imply that they are main-sequence star-forming galaxies.

preprint2016arXiv

SXDF-ALMA 2 arcmin$^2$ Deep Survey: 1.1-mm Number Counts

We report 1.1 mm number counts revealed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey Field (SXDF). The advent of ALMA enables us to reveal millimeter-wavelength number counts down to the faint end without source confusion. However, previous studies are based on the ensemble of serendipitously-detected sources in fields originally targeting different sources and could be biased due to the clustering of sources around the targets. We derive number counts in the flux range of 0.2-2 mJy by using 23 (>=4sigma) sources detected in a continuous 2.0 arcmin$^2$ area of the SXDF. The number counts are consistent with previous results within errors, suggesting that the counts derived from serendipitously-detected sources are not significantly biased, although there could be field-to-field variation due to the small survey area. By using the best-fit function of the number counts, we find that ~40% of the extragalactic background light at 1.1 mm is resolved at S(1.1mm) > 0.2 mJy.

preprint2016arXiv

SXDF-ALMA 2-arcmin^2 Deep Survey: Stacking of Rest-Frame Near-Infrared Selected Objects

We present stacking analyses on our ALMA deep 1.1 mm imaging in the SXDF using 1.6 μm and 3.6 μm selected galaxies in the CANDELS WFC3 catalog. We detect a stacked flux of ~0.03-0.05 mJy, corresponding to LIR < 10^11 Lsun and a star formation rate (SFR) of ~ 15 Msun/yr at z = 2. We find that galaxies brighter in the rest-frame near-infrared tend to be also brighter at 1.1 mm, and galaxies fainter than m[3.6um] = 23 do not produce detectable 1.1 mm emission. This suggests a correlation between stellar mass and SFR, but outliers to this correlation are also observed, suggesting strongly boosted star formation or extremely large extinction. We also find tendencies that redder galaxies and galaxies at higher redshifts are brighter at 1.1 mm. Our field contains z ~ 2.5 H-alpha emitters and a bright single-dish source. However, we do not find evidence of bias in our results caused by the bright source. By combining the fluxes of sources detected by ALMA and fluxes of faint sources detected with stacking, we recover a 1.1 mm surface brightness of up to 20.3 +/- 1.2 Jy/deg, comparable to the extragalactic background light measured by COBE. Based on the fractions of optically faint sources in our and previous ALMA studies and the COBE measurements, we find that approximately half of the cosmic star formation may be obscured by dust and missed by deep optical surveys, Much deeper and wider ALMA imaging is therefore needed to better constrain the obscured cosmic star formation history.

preprint2015arXiv

High-resolution ALMA observations of SDP.81. I. The innermost mass profile of the lensing elliptical galaxy probed by 30 milli-arcsecond images

We report a detailed modeling of a mass profile of a $z = 0.2999$ massive elliptical galaxy using 30 milli-arcsecond resolution 1-mm Atacama Large Millimeter/submillimeter Array (ALMA) images of the galaxy-galaxy lensing system SDP.81. The detailed morphology of the lensed multiple images of the $z = 3.042$ infrared-luminous galaxy, which is found to consist of tens of $\lesssim 100$-pc-sized star-forming clumps embedded in a $\sim 2$ kpc disk, are well reproduced by a lensing galaxy modeled by an isothermal ellipsoid with a 400 pc core. The core radius is consistent with that of the visible stellar light, and the mass-to-light ratio of $\sim 2\,M_{\odot}\,L_{\odot}^{-1}$ is comparable to the locally measured value, suggesting that the inner 1 kpc region is dominated by luminous matter. The position of the predicted mass centroid is consistent to within $\simeq 30$ mas with that of a non-thermal source detected with ALMA, which likely traces an active galactic nucleus of the foreground elliptical galaxy. A point source mass of $> 3 \times 10^8\,M_{\odot}$ mimicking a supermassive black hole is required to explain the non-detection of a central image of the background galaxy, although the black hole mass degenerates with the core radius of the elliptical galaxy. The required mass is consistent with that predicted from the well-known correlation between black hole mass and host velocity dispersion. Our analysis demonstrates the power of ALMA imaging of strong gravitational lensing for studying the innermost mass profiles and the central supermassive black hole of distant elliptical galaxies.

preprint2015arXiv

High-resolution ALMA Observations of SDP.81. II. Molecular Clump Properties of a Lensed Submillimeter Galaxy at z=3.042

We present spatially-resolved properties of molecular gas and dust in a gravitationally-lensed submillimeter galaxy H-ATLAS J090311.6+003906 (SDP.81) at $z=3.042$ revealed by the Atacama Large Millimeter/submillimeter Array (ALMA). We identified 14 molecular clumps in the CO(5-4) line data, all with a spatial scale of $\sim$50-300 pc in the source plane. The surface density of molecular gas ($Σ_{\rm H_2}$) and star-formation rate ($Σ_{\rm SFR}$) of the clumps are more than three orders of magnitude higher than those found in local spiral galaxies. The clumps are placed in the `burst' sequence in the $Σ_{\rm H_2}$-$Σ_{\rm SFR}$ plane, suggesting that $z \sim 3$ molecular clumps follow the star-formation law derived for local starburst galaxies. With our gravitational lens model, the positions in the source plane are derived for the molecular clumps, dust clumps, and stellar components identified in the {\sl Hubble Space Telescope} image. The molecular and dust clumps coexist in a similar region over $\sim$2 kpc, while the stellar components are offset at most by $\sim$5 kpc. The molecular clumps have a systematic velocity gradient in the north-south direction, which may indicate a rotating gas disk. One possible scenario is that the components of molecular gas, dust, and stars are distributed in a several-kpc scale rotating disk, and the stellar emission is heavily obscured by dust in the central star-forming region. Alternatively, SDP.81 can be explained by a merging system, where dusty starbursts occur in the region where the two galaxies collide, surrounded by tidal features traced in the stellar components.

preprint2015arXiv

SXDF-ALMA 1.5 arcmin^2 deep survey. A compact dusty star-forming galaxy at z=2.5

We present first results from the SXDF-ALMA 1.5 arcmin^2 deep survey at 1.1 mm using Atacama Large Millimeter Array (ALMA). The map reaches a 1sigma depth of 55 uJy/beam and covers 12 Halpha-selected star-forming galaxies at z = 2.19 or z=2.53. We have detected continuum emission from three of our Halpha-selected sample, including one compact star-forming galaxy with high stellar surface density, NB2315-07. They are all red in the rest-frame optical and have stellar masses of log (M*/Msun)>10.9 whereas the other blue, main-sequence galaxies with log(M*/Msun)=10.0-10.8 are exceedingly faint, <290 uJy (2sigma upper limit). We also find the 1.1 mm-brightest galaxy, NB2315-02, to be associated with a compact (R_e=0.7+-0.1 kpc), dusty star-forming component. Given high gas fraction (44^{+20}_{-8}% or 37^{+25}_{-3}%) and high star formation rate surface density (126^{+27}_{-30} Msun yr^{-1}kpc^{-2}), the concentrated starburst can within less than 50^{+12}_{-11} Myr build up a stellar surface density matching that of massive compact galaxies at z~2, provided at least 19+-3% of the total gas is converted into stars in the galaxy centre. On the other hand, NB2315-07, which already has such a high stellar surface density core, shows a gas fraction (23+-8%) and is located in the lower envelope of the star formation main-sequence. This compact less star-forming galaxy is likely to be in an intermediate phase between compact dusty star-forming and quiescent galaxies.

preprint2014arXiv

A Deep Narrowband Imaging Search for CIV and He II Emission from Ly$α$ Blobs

We conduct a deep narrow-band imaging survey of 13 Ly$α$ blobs (LABs) located in the SSA22 proto-cluster at z~3.1 in the CIV and HeII emission lines in an effort to constrain the physical process powering the Ly$α$ emission in LABs. Our observations probe down to unprecedented surface brightness limits of 2.1 $-$ 3.4 $\times$ 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ per 1 arcsec$^2$ aperture (5$σ$) for the HeII$λ$1640 and CIV$λ$1549 lines, respectively. We do not detect extended HeII and CIV emission in any of the LABs, placing strong upper limits on the HeII/Ly$α$ and CIV/Ly$α$ line ratios, of 0.11 and 0.16, for the brightest two LABs in the field. We conduct detailed photoionization modeling of the expected line ratios and find that, although our data constitute the deepest ever observations of these lines, they are still not deep enough to rule out a scenario where the Ly$α$ emission is powered by the ionizing luminosity of an obscured AGN. Our models can accommodate HeII/Ly$α$ and CIV/Ly$α$ ratios as low as $\simeq$0.05 and $\simeq$0.07 respectively, implying that one needs to reach surface brightness as low as 1 $-$ 1.5 $\times$ 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ (at 5$σ$) in order to rule out a photoionization scenario. These depths will be achievable with the new generation of image-slicing integral field units such as VLT/MUSE or Keck/KCWI. We also model the expected HeII/Ly$α$ and CIV/Ly$α$ in a different scenario, where Ly$α$ emission is powered by shocks generated in a large-scale superwind, but find that our observational constraints can only be met for shock velocities $v_{\rm s} \gtrsim$ 250 km s$^{-1}$, which appear to be in conflict with recent observations of quiescent kinematics in LABs.

preprint2014arXiv

Cold Molecular Gas in Merger Remnants: I. Formation of Molecular Gas Disks

We present < 1 kpc resolution CO imaging study of 37 optically-selected local merger remnants using new and archival interferometric maps obtained with ALMA, CARMA, SMA and PdBI. We supplement a sub-sample with single-dish measurements obtained at the NRO 45 m telescope for estimating the molecular gas mass (10^7 - 10^11 M_sun), and evaluating the missing flux of the interferometric measurements. Among the sources with robust CO detections, we find that 80 % (24/30) of the sample show kinematical signatures of rotating molecular gas disks (including nuclear rings) in their velocity fields, and the sizes of these disks vary significantly from 1.1 kpc to 9.3 kpc. The size of the molecular gas disks in 54 % of the sources is more compact than the K-band effective radius. These small gas disks may have formed from a past gas inflow that was triggered by a dynamical instability during a potential merging event. On the other hand, the rest (46 %) of the sources have gas disks which are extended relative to the stellar component, possibly forming a late-type galaxy with a central stellar bulge. Our new compilation of observational data suggests that nuclear and extended molecular gas disks are common in the final stages of mergers. This finding is consistent with recent major-merger simulations of gas rich progenitor disks. Finally, we suggest that some of the rotation-supported turbulent disks observed at high redshifts may result from galaxies that have experienced a recent major merger.

preprint2014arXiv

IGM Emission Observations with the Cosmic Web Imager: II. Discovery of Extended, Kinematically-Linked Emission around SSA22 Lyman-alpha Blob 2

The intergalactic medium (IGM) is the dominant reservoir of baryons, delineates the large scale structure at low to moderate overdensities, and provides gas from which galaxies form and evolve. Simulations of a Cold Dark Matter (CDM) dominated universe predict that the IGM is distributed in a cosmic web of filaments, and that galaxies should form along and at the intersections of these filaments (Bond et al. 1994; Miralda-Escude et al. 1996). While observations of QSO absorption lines and the large scale distribution of galaxies have confirmed the CDM paradigm, the cosmic web has never been confirmed by direct imaging. Here we report the Lyman-alpha blob 2 (LAB2) in SSA22, with the Cosmic Web Imager. This is an integral field spectrograph optimized for low surface brightness, extended emission. With 22 hours of total source exposure, CWI has revealed that LAB2 has extended Lyman-alpha emission which is consistent with filaments. We perform tests to secure the robustness of this result, which relies on data with modest signal-to-noise ratio. We have developed a smoothing algorithm that permits visualization of data cube slices along image or spectral-image planes. With both raw and smoothed data cubes we demonstrate that the filaments are kinematically associated with LAB2 and display double-peaked profiles of optically thick Lyman-alpha emission. The flux is 10 to 20x brighter than expected for the emission from the IGM but is consistent with boosted fluorescence from a buried QSO or gravitation cooling radiation. Using emission models we infer a baryon mass in the filaments of at least 1-4x10e11 Solar Mass, and the dark halo mass is at least 2x10e12 Solar Mass. The spatial kinematic morphology is more consistent with inflow from the cosmic web than outflow. LAB2 and the surrounding gas have significant and coaligned angular momentum, strengthening the case for their association.

preprint2014arXiv

Mapping the large scale structure around a z=1.46 galaxy cluster in 3-D using two adjacent narrow-band filters

We present a novel method to estimate accurate redshifts of star-forming galaxies by measuring the flux ratio of the same emission line observed through two adjacent narrow-band filters. We apply this method to our NB912 and new NB921 data taken with Suprime-Cam on the Subaru Telescope of a galaxy cluster, XMMXCS J2215.9-1738, at z=1.46 and its surrounding structures. We obtain redshifts for 170 [OII] emission line galaxies at z~1.46, among which 41 galaxies are spectroscopically confirmed with MOIRCS and FMOS on the Subaru mainly, showing an accuracy of $σ$((z-z_spec})/(1+z_spec))=0.002. This allows us to reveal filamentary structures that penetrate towards the centre of the galaxy cluster and intersect with other structures, consistent with the picture of hierarchical cluster formation. We also find that the projected celestial distribution does not precisely trace the real distribution of galaxies, indicating the importance of the three dimensional view of structures to properly identify and quantify galaxy environments. We investigate the environmental dependence of galaxy properties with local density, confirming that the median colour of galaxies becomes redder in higher density region while the star-formation rate of star-forming galaxies does not depend strongly on local environment in this structure. This implies that the star-forming activity in galaxies is truncated on a relatively short time scale in the cluster centre.

preprint2014arXiv

NIR spectroscopic observation of massive galaxies in the protocluster at z = 3.09

We present the results of near-infrared spectroscopic observations of the $K$-band selected candidate galaxies in the protocluster at $z=3.09$ in the SSA22 field. We observed 67 candidates with $K_{\rm AB}<24$ and confirmed redshifts of the 39 galaxies at $2.0< z_{\rm spec}< 3.4$. Of the 67 candidates, 24 are certainly protocluster members with $3.04\leq z_{\rm spec}\leq 3.12$, which are massive red galaxies those have been unidentified in previous optical observations of the SSA22 protocluster. Many distant red galaxies (DRGs; $J-K_{\rm AB}>1.4$), hyper extremely red objects (HEROs; $J-K_{\rm AB}>2.1$), {\it Spitzer} MIPS 24 $μ$m sources, active galactic nuclei (AGNs) as well as the counterparts of Ly$α$ blobs and the AzTEC/ASTE 1.1-mm sources in the SSA22 field are also found to be the protocluster members. The mass of the SSA22 protocluster is estimated to be $\sim2-5\times10^{14}~M_{\odot}$ and this system is plausibly a progenitor of the most massive clusters of galaxies in the current Universe. The reddest ($J-K_{\rm AB}\geq 2.4$) protocluster galaxies are massive galaxies with $M_{\rm star}\sim10^{11}~M_{\odot}$ showing quiescent star formation activities and plausibly dominated by old stellar populations. Most of these massive quiescent galaxies host moderately luminous AGNs detected by X-ray. There are no significant differences in the [O{\footnotesize III}] $λ$5007/H$β$ emission line ratios, and [O{\footnotesize III}] $λ$5007 line widths and spatial extents of the protocluster galaxies from those of massive galaxies at $z\sim2-3$ in the general field.

preprint2014arXiv

The environments of Ly$α$ blobs I: Wide-field Ly$α$ imaging of TN J1338-1942, a powerful radio galaxy at $z\simeq 4.1$ associated with a giant Ly$α$ nebula

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$α$ halo at z=4.11. We used a sample of Ly$α$ emitters (LAEs) down to $\log(L_{\rm Lyα} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $δ_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$α$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Lyα}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$α$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings.

preprint2013arXiv

A fundamental metallicity relation for galaxies at z = 0.84 - 1.47 from HiZELS

We obtained Subaru FMOS observations of Halpha emitting galaxies selected from the HiZELS narrow-band survey, to investigate the relationship between stellar mass, metallicity and star-formation rate at z = 0.84 - 1.47, for comparison with the Fundamental Metallicity Relation seen at low redshift. Our findings demonstrate, for the first time with a homogeneously selected sample, that a relationship exists for typical star-forming galaxies at z = 1 - 1.5 and that it is surprisingly similar to that seen locally. Therefore, star-forming galaxies at z = 1 - 1.5 are no less metal abundant than galaxies of similar mass and star formation rate (SFR) at z = 0.1, contrary to claims from some earlier studies. We conclude that the bulk of the metal enrichment for this star-forming galaxy population takes place in the 4 Gyr before z = 1.5. We fit a new mass-metallicity-SFR plane to our data which is consistent with other high redshift studies. However, there is some evidence that the mass-metallicity component of this high redshift plane is flattened, at all SFR, compared with z = 0.1, suggesting that processes such as star-formation driven winds, thought to remove enriched gas from low mass halos, are yet to have as large an impact at this early epoch. The negative slope of the SFR-metallicity relation from this new plane is consistent with the picture that the elevation in the SFR of typical galaxies at z > 1 is fuelled by the inflow of metal-poor gas and not major merging.

preprint2013arXiv

Obscured star formation in Ly-alpha blobs at z = 3.1

We present results from the AzTEC/ASTE 1.1-mm imaging survey of 35 Ly-alpha blobs (LABs) found in the SSA22 protocluster at z = 3.1. These 1.1-mm data reach an r.m.s. noise level of 0.7-1 mJy/beam, making this the largest millimetre-wave survey of LABs to date. No significant (> 3.5-sigma) emission is found in any of individual 35 LABs, and from this, we estimate 3-sigma upper limits on the far-infrared luminosity of L_FIR < 2e+12 Lsun. Stacking analysis reveals that the 1.1-mm flux density averaged over the LABs is S(1.1mm) < 0.40 mJy (3-sigma), which places a constraint of L_FIR < 4.5e+11 Lsun. This indicates that earlier 850-um measurements of the LABs may have overestimated their flux densities. Our results suggest that LABs on average have little ultra-luminous obscured star-formation, in contrast to a long-believed picture that LABs undergo an intense episode of dusty star-formation activities with star-formation rates of ~ 10^3 Msun/yr. Observations with ALMA are needed to directly study the obscured part of star-formation activity in the LABs.

preprint2013arXiv

The Formation of the massive galaxies in the SSA22 z=3.1 protocluster

The properties of K-band selected galaxies (K_AB<24) in the z = 3.09 SSA22 protocluster field are studied. 430 galaxies at 2.6 < z_phot < 3.6 are selected as potential protocluster members in a 112 arcmin^2 area based on their photometric redshifts. We find that \approx 20% of the massive galaxies with stellar masses >10^11 M_sun at z_phot \sim 3.1 have colors consistent with those of quiescent galaxies with ages > 0.5 Gyr. This fraction increases to \approx 50% after correcting for unrelated foreground/background objects. We also find that 30% of the massive galaxies are heavily reddened dusty star-forming galaxies. Few such quiescent galaxies at similar redshifts are seen in typical survey fields. An excess surface density of 24μm sources at z_phot \sim 3.1 is also observed, implying the presence of dusty star-formation activity in the protocluster. Cross-correlation with the X-ray data indicates that the fraction of K-band selected protocluster galaxies hosting active galactic nuclei (AGN) is also high compared with the field. The sky distribution of the quiescent galaxies, the 24μm sources, and the X-ray AGNs show clustering around a density peak of z=3.1 Lyαemitters (LAEs). A significant fraction of the massive galaxies have already become quiescent, while the dusty star-formation is still active in the SSA22 protocluster. These findings indicate that we are witnessing the formation epoch of massive early-type galaxies at the center of predecessors to present-day rich galaxy clusters.

preprint2012arXiv

A large Hα survey at z=2.23, 1.47, 0.84 & 0.40: the 11 Gyr evolution of star-forming galaxies from HiZELS

This paper presents new deep and wide narrow-band surveys undertaken with UKIRT, Subaru and the VLT; a combined effort to select large, robust samples of H-alpha (Ha) emitters at z=0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyrs) in a uniform manner over ~2 deg^2 in the COSMOS and UDS fields. The deep Ha surveys reach ~3M_sun/yr out to z=2.2 for the first time, while the wide area and the coverage over two independent fields allow to greatly overcome cosmic variance. A total of 1742, 637, 515 and 807 Ha emitters are homogeneously selected at z=0.40, 0.84, 1.47 and 2.23, respectively, and used to determine the Ha luminosity function and its evolution. The faint-end slope is found to be -1.60+-0.08 over z=0-2.23, showing no evolution. The characteristic luminosity of SF galaxies, L*, evolves significantly as log[L*(z)]=0.45z+log[L*(z=0)]. This is the first time Ha has been used to trace SF activity with a single homogeneous survey at z=0.4-2.23. Overall, the evolution seen in Ha is in good agreement with the evolution seen using inhomogeneous compilations of other tracers of star formation, such as FIR and UV, jointly pointing towards the bulk of the evolution in the last 11 Gyrs being driven by a strong luminosity increase from z~0 to z~2.2. Our uniform analysis reveals an Ha star formation history of the Universe which can be simply parameterised by log(SFRD)=-2.1/(1+z) for z<2.23. Both the shape and normalisation of the Ha star formation history are consistent with the measurements of the stellar mass density growth, confirming that our Ha analysis traces the bulk of the formation of stars in the Universe up to z~2.2. Star formation activity since z~2.2 is responsible for ~95% of the total stellar mass density observed locally today (with half being assembled between z=1.2-2.2, and the rest since z<1.2).

preprint2012arXiv

A Ly-alpha Emitter with an Extremely Large Rest-frame Equivalent Width of ~900A at z=6.5: A Candidate of Population III-dominated Galaxy?

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

preprint2012arXiv

Assembly of Massive Galaxies in a High-z Protocluster

We present the results of wide-field deep JHK imaging of the SSA22 field using MOIRCS instrument equipped with Subaru telescope. The observed field is 112 arcmin^2 in area, which covers the z=3.1 protocluster characterized by the overdensities of Ly Alpha emitters (LAEs) and Ly Alpha Blobs (LABs). The 5 sigma limiting magnitude is K_{AB} = 24.3. We extract the potential protocluster members from the K-selected sample by using the multi-band photometric-redshift selection as well as the simple color cut for distant red galaxies (DRGs; J-K_{AB}>1.4). The surface number density of DRGs in our observed fields shows clear excess compared with those in the blank fields, and the location of the densest area whose projected overdensity is twice the average coincides with the large-scale density peak of LAEs. We also found that K-band counterparts with z_{phot} = 3.1 are detected for 75% (15/20) of the LABs within their Ly Alpha halo, and the 40 % (8/20) of LABs have multiple components, which gives a direct evidence of the hierarchical multiple merging in galaxy formation. The stellar mass ofLABs correlates with their luminosity, isophotal area, and the Ly Alpha velocity widths, implying that the physical scale and the dynamical motion of Ly Alpha emission are closely related to their previous star-formation activities. Highly dust-obscured galaxies such as hyper extremely red objects (HEROs; J-K_{AB}>2.1) and plausible K-band counterparts of submillimeter sources are also populated in the high density region.

preprint2012arXiv

Characterization of the distribution of the Llyαemitters in the 53W002 field at z = 2.4

We present the results of our wide-field narrow band imaging of the field around the radio galaxy 53W002 at z = 2.390 with Subaru/Suprime-Cam. A custom made filter NB413 centered at 4140 Å with the width of 83 Å is used to observe the 31' x 24' area around the radio galaxy. We detected 204 Lyαemitters (LAEs) at z = 2.4 with a rest frame equivalent width larger than 25 Å to the depth of 26 AB mag (in NB413). The entire LAE population in the 53W002 field has an average number density and distributions of equivalent width and size that are similar to those of other fields at z ~ 2. We identify a significant high density region (53W002F-HDR) that spreads over ~ 5' x 4' near 53W002 where the LAE number density is nearly four times as large as the average of the entire field. Using the probability distribution function of density fluctuation, we evaluate the rareness probability of the 53W002F-HDR to be 0.9^{+2.4}_{-0.62}%, which corresponds to a moderately rich structure. No notable environmental dependency at the comoving scale of 10 Mpc is found for the distributions of the Lyαequivalent width and luminosity in the field. We also detected 4 Lyαblobs (LABs), one of which is newly discovered. They are all found to be located in the rims of high density regions. The biased location and unique morphologies in Lyαsuggest that galaxy interaction play a key role in their formation.

preprint2012arXiv

Resolving the Galaxies within a Giant Lya Nebula: Witnessing the Formation of a Galaxy Group?

Detailed analysis of the substructure of Lya nebulae can put important constraints on the physical mechanisms at work and the properties of galaxies forming within them. Using high resolution HST imaging of a Lya nebula at z~2.656, we have taken a census of the compact galaxies in the vicinity, used optical/near-infrared colors to select system members, and put constraints on the morphology of the spatially-extended emission. The system is characterized by (a) a population of compact, low luminosity (~0.1 L*) sources --- 17 primarily young, small (Re~1-2 kpc), disky galaxies including an obscured AGN --- that are all substantially offset (>20 kpc) from the line-emitting nebula; (b) the lack of a central galaxy at or near the peak of the Lya emission; and (c) several nearly coincident, spatially extended emission components --- Lya, HeII, and UV continuum --- that are extremely smooth. These morphological findings are difficult to reconcile with theoretical models that invoke outflows, cold flows, or resonant scattering, suggesting that while all of these physical phenomena may be occurring, they are not sufficient to explain the powering and large extent of Lya nebulae. In addition, although the compact galaxies within the system are irrelevant as power sources, the region is significantly overdense relative to the field galaxy population (by at least a factor of 4). These observations provide the first estimate of the luminosity function of galaxies within an individual Lya nebula system, and suggest that large Lya nebulae may be the seeds of galaxy groups or low-mass clusters.

preprint2011arXiv

Completing the Census of Ly-alpha Emitters at the Reionization Epoch

We carried out extended spectroscopic confirmations of Ly-alpha emitters (LAEs) at z=6.5 and 5.7 in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 at z=6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z=6.5 (5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Ly-alpha luminosity, both photometrically and spectroscopically, to accurately determine the Ly-alpha and rest-UV luminosity functions (LFs). The substantially improved evaluation of the Ly-alpha LF at z=6.5 shows an apparent deficit from z=5.7 at least at the bright end, and a possible decline even at the faint end, though small uncertainties remain. The rest-UV LFs at z=6.5 and 5.7 are in good agreement, at least at the bright end, in clear contrast to the differences seen in the Ly-alpha LF. These results imply an increase in the neutral fraction of the intergalactic medium from z=5.7 to 6.5. The rest-frame equivalent width (EW_0) distribution at z=6.5 seems to be systematically smaller than z=5.7, and it shows an extended tail toward larger EW_0. The bright end of the rest-UV LF can be reproduced from the observed Ly-alpha LF and a reasonable EW_0-UV luminosity relation. Integrating this rest-UV LF provides the first measurement of the contribution of LAEs to the photon budget required for reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon budget among Lyman break galaxies significantly increases towards faint magnitudes. Low-luminosity LAEs could dominate the ionizing photon budget, though this inference depends strongly on the uncertain faint-end slope of the Ly-alpha LF.

preprint2011arXiv

Constraining Dust and Molecular Gas Properties in Lyman Alpha Blobs at z~3

In order to constrain the bolometric luminosities, dust properties and molecular gas content of giant Lyman alpha nebulae, the so-called Lyman alpha blobs, we have carried out a study of dust continuum and CO line emission in two well-studied representatives of this population at z ~ 3: a Lya blob discovered by its strong Spitzer MIPS 24um detection (LABd05; Dey et al. 2005) and the Steidel blob 1 (SSA22-LAB01; Steidel et al. 2000). We find that the spectral energy distribution of LABd05 is well described by an AGN-starburst composite template with L(FIR) = (4.0 +/- 0.5) x 10^12 Lsun, comparable to high-z sub-millimeter galaxies and ultraluminous infrared galaxies. New APEX/LABOCA 870um measurements rule out the reported SCUBA detection of the SSA22-LAB01 (S[850um] = 16.8 mJy) at the > 4sigma level. Consistent with this, ultra-deep Plateau de Bure Interferometer (PdBI) observations with ~2arcsec spatial resolution also fail to detect any 1.2mm continuum source down to ~0.45mJy per beam (3sigma). Combined with the existing (sub)mm observations in the literature, we conclude that the FIR luminosity of SSA22-LAB01 remains uncertain. No CO line is detected in either case down to integrated flux limits of (Snu dV) < 0.25--1.0 Jy km/s, indicating a modest molecular gas reservoir, M(H_2) < 1--3 x 10^10 Msun. The non-detections exclude, with high significance (12 sigma), the previous tentative detection of a CO(4-3) line in the SSA22-LAB01. The increased sensitivity afforded by ALMA will be critical in studying molecular gas and dust in these interesting systems.

preprint2011arXiv

Keck Spectroscopy of Lyman-break Galaxies and its Implications for the UV-continuum and Ly_alpha Luminosity Functions at z>6

We present Keck spectroscopic observations of z>6 Lyman-break galaxy (LBG) candidates in the Subaru Deep Field (SDF). The candidates were selected as i'-dropout objects down to z'=27 AB magnitudes from an ultra-deep SDF z'-band image. With the Keck spectroscopy we identified 19 LBGs with prominent Ly_alpha emission lines at 6< z < 6.4. The median value of the Ly_alpha rest-frame equivalent widths (EWs) is ~50 A, with four EWs >100 A. This well-defined spectroscopic sample spans a UV-continuum luminosity range of -21.8< M_{UV}<-19.5 (0.6~5 L*_{UV}) and a Ly_alpha luminosity range of (0.3~3) x 10^{43} erg s^{-1} (0.3~3 L*_ {Ly_alpha}). We derive the UV and Ly_alpha luminosity functions (LFs) from our sample at <z>~6.2 after we correct for sample incompleteness. We find that our measurement of the UV LF is consistent with the results of previous studies based on photometric LBG samples at 5<z<7. Our Ly_alpha LF is also generally in agreement with the results of Ly_alpha-emitter surveys at z~5.7 and 6.6. This study shows that deep spectroscopic observations of LBGs can provide unique constraints on both the UV and Ly_alpha LFs at z>6.

preprint2011arXiv

Star formation at z=1.47 from HiZELS: An Hα+[OII] double-blind study

This paper presents the results from the first wide and deep dual narrow-band survey to select H-alpha (Ha) and [OII] line emitters at z=1.47+-0.02 (using matched narrow-band filters in the H and z' bands), exploiting synergies between the UKIRT and Subaru telescopes. The Ha survey at z=1.47 reaches a flux limit of ~7x10^-17 erg/s/cm^2 and detects ~200 Ha emitters over 0.7deg^2, while the much deeper [OII] survey reaches an effective flux of ~7x10^-18 erg/s/cm^2, detecting ~1400 z=1.47 [OII] emitters in a matched co-moving volume of ~2.5x10^5 Mpc^3. The combined survey results in the identification of 190 simultaneous Ha and [OII] emitters at z=1.47. Ha and [OII] luminosity functions are derived and both are shown to evolve significantly from z~0 in a consistent way. The star formation rate density of the Universe at z=1.47 is evaluated, with the Ha analysis yielding 0.16+-0.05 M_sun/yr/Mpc^3 and the [OII] analysis 0.17+-0.04 M_sun/yr/Mpc^3. The measurements are combined with other studies, providing a self-consistent measurement of the star formation history of the Universe over the last ~11Gyrs. By using a large comparison sample at z~0.1 (from the SDSS), [OII]/Ha line ratios are calibrated as probes of dust-extinction. Ha emitters at z~1.47 show on average 1 mag of extinction at Ha, similar to the SDSS sources at z~0. Although we find that dust extinction correlates with SFR, the relation evolves by about ~0.5 mag from z~1.5 to z~0, with z~0 relations over-predicting the dust extinction corrections at high-z by that amount. Stellar mass is found to be a much more fundamental extinction predictor, with the relation between mass and extinction being valid at both z~0 and z~1.5. Dust extinction corrections as a function of optical colours are also derived, offering simpler mechanisms for estimating extinction in moderately star-forming systems over the last ~9Gyrs [Abridged].

preprint2010arXiv

Discovery of an Excess of Halpha Emitters around 4C 23.56 at z=2.48

We report the discovery of a significant excess of candidate Halpha emitters (HAEs) in the field of the radio galaxy 4C 23.56 at z=2.483. Using the MOIRCS near-infrared imager on the Subaru Telescope we found 11 candidate emission-line galaxies to a flux limit of ~7.5 10^-17 erg s-1 cm-2, which is about 5 times excess from the expected field counts with ~3-sigma significance. Three of these are spectroscopically confirmed as redshifted Halpha at z=2.49. The distribution of candidate emitters on the sky is tightly confined to a 1.2-Mpc-radius area at z=2.49, locating 4C 23.56 at the western edge of the distribution. Analysis of the deep Spitzer MIPS 24 mu m imaging shows that there is also an excess of faint MIPS sources. All but two of the 11 HAEs are also found in the MIPS data. The inferred star-formation rate (SFR) of the HAEs based on the extinction-corrected Halpha luminosity (median SFR >~100 M_solar yr-1) is similar to those of HAEs in random fields at z~2. On the other hand, the MIPS-based SFR for the HAEs is on average 3.6 times larger, suggesting the existence of the star-formation significanly obscured by dust. The comparison of the Halpha-based star-formation activities of the HAEs in the 4C 23.56 field to those in another proto-cluster around PKS 1138-262 at z=2.16 reveals that the latter tend to have fainter Halpha emission despite similar K-band magnitudes. This suggests that star-formation may be suppressed in the PKS 1138-262 protocluster relative to the 4C 23.56 protocluster. This difference among the HAEs in the two proto-clusters at z > 2 may imply that some massive cluster galaxies are just forming at these epochs with some variation among clusters.

preprint2008arXiv

Subaru/MOIRCS Near-Infrared Imaging in the Proto-Cluster Region at z=3.1

We present the results of deep near-infrared imaging observations of the z=3.1 proto-cluster region in the SSA22a field taken by MOIRCS mounted on the Subaru Telescope. We observed a 21.7 arcmin^2 field to the depths of J=24.5, H=24.3, and K=23.9 (5 sigma). We examine the distribution of the K-selected galaxies at z~3 by using the simple color cut for distant red galaxies (DRGs) as well as the photometric-redshift selection technique. The marginal density excess of DRGs and the photo-z selected objects are found around the two most luminous Ly alpha blobs (LABs). We investigate the correlation between the K-selected objects and the LABs, and find that several galaxies with stellar mass M_* = 10^9-10^11 M_solar exist in vicinity of LABs, especially around the two most luminous ones. We also find that 7 of the 8 LABs in the field have plausible K_s-band counterparts and the sum of the stellar mass possibly associated with LABs correlates with the luminosity and surface brightness of them, which implies that the origin of Ly alpha emission may be closely correlated with their stellar mass or their previous star formation phenomena.

preprint2008arXiv

The Reionization and Galaxy Evolution Probed by z=7 Lyman Alpha Emitters

We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7.