Source author record

Nobunari Kashikawa

Nobunari Kashikawa 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

82works
5topics
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

82 published item(s)

preprint2023arXiv

MAMMOTH-Subaru III. Ly$α$ Halo Extended to $\sim200$ kpc Identified by Stacking $\sim 3300$ Ly$α$ Emitters at $z=2.2-2.3$

In this paper, we present a Ly$α$ halo extended to $\sim200$ kpc identified by stacking $\sim 3300$ Ly$α$ emitters at $z=2.2-2.3$. We carry out imaging observations and data reduction with Subaru/Hyper Suprime-Cam (HSC). Our total survey area is $\sim12$ deg$^2$ and imaging depths are $25.5-27.0$ mag. Using the imaging data, we select 1,240 and 2,101 LAE candidates at $z=2.2$ and 2.3, respectively. We carry out spectroscopic observations of our LAE candidates and data reduction with Magellan/IMACS to estimate the contamination rate of our LAE candidates. We find that the contamination rate of our sample is low (8%). We stack our LAE candidates with a median stacking method to identify the Ly$α$ halo at $z=2$. We show that the Ly$α$ halo is extended to $\sim200$ kpc at a surface brightness level of $10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Comparing to previous studies, our Ly$α$ halo is more extended at radii of $\sim25-100$ kpc, which is not likely caused by the contamination in our sample but by different redshifts and fields instead. To investigate how central galaxies affect surrounding LAHs, we divide our LAEs into subsamples based on the Ly$α$ luminosity ($L_{\rm Lyα}$), rest-frame Ly$α$ equivalent width (EW$_0$), and UV magnitude (M$_{\rm uv}$). We stack the subsamples and find that higher $L_{\rm Lyα}$, lower EW$_0$, and brighter M$_{\rm uv}$ cause more extended halos. Our results suggest that more massive LAEs generally have more extended Ly$α$ halos.

preprint2023arXiv

MAMMOTH-Subaru V. Effects of Cosmic Variance on Ly$α$ Luminosity Functions at $z=2.2-2.3$

Cosmic variance introduces significant uncertainties into galaxy number density properties when surveying the high-z Universe with a small volume, such uncertainties produce the field-to-field variance of galaxy number $σ_{g}$ in observational astronomy. This uncertainty significantly affects the Luminosity Functions (LF) measurement of Lya Emitters (LAEs). For most previous Lya LF studies, $σ_{g}$ is often adopted from predictions by cosmological simulations, but barely confirmed by observations. Measuring cosmic variance requires a huge sample over a large volume, exceeding the capabilities of most astronomical instruments. In this study, we demonstrate an observational approach for measuring the cosmic variance contribution for $z\approx2.2$ Lya LFs. The LAE candidates are observed using narrowband and broadband of the Subaru/Hyper Suprime-Cam (HSC), with 8 independent fields, making the total survey area $\simeq11.62$deg$^2$ and a comoving volume of $\simeq8.71\times10^6$Mpc$^3$. These eight fields are selected using the project of MAMMOTH. We report a best-fit Schechter function with parameters $α=-1.75$ (fixed), $L_{Lyα}^{*}=5.18_{-0.40}^{+0.43} \times 10^{42}$erg s$^{-1}$ and $ϕ_{Lya}^{*}=4.87_{-0.55}^{+0.54}\times10^{-4}$Mpc$^{-3}$ for the overall Lya LFs. After clipping out the regions that can bias the cosmic variance measurements, we calculate $σ_{g}$, by sampling LAEs within multiple pointings assigned on the field image. We investigate the relation between $σ_{g}$ and survey volume $V$, and fit a simple power law: $σ_g=k\times(\frac{V_{\rm eff}}{10^5 {\rm Mpc}^3})^β$. We find best-fit values of $-1.209_{-0.106}^{+0.106}$ for $β$ and $0.986_{-0.100}^{+0.108}$ for k. We compare our measurements with predictions from simulations and find that the cosmic variance of LAEs might be larger than that of general star-forming galaxies.

preprint2022arXiv

CHORUS IV: Mapping the Spatially Inhomogeneous Cosmic Reionization with Subaru HSC

The spatial inhomogeneity is one of the important features for understanding the reionization process; however, it has not yet been fully quantified. To map this inhomogeneous distribution, we simultaneously detect Ly$α$ emitters (LAEs) and Lyman break galaxies (LBGs) at $z \sim 6.6$ from the Subaru/Hyper Suprime-Cam (HSC) large-area ($\sim1.5\,\mathrm{ deg}^2 = 34000\,\mathrm{cMpc}^2$) deep survey. We estimate the neutral fraction, $x_\mathrm{HI}$, from the observed number density ratio of LAEs to LBGs, $n(\mathrm{LAE})/n(\mathrm{LBG})$ based on numerical radiative transfer simulation, in which model galaxies are selected to satisfy the observed selection function. While the average $x_\mathrm{HI}$ within the field of view is found to be $x_\mathrm{HI} < 0.4$, which is consistent with previous studies, the variation of $n(\mathrm{LAE})/n(\mathrm{LBG})$ within the field of view for each $140\,\mathrm{pMpc}^2$ is found to be as large as a factor of three. This may suggest a spatially inhomogeneous topology of reionization, but it also leaves open the possibility that the variation is based on the inherent large-scale structure of the galaxy distribution. Based on the simulations, it may be difficult to distinguish between the two from the current survey. We also find that LAEs in the high LAE density region are more populate high $\mathrm{EW}_0$, supporting that the observed $n(\mathrm{LAE})/n(\mathrm{LBG})$ is more or less driven by the neutral fraction, though the statistical significance is not high.

preprint2022arXiv

Exploratory Study of Transverse Proximity Effect around BAL Quasars

We aim to find out the reason why there exists an anisotropic HI absorption around quasars; i.e., the environments around quasars are highly biased toward producing strong HI absorption in transverse direction while there exists a significant deficit of HI absorption within a few Mpc of quasars along line-of-sight. The most plausible explanation for this opposite trend is that the transverse direction is shadowed from the quasar UV radiation due to dust torus. However, a critical weakness of this idea is that we have no information on inclination angle of our sightline relative to the torus. In this study, we examine environments of quasars with broad absorption troughs in their spectra (i.e., BAL quasars) because it is widely believed that BAL troughs are observed if the central continuum is viewed from the side through their powerful outflows near the dust torus. With closely separated 12 projected quasar pairs at different redshift with separation angle of $θ$$<$120$^{\prime\prime}$, we examine HI absorption at foreground BAL quasars in spectra of background quasars. We confirm there exist optically thick gas around two of 12 BAL quasars, and that the mean HI absorption strength is EW$_{\rm rest}$$\sim$1A. These are consistent to the past results around non-BAL quasars, although not statistically significant. However, the origins of optically thick HI absorbers around BAL and non-BAL quasars could be different since their column densities are different by $\sim$3 orders of magnitude. The larger sample would be required for narrowing down possible scenarios for the anisotropic HI absorption around quasars.

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

The physical origin for spatially large scatter of IGM opacity at the end of reionization: the IGM Ly$α$ opacity-galaxy density relation

The large opacity fluctuations in the $z > 5.5$ Ly$α$ forest may indicate inhomogeneous progress of reionization. To explain the observed large scatter of the effective Ly$α$ optical depth ($τ_{\rm eff}$) of the intergalactic medium (IGM), fluctuation of UV background ($Γ$ model) or the IGM gas temperature ($T$ model) have been proposed, which predict opposite correlations between $τ_{\rm eff}$ and galaxy density. In order to address which model can explain the large scatter of $τ_{\rm eff}$, we search for Ly$α$ emitters (LAEs) around two (J1137+3549 and J1602+4228) quasar sightlines with $τ_{\rm eff}\sim3$ and J1630+4012 sightline with $τ_{\rm eff}\sim5.5$. Using a narrowband imaging with Subaru/Hyper Suprime-Cam, we draw LAE density maps to explore their spatial distributions. Overdensities are found within 20 $h^{-1}$Mpc of the quasar sightlines in the low $τ_{\rm eff}$ regions, while a deficit of LAEs is found in the high $τ_{\rm eff}$ region. Although the $τ_{\rm eff}$ of the three quasar sightlines are neither high nor low enough to clearly distinguish the two models, these observed $τ_{\rm eff}$-galaxy density relations all consistently support the $Γ$ model rather than the $T$ model in the three fields, along with the previous studies. The observed overdensities near the low $τ_{\rm eff}$ sightlines may suggest that the relic temperature fluctuation does not affect reionization that much. Otherwise, these overdensities could be attributed to other factors besides the reionization process, such as the nature of LAEs as poor tracers of underlying large-scale structures.

preprint2022arXiv

Where's Swimmy?: Mining unique color features buried in galaxies by deep anomaly detection using Subaru Hyper Suprime-Cam data

We present the Swimmy (Subaru WIde-field Machine-learning anoMalY) survey program, a deep-learning-based search for unique sources using multicolored ($grizy$) imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). This program aims to detect unexpected, novel, and rare populations and phenomena, by utilizing the deep imaging data acquired from the wide-field coverage of the HSC-SSP. This article, as the first paper in the Swimmy series, describes an anomaly detection technique to select unique populations as "outliers" from the data-set. The model was tested with known extreme emission-line galaxies (XELGs) and quasars, which consequently confirmed that the proposed method successfully selected 60-70% of the quasars and 60% of the XELGs without labeled training data. In reference to the spectral information of local galaxies at $z=$0.05-0.2 obtained from the Sloan Digital Sky Survey, we investigated the physical properties of the selected anomalies and compared them based on the significance of their outlier values. The results revealed that XELGs constitute notable fractions of the most anomalous galaxies, and certain galaxies manifest unique morphological features. In summary, a deep anomaly detection is an effective tool that can search rare objects, and ultimately, unknown unknowns with large data-sets. Further development of the proposed model and selection process can promote the practical applications required to achieve specific scientific goals.

preprint2021arXiv

"A more probable explanation" is still impossible to explain GN-z11-flash: in response to Steinhardt et al. (arXiv:2101.12738)

In Jiang et al. (2020), we reported a possible bright flash (hereafter GN-z11-flash) from a galaxy GN-z11 at z ~ 11. Recently, Steinhardt et al. (2021; arXiv:2101.12738) found 27 images with transient signals in Keck MOSFIRE archival data and claimed that GN-z11-flash was more likely from a moving object in our Solar system. We show that the Steinhardt et al.'s definition of the chance probability and their methodology of finding GN-z11-flash-like transients are problematic in several aspects. In particular, none of their transients is analogous to GN-z11-flash, and none of them is positionally coincident with a known object in their imaging data. In Jiang et al., we performed a comprehensive analysis of the origin of GN-z11-flash and ruled out, to the best of our knowledge, the possibility of known man-made objects or moving objects in the Solar system, based on all available information and our current understanding of these objects. Steinhardt et al. did not use such information and did not analyse the GN-z11-flash event itself. The majority of their transients are apparently low-Earth orbit satellites or aircrafts. Therefore, their analysis can neither prove nor disprove our results. Finally, we present a method to estimate the chance probability of finding GN-z11-flash-like transients in archival data. Based on this method and the archival data used by Steinhardt et al., we obtain a loose upper limit of the probability that actually support the original results of Jiang et al. (2020).

preprint2021arXiv

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VI. Distant Filamentary Structures Pointed by High-z Radio Galaxies at z~4

We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galaxy overdensities around 21 HzRGs whose rest-frame 1.4 GHz radio luminosities ($L_{1.4\mathrm{GHz}}$) are $10^{26-27}$ W Hz$^{-1}$. We find that the overdensities around the faint HzRGs with $L_{1.4\mathrm{GHz}}\sim10^{26.0-26.5}$ W Hz$^{-1}$ tend to be higher than that of the $g$-dropout galaxies. On the other hand, no significant difference of density environments is found between the luminous HzRGs with $ L_{1.4\mathrm{GHz}}\sim10^{26.5-27.0} $ W Hz$^{-1}$ and the $g$-dropout galaxies. The HzRGs are found to occupy more massive halos than $g$-dropout galaxies through a cross-correlation between the HzRGs and $g$-dropout galaxies. This trend is more pronounced in the faint HzRGs. These results are consistent with a scenario where HzRGs get older and more massive as the radio-luminosity decreases. The HzRGs are expected to trace the progenitors of local cluster halos from their calculated halo mass. In addition, we find that surrounding galaxies tend to distribute along the radio-jet major axis of the HzRGs at angular distances less than $\lesssim500$ physical kpc. Our findings imply the onset of the filamentary structures around the HzRGs at $z\sim4$.

preprint2021arXiv

GOLDRUSH. IV. Luminosity Functions and Clustering Revealed with ~4,000,000 Galaxies at z~2-7: Galaxy-AGN Transition, Star Formation Efficiency, and Implication for Evolution at z>10

We present new measurements of rest-UV luminosity functions and angular correlation functions from 4,100,221 galaxies at z~2-7 identified in the Subaru/Hyper Suprime-Cam survey and CFHT Large-Area U-band Survey. The obtained luminosity functions at z~4-7 cover a very wide UV luminosity range of ~0.002-2000L*uv combined with previous studies, revealing that the dropout luminosity function is a superposition of the AGN luminosity function dominant at Muv<-24 mag and the galaxy luminosity function dominant at Muv>-22 mag, consistent with galaxy fractions based on 1037 spectroscopically-identified sources. Galaxy luminosity functions estimated from the spectroscopic galaxy fractions show the bright end excess beyond the Schechter function at >2sigma levels, which is possibly made by inefficient mass quenching, low dust obscuration, and/or hidden AGN activity. By analyzing the correlation functions at z~2-6 with halo occupation distribution models, we find a weak redshift evolution (within 0.3 dex) of the ratio of the star formation rate (SFR) to the dark matter accretion rate, SFR/(dMh/dt), indicating the almost constant star formation efficiency at z~2-6, as suggested by our earlier work at z~4-7. Meanwhile, the ratio gradually increases with decreasing redshift at z<5 within 0.3 dex, which quantitatively reproduces the redshift evolution of the cosmic SFR density, suggesting that the evolution is primarily driven by the increase of the halo number density due to the structure formation, and the decrease of the accretion rate due to the cosmic expansion. Extrapolating this calculation to higher redshifts assuming the constant efficiency suggests a rapid decrease of the SFR density at z>10 with $\propto10^{-0.5(1+z)}$, which will be directly tested with JWST.

preprint2021arXiv

SILVERRUSH. XI. Intensity Mapping for Lya Emission Extending over $100-1000$ comoving kpc around $z\sim2-7$ LAEs with Subaru HSC-SSP and CHORUS Data

We conduct intensity mapping to probe for extended diffuse Ly$α$ emission around Ly$α$ emitters (LAEs) at $z\sim2-7$, exploiting very deep ($\sim26$ mag at $5σ$) and large-area ($\sim4.5$ deg$^2$) Subaru/Hyper Suprime-Cam narrow-band (NB) images and large LAE catalogs consisting of a total of 1781 LAEs at $z=2.2$, $3.3$, $5.7$, and $6.6$ obtained by the HSC-SSP SILVERRUSH and CHORUS projects. We calculate the spatial correlations of these LAEs with $\sim1-2$ billion pixel flux values of the NB images, deriving the average Ly$α$ surface brightness (${\rm SB_{Lyα}}$) radial profiles around the LAEs. By carefully estimating systematics such as fluctuations of sky background and point spread functions, we detect diffuse Ly$α$ emission ($\sim10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$) at $100-1000$ comoving kpc around $z=3.3$ LAEs at the $4.1σ$ level and tentatively ($\sim2σ$) at the other redshifts, beyond the virial radius of a dark-matter halo with a mass of $10^{11}\ M_\odot$. While the observed ${\rm SB_{Lyα}}$ profiles have similar amplitudes at $z=2.2-6.6$ within the uncertainties, the intrinsic ${\rm SB_{Lyα}}$ profiles (corrected for the cosmological dimming effect) increase toward high redshifts. This trend may be explained by increasing hydrogen gas density due to the evolution of the cosmic volume. Comparisons with theoretical models suggest that extended Ly$α$ emission around a LAE is powered by resonantly scattered Ly$α$ photons in the CGM and IGM that originates from the inner part of the LAE, and/or neighboring galaxies around the LAE.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) XII. Extended [C II] Structure (Merger or Outflow) in a z = 6.72 Red Quasar

We present ALMA [C II] 158 $μ$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$ with the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5σ} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag.

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

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

Evidence for a highly opaque large-scale galaxy void at the end of reionization

We present evidence that a region of high effective Ly$α$ optical depth at $z\sim5.7$ is associated with an underdense region at the tail end of cosmic reionization. We carried out a survey of Lyman-break Galaxies (LBGs) using Subaru Hyper Suprime-Cam in the field of the $z=5.98$ quasar J0148+0600, whose spectrum presents an unusually long ($\sim160 \mathrm{cMpc}$) and opaque ($τ\gtrsim7$) Ly$α$ trough at $5.5\le z\le 5.9$. LBG candidates were selected to lie within the redshift range of the trough, and the projected number densities were measured within 90~cMpc of the quasar sightline. The region within $8'$ (or $\approx 19~\mathrm{cMpc}$) of the quasar position is the most underdense of the whole field. The significance of the presence of the void is estimated to be $99\%$. This is consistent with the significant deficit of Ly$α$ emitters (LAEs) at $z=5.72$ reported by Becker et al. and suggests that the paucity of LAEs is not purely due to the removal of the Ly$α$ emission by the high opacity but reflects a real coherent underdensity of galaxies across the entire redshift range of the trough. These observations are consistent with scenarios in which large optical depth fluctuations arise due to fluctuations in the galaxy-dominant UV background or due to residual neutral islands that are expected from reionization that is completed at redshifts as low as $z\lesssim5.5$.

preprint2020arXiv

Faint Quasars Live in the Same Number Density Environments as Lyman Break Galaxies at z ~ 4

Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies around the quasars with those around g-dropout galaxies, and find no significant difference between their distributions. A total of 4 (22) out of the 570 faint quasars, 0.7_{-0.4}^{+0.4} (3.9_{-0.8}^{+0.8}) %, are found to be associated with the > 4 sigma overdense regions within an angular separation of 1.8 (3.0) arcmin, which is the typical size of protoclusters at this epoch. This is similar to the fraction of g-dropout galaxies associated with the > 4 sigma overdense regions. This result is consistent with our previous work that 1.3_{-0.9}^{+0.9} % and 2.0_{-1.1}^{+1.1} % of luminous quasars detected in the Sloan Digital Sky Survey exist in the > 4 sigma overdense regions within 1.8 and 3.0 arcmin separations, respectively. Therefore, we suggest that the galaxy number densities around quasars are independent of their luminosity, and most quasars do not preferentially appear in the richest protocluster regions at z ~ 4. The lack of an apparent positive correlation between the quasars and the protoclusters implies that: i) the gas-rich major merger rate is relatively low in the protocluster regions, ii) most high-z quasars may appear through secular processes, or iii) some dust-obscured quasars exist in the protocluster regions.

preprint2020arXiv

Possible evolution of the circum-galactic medium around QSOs with QSO age and cosmic time revealed by Ly$α$ halos

We first present new Subaru narrow-band observations of the Ly$α$ halo around the quasi-stellar object (QSO) CFHQ J232908$-$030158 at $z=6.42$, which appears the most luminous and extended halo at $z>5$ ($L_{Lyα}=9.8\times10^{43}$ erg s$^{-1}$ within $37$ pkpc diameter). Then, combining these measurements with available data in the literature, we find two different evolutions of QSOs' Ly$α$ halos. First is a possible short-term evolution with QSO age seen in four $z>6$ QSOs. We find the anti-correlation between the Ly$α$ halo scales with QSOs' IR luminosity, with J2329-0301's halo being the brightest and largest. It indicates that ionizing photons escape more easily out to circum-galactic regions when host galaxies are less dusty. We also find a positive correlation between IR luminosity and black hole mass ($M_\text{BH}$). Given $M_\text{BH}$ as an indicator of QSO age, we propose a hypothesis that a large Ly$α$ halo mainly exists around QSOs in the young phase of their activity due to a small amount of dust. The second is an evolution with cosmic time seen over $z\sim2-5$. We find the increase of surface brightness toward lower-redshift with a similar growth rate to that of dark matter halos (DHs) which evolve to $M_\text{DH}=10^{12}-10^{13}$ M$_\odot$ at $z=2$. The extent of Ly$α$ halos is also found to increase at a rate scaling with the virial radius of growing DHs, $r_\text{vir} \propto M_\text{DH}^{1/3}$($1+z$)$^{-1}$. These increases are consistent with a scenario that the CGM around QSOs evolves in mass and size keeping pace with hosting DHs.

preprint2020arXiv

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$.

preprint2020arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). IX. Identification of Two Red Quasars at z > 5.6

We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z >5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the spectroscopically confirmed 93 high-z quasars in the SHELLQs sample, based on detections in the Wide-field Infrared Survey Explorer (WISE) data at 3.4 and 4.6 um (rest-frame ~5000-6500 A). The amount of dust reddening was estimated with spectral energy distribution (SED) fits over optical and mid-infrared wavelengths. Two of the four candidates were found to be red quasars with dust reddening of E(B-V) > 0.1. The remaining SHELLQs quasars without individual WISE detections are significantly fainter in the WISE bands and bluer than the red quasars, although we did detect them in the W1 band in a stacked image. We also conducted the same SED fits for high-z optically-luminous quasars, but no red quasar was found. This demonstrates the power of Subaru HSC to discover high-z red quasars, which are fainter than the limiting magnitudes of past surveys in the rest-frame ultraviolet, due to dust extinction.

preprint2020arXiv

The UV Luminosity Function of Protocluster Galaxies at $z\sim4$: the Bright-end Excess and the Enhanced Star Formation Rate Density

We report the rest-frame ultraviolet luminosity function of $g$-dropout galaxies in 177 protocluster candidates (PC UVLF) at $z\sim4$ selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess towards the bright-end. This excess can not be explained by the contribution of only active galactic nuclei, and we also find that this is more significant in higher dense regions. Assuming that all protocluster members are located on the star formation main sequence, the PC UVLF can be converted into a stellar mass function. Consequently, our protocluster members are inferred to have a 2.8 times more massive characteristic stellar mass than that of the field Lyman break galaxies at the same redshift. This study, for the first time, clearly shows that the enhancement in star formation or stellar mass in overdense regions can generally be seen as early as at $z\sim4$. We also estimate the star formation rate density (SFRD) in protocluster regions as $\simeq 6-20\%$ of the cosmic SFRD, based on the measured PC UVLF after correcting for the selection incompleteness in our protocluster sample. This high value suggests that protoclusters make a non-negligible contribution to the cosmic SFRD at $z\sim4$, as previously suggested by simulations. Our results suggest that protoclusters are essential components for the galaxy evolution at $z\sim4$.

preprint2019arXiv

CHORUS. III. Photometric and Spectroscopic Properties of Ly$α$ Blobs at $z=4.9-7.0$

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$α$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$α$ emitters with a bright ($\log L_{\rm Lyα}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$α$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$α$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$λ$1240 emission, or Ly$α$ line broadening, while z49-1 has a strong {\sc Civ}$λ$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$α$ emission of the LABs. The Ly$α$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$α$ luminosity relation of the Ly$α$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

preprint2019arXiv

Discovery of Protoclusters at z~3.7 & 4.9: Embedded in Primordial Superclusters

We have carried out follow-up spectroscopy on three overdense regions of $g$- and $r$-dropout galaxies in the Canada-France-Hawaii Telescope Legacy Survey Deep Fields, finding two new protoclusters at $z=4.898$, 3.721 and a possible protocluster at $z=3.834$. The $z=3.721$ protocluster overlaps with a previously identified protocluster at $z=3.675$. The redshift separation between these two protoclusters is $Δz=0.05$, which is slightly larger than the size of typical protoclusters. Therefore, if they are not the progenitors of a $>10^{15}\,\mathrm{M_\odot}$ halo, they would grow into closely-located independent halos like a supercluster. The other protocluster at $z=4.898$ is also surrounded by smaller galaxy groups. These systems including protoclusters and neighboring groups are regarded as the early phase of superclusters. We quantify the spatial distribution of member galaxies of the protoclusters at $z=3.675$ and 3.721 by fitting triaxial ellipsoids, finding a tentative difference: one has a pancake-like shape while the other is filamentary. This could indicate that these two protoclusters are in different stages of formation. We investigate the relation between redshift and the velocity dispersion of protoclusters, including other protoclusters from the literature, in order to compare their dynamical states. Although there is no significant systematic trend in the velocity dispersions of protoclusters with redshift, the distribution is skewed to higher velocity dispersion over the redshift range of $z=2\mathrm{-}6$. This could be interpreted as two phases of cluster formation, one dominated by the steady accretion of galaxies, and the other by the merging between group-size halos, perhaps depending on the surrounding large-scale environments.

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.

preprint2019arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) VIII. A less biased view of the early co-evolution of black holes and host galaxies

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$σ$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

preprint2016arXiv

A Systematic Survey of Protoclusters at $z\sim3\mathrm{-}6$ in the CFHTLS Deep Fields

We present the discovery of three protoclusters at $z\sim3\mathrm{-}4$ with spectroscopic confirmation in the Canada-France-Hawaii Telescope (CFHT) Legacy Survey Deep Fields. In these fields, we investigate the large-scale projected sky distribution of $z\sim3\mathrm{-}6$ Lyman break galaxies and identify 21 protocluster candidates from regions that are overdense at more than $4σ$ overdensity significance. Based on cosmological simulations, it is expected that more than $76\%$ of these candidates will evolve into a galaxy cluster of at least a halo mass of $10^{14}\,\mathrm{M_\odot}$ at $z=0$. We perform follow-up spectroscopy for eight of the candidates using Subaru/FOCAS, KeckII/DEIMOS, and Gemini-N/GMOS. In total we target 462 dropout candidates and obtain 138 spectroscopic redshifts. We confirm three real protoclusters at $z=3\mathrm{-}4$ with more than five members spectroscopically identified, and find one to be an incidental overdense region by mere chance alignment. The other four candidate regions at $z\sim5\mathrm{-}6$ require more spectroscopic follow-up in order to be conclusive. A $z=3.67$ protocluster, which has eleven spectroscopically confirmed members, shows a remarkable core-like structure composed of a central small region ($<0.5\,\mathrm{physical\>Mpc}$) and an outskirts region ($\sim1.0\,\mathrm{physical\>Mpc}$). The Ly$α$ equivalent widths of members of the protocluster are significantly smaller than those of field galaxies at the same redshift while there is no difference in the UV luminosity distributions. These results imply that some environmental effects start operating as early as at $z\sim4$ along with the growth of the protocluster structure.

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

Detection of emission lines from z ~ 3 DLAs towards the QSO J2358+0149

Using VLT/X-shooter we searched for emission line galaxies associated to four damped Lyman-$α$ systems (DLAs) and one sub-DLA at 2.73<=z<=3.25 towards QSO J2358+0149. We detect [O III] emission from a "low-cool" DLA at z_abs = 2.9791 (having log N(HI)=21.69+\-0.10, [Zn/H] = -1.83+\-0.18) at an impact parameter of, $ρ$ ~12 kpc. The associated galaxy is compact with a dynamical mass of (1-6)x10^9 M_solar, very high excitation ([O III]/[O II] and [O III]/[H$β$] both greater than 10), 12+[O/H]<=8.5 and moderate star formation rate (SFR <=2 M_solar yr^{-1}). Such properties are typically seen in the low-z extreme blue compact dwarf galaxies. The kinematics of the gas is inconsistent with that of an extended disk and the gas is part of either a large scale wind or cold accretion. We detect Ly$α$ emission from the z_abs = 3.2477 DLA (having log N(HI)=21.12+\-0.10 and [Zn/H]=-0.97+\-0.13).The Ly$α$ emission is redshifted with respect to the metal absorption lines by 320 km s^{-1}, consistent with the location of the red hump expected in radiative transport models. We derive SFR ~0.2-1.7 M_solar yr^{-1} and Ly$α$ escape fraction of >=10 per cent. No other emission line is detected from this system. Because the DLA has a small velocity separation from the quasar (~500 km s^{-1}) and the DLA emission is located within a small projected distance ($ρ<5$ kpc), we also explore the possibility that the Ly$α$ emission is being induced by the QSO itself. QSO induced Ly$α$ fluorescence is possible if the DLA is within a physical separation of 340 kpc to the QSO. Detection of stellar continuum light and/or the oxygen emission lines would disfavor this possibility. We do not detect any emission line from the remaining three systems.

preprint2016arXiv

Subaru high-z exploration of low-luminosity quasars (SHELLQs). I. Discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

preprint2016arXiv

The Metal Abundances across Cosmic Time ($\mathcal{MACT}$) Survey. I. Optical Spectroscopy in the Subaru Deep Field

Deep rest-frame optical spectroscopy is critical for characterizing and understanding the physical conditions and properties of the ionized gas in galaxies. Here, we present a new spectroscopic survey called "Metal Abundances across Cosmic Time" or $\mathcal{MACT}$, which will obtain rest-frame optical spectra for $\sim$3000 emission-line galaxies. This paper describes the optical spectroscopy that has been conducted with MMT/Hectospec and Keck/DEIMOS for $\approx$1900 $z=0.1-1$ emission-line galaxies selected from our narrowband and intermediate-band imaging in the Subaru Deep Field. In addition, we present a sample of 164 galaxies for which we have measured the weak [OIII]$λ$4363 line (66 with at least 3$σ$ detections and 98 with significant upper limits). This nebular emission line determines the gas-phase metallicity by measuring the electron temperature of the ionized gas. This paper presents the optical spectra, emission-line measurements, interstellar properties (e.g., metallicity, gas density), and stellar properties (e.g., star formation rates, stellar mass). Paper II of the $\mathcal{MACT}$ survey (Ly et al.) presents the first results on the stellar mass--gas metallicity relation at $z\lesssim1$ using the sample with [OIII]$λ$4363 measurements.

preprint2016arXiv

The very wide-field $gzK$ galaxy survey -- I. Details of the clustering properties of star-forming galaxies at $z \sim 2$

We present the results of clustering analysis on the $z \sim 2$ star-forming galaxies. By combining our data with data from publicly available archives, we collect $g$-, $\zb / z$-, and $K$-band imaging data over 5.2 deg$^{2}$, which represents the largest area BzK/gzK survey. We apply colour corrections to translate our filter-set to those used in the original BzK selection for the gzK selection. Because of the wide survey area, we obtain a sample of 41,112 star-forming gzK galaxies at $z \sim 2$ (sgzKs) down to $\KAB < 23.0$, and determine high-quality two-point angular correlation functions (ACFs). Our ACFs show an apparent excess from power-law behaviour at small angular scale $(θ\la 0.01^{\circ})$, which corresponds the virial radius of a dark halo at $z \sim 2$ with a mass of $\sim 10^{13} \Msun$. We find that the correlation lengths are consistent with the previous estimates over all magnitude range; however, our results are evaluated with a smaller margin of error than that in previous studies. The large amount of data enables us to determine ACFs differentially depending on the luminosity of the subset of the data. The mean halo mass of faint sgzKs $(22.0 < K \leq 23.0)$ was found to be $\Mh = (1.32^{+0.09}_{-0.12}) \times 10^{12} h^{-1} \Msun$, whereas bright sgzKs ($18.0 \leq K \leq 21.0)$ were found to reside in dark haloes with a mass of $\Mh = (3.26^{+1.23}_{-1.02}) \times 10^{13} h^{-1} \Msun$.

preprint2016arXiv

The very wide-field $gzK$ galaxy survey -- II. The relationship between star-forming galaxies at $z \sim 2$ and their host haloes based upon HOD modelling

We present the results of an halo occupation distribution (HOD) analysis of star-forming galaxies at $z \sim 2$. We obtained high-quality angular correlation functions based on a large sgzK sample, which enabled us to carry out the HOD analysis. The mean halo mass and the HOD mass parameters are found to increase monotonically with increasing $K$-band magnitude, suggesting that more luminous galaxies reside in more massive dark haloes. The luminosity dependence of the HOD mass parameters was found to be the same as in the local Universe; however, the masses were larger than in the local Universe over all ranges of magnitude. This implies that galaxies at $z \sim 2$ tend to form in more massive dark haloes than in the local Universe, a process known as downsizing. By analysing the dark halo mass evolution using the extended Press--Schechter formalism and the number evolution of satellite galaxies in a dark halo, we find that faint Lyman break galaxies at $z \sim 4$ could evolve into the faintest sgzKs $(22.0 < K \leq 23.0)$ at $z \sim 2$ and into the Milky-Way-like galaxies or elliptical galaxies in the local Universe, whereas the most luminous sgzKs $(18.0 \leq K \leq 21.0)$ could evolve into the most massive systems in the local Universe. The stellar-to-halo mass ratio (SHMR) of the sgzKs was found to be consistent with the prediction of the model, except that the SHMR of the faintest sgzKs was smaller than the prediction at $z \sim 2$. This discrepancy may be attributed that our samples are confined to star-forming galaxies.

preprint2015arXiv

Giant H$α$ nebula surrounding the starburst merger NGC 6240

We revealed the detailed structure of a vastly extended H$α$-emitting nebula ("H$α$ nebula") surrounding the starburst/merging galaxy NGC 6240 by deep narrow-band imaging observations with the Subaru Suprime-Cam. The extent of the nebula is $\sim$90 kpc in diameter and the total H$α$ luminosity amounts to $\approx 1.6 \times 10^{42}$ erg s$^{-1}$. The volume filling factor and the mass of the warm ionized gas are $\sim$10$^{-4}$--10$^{-5}$ and $\sim$$5 \times 10^8$ $M_\odot$, respectively. The nebula has a complicated structure, which includes numerous filaments, loops, bubbles, and knots. We found that there is a tight spatial correlation between the H$α$ nebula and the extended soft X-ray-emitting gas, both in large and small scales. The overall morphology of the nebula is dominated by filamentary structures radially extending from the center of the galaxy. A large-scale bi-polar bubble extends along the minor axis of the main stellar disk. The morphology strongly suggests that the nebula was formed by intense outflows - superwinds - driven by starbursts. We also found three bright knots embedded in a looped filament of ionized gas that show head-tail morphologies in both emission-line and continuum, suggesting close interactions between the outflows and star forming regions. Based on the morphology and surface brightness distribution of the H$α$ nebula, we propose the scenario that three major episodes of starburst/superwind activities which were initiated $\sim$10$^2$ Myr ago formed the extended ionized gas nebula of NGC 6240.

preprint2015arXiv

Physical conditions of the interstellar medium in star-forming galaxies at z~1.5

We present results from Subaru/FMOS near-infrared (NIR) spectroscopy of 118 star-forming galaxies at $z\sim1.5$ in the Subaru Deep Field. These galaxies are selected as [OII]$λ$3727 emitters at $z\approx$ 1.47 and 1.62 from narrow-band imaging. We detect H$α$ emission line in 115 galaxies, [OIII]$λ$5007 emission line in 45 galaxies, and H$β$, [NII]$λ$6584, and [SII]$λλ$6716,6731 in 13, 16, and 6 galaxies, respectively. Including the [OII] emission line, we use the six strong nebular emission lines in the individual and composite rest-frame optical spectra to investigate physical conditions of the interstellar medium in star-forming galaxies at $z\sim$1.5. We find a tight correlation between H$α$ and [OII], which suggests that [OII] can be a good star formation rate (SFR) indicator for galaxies at $z\sim1.5$. The line ratios of H$α$/[OII] are consistent with those of local galaxies. We also find that [OII] emitters have strong [OIII] emission lines. The [OIII]/[OII] ratios are larger than normal star-forming galaxies in the local Universe, suggesting a higher ionization parameter. Less massive galaxies have larger [OIII]/[OII] ratios. With evidence that the electron density is consistent with local galaxies, the high ionization of galaxies at high redshifts may be attributed to a harder radiation field by a young stellar population and/or an increase in the number of ionizing photons from each massive star.

preprint2015arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at $z\ge6$. III. Stellar Populations from SED Modeling with Secure Ly$α$ Emission and Redshifts

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$α$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$α$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$α$ emission are indistinguishable in terms of age, stellar mass, and SFR.

preprint2014arXiv

A First Site of Galaxy Cluster Formation: Complete Spectroscopy of a Protocluster at $z=6.01$

We performed a systematic spectroscopic observation of a protocluster at $z=6.01$ in the Subaru Deep Field. We took spectroscopy for all 53 $i'$-dropout galaxies down to $z'=27.09\,\mathrm{mag}$ in/around the protocluster region. From these observations, we confirmed that 28 galaxies are at $z\sim6$, of which ten are clustered in a narrow redshift range of $Δz<0.06$. To trace the evolution of this primordial structure, we applied the same $i'$-dropout selection and the same overdensity measurements used in the observations to a semi-analytic model built upon the Millennium Simulation. We obtain a relation between the significance of overdensities observed at $z\sim6$ and the predicted dark matter halo mass at $z=0$. This protocluster with $6σ$ overdensity is expected to grow into a galaxy cluster with a mass of $\sim5\times10^{14}\,\mathrm{M_\odot}$ at $z=0$. Ten galaxies within $10\,\mathrm{comoving\>Mpc}$ of the overdense region can, with more than an 80% probability, merge into a single dark matter halo by $z=0$. No significant differences appeared in UV and Ly$α$ luminosities between the protocluster and field galaxies, suggesting that this protocluster is still in the early phase of cluster formation before the onset of any obvious environmental effects. However, further observations are required to study other properties, such as stellar mass, dust, and age. We do find that galaxies tend to be in close pairs in this protocluster. These pair-like subgroups will coalesce into a single halo and grow into a more massive structure. We may witness an onset of cluster formation at $z\sim6$ toward a cluster as seen in local universe.

preprint2014arXiv

ALMA Observation of 158 micron [CII] Line and Dust Continuum of a z=7 Normally Star-forming Galaxy in the Epoch of Reionization

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

preprint2014arXiv

The Subaru high-z quasar survey: discovery of faint z~6 quasars

We present the discovery of one or two extremely faint z~6 quasars in 6.5 deg^2 utilizing a unique capability of the wide-field imaging of the Subaru/Suprime-Cam. The quasar selection was made in (i'-z_B) and (z_B-z_R) colors, where z_B and z_R are bandpasses with central wavelengths of 8842A and 9841A, respectively. The color selection can effectively isolate quasars at z~6 from M/L/T dwarfs without the J-band photometry down to z_R<24.0, which is 3.5 mag. deeper than SDSS. We have selected 17 promising quasar candidates. The follow-up spectroscopy for seven targets identified one apparent quasar at z=6.156 with M_1450=-23.10. We also identified one possible quasar at z=6.041 with a faint continuum of M_1450=-22.58 and a narrow Lyman-alpha emission with HWHM=427 km/s, which cannot be distinguished from Lyman-alpha emitters. We derive the quasar luminosity function at z~6 by combining our faint quasar sample with the bright quasar samples by SDSS and CFHQS. Including our data points invokes a higher number density in the faintest bin of the quasar luminosity function than the previous estimate employed. This suggests a steeper faint-end slope than lower-z, though it is yet uncertain based on a small number of spectroscopically identified faint quasars and several quasar candidates are still remain to be diagnosed. The steepening of the quasar luminosity function at the faint-end does increase the expected emission rate of the ionizing photon, however, it only changes by a factor of ~2-6. This was found to be still insufficient for the required photon budget of reionization at z~6.

preprint2013arXiv

"Direct" Gas-phase Metallicities, Stellar Properties, and Local Environments of Emission-line Galaxies at Redshift below 0.90

Using deep narrow-band (NB) imaging and optical spectroscopy from the Keck telescope and MMT, we identify a sample of 20 emission-line galaxies (ELGs) at z=0.065-0.90 where the weak auroral emission line, [OIII]4363, is detected at >3σ. These detections allow us to determine the gas-phase metallicity using the "direct'' method. With electron temperature measurements and dust attenuation corrections from Balmer decrements, we find that 4 of these low-mass galaxies are extremely metal-poor with 12+log(O/H) <= 7.65 or one-tenth solar. Our most metal-deficient galaxy has 12+log(O/H) = 7.24^{+0.45}_{-0.30} (95% confidence), similar to some of the lowest metallicity galaxies identified in the local universe. We find that our galaxies are all undergoing significant star formation with average specific star formation rate (SFR) of (100 Myr)^{-1}, and that they have high central SFR surface densities (average of 0.5 Msun/yr/kpc^2. In addition, more than two-thirds of our galaxies have between one and four nearby companions within a projected radius of 100 kpc, which we find is an excess among star-forming galaxies at z=0.4-0.85. We also find that the gas-phase metallicities for a given stellar mass and SFR lie systematically below the local M-Z-(SFR) relation by \approx0.2 dex (2σ significance). These results are partly due to selection effects, since galaxies with strong star formation and low metallicity are more likely to yield [OIII]4363 detections. Finally, the observed higher ionization parameter and electron density suggest that they are lower redshift analogs to typical z>1 galaxies.

preprint2013arXiv

Extended Lyman-alpha emission from a high-z DLA at z=3.115

We searched for star formation activity associated with high-z Damped Lyman-alpha systems (DLAs) with Subaru telescope. We used a set of narrow-band (NB) filters whose central wavelengths correspond to the redshifted Lyman-alpha emission lines of targeted DLA absorbers at 3<z<4.5. We detected one apparent NB-excess object located 3.80 arcsec (~28kpc) away from the quasar SDSS J031036.84+005521.7. Follow-up spectroscopy revealed an asymmetric Lyman-alpha emission at z_em=3.115+/-0.003, which perfectly matches the sub-DLA trough at z_abs=3.1150 with logN(HI)/cm^-2=20.05. The Lyman-alpha luminosity is estimated to be L(LyA)=1.07x10^42 erg s^-1, which corresponds to a star formation rate of 0.97 M_\odot yr^-1. Interestingly, the detected Lyman-alpha emission is spatially extended with a sharp peak. The large extent of the Lyman-alpha emission is remarkably one-sided toward the quasar line-of-sight, and is redshifted. The observed spatially asymmetric surface brightness profile can be qualitatively explained by a model of a DLA host galaxy, assuming a galactic outflow and a clumpy distribution of HI clouds in the circumgalactic medium. This large Lyman-alpha extension, which is similar to those found in Rauch et al. (2008), could be the result of complicated anisotropic radiative transfer through the surrounding neutral gas embedded in the DLA.

preprint2013arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. I. Basic Characteristics of the Rest-Frame UV Continuum and Lyman-alpha emission

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 Å. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 Å) that could be the dominant contribution to the total ionizing flux at z >= 6.

preprint2013arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. II. Morphology of the Rest-Frame UV Continuum and Lyman-alpha Emission

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

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

A Search for z=7.3 Lyα Emitters behind Gravitationally Lensing Clusters

We searched for z=7.3 Lya emitters (LAEs) behind two lensing clusters, Abell 2390 and CL 0024, with the Subaru Telescope Suprime-Cam and a narrowband NB1006 (FWHM ~ 21 nm centered at 1005 nm). We investigated if there exist objects consistent with the color of z=7.3 LAEs behind the clusters but could not detect any LAEs to the unlensed line limit F(Lya) ~ 6.9 x 10^{-18} erg/s/cm^2. Using several z=7 Lya luminosity functions (LFs) from the literature, we estimated and compared the expected detection numbers of z ~ 7 LAEs in lensing and blank field surveys in the case of using an 8m class ground based telescope. Given the steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient in detecting z ~ 7 LAEs than imaging a blank field. However, the gain is expected to be modest, a factor of two at most and likely much less depending on the adopted LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain in depth and not necessarily in detection efficiency. For much larger detectors, the lensing effect becomes negligible and the efficiency of LAE detection is proportional to the instrumental FOV. We also inspected NB1006 images of three z ~ 7 z-dropouts previously detected in Abell 2390 and found that none of them are detected in NB1006. Two of them are consistent with predictions from the previous studies that they would be at lower redshifts. The other one has a photometric redshift of z ~ 7.3, and if it is at z=7.3, its unlensed Lya line flux would be very faint: F(Lya) < 4.4 x 10^{-18} erg/s/cm^2 (1 sigma upper limit) or rest frame equivalent width of W(Lya) < 26A. Its Lya emission might be attenuated by neutral hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying strong Lya emission is lower at z ~ 7 than at z <~ 6.

preprint2012arXiv

Discovery of a protocluster at z~6

We report the discovery of a protocluster at z~6 containing at least eight cluster member galaxies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The overdensity of the protocluster is significant at the 6 sigma level, based on the surface number density of i'-dropout galaxies. The overdense region covers ~36 sq. arcmin, and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real z~6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range centered at z=6.01, corresponding to a seven-fold increase in number density over the average in redshift space. We found no significant difference in the observed properties, such as Ly-alpha luminosities and UV continuum magnitudes, between the eight protocluster members and the seven non-members. The velocity dispersion of the eight protocluster members is 647 km/s, which is about three times higher than that predicted by the standard cold dark matter model. This discrepancy could be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. We discuss two possible explanations for this discrepancy: either the protocluster is already mature, with old galaxies at the center, or it is still immature and composed of three subgroups merging to become a larger cluster. In either case, this concentration of z=6.01 galaxies in the SDF may be one of the first sites of formation of a galaxy cluster in the universe.

preprint2012arXiv

Dust Attenuation and H-alpha Star Formation Rates of z~0.5 Galaxies

Using deep narrow-band and broad-band imaging, we identify 401 z~0.40 and 249 z~0.49 H-alpha line-emitting galaxies in the Subaru Deep Field. Compared to other H-alpha surveys at similar redshifts, our samples are unique since they probe lower H-alpha luminosities, are augmented with multi-wavelength (rest-frame 1000AA--1.5 microns) coverage, and a large fraction (20%) of our samples has already been spectroscopically confirmed. Our spectra allow us to measure the Balmer decrement for nearly 60 galaxies with H-beta detected above 5-sigma. The Balmer decrements indicate an average extinction of A(H-alpha)=0.7^{+1.4}_{-0.7} mag. We find that the Balmer decrement systematically increases with higher H-alpha luminosities and with larger stellar masses, in agreement with previous studies with sparser samples. We find that the SFRs estimated from modeling the spectral energy distribution (SED) is reliable---we derived an "intrinsic" H-alpha luminosity which is then reddened assuming the color excess from SED modeling. The SED-predicted H-alpha luminosity agrees with H-alpha narrow-band measurements over 3 dex (rms of 0.25 dex). We then use the SED SFRs to test different statistically-based dust corrections for H-alpha and find that adopting one magnitude of extinction is inappropriate: galaxies with lower luminosities are less reddened. We find that the luminosity-dependent dust correction of Hopkins et al. yields consistent results over 3 dex (rms of 0.3 dex). Our comparisons are only possible by assuming that stellar reddening is roughly half of nebular reddening. The strong correspondence argue that with SED modeling, we can derive reliable intrinsic SFRs even in the absence of H-alpha measurements at z~0.5.

preprint2012arXiv

Kinematics and excitation of the ram pressure stripped ionized gas filaments in the Coma cluster of galaxies

We present the results of deep imaging and spectroscopic observations of very extended ionized gas (EIG) around four member galaxies of the Coma cluster of galaxies: RB199, IC4040, GMP2923 and GMP3071. The EIGs were serendipitously found in an H-alpha narrow band imaging survey of the central region of the Coma cluster. The relative radial velocities of the EIGs with respect to the systemic velocities of the parent galaxies from which they emanate increase almost monotonically with the distance from the nucleus of the respective galaxies, reaching -400 - -800 km/s at around 40 - 80 kpc from the galaxies. The one-sided morphologies and the velocity fields of the EIGs are consistent with the predictions of numerical simulations of ram pressure stripping. We found a very low-velocity filament (v_rel = -1300 km/s) at the southeastern edge of the disk of IC4040. Some bright compact knots in the EIGs of RB199 and IC4040 exhibit blue continuum and strong H-alpha emission. The equivalent widths of the H-alpha emission exceed 200 A, and are greater than 1000 A for some knots. The emission line intensity ratios of the knots are basically consistent with those of sub-solar abundance HII regions. These facts indicate that intensive star formation occurs in the knots. Some filaments, including the low velocity filament of the IC4040 EIG, exhibit shock-like emission line spectra, suggesting that shock heating plays an important role in ionization and excitation of the EIGs.

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.

preprint2012arXiv

The First Systematic Survey for Lyman Alpha Emitters at z=7.3 with Red-sensitive Subaru/Suprime-Cam

We have performed deep imaging surveys for LyA emitters (LAEs) at redshift ~7.3 in two blank fields, the Subaru Deep Field (SDF) and the Subaru/XMM-Newton Deep survey Field (SXDF), using the Subaru/Suprime-Cam equipped with new red-sensitive CCDs and a new narrow-band filter, NB1006 (lambda_c=10052 Ang, FWHM=214 Ang). We identified four objects as LAE candidates that exhibit luminosity excess in NB1006. By carrying out deep follow-up spectroscopy for three of them using Subaru/FOCAS and Keck/DEIMOS, a definitively asymmetric emission line is detected for one of them, SXDF-NB1006-2. Assuming this line is LyA, this object is a LAE at z=7.215 which has luminosity of 1.2^{+1.5}_{-0.6} x 10^43 [erg s-1] and a weighted skewness S_w=4.90+-0.86. Another object, SDF-NB1006-2, shows variable photometry and is thus probably a quasar (QSO) or an active galactic nucleus (AGN). It shows an asymmetric emission line at 10076 Ang, which may be due to either LyA at z=7.288 or [OII] at z=1.703. The third object, SDF-NB1006-1, is likely a galaxy with temporal luminosity enhancement associated with a supernova explosion, as the brightness of this object varies between the observed epochs. Its spectrum does not show any emission lines. The inferred decrease in the number density of LAEs toward higher redshift is n_LyA(z=7.3)/n_LyA(z=5.7) = 0.05^+0.11_-0.05 from z=5.7 to 7.3 down to L(LyA)=1.0 x 10^43 [erg s-1]. The present result is consistent with the interpretation in previous studies that the neutral hydrogen fraction is rapidly increasing from z=5.7 to 7.3.

preprint2012arXiv

The Galaxy Luminosity Functions down to M_R=-10 in the Coma Cluster

We derived the luminosity function (LF) of dwarf galaxies in the Coma Cluster down to M_R=-10 at three fields located at the center, intermediate, and outskirt. The LF (-19<M_R<-10) shows no significant differences among the three fields. It shows a clear dip at M_R\sim-13, and is composed of two distinct components of different slopes; the bright component with -19<M_R<-13 has a flatter slope than the faint component with -13<M_R<-10 which has a steep slope. The bright component (-19<M_R<-13) consists of mostly red extended galaxies including few blue galaxies whose colors are typical of late-type galaxies. On the other hand, the faint component (-13<M_R<-10) consists of largely PSF-like compact galaxies. We found that both these compact galaxies and some extended galaxies are present in the center while only compact galaxies are seen in the outskirt. In the faint component, the fraction of blue galaxies is larger in the outskirt than in the center. We suggest that the dwarf galaxies in the Coma Cluster, which make up the two components in the LF, are heterogeneous with some different origins.

preprint2012arXiv

The Stellar Population and Star Formation Rates of z~1.5-1.6 [O II] Emitting Galaxies Selected from Narrow-Band Emission-Line Surveys

We present the first detailed study of the stellar populations of star-forming galaxies at z~1.5, which are selected by their [O II] emission line, detected in narrow-band surveys. We identified ~1,300 [O II] emitters at z=1.47 and z=1.62 in the Subaru Deep Field with rest-frame EWs above 13Å. Optical and near-infrared spectroscopic observations for ~10% of our samples show that our separation of [O II] from [O III] emission-line galaxies in two-color space is 99% successful. We analyze the multi-wavelength properties of a subset of ~1,200 galaxies with the best photometry. They have average rest-frame EW of 45Å, stellar mass of 3 x 10^9 M_sun, and stellar age of 100 Myr. In addition, our SED fitting and broad-band colors indicate that [O II] emitters span the full range of galaxy populations at z~1.5. We also find that 80% of [O II] emitters are also photometrically classified as "BX/BM" (UV) galaxies and/or the star-forming "BzK" (near-IR) galaxies. Our [O II] emission line survey produces a far more complete, and somewhat deeper sample of z~1.5 galaxies than either the BX/BM or sBzK selection alone. We constructed average SEDs and find that higher [O II] EW galaxies have somewhat bluer continua. SED model-fitting shows that they have on average half the stellar mass of galaxies with lower [O II] EW. The observed [O II] luminosity is well-correlated with the far-UV continuum with a logarithmic slope slightly 0f 0.89\pm0.22. The scatter of the [O II] luminosity against the far-UV continuum suggests that [O II] can be used as a SFR indicator with a reliability of 0.23 dex.

preprint2012arXiv

Unusual Long and Luminous Optical Transient in the Subaru Deep Field

We present observations of SDF-05M05, an unusual optical transient discovered in the Subaru Deep Field (SDF). The duration of the transient is > ~800 d in the observer frame, and the maximum brightness during observation reached approximately 23 mag in the i' and z' bands. The faint host galaxy is clearly identified in all 5 optical bands of the deep SDF images. The photometric redshift of the host yields z~0.6 and the corresponding absolute magnitude at maximum is ~-20. This implies that this event shone with an absolute magnitude brighter than -19 mag for approximately 300 d in the rest frame, which is significantly longer than a typical supernova and ultra-luminous supernova. The total radiated energy during our observation was 1x10^51 erg. The light curves and color evolution are marginally consistent with some of luminous IIn supernova. We suggest that the transient may be a unique and peculiar supernova at intermediate redshift.

preprint2011arXiv

A Census of Star-Forming Galaxies at z = 1-3 in the Subaru Deep Field

Several UV and near-infrared color selection methods have identified galaxies at z = 1-3. Since each method suffers from selection biases, we have applied three leading techniques (Lyman break, BX/BM, and BzK selection) simultaneously in the Subaru Deep Field. This field has reliable (Δz/(1 + z) = 0.02--0.09) photometric redshifts for ~53,000 galaxies from 20 bands (1500Å--2.2μm). The BzK, LBG, and BX/BM samples suffer contamination from z<1 interlopers of 6%, 8%, and 20%, respectively. Around the redshifts where it is most sensitive (z~1.9 for star-forming BzK, z~1.8 for z~2 LBGs, z~1.6 for BM, and z~2.3 for BX), each technique finds 60-80% of the census of the three methods. In addition, each of the color techniques shares 75-96% of its galaxies with another method, which is consistent with previous studies that adopt identical criteria on magnitudes and colors. Combining the three samples gives a comprehensive census that includes ~90% of z-phot = 1-3 galaxies, using standard magnitude limits similar to previous studies. In fact, we find that among z = 1-2.5 galaxies in the color selection census, 81-90% of them can be selected by just combining the BzK selection with one of the UV techniques (z~2 LBG or BX and BM). The average galaxy stellar mass, reddening and SFRs all decrease systematically from the sBzK population to the LBGs, and to the BX/BMs. The combined color selections yield a total cosmic SFR density of 0.18 $\pm$ 0.03 M_sun yr^{-1} Mpc^{-3} for K_AB <= 24. We find that 65% of the star formation is in galaxies with E(B-V) > 0.25 mag, even though they are only one-fourth of the census by number.

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

Cryogenic Volume-Phase Holograpic Grisms for MOIRCS

We have developed high dispersion VPH (volume phase holographic) grisms with zinc selenide (ZnSe) prisms for the cryogenic optical system of MOIRCS (Multi-Object near InfraRed Camera and Spectrograph) for Y-, J-, H- and K- band observations. We fabricated the VPH gratings using a hologram resin. After several heat cycles at between room temperature and 120 K, the VPH gratings were assembled to grisms by gluing with two ZnSe prisms. Several heat cycles were also carried out for the grisms before being installed into MOIRCS. We measured the efficiencies of the VPH grisms in a laboratory, and found them to be 70% - 82%. The performances obtained by observations of MOIRCS with the 8.2 m Subaru Telescope have been found to be very consistent with the results in the laboratory test. This is the first astronomical application of cryogenic VPH grisms.

preprint2011arXiv

Grisms Developed for FOCAS

Faint Object Camera and Spectrograph (FOCAS) is a versatile common-use optical instrument for the 8.2m Subaru Telescope, offering imaging and spectroscopic observations. FOCAS employs grisms with resolving powers ranging from 280 to 8200 as dispersive optical elements. A grism is a direct-vision grating composed of a transmission grating and prism(s). FOCAS has five grisms with replica surface-relief gratings including an echelle-type grism, and eight grisms with volume-phase holographic (VPH) gratings. The size of these grisms is 110 mm X 106 mm in aperture with a maximum thickness of 110 mm. We employ not only the dichromated gelatin, but also the hologram resin as a recording material for VPH gratings. We discuss the performance of these FOCAS grisms measured in the laboratory, and verify it by test observations, and show examples of astronomical spectroscopic observations.

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

Search for z~7 Ly-alpha emitters with Suprime-Cam at the Subaru Telescope

We report a search for z=7 Ly-alpha emitters (LAEs) using a custom-made Narrow-Band filter, centered at 9755 Angstroms, with the instrument Suprime-Cam installed at the Subaru telescope. We observed two different fields and obtained two sample of 7 Ly-alpha emitters of which 4 are robust in each field. We are covering the luminosity range of 9.10^{42} - 2.10^{43} erg/s in comoving volumes of ~ 4.10^{5} and 4.3.10^{5} Mpc^{3}. From this result, we derived possible z~7 Ly-alpha luminosity functions for the full samples and for a subsample of 4 objects in each field. We do not observe, in each case, any strong evolution between the z=6.5 and z~7 Ly-alpha luminosity functions. Spectroscopic confirmation for these candidate samples is required to establish a definitive measure of the luminosity function at z~7.

preprint2011arXiv

Spectroscopic Confirmation of Three z-Dropout Galaxies at z = 6.844 - 7.213: Demographics of Lyman-Alpha Emission in z ~ 7 Galaxies

We present the results of our ultra-deep Keck/DEIMOS spectroscopy of z-dropout galaxies in the SDF and GOODS-N. For 3 out of 11 objects, we detect an emission line at ~ 1um with a signal-to-noise ratio of ~ 10. The lines show asymmetric profiles with high weighted skewness values, consistent with being Lya, yielding redshifts of z=7.213, 6.965, and 6.844. Specifically, we confirm the z=7.213 object in two independent DEIMOS runs with different spectroscopic configurations. The z=6.965 object is a known Lya emitter, IOK-1, for which our improved spectrum at a higher resolution yields a robust skewness measurement. The three z-dropouts have Lya fluxes of 3 x 10^-17 erg s^-1 cm^-2 and rest-frame equivalent widths EW_0^Lya = 33-43A. Based on the largest spectroscopic sample of 43 z-dropouts that is the combination of our and previous data, we find that the fraction of Lya-emitting galaxies (EW_0^Lya > 25A) is low at z ~ 7; 17 +- 10% and 24 +- 12% for bright (Muv ~= -21) and faint (Muv ~= -19.5) galaxies, respectively. The fractions of Lya-emitting galaxies drop from z ~ 6 to 7 and the amplitude of the drop is larger for faint galaxies than for bright galaxies. These two pieces of evidence would indicate that the neutral hydrogen fraction of the IGM increases from z ~ 6 to 7, and that the reionization proceeds from high- to low-density environments, as suggested by an inside-out reionization model.

preprint2010arXiv

A dozen new galaxies caught in the act: Gas stripping and extended emission line regions in the Coma cluster

We present images of extended H-alpha clouds associated with 14 member galaxies in the Coma cluster obtained from deep narrow band imaging observations with Suprime-Cam at the Subaru Telescope. The parent galaxies of the extended H-alpha clouds are distributed farther than 0.2 Mpc from the peak of X-ray emission of the cluster. Most of the galaxies have colors bluer than g-r approx 0.5 and they account for 57% of the blue (g-r<0.5) bright (r<17.8 mag) galaxies in the central region of the Coma cluster. They reside near the red- and blue-shifted edges of the Coma cluster's radial velocity distribution. These findings suggest that the most of the parent galaxies were recently captured by the Coma cluster potential and are now infalling toward the cluster center with their disk gas being stripped off and producing the observed H-alpha clouds.

preprint2010arXiv

A Large Number of z > 6 Galaxies around a QSO at z = 6.43: Evidence for a Protocluster?

QSOs have been thought to be important for tracing highly biased regions in the early universe, from which the present-day massive galaxies and galaxy clusters formed. While overdensities of star-forming galaxies have been found around QSOs at 2<z<5, the case for excess galaxy clustering around QSOs at z>6 is less clear. Previous studies with HST have reported the detection of small excesses of faint dropout galaxies in some QSO fields, but these surveys probed a relatively small region surrounding the QSOs. To overcome this problem, we have observed the most distant QSO at z=6.4 using the large field of view of the Suprime-Cam (34' x 27'). Newly-installed CCDs allowed us to select Lyman break galaxies (LBG) at z~6.4 more efficiently. We found seven LBGs in the QSO field, whereas only one exists in a comparison field. The significance of this apparent excess is difficult to quantify without spectroscopic confirmation and additional control fields. The Poisson probability to find seven objects when one expects four is ~10%, while the probability to find seven objects in one field and only one in the other is less than 0.4%, suggesting that the QSO field is significantly overdense relative to the control field. We find some evidence that the LBGs are distributed in a ring-like shape centered on the QSO with a radius of ~3 Mpc. There are no candidate LBGs within 2 Mpc from the QSO, i.e., galaxies are clustered around the QSO but appear to avoid the very center. These results suggest that the QSO is embedded in an overdense region when defined on a sufficiently large scale. This suggests that the QSO was indeed born in a massive halo. The central deficit of galaxies may indicate that (1) the strong UV radiation from the QSO suppressed galaxy formation in its vicinity, or (2) that star-formation closest to the QSO occurs mostly in an obscured mode that is missed by our UV selection.

preprint2010arXiv

Lyman Alpha Emitters at z=7 in the Subaru/XMM-Newton Deep Survey Field: Photometric Candidates and Luminosity Function

We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts.

preprint2010arXiv

Overdensity of i'-Dropout Galaxies in the Subaru Deep Field: A Candidate Protocluster at z ~ 6

We investigate the sky distribution of z ~ 6 Lyman break galaxies selected as i'-dropouts having i' - z' > 1.45 down to z' < 26.5 in the Subaru Deep Field (SDF). We discover 37 i'-dropouts clustered in a projected comoving 21.6 x 21.6 Mpc^2 region at z = 6, showing a local density excess. Carrying out follow-up spectroscopy, we identify four of them as Lyman-alpha emitters at z = 5.92, 6.01, 6.03 and 6.03 (spread over a distance of 46.6 Mpc). The number density of the cluster itself in SDF is ~ 2.2 x 10^{-7} Mpc^{-3}, smaller than those of protoclusters (i.e., forming galaxy clusters) at z ~ 2-5.7. Also, the structure shows ~4-21 times larger galaxy number density than those of z ~ 6 galaxies in a general field. It has a mass of M ~ 1.5^{+1.8}_{-0.5} x 10^{15}M_sun, comparable to those of z ~ 0-5 protoclusters. Since the contamination of our sample by interlopers is estimated to be quite low, 5.9%, most of the other unconfirmed i'-dropouts in the overdense region can be also z ~ 6 galaxies. Hence, it could be a candidate forming cluster at z ~ 6, representing a progenitor of galaxy clusters seen in the recent-day Universe.

preprint2010arXiv

Spitzer Space Telescope Constraint on the Stellar Mass of a z = 6.96 Lyman Alpha Emitter

We obtained mid-infrared 3.6 and 4.5 micron imaging of a z=6.96 Lyman alpha emitter (LAE) IOK-1 discovered in the Subaru Deep Field, using Spitzer Space Telescope Infrared Array Camera observations. After removal of a nearby bright source, we find that IOK-1 is not significantly detected in any of these infrared bands to m_3.6 ~ 24.00 and m_4.5 ~ 23.54 at 3 sigma. Fitting population synthesis models to the spectral energy distribution consisting of the upper limit fluxes of the optical to infrared non-detection images and fluxes in detection images, we constrain the stellar mass M* of IOK-1. This LAE could have either a mass as low as M* <~ 2-9 x 10^8 Msun for the young age (<~ 10 Myr) and the low dust reddening (A_V ~ 0) or a mass as large as M* <~ 1-4 x 10^{10} Msun for either the old age (> 100 Myr) or the high dust reddening (A_V ~ 1.5). This would be within the range of masses of z ~ 3-6.6 LAEs studied to date, ~ 10^6-10^{10} Msun. Hence, IOK-1 is not a particularly unique galaxy with extremely high mass or low mass but is similar to one of the LAEs seen at the later epochs.

preprint2009arXiv

Discovery of a Giant Lya Emitter Near the Reionization Epoch

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

preprint2009arXiv

Large Area Survey for z=7 Galaxies in SDF and GOODS-N: Implications for Galaxy Formation and Cosmic Reionization

We present results of our large-area survey for z'-band dropout galaxies at z=7 in a 1568 arcmin^2 sky area covering the SDF and GOODS-N fields. Combining our ultra-deep Subaru/Suprime-Cam z'- and y-band (lambda_eff=1um) images with legacy data of Subaru and HST, we have identified 22 bright z-dropout galaxies down to y=26, one of which has a spectroscopic redshift of z=6.96 determined from Lya emission. The z=7 luminosity function (LF) yields the best-fit Schechter parameters of phi*=0.69 +2.62/-0.55 x10^(-3) Mpc^(-3), Muv*=-20.10 +/-0.76 mag, and alpha=-1.72 +/-0.65, and indicates a decrease from z=6 at a >95% confidence level. This decrease is beyond the cosmic variance in our two fields, which is estimated to be a factor of <~2. We have found that the cosmic star formation rate density drops from the peak at z=2-3 to z=7 roughly by a factor of ~10 but not larger than ~100. A comparison with the reionization models suggests either that the Universe could not be totally ionized by only galaxies at z=7, or more likely that properties of galaxies at z=7 are different from those at low redshifts having, e.g., a larger escape fraction (>~0.2), and/or a flatter IMF. Our SDF z-dropout galaxies appear to form 60-Mpc long filamentary structures, and the z=6.96 galaxy with Lya emission is located at the center of an overdense region consisting of four UV bright dropout candidates, which might suggest an existence of a well-developed ionized bubble at z=7.

preprint2009arXiv

Lyman Break Galaxies at z~1.8-2.8: GALEX/NUV Imaging of the Subaru Deep Field

Abridged: A photometric sample of ~7100 V<25.3 Lyman break galaxies (LBGs) has been selected by combining Subaru/Suprime-Cam BVRci'z' data with deep GALEX/NUV imaging of the Subaru Deep Field. Follow-up spectroscopy confirmed 24 LBGs at 1.5<z<2.7. Among the optical spectra, 12 have Ly-alpha emission with rest-frame equivalent widths of ~5-60AA. The success rate for identifying LBGs as NUV-dropouts at 1.5<z<2.7 is 86%. The rest-frame UV (1700AA) luminosity function (LF) is constructed from the photometric sample with corrections for stellar contamination and z<1.5 interlopers. The LF is 1.7+/-0.1 times higher than those of z~2 BXs and z~3 LBGs. Three explanations were considered, and it is argued that significantly underestimating low-z contamination or effective comoving volume is unlikely: the former would be inconsistent with the spectroscopic sample at 93% confidence, and the second explanation would not resolve the discrepancy. The third scenario is that different photometric selection of the samples yields non-identical galaxy populations, such that some BX galaxies are LBGs and vice versa. This argument is supported by a higher surface density of LBGs at all magnitudes while the redshift distribution of the two populations is nearly identical. This study, when combined with other star-formation rate (SFR) density UV measurements from LBG surveys, indicates that there is a rise in the SFR density: a factor of 3-6 (3-10) increase from z~5 (z~6) to z~2, followed by a decrease to z~0. This result, along with past sub-mm studies that find a peak at z~2 in their redshift distribution, suggest that z~2 is the epoch of peak star-formation. Additional spectroscopy is required to characterize the complete shape of the z~2 LBG UV LF via measurements of contamination and accurate distances.

preprint2009arXiv

Subaru FOCAS Spectroscopic Observations of High-Redshift Supernovae

We present spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

preprint2009arXiv

Subaru Mid-infrared Imaging of the Quadruple Lenses. II. Unveiling Lens Structure of MG0414+0534 and Q2237+030

We present mid-infrared imaging at 11.7 $μ$m for the quadruple lens systems, MG0414+0534 and Q2237+030, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. MG0414+0534 is characterized by a bright pair of lensed images (A1, A2) and their optical flux ratio A2/A1 deviates significantly from the prediction of a smooth lens model. Q2237+030 is comprised of four lensed images, which are significantly affected by microlensing in a foreground lensing galaxy. Our mid-infrared observations of these lensed images have revealed that the mid-infrared flux ratio for A2/A1 of MG0414+0534 is nearly unity (0.90 +- 0.04). We find that this flux ratio is systematically small, at 4 to 5 sigma level, compared with the prediction of a best smooth lens model (1.09) represented by a singular isothermal ellipsoid and external shear, thereby suggesting the presence of substructures to explain our observational result. In contrast, for Q2237+030, our high signal-to-noise observation indicates that the mid-infrared flux ratios between all the four images of Q2237+030 are virtually consistent with the prediction of a smooth lens model. Based on the size estimate of the dust torus surrounding the nuclei of these QSOs, we set limits on the mass of a substructure in these lens systems, which can cause anomalies in the flux ratios. For MG0414+0534, since the required mass of a substructure inside its Einstein radius is > 360 M_sun, millilensing by a CDM substructure is most likely. If it is modeled as a singular isothermal sphere, the mass inside radius of 100 pc is given as > 1.0x10^5 M_sun. For Q2237+030, there is no significant evidence of millilensing, so the reported anomalous flux ratios in shorter wavelengths are entirely caused by microlensing by stars (abridged).

preprint2008arXiv

Star Formation Rates and Metallicities of K-selected Star Forming Galaxies at z~2

We present spectroscopy of 15 star-forming BzK galaxies (sBzKs) with K(AB)<23 in the Subaru Deep Field, for which Halpha and some other emission lines are detected in 0.9 to 2.3 micron spectra with a resolution of R=500. Using Halpha luminosities, we obtain star formation rates (SFRs), and then specific SFRs (SSFRs) dividing SFRs by stellar masses, which are derived from SED fitting to BVRi'z'K photometry. It is found that sBzKs with higher stellar masses have larger SFRs. A negative correlation is seen between stellar mass and SSFR, which is consistent with the previous results for z~2 galaxies. This implies that a larger growth of stellar mass occurs in less massive galaxies. In addition, gas-phase oxygen abundances, 12+log(O/H), are derived from the ratio of NII(lambda 6584) to Halpha using the N2 index method. We have found a correlation between stellar mass and oxygen abundance in the sense that more massive sBzKs tend to be more metal rich, which is qualitatively consistent with the relation for UV-selected z~2 galaxies. However, the metallicity of the sBzKs is ~0.2 dex higher than that of UV-selected galaxies with similar stellar masses, which is significant considering the small uncertainties. The sBzKs in our sample have redder R-K colors than the UV-selected galaxies. This galaxy color-dependence in the oxygen abundance may be caused by older or dustier galaxies having higher metallicities at z~2.

preprint2008arXiv

Subaru Weak Lensing survey -- II: Multi-object Spectroscopy and Cluster Masses

We present the first results of a MOS campaign to follow up cluster candidates located via weak lensing. Our main goals are to search for spatial concentrations of galaxies that are plausible optical counterparts of the weak lensing signals, and to determine the cluster redshifts from those of member galaxies. Around each of 36 targeted cluster candidates, we obtain 15-32 galaxy redshifts. For 28 of these targets, we confirm a secure cluster identification. This includes three cases where two clusters at different redshifts are projected along the same line-of-sight. In 6 of the 8 unconfirmed targets, we find multiple small galaxy concentrations at different redshifts. In both the remaining two targets, a single small galaxy concentration is found. We evaluate the weak lensing mass of confirmed clusters. For a subsample of our most cleanly measured clusters, we investigate the statistical relation between their weak lensing mass and the velocity dispersion of their member galaxies, comparing our sample with optically and X-ray selected samples from the literature. Our lensing-selected clusters are consistent with sigma_v=sigma_sis, with a similar scatter to the optically and X-ray selected clusters. We thus find no evidence of selection bias compared to these other techniques. We also derive an empirical relation between the cluster mass and the galaxy velocity dispersion, which is in reasonable agreement with the prediction of N-body simulations in the LCDM cosmology.

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.

preprint2007arXiv

Luminosity dependent clustering of star-forming BzK galaxies at redshift 2

We use the BzK color selection proposed by Daddi et al. (2004) to obtain a sample of 1092 faint star-forming galaxies (hereafter sBzKs) from 180 arcmin^2 in the Subaru Deep Field. This sample represents star-forming galaxies at 1.4 < z < 2.5 down to K(AB)=23.2, which roughly corresponds to a stellar-mass limit of ~ 1 x 10^{10} Msun. We measure the angular correlation function (ACF) of these sBzKs to be w(theta) = (0.58 +- 0.13) x theta["]^{-0.8} and translate the amplitude into the correlation length assuming a reasonable redshift distribution. The resulting value, r0 = 3.2^{+0.6}_{-0.7} h^{-1} Mpc, suggests that our sBzKs reside in haloes with a typical mass of 2.8 x 10^{11} Msun. Combining this halo mass estimate with those for brighter samples of Kong et al. (2006), we find that the mass of dark haloes largely increases with K brightness, a measure of the stellar mass. Comparison with other galaxy populations suggests that faint sBzKs (K(AB)<23.2) and Lyman Break Galaxies at z ~ 2 are similar populations hosted by relatively low-mass haloes, while bright sBzKs (K(AB)<21) reside in haloes comparable to or more massive than those of Distant Red Galaxies and Extremely Red Objects. Using the extended Press-Schechter formalism, we predict that present-day descendants of haloes hosting sBzKs span a wide mass range depending on K brightness, from lower than that of the Milky Way up to those of richest clusters.

preprint2007arXiv

Results of Monitoring the Dramatically Variable C IV Mini-BAL System in the Quasar HS 1603+3820

We present six new and two previously published high-resolution spectra of the quasar HS 1603+3820 (z_em = 2.542) taken over an interval of 4.2 years (1.2 years in the quasar rest frame). The observations were made with the High-Dispersion Spectrograph on the Subaru telescope and Medium-Resolution Spectrograph on the Hobby-Eberly Telescope. The purpose was to study the narrow absorption lines (NALs). We use time variability as well as coverage fraction analysis to separate intrinsic absorption lines, which are physically related to the quasar, from intervening absorption lines. By fitting models to the line profiles, we derive the parameters of the respective absorbers as a function of time. Only the mini-BAL system at z_abs ~ 2.43 (v_off ~ 9,500 km/s) shows both partial coverage and time variability, although two NAL systems possibly show evidence of partial coverage. We find that all the troughs of the mini-BAL system vary in concert and its total equivalent width variations resemble those of the coverage fraction. However, no other correlations are seen between the variations of different model parameters. Thus, the observed variations cannot be reproduced by a simple change of ionization state nor by motion of a homogeneous parcel of gas across the cylinder of sight. We propose that the observed variations are a result of rapid continuum fluctuations, coupled with coverage fraction fluctuations caused by a clumpy screen of variable optical depth located between the continuum source and the mini-BAL gas. An alternative explanation is that the observed partial coverage signature is the result of scattering of continuum photons around the absorber, thus the equivalent width of the mini-BAL can vary as the intensity of the scattered continuum changes.

preprint2007arXiv

Searching for a Companion Star of Tycho's Type Ia Supernova with Optical Spectroscopic Observations

We report our first results of photometric and spectroscopic observations for Tycho's supernova remnant (SNR Tycho) to search for the companion star of a type Ia supernova (SN Ia). From photometric observations using Suprime-Cam on the Subaru Telescope, we have picked up stars brighter than 22 mag (in $V$-band) for spectroscopy, which are located within a circular region with the radius of 30" around the center of SNR Tycho. If the ejecta of young supernova remnants, such as SNR Tycho, have a sufficient amount of Fe I, we should be able to detect absorption lines at 3720 Åand 3860 Åassociated with transitions from the ground state of Fe I in the spectrum of the companion star. To identify the companion star of a SN Ia using these characteristic absorption lines of Fe I, we made optical low-resolution spectroscopy of their targets using FOCAS on the Subaru Telescope. In our spectroscopic observations, we obtained spectra of 17 stars in the SNR Tycho region and compare them with template stellar spectra. We detect significant absorption lines from two stars at 3720 Å. Since widths of their absorption lines are broad, it is likely that the detected absorptions are due to Fe I in the expanding ejecta of SNR Tycho. However, none of stars exhibits a clear red wing in the observed profiles of the absorption, though a star in the background of the SNR should show it. Hence, we suggest another interpretation that the detected absorption lines might be caused by the peculiarity of stars. A star named Tycho(E) has the absorption line at 3720 Åand its projected position is close to the center of SNR Tycho. Based on our observations, Tycho(E) is a new candidate as the companion star of Tycho's supernova.

preprint2007arXiv

The Remarkable 60 x 2 kpc Optical Filament Associated with a Poststarburst Galaxy in Coma Cluster

In the deep narrow band image of the Coma Cluster taken with Suprime-Cam of the Subaru telescope, we found an extremely long and narrow (~ 60 kpc x 2 kpc) Halpha emitting region associated with a poststarburst galaxy (D100). Follow up spectroscopy shows that the region has the same redshift as D100. The surface brightness of the region is typically 25 mag(AB) arcsec^-2 in Halpha, which corresponds to 0.5--4x10^-17 erg s^-1 cm^-2 arcsec^-2. We set two possible explanations for the origin of the region; gas stripped off from a merged dwarf, or gas stripped off from D100 by ram pressure. Either scenario has difficulty to fully explain all the observed characteristics of the region.

preprint2005arXiv

A survey of NB921 dropouts in the Subaru Deep Field

In order to search for high-redshift galaxies beyond $z = 6.6$ in the Subaru Deep Field, we have investigated NB921-dropout galaxies where NB921 is the narrowband filter centered at 919.6 nm with FWHM of 13.2 nm for the Suprime-Cam on the Subaru Telescope. There are no secure NB921-dropout candidates brighter than $z^\prime = 25.5$. Based on this result, we discuss the UV luminosity function of star-forming galaxies at $z > 6.6$.

preprint2005arXiv

Discovery of a Strongly Lensed Galaxy at z=3.9 behind a z=0.83 Galaxy Cluster

We report the discovery and spectroscopic confirmation of three gravitationally-lensed images of a galaxy at z=3.9 in the background of a distant, rich cluster of galaxies at z=0.83, on the basis of observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estimates from deep, high-resolution images by Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The lens model can account simultaneously for the observed image configuration and the flux ratio of the closer pair located close to the critical curve. The parities of the three images are also consistent with the lensing hypothesis. Since this galaxy is apparently bright (i'_{AB} \sim 23.7) for its redshift due to the magnification, it serves as a good high redshift target on which we can make extensive and detailed studies based on multi-wavelength imaging and spectroscopy.

preprint2005arXiv

New High-Redshift Galaxies at $z$ = 5.8 - 6.5 in the Subaru Deep Field

In order to search for high-redshift galaxies around $z \sim 6$ in the Subaru Deep Field, we have investigated NB816-dropout galaxies where NB816 is the narrowband filter centered at 815nm with FWHM of 12.5 nm for the Suprime-Cam on the Subaru Telescope. Since the NB816 imaging is so deep we can detect 10 well-defined NB816-dropout galaxies that we identify to lie at $z$ = 5.8 -- 6.5. We discuss their observational properties.

preprint2003arXiv

Are Two z~6 Quasars Gravitationally Lensed ?

Several high-z (z > 5.7) quasars have been found in the course of Sloan Digital Sky Survey. The presence of such very high-z quasars is expected to give constraints on early structure formation. On one hand, it is suggested that these most luminous objects at high redshift are biased toward the highly magnified objects by gravitational lensing. To clarify the effect of gravitational lensing on the high-z quasars, we began the imaging survey of intervening lensing galaxies. Indeed our previous optical image showed that SDSSp J104433.04+012502.2 at z=5.74 is gravitationally magnified by a factor 2. In this paper, we report our new optical imaging of other two high-z quasars, SDSSp J103027.10+052455.0 at z=6.28 and SDSSp J130608.26+035626.3 at z=5.99. Since we find neither intervening galaxy nor counter image with i^{\prime} < 25.4-25.8 around each quasar, we conclude that they are not strongly magnified regardless that a lens galaxy is dusty.

preprint2003arXiv

Dusty ERO Search behind Two Massive Clusters

We performed deep K'-band imaging observations of 2 massive clusters, MS 0451.6-0305 at z = 0.55 and MS 0440.5+0204 at z = 0.19, for searching counterparts of the faint sub-mm sources behind these clusters, which would provide one of the deepest extremely red object(ERO) samples. Comparing our near-infrared images with optical images taken by the Hubble Space Telescope and by the Subaru Telescope, we identified 13 EROs in these fields. The sky distributions of EROs are consistent with the previous results, that there is a sign of strong clustering among detected EROs. Also, the surface density with corrected lensing amplification factors in both clusters are in good agreement with that derived from previous surveys. We found 7 EROs and 3 additional very red objects in a small area (\sim 0.6 arcmin^2) of the MS 0451.6-0305 field around an extended SCUBA source. Many of their optical and near-infrared colors are consistent with dusty star-forming galaxies at high redshifts(z \sim 1.0-4.0), and they may be constituting a cluster of dusty starburst galaxies and/or lensed star-forming galaxies at high redshift. Their red J-K' colors and faint optical magnitudes suggest they are relatively old massive stellar systems with ages(>300 Mega years) suffering from dust obscuration. We also found a surface-density enhancement of EROs around the SCUBA source in the MS 0440.5+0204 field.

preprint2002arXiv

Decomposition of the Superwind in M82

We present new optical images ($B$, $V$, and H$α$) of the archetypical starburst/superwind galaxy M82 obtained with the 8.2 m Subaru Telescope to reveal new detailed structures of the superwind-driven nebula and the high-latitude dark lanes. The emission-line nebula is decomposed into (1) a ridge-dominated component comprising numerous filament/loop sub-structures whose overall morphology appears as a pair of narrow cylinders, and (2) a diffuse component extended over much wider opening angle from the nucleus. We suggest that these two components have different origins. The ridge-dominated component appears as a pair of cylinders rather than a pair of cones. Since this morphological property is similar to that of hot plasma probed by soft X-ray, this component seems to surround the hot plasma. On the other hand, the diffuse component may arise from dust grains which scatter stellar light from the galaxy. Since inner region of this component is seen over the prominent ^^ ^^ X"-shaped dark lanes streaming out from the nuclear region and they can be reproduced as a conical distribution of dust grains, there seems to be a dusty cold outflow as well as the hot one probed by soft X-ray and shock-excited optical emission lines. If this is the case, the presence of such high-latitude dust grains implies that neutral gaseous matter is also blown out during the course of the superwind activity.