Source author record

Kazuhiro Shimasaku

Kazuhiro Shimasaku 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

72works
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

72 published item(s)

preprint2022arXiv

A systematic search for galaxy protocluster cores at the transition epoch of their star formation activity

The redshift of $z\sim1.5$ is the transition epoch of protoclusters (PCs) from the star-forming phase into the quenching phase, and hence an appropriate era to investigate the build up of the quenched population. We define a `core' as the most massive halo in a given PC, where environmental effects are likely to work most effectively, and search for cores at $1<z<1.5$. We use a photometric redshift catalogue of a wide (effective area of $\sim22.2\,\mathrm{deg}^{2}$) and deep ($i\sim26.8\,\mathrm{mag}$) optical survey with Subaru Hyper-Suprime Cam. Regarding galaxies with $\log(M_{*}/M_{\odot})>11.3$ as the central galaxies of PC cores, we estimate their average halo mass by clustering analysis and find it to be $\log(M_\mathrm{h}/M_{\odot})\sim13.7$. An expected mass growth by the IllustrisTNG simulation and the observed overdensities around them suggest that the PC cores we find are progenitors of present-day clusters. Classifying our galaxy sample into red and blue galaxies, we calculate the stellar mass function (SMF) and the red galaxy fraction. The SMFs in the PC cores are more-top heavy than field, implying early high-mass galaxy formation and disruption of low-mass galaxies. We also find that the red fraction increases with stellar mass, consistent with stellar-mass dependent environmental quenching recently found at $z>1$. Interestingly, although the cores with red and blue centrals have similar halo masses, only those with red centrals show a significant red fraction excess compared to the field, suggesting a conformity effect. Some observational features of PC cores may imply that the conformity is caused by assembly bias.

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

Detection of anisotropic satellite quenching in galaxy clusters up to $z\sim1$

Satellite galaxies in the cluster environment are more likely to be quenched than galaxies in the general field. Recently, it has been reported that satellite galaxy quenching depends on the orientation relative to their central galaxies: satellites along the major axis of centrals are more likely to be quenched than those along the minor axis. In this paper, we report a detection of such anisotropic quenching up to $z\sim1$ based on a large optically-selected cluster catalogue constructed from the Hyper Suprime-Cam Subaru Strategic Program. We calculate the quiescent satellite galaxy fraction as a function of orientation angle measured from the major axis of central galaxies and find that the quiescent fractions at $0.25<z<1$ are reasonably fitted by sinusoidal functions with amplitudes of a few percent. Anisotropy is clearer in inner regions ($<r_\mathrm{200m}$) of clusters and not significant in cluster outskirts ($>r_\mathrm{200m}$). We also confirm that the observed anisotropy cannot be explained by differences in local galaxy density or stellar mass distribution along the two axes. Quiescent fraction excesses between the two axes suggest that the quenching efficiency contributing to the anisotropy is almost independent of stellar mass, at least down to our stellar mass limit of $M_{*}=1\times10^{10}\,M_{\odot}$. Finally, we argue that the physical origins of the observed anisotropy should have shorter quenching timescales than $\sim1\,\mathrm{Gyr}$, like ram-pressure stripping, because, for anisotropic quenching to be observed, satellites must be quenched before their initial orientation angles are significantly changed.

preprint2022arXiv

MUSUBI (MegaCam Ultra-deep Survey: $u^\ast$-Band Imaging)-Data for the COSMOS and SXDS Fields

The Subaru Hyper Suprime-Cam (HSC) Strategic Survey is the latest-generation multi-band optical imaging survey for galaxy evolution and structure formation. The "Ultra-Deep" component of the HSC survey provides $grizy$ broad-band images over $\sim3.4$ deg$^2$ to detection limits of $\sim26$-28 AB, along with narrow-band images, in the COSMOS and the SXDS fields. These images provide an unprecedented combination of depths and area coverage, for the studies galaxies up to $z\sim7$. However, the lack of coverage at $<4000$ Ang implies incomplete sampling of the rest-frame UV at $z\lesssim 3$, which is critically needed for understanding the buildup of stellar mass in the later cosmic time. We conducted a multi-year CFHT $u^\ast$-band imaging campaign in the two HSC Ultra-Deep fields with CFHT MegaCam. By including shallower archival data, we reach 5-$σ$ depths of $u^\ast=28.1$ and 28.4 (AB) at the centers of the COSMOS and SXDS fields, respectively, and $u^\ast=27.7$ and 27.8 in the central 1 deg$^2$ fields. The image quality is $\gtrsim0.90$ arcsec, fairly good for the $u^\ast$ band. Both the photometric and astrometric quality of our data are excellent. We show that the combination of our $u^\ast$-band and HSC data can lead to high-quality photometric redshifts at $z=0$-3, and robust measurements of rest-frame UV on galaxies at $0.4<z<0.6$ for distinguishing green-valley galaxies from star-forming and quiescent galaxies. We publicly release our reduced $u^\ast$-band images and reference catalogs that can be used readily for scientific studies.

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

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: Constraints on the Ly$α$ luminosity function and cosmic reionization at $z=7.3$ with Subaru/Hyper Suprime-Cam

The Ly$α$ luminosity function (LF) of Ly$α$ emitters (LAEs) has been used to constrain the neutral hydrogen fraction in the intergalactic medium (IGM) and thus the timeline of cosmic reionization. Here we present the results of a new narrow-band imaging survey for $z=7.3$ LAEs in a large area of $\sim 3\ \mathrm{deg}^2$ with Subaru/Hyper Suprime-Cam. No LAEs are detected down to $L_{\mathrm{Ly}α}\simeq 10^{43.2}\ \mathrm{erg\ s^{-1}}$ in an effective cosmic volume of $\sim 2\times 10^6$ Mpc$^3$, placing an upper limit to the bright part of the $z=7.3$ Ly$α$ LF for the first time and confirming a decrease in bright LAEs from $z=7.0$. By comparing this upper limit with the Ly$α$ LF in the case of the fully ionized IGM, which is predicted using an observed $z=5.7$ Ly$α$ LF on the assumption that the intrinsic Ly$α$ LF evolves in the same way as the UV LF, we obtain the relative IGM transmission $T^\mathrm{IGM}_{\mathrm{Ly}α}(7.3)/T^\mathrm{IGM}_{\mathrm{Ly}α}(5.7)<0.77$, and then the volume-averaged neutral fraction $x_\mathrm{HI}(7.3)>0.28$. Cosmic reionization is thus still ongoing at $z=7.3$, being consistent with results from other $x_\mathrm{HI}$ estimation methods. A similar analysis using literature Ly$α$ LFs finds that at $z=6.6$ and 7.0 the observed Ly$α$ LF agrees with the predicted one, consistent with full ionization.

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.

preprint2020arXiv

A Systematic Search for Galaxy Proto-Cluster Cores at $z\sim 2$

A proto-cluster core is the most massive dark matter halo (DMH) in a given proto-cluster. To reveal the galaxy formation in core regions, we search for proto-cluster cores at $z\sim 2$ in $\sim 1.5\, \mathrm{deg}^{2}$ of the COSMOS field. Using pairs of massive galaxies ($\log(M_{*}/M_{\odot})\geq11$) as tracers of cores, we find 75 candidate cores, among which 54\% are estimated to be real. A clustering analysis finds that these cores have an average DMH mass of $2.6_{-0.8}^{+0.9}\times 10^{13}\, M_{\odot}$, or $4.0_{-1.5}^{+1.8}\, \times 10^{13} M_{\odot}$ after contamination correction. The extended Press-Schechter model shows that their descendant mass at $z=0$ is consistent with Fornax-like or Virgo-like clusters. Moreover, using the IllustrisTNG simulation, we confirm that pairs of massive galaxies are good tracers of DMHs massive enough to be regarded as proto-cluster cores. We then derive the stellar mass function (SMF) and the quiescent fraction for member galaxies of the 75 candidate cores. We find that the core galaxies have a more top-heavy SMF than field galaxies at the same redshift, showing an excess at $\log(M_{*}/M_{\odot})\gtrsim 10.5$. The quiescent fraction, $0.17_{-0.04}^{+0.04}$ in the mass range $9.0\leq \log(M_{*}/M_{\odot})\leq 11.0$, is about three times higher than that of field counterparts, giving an environmental quenching efficiency of $0.13_{-0.04}^{+0.04}$. These results suggest that stellar mass assembly and quenching are accelerated as early as at $z\sim 2$ in proto-cluster cores.

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

Implications of the mild gas motion found with Hitomi in the core of the Perseus cluster

Based mainly on X-ray observations, studies are made on interactions between the intra-cluster medium (ICM) in clusters of galaxies and their member galaxies. Through (magneto)hydrodynamic and gravitational channels, the moving galaxies are expected to drag the ICM around them, and transfer to the ICM some fraction of their dynamical energies on cosmological time scales. This hypothesis is in line with several observations, including the possible cosmological infall of galaxies towards the cluster center, found over redshifts of z~1 to z~0. Further assuming that the energy lost by the galaxies is first converted into ICM turbulence and then dissipated, this picture can explain the subsonic and uniform ICM turbulence, measured with Hitomi in the core region of the Perseus cluster. The scenario may also explain several other unanswered problems regarding clusters of galaxies, including what prevents the ICM from the expected radiative cooling, how the various mass components in nearby clusters have attained different radial distributions, and how a thermal stability is realized between hot and cool ICM components that co-exist around cD galaxies. This view is also considered to pertain to the general scenario of galaxy evolution, including their environmental effects.

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$.

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.

preprint2016arXiv

A Close Comparison between Observed and Modeled Lyα Lines for z ~ 2.2 Lyman Alpha Emitters

We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z = 2.2 with high-resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if any) with respect to the systemic redshift is Delta_v_Lya,r = 174+/- 19 km s-1 (Delta_v_Lya,b = -316+/-45 km s-1), which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities inferred from metal absorption lines in three individual and one stacked spectra are comparable to those of LBGs. The expanding shell model constructed by Verhamme et al. (2006) reproduces not only the Lya profiles but also other observed quantities including the outflow velocity and the FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear if we introduce additional Lya photons produced by gravitational cooling. We show that the small Delta_v_Lya,r values of our sample can be explained by low neutral-hydrogen column densities of log(NHI) = 18.9 cm-2 on average. This value is more than one order of magnitude lower than those of LBGs but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas masses. This result suggests that low NHI values, giving reduced numbers of resonant scattering of Lya photons, are the key to the strong Lya emission of LAEs.

preprint2016arXiv

Bright and Faint Ends of Ly$α$ Luminosity Functions at $\textit{z} = 2$ Determined by the Subaru Survey: Implications for AGN, Magnification Bias, and ISM HI Evolution

We present the Lya luminosity functions (LFs) derived by our deep Subaru narrowband survey that identifies a total of 3,137 Lya emitters (LAEs) at $z = 2.2$ in five independent blank fields. The sample of these LAEs is the largest, to date, and covers a very wide Lya luminosity range of $\log L_{Lyα} = 41.7-44.4$ erg s$^{-1}$. We determine the Lya LF at $z = 2.2$ with unprecedented accuracies, and obtain the best-fit Schechter parameters of $L^{*}_{Lyα} = 5.29^{+1.67}_{-1.13} \times 10^{42}$ erg s$^{-1}$, $ϕ^{*}_{Lyα} = 6.32^{+3.08}_{-2.31} \times 10^{-4}$ Mpc$^{-3}$, and $α= -1.75^{+0.10}_{-0.09}$ showing a steep faint-end slope. We identify a significant hump at the LF bright end ($\log L_{Lyα} > 43.4$ erg s$^{-1}$). Because all of the LAEs in the bright-end hump have (a) bright counterpart(s) either in the X-ray, UV, or radio data, this bright-end hump is not made by gravitational lensing magnification bias but AGNs. These AGNs allow us to derive the AGN UV LF at $z \sim 2$ down to the faint magnitude limit of $M_{UV} \simeq -22.5$, and to constrain the faint-end slope of AGN UV LF, $α_{AGN}=-1.2 \pm 0.1$, that is flatter than those at $z > 4$. Based on the Lya and UV LFs from our and previous studies, we find the increase of Lya escape fraction $f^{Lyα}_{esc}$ from $z \sim 0$ to $6$ by two orders of magnitude. This large $f^{Lyα}_{esc}$ increase can be explained neither by the evolution of stellar population nor outflow alone, but the evolution of neutral hydrogen HI density in inter-stellar medium that enhances dust attenuation for Lya by resonance scattering. Our uniform expanding shell models suggest that the typical HI column density decreases from $N_{HI} \sim 7 \times 10^{19}$ ($z \sim 0$) to $\sim 1 \times 10^{18}$ cm$^{-2}$ ($z \sim 6$) to explain the large $f^{Lyα}_{esc}$ increase.

preprint2016arXiv

Evolution of Stellar-to-Halo Mass Ratio at z=0-7 Identified by Clustering Analysis with the Hubble Legacy Imaging and Early Subaru/Hyper Suprime-Cam Survey Data

We present clustering analysis results from 10,381 Lyman break galaxies (LBGs) at z~ 4-7, identified in the Hubble legacy deep imaging and new complimentary large-area Subaru/Hyper Suprime-Cam data. We measure the angular correlation functions (ACFs) of these LBGs at z~4, 5, 6, and 7, and fit these measurements using halo occupation distribution (HOD) models that provide an estimate of halo masses, M_h~(1-20)x10^11 Msun. Our M_h estimates agree with those obtained by previous clustering studies in a UV-magnitude vs. M_h plane, and allow us to calculate stellar-to-halo mass ratios (SHMRs) of LBGs. By comparison with the z~0 SHMR, we identify evolution of the SHMR from z~0 to z~4, and z~4 to z~7 at the >98% confidence levels. The SHMR decreases by a factor of ~2 from z~0 to 4, and increases by a factor of ~4 from z~4 to 7. We compare our SHMRs with results of a hydrodynamic simulation and a semi-analytic model, and find that these theoretical studies do not predict the SHMR increase from z~4 to 7. We obtain the baryon conversion efficiency (BCE) of LBGs at z~4, and find that the BCE increases with increasing dark matter halo mass. Finally, we compare our clustering+HOD estimates with results from abundance matching techniques, and conclude that the M_h estimates of the clustering+HOD analyses agree with those of the simple abundance matching within a factor of 3, and that the agreement improves when using more sophisticated abundance matching techniques that include subhalos, incompleteness, and/or evolution in the star formation and stellar mass functions.

preprint2016arXiv

Ly$α$ Emitters with Very Large Ly$α$ Equivalent Widths, EW$_{\rm 0}$(Ly$α$) $\simeq 200-400$ Å, at $z\sim 2$

We present physical properties of spectroscopically confirmed Ly$α$ emitters (LAEs) with very large rest-frame Ly$α$ equivalent widths EW$_{\rm 0}$(Ly$α$). Although the definition of large EW$_{\rm 0}$(Ly$α$) LAEs is usually difficult due to limited statistical and systematic uncertainties, we identify six LAEs selected from $\sim 3000$ LAEs at $z\sim 2$ with reliable measurements of EW$_{\rm 0}$ (Ly$α$) $\simeq 200-400$ Å given by careful continuum determinations with our deep photometric and spectroscopic data. These large EW$_{\rm 0}$(Ly$α$) LAEs do not have signatures of AGN, but notably small stellar masses of $M_{\rm *} = 10^{7-8}$ $M_{\rm \odot}$ and high specific star-formation rates (star formation rate per unit galaxy stellar mass) of $\sim 100$ Gyr$^{-1}$. These LAEs are characterized by the median values of $L({\rm Lyα})=3.7\times 10^{42}$ erg s$^{-1}$ and $M_{\rm UV}=-18.0$ as well as the blue UV continuum slope of $β= -2.5\pm0.2$ and the low dust extinction $E(B-V)_{\rm *} = 0.02^{+0.04}_{-0.02}$, which indicate a high median Ly$α$ escape fraction of $f_{\rm esc}^{\rm Lyα}=0.68\pm0.30$. This large $f_{\rm esc}^{\rm Lyα}$ value is explained by the low {\sc Hi} column density in the ISM that is consistent with FWHM of the Ly$α$ line, ${\rm FWHM (Lyα)}=212\pm32$ km s$^{-1}$, significantly narrower than those of small EW$_{\rm 0}$(Ly$α$) LAEs. Based on the stellar evolution models, our observational constraints of the large EW$_{\rm 0}$ (Ly$α$), the small $β$, and the rest-frame He{\sc ii} equivalent width imply that at least a half of our large EW$_{\rm 0}$(Ly$α$) LAEs would have young stellar ages of $\lesssim 20$ Myr and very low metallicities of $Z<0.02 Z_\odot$ regardless of the star-formation history.

preprint2016arXiv

Precise Strong Lensing Mass Modeling of Four Hubble Frontier Fields Clusters and a Sample of Magnified High-Redshift Galaxies

We conduct precise strong lensing mass modeling of four ${\it Hubble}$ Frontier Fields (HFF) clusters, Abell$~$2744, MACS$~$J0416.1$-$2403, MACS$~$J0717.5$+$3745, and MACS$~$J1149.6$+$2223, for which HFF imaging observations are completed. We construct a refined sample of more than 100 multiple images for each cluster by taking advantage of the full depth HFF images, and conduct mass modeling using the ${\small \rm GLAFIC}$ software, which assumes simply parametrized mass distributions. Our mass modeling also exploits a magnification constraint from the lensed Type Ia supernova HFF14Tom for Abell$~$2744 and positional constraints from the multiple images S1$-$S4 of the lensed supernova SN Refsdal for MACS$~$J1149.6$+$2223. We find that our best-fitting mass models reproduce the observed image positions with RMS errors of $\sim 0.4$ arcsec, which are smaller than RMS errors in previous mass modeling that adopted similar numbers of multiple images. Our model predicts a new image of SN Refsdal with a relative time delay and magnification that are fully consistent with a recent detection of reappearance. We then construct catalogs of $z\sim 6-9$ dropout galaxies behind the four clusters and estimate magnification factors for these dropout galaxies with our best-fitting mass models. The dropout sample from the four cluster fields contains $\sim 120$ galaxies at $z\gtrsim 6$, about 20 of which are predicted to be magnified by a factor of more than 10. Some of the high-redshift galaxies detected in the HFF have lensing-corrected magnitudes of $M_{\rm UV}\sim -15$ to $-14$. Our analysis demonstrates that the HFF data indeed offer an ideal opportunity to study faint high-redshift galaxies. All lensing maps produced from our mass modeling will be made available on the STScI website.

preprint2015arXiv

First Infrared-based implications for the dust attenuation and star formation of typical Ly$α$ emitters

By stacking publicly available deep Spitzer/MIPS 24 $μ$m and Herschel/PACS images for 213 $z \simeq 2.18$ Ly$α$ Emitters (LAEs) in GOODS-South, we obtain a strong upper limit to the IR luminosity of typical LAEs and discuss their attenuation curve for the first time. The $3σ$ upper limit $L_{\rm TIR}^{3σ}= 1.1 \times 10^{10} L_\odot$, determined from the MIPS data providing the lowest limit, gives $IRX \equiv L_{\rm TIR}/L_{\rm UV} \leq 2.2$. Here we assume that the local calibration between the 8 $μ$m emission and the dust SED shape and metallicity applies at high redshifts and that our LAEs have low metallicities as suggested by previous studies. The inferred escape fractions of Ly$α$, $16$--$37%$, and UV continuum, $\ge 44%$, are higher than the cosmic averages at the same epoch. The SMC attenuation curve is consistent with the IRX and the UV slope $β= -1.4^{+0.2}_{-0.2}$ of our stacked LAE, while the Meurer's relation (Calzetti curve) predicts a 3.8 times higher IRX; we also discuss the validity of PACS-based $L_{\rm TIR}^{3σ}$ allowing the Meurer's relation. SED fitting using the Calzetti curve also gives a $\sim 10$ times higher SFR than from the $L_{\rm TIR}^{3σ}$ and $L_{\rm UV}$. With $M_{\star}=6.3^{+0.8}_{-2.0} \times10^8 \mathrm{M_{\odot}}$, our LAEs lie on a lower-mass extrapolation of the star formation main sequence at $z \sim 2$, suggesting that the majority of $z \sim 2$ LAEs are mildly star forming with relatively old ages of $\sim 200$ Myr. The faint $L_{\rm TIR}^{3σ}$ implies that LAEs contribute little to the faint ($\gtrsim 100 μ$Jy) submm number counts by ALMA.

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.

preprint2015arXiv

Statistical Properties of Diffuse Lyman-alpha Halos around Star-forming Galaxies at z~2

We present statistical properties of diffuse Lyman-alpha halos (LAHs) around high-$z$ star-forming galaxies with large Subaru samples of Lyman-alpha emitters (LAEs) at $z=2.2$. We make subsamples defined by the physical quantities of LAEs' central Lyman-alpha luminosities, UV magnitudes, Lyman-alpha equivalent widths, and UV slopes, and investigate LAHs' radial surface brightness (SB) profiles and scale lengths $r_n$ as a function of these physical quantities. We find that there exist prominent LAHs around LAEs with faint Lyman-alpha luminosities, bright UV luminosities, and small Lyman-alpha equivalent widths in cumulative radial Lyman-alpha SB profiles. We confirm this trend with the anti-correlation between $r_n$ and Lyman-alpha luminosities (equivalent widths) based on the Spearman's rank correlation coefficient that is $ρ=-0.9$ ($-0.7$) corresponding to the $96\%$ ($93\%$) confidence level, although the correlation between $r_n$ and UV magnitudes is not clearly found in the rank correlation coefficient. Our results suggest that LAEs with properties similar to typical Lyman-break galaxies (with faint Lyman-alpha luminosities and small equivalent widths) possess more prominent LAHs. We investigate scenarios for the major physical origins of LAHs with our results, and find that the cold stream scenario is not preferred, due to the relatively small equivalent widths up to $77$Å in LAHs that include LAEs' central components. There remain two possible scenarios of Lyman-alpha scattering in circum-galactic medium and satellite galaxies that cannot be tested with our observational data.

preprint2015arXiv

The sizes of $z\sim6-8$ lensed galaxies from the Hubble Frontier Fields Abell 2744 data

We investigate sizes of $z\sim6-8$ dropout galaxies using the complete data of the Abell 2744 cluster and parallel fields in the Hubble Frontier Fields program. By directly fitting light profiles of observed galaxies with lensing-distorted Sérsic profiles on the image plane with the \texttt{glafic} software, we accurately measure intrinsic sizes of 31 $z\sim6-7$ and eight $z\sim8$ galaxies, including those as faint as $M_{\mathrm{UV}}\simeq-16.6$. We find that half-light radii $r_\mathrm{e}$ positively correlates with UV luminosity at each redshift, although the correlation is not very tight. Largest ($r_\mathrm{e}>0.8$ kpc) galaxies are mostly red in UV color while smallest ($r_\mathrm{e} < 0.08$ kpc) ones tend to be blue. We also find that galaxies with multiple cores tend to be brighter. Combined with previous results at $2.5\lesssim z\lesssim12$, our result confirms that the average $r_{\mathrm{e}}$ of bright ($(0.3-1)L^\star_{z=3}$) galaxies scales as $r_{\mathrm{e}}\propto(1+z)^{-m}$ with $m=1.24\pm0.1$. We find that the ratio of $r_\mathrm{e}$ to virial radius is virtually constant at $3.3\pm0.1\%$ over a wide redshift range, where the virial radii of hosting dark matter halos are derived based on the abundance matching. This constant ratio is consistent with the disk formation model by Mo et al. (1998) with $j_\mathrm{d}\sim m_\mathrm{d}$, where $j_\mathrm{d}$ and $m_\mathrm{d}$ are the fractions of the angular momentum and mass within halos confined in the disks. A comparison with various types of local galaxies indicates that our galaxies are most similar to circumnuclear star-forming regions of barred galaxies in the sense that a sizable amount of stars are forming in a very small area.

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

Accelerated Evolution of Ly$α$ Luminosity Function at $\textit{z} \gtrsim 7$ Revealed by the Subaru Ultra-Deep Survey for Ly$α$ Emitters at $\textit{z}=7.3$

We present the ultra-deep Subaru narrowband imaging survey for Lya emitters (LAEs) at $z=7.3$ in SXDS and COSMOS fields with a total integration time of 106 hours. Exploiting our new sharp bandwidth filter, NB101, installed on Suprime-Cam, we have reached $L(Lya)=2.4\times10^{42} \ erg \ s^{-1}$ ($5σ$) for $z=7.3$ LAEs, about 4 times deeper than previous Subaru $z \gtrsim 7$ studies, which allows us to reliably investigate evolution of Lya luminosity function (LF), for the first time, down to the luminosity limit same as those of Subaru $z=3.1-6.6$ LAE samples. Surprisingly, we only find three and four LAEs in SXDS and COSMOS fields, respectively, while one expects a total of $\sim 65$ LAEs by our survey in the case of no Lya LF evolution from $z=6.6$ to $7.3$.We identify a decrease of Lya LF from $z=6.6$ to $7.3$ at the $>90\%$ confidence level from our $z=7.3$ Lya LF.Moreover, the evolution of Lya LF is clearly accelerated at $z>6.6$ beyond the measurement uncertainties including cosmic variance. Because no such accelerated evolution of UV-continuum LF or cosmic star-formation rate (SFR) is found at $z\sim 7$, but suggested only at $z>8$ (Oesch et al. 2013, Bouwens et al. 2014), this accelerated Lya LF evolution is explained by physical mechanisms different from a pure SFR decrease but related to Lya production and escape in the process of cosmic reionization. Because a simple accelerating increase of IGM neutral hydrogen absorbing Lya would not reconcile with Thomson scattering optical depth measurements from WMAP and Planck, our findings may support new physical pictures suggested by recent theoretical studies, such as the existence of HI clumpy clouds within cosmic ionized bubbles selectively absorbing Lya and the large ionizing photon escape fraction of galaxies making weak Lya emission.

preprint2014arXiv

Diffuse Lyman-alpha Halos around Galaxies at z=2.2-6.6: Implications for Galaxy Formation and Cosmic Reionization

We present diffuse Lyman-alpha halos (LAHs) identified in the composite Subaru narrowband images of 100-3600 Lyman-alpha emitters (LAEs) at z=2.2, 3.1, 3.7, 5.7, and 6.6. First, we carefully examine potential artifacts mimicking LAHs that include a large-scale point-spread function (PSF) made by instrumental and atmospheric effects. Based on our critical test with composite images of non-LAE samples whose narrowband-magnitude and source-size distributions are the same as our LAE samples, we confirm that no artifacts can produce a diffuse extended feature similar to our LAHs. After this test, we measure the scale lengths of exponential profile for the LAHs estimated from our z=2.2-6.6 LAE samples of L(Lyman-alpha) > 2 x 10^42 erg s^-1. We obtain the scale lengths of ~ 5-10 kpc at z=2.2-5.7, and find no evolution of scale lengths in this redshift range beyond our measurement uncertainties. Combining this result and the previously-known UV-continuum size evolution, we infer that the ratio of LAH to UV-continuum sizes is nearly constant at z=2.2-5.7. The scale length of our z=6.6 LAH is larger than 5-10 kpc just beyond the error bar, which is a hint that the scale lengths of LAHs would increase from z=5.7 to 6.6. If this increase is confirmed by future large surveys with significant improvements of statistical and systematical errors, this scale length change at z > 6 would be a signature of increasing fraction of neutral hydrogen scattering Lyman-alpha photons, due to cosmic reionization.

preprint2014arXiv

Hubble Frontier Fields First Complete Cluster Data: Faint Galaxies at $z\sim 5-10$ for UV Luminosity Functions and Cosmic Reionization

We present the comprehensive analyses of faint dropout galaxies up to $z\sim10$ with the first full-depth data set of Abell 2744 lensing cluster and parallel fields observed by the Hubble Frontier Fields (HFF) program. We identify $54$ dropouts at $z\sim5-10$ in the HFF fields, and enlarge the size of $z\sim9$ galaxy sample obtained to date. Although the number of highly magnified ($μ\sim10$) galaxies is small due to the tiny survey volume of strong lensing, our study reaches the galaxies' intrinsic luminosities comparable to the deepest-field HUDF studies. We derive UV luminosity functions with these faint dropouts, carefully evaluating the combination of observational incompleteness and lensing effects in the image plane by intensive simulations including magnification, distortion, and multiplication of images, with the evaluations of mass model dependences. Our results confirm that the faint-end slope, $α$, is as steep as $-2$ at $z\sim6-8$, and strengthen the evidence of the rapid decrease of UV luminosity densities, $ρ_\mathrm{UV}$, at $z>8$ from the large $z\sim9$ sample. We examine whether the rapid $ρ_\mathrm{UV}$ decrease trend can reconcile with the large Thomson scattering optical depth, $τ_\mathrm{e}$, measured by CMB experiments allowing a large space of free parameters such as average ionizing photon escape fraction and stellar-population dependent conversion factor. No parameter set can reproduce both the rapid $ρ_\mathrm{UV}$ decrease and the large $τ_\mathrm{e}$. It is possible that the $ρ_\mathrm{UV}$ decrease moderates at $z\gtrsim11$, that the free parameters significantly evolve towards high-$z$, or that there exist additional sources of reionization such as X-ray binaries and faint AGNs.

preprint2014arXiv

Large-scale environment of z~5.7 CIV absorption systems I: projected distribution of galaxies

Metal absorption systems are products of star formation. They are believed to be associated with massive star forming galaxies, which have significantly enriched their surroundings. To test this idea with high column density CIV absorption systems at z~5.7, we study the projected distribution of galaxies and characterise the environment of CIV systems in two independent quasar lines-of-sight: J103027.01+052455.0 and J113717.73+354956.9. Using wide field photometry (~80x60h$^{-1}$ comoving Mpc), we select bright (Muv(1350Å)<-21.0 mag.) Lyman break galaxies (LBGs) at z~5.7 in a redshift slice Δz~0.2 and we compare their projected distribution with z~5.7 narrow-band selected Lyman alpha emitters (LAEs, Δz~0.08). We find that the CIV systems are located more than 10h$^{-1}$ projected comoving Mpc from the main concentrations of LBGs and no candidate is closer than ~5h$^{-1}$ projected comoving Mpc. In contrast, an excess of LAEs -lower mass galaxies- is found on scales of ~10h$^{-1}$ comoving Mpc, suggesting that LAEs are the primary candidates for the source of the CIV systems. Furthermore, the closest object to the system in the field J1030+0524 is a faint LAE at a projected distance of 212h$^{-1}$ physical kpc. However, this work cannot rule out undiscovered lower mass galaxies as the origin of these absorption systems. We conclude that, in contrast with lower redshift examples (z<3.5), strong CIV absorption systems at z~5.7 trace low-to-intermediate density environments dominated by low-mass galaxies. Moreover, the excess of LAEs associated with high levels of ionizing flux agrees with the idea that faint galaxies dominate the ionizing photon budget at this redshift.

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.

preprint2014arXiv

What is the Physical Origin of Strong Lya Emission? I. Demographics of Lya Emitter Structures

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

preprint2014arXiv

What is the physical origin of strong Lya emission? II. Gas Kinematics and Distribution of Lya Emitters

We present a statistical study of velocities of Lya, interstellar (IS) absorption, and nebular lines and gas covering fraction for Lya emitters (LAEs) at z~2. We make a sample of 22 LAEs with a large Lya equivalent width (EW) of > 50A based on our deep Keck/LRIS observations, in conjunction with spectroscopic data from the Subaru/FMOS program and the literature. We estimate the average velocity offset of Lya from a systemic redshift determined with nebular lines to be dv_Lya=234+-9 km s-1. Using a Kolmogorv-Smirnov test, we confirm the previous claim of Hashimoto et al. (2013) that the average dv_Lya of LAEs is smaller than that of LBGs. Our LRIS data successfully identify blue-shifted multiple IS absorption lines in the UV continua of four LAEs on an individual basis. The average velocity offset of IS absorption lines from a systemic redshift is dv_IS=204+-27 km s-1, indicating LAE's gas outflow with a velocity comparable to typical LBGs. Thus, the ratio, R^Lya_ IS = dv_Lya/dv_IS of LAEs, is around unity, suggestive of low impacts on Lya transmission by resonant scattering of neutral hydrogen in the IS medium. We find an anti-correlation between Lya EW and the covering fraction, f_c, estimated from the depth of absorption lines, where f_c is an indicator of average neutral hydrogen column density, N_HI. The results of our study support the idea that N_HI is a key quantity determining Lya emissivity.

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

An Intensely Star-Forming Galaxy at z~7 with Low Dust and Metal Content Revealed by Deep ALMA and HST Observations

We report deep ALMA observations complemented with associated HST imaging for a luminous (m_uv=25) galaxy, `Himiko', at a redshift z=6.595. The galaxy is remarkable for its high star formation rate, 100 Mo/yr, securely estimated from our deep HST and Spitzer photometry, and the absence of any evidence for strong AGN activity or gravitational lensing magnification. Our ALMA observations probe an order of magnitude deeper than previous IRAM observations, yet fail to detect a 1.2mm dust continuum, indicating a flux <52uJy comparable with or weaker than that of local dwarf irregulars with much lower star formation rates. We likewise provide a strong upper limit for the flux of [CII] 158um, L([CII]) < 5.4x10^7 Lo, a diagnostic of the hot interstellar gas often described as a valuable probe for early galaxies. In fact, our observations indicate Himiko lies off the local L([CII]) - star formation rate scaling relation by a factor of more than 30. Both aspects of our ALMA observations suggest Himiko is an unique object with a very low dust content and perhaps nearly primordial interstellar gas. Our HST images provide unique insight into the morphology of this remarkable source, highlighting an extremely blue core of activity and two less extreme associated clumps. Himiko is undergoing a triple major merger event whose extensive ionized nebula of Lyman alpha emitting gas, discovered in our earlier work with Subaru, is powered by star formation and the dense circum-galactic gas. We are likely witnessing an early massive galaxy during a key period of its mass assembly close to the end of the reionization era.

preprint2013arXiv

First Spectroscopic Evidence for High Ionization State and Low Oxygen Abundance in Lya Emitters

We present results from Keck/NIRSPEC and Magellan/MMIRS follow-up spectroscopy of Lya emitters (LAEs) at z=2.2 identified in our Subaru narrowband survey. We successfully detect Ha emission from seven LAEs, and perform a detailed analysis of six LAEs free from AGN activity, two out of which, CDFS-3865 and COSMOS-30679, have [OII] and [OIII] line detections. They are the first [OII]-detected LAEs at high-z, and their [OIII]/[OII] ratios and R23-indices provide the first simultaneous determinations of ionization parameter and oxygen abundance for LAEs. CDFS-3865 has a very high ionization parameter (q_{ion}=2.5^{+1.7}_{-0.8}x10^8 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=7.84^{+0.24}_{-0.25}) in contrast with moderate values of other high-z galaxies such as Lyman-break galaxies (LBGs). COSMOS-30679 also possesses a relatively high ionization parameter (q_{ion}=8^{+10}_{-4}x10^7 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=8.18^{+0.28}_{-0.28}). Both LAEs appear to fall below the mass-metallicity relation of z~2 LBGs. Similarly, a low metallicity of 12+log(O/H)<8.4 is independently indicated for typical LAEs from a composite spectrum and the [NII]/Ha index. Such high ionization parameters and low oxygen abundances can be found in local star-forming galaxies, but this extreme local population occupies only ~0.06% of the SDSS spectroscopic galaxy sample with a number density ~100 times smaller than that of LAEs. With their high ionization parameters and low oxygen abundances, LAEs would represent an early stage of galaxy formation dominated by massive stars in compact star-forming regions. High-q_{ion} galaxies like LAEs would produce ionizing photons efficiently with a high escape fraction achieved by density-bounded HII regions, which would significantly contribute to cosmic reionization at z>6.

preprint2013arXiv

First Systematic Search for Oxygen-Line Blobs at High Redshift: Uncovering AGN Feedback and Star-Formation Quenching

We present the first systematic search for extended metal-line [OII]λλ3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]λ3426 line at the 6σ level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeIIλ2587 absorption and FeII*λ2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

preprint2013arXiv

Gas Motion Study of Lya Emitters at z~2 Using UV and Optical Spectral Lines

We present the results of Magellan/MMIRS and Keck/NIRSPEC spectroscopy for five Lya emitters (LAEs) at z=2.2 for which high-resolution FUV spectra from Magellan/MagE are available. We detect nebular emission lines including Ha on the individual basis and low-ionization interstellar (LIS) absorption lines in a stacked FUV spectrum, and measure average offset velocities of the Lya line, Delta_v_Lya, and LIS absorption lines, Delta_v_abs, with respect to the systemic velocity defined by the nebular lines. For a sample of eight z~2-3 LAEs without AGN from our study and the literature, we obtain Delta_v_Lya = 175+/-35 km s^{-1}, which is significantly smaller than that of Lyman-break Galaxies (LBGs), Delta_v_Lya=400 km s^{-1}. The stacked FUV spectrum gives Delta_v_abs = -179 +/- 73 km s^{-1}, comparable to that of LBGs. These positive Delta_v_Lya and negative Delta_v_abs suggest that LAEs also have outflows. In contrast to LBGs, however, the LAEs' Delta_v_Lya is as small as |Delta_v_abs|, suggesting low neutral hydrogen column densities. Such a low column density with a small number of resonant scattering may cause the observed strong Lya emission of LAEs. We find an anti-correlation between Lya equivalent width (EW) and Delta_v_Lya in a compilation of LAE and LBG samples. Although its physical origin is not clear, this anti-correlation result appears to challenge the hypothesis that a strong outflow, by means of a reduced number of resonant scattering, produces a large EW. If LAEs at z>6 have similarly small Delta_v_Lya values, constraints on the reionization history derived from the Lya transmissivity may need to be revised.

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.

preprint2013arXiv

Probing of the Interactions Between the Hot Plasmas and Galaxies in Clusters from z=0.1 to 0.9

Based on optical and X-ray data for a sample of 34 relaxed rich clusters of galaxies with redshifts of 0.1-0.9, we studied relative spatial distributions of the two major baryon contents, the cluster galaxies and the hot plasmas. Using multi-band photometric data taken with the UH88 telescope, we determined the integrated (two dimensional) radial light profiles of member galaxies in each cluster using two independent approaches, i.e., the background subtraction and the color-magnitude filtering. The ICM mass profile of each cluster in our sample, also integrated in two dimensions, was derived from a spatially-resolved spectral analysis using XMM-Newton and Chandra data. Then, the radially-integrated light profile of each cluster was divided by its ICM mass profile, to obtain a profile of "galaxy light vs. ICM mass ratio". The ratio profiles over the central 0.65 R500 regions were found to steepen from the higher- to lower- redshift subsamples, meaning that the galaxies become more concentrated in the ICM sphere towards lower redshifts. The evolution is also seen in galaxy number vs. ICM mass ratio profiles. A range of systematic uncertainties in the galaxy light measurements, as well as many radius-/redshift- dependent biases to the galaxy vs. ICM profiles have been assessed, but none of them is significant against the observed evolution. Besides, the galaxy light vs. total mass ratio profiles also exhibit gradual concentration towards lower redshift. We interpret in the context that the galaxies, the ICM, and the dark matter components followed a similar spatial distribution in the early phase (z>0.5), while the galaxies have fallen towards the center relative to the others at a later phase.

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

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

Evolution of the Sizes of Galaxies over 7<z<12 Revealed by the 2012 Hubble Ultra Deep Field Campaign

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

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 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.

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

Average Metallicity and Star Formation Rate of Lya Emitters Probed by a Triple Narrow-Band Survey

We present the average metallicity and star-formation rate of Lya emitters (LAEs) measured from our large-area survey with three narrow-band (NB) filters covering the Lya, [OII]3727, and Ha+[NII] lines of LAEs at z=2.2. We select 919 z=2.2 LAEs from Subaru/Suprime-Cam NB data in conjunction with Magellan/IMACS spectroscopy. Of these LAEs, 561 and 105 are observed with KPNO/NEWFIRM near-infrared NB filters whose central wavelengths are matched to redshifted [OII] and Ha nebular lines, respectively. By stacking the near-infrared images of the LAEs, we successfully obtain average nebular-line fluxes of LAEs, the majority of which are too faint to be identified individually by narrow-band imaging or deep spectroscopy. The stacked object has an Ha luminosity of 1.7x10^{42} erg s^{-1} corresponding to a star formation rate (SFR) of 14 M_{sun} yr^{-1}. We place, for the first time, a firm lower limit to the average metallicity of LAEs of Z>~0.09 Z_{sun} (2sigma) based on the [OII]/(Ha+[NII]) index together with photo-ionization models and empirical relations. This lower limit of metallicity rules out the hypothesis that LAEs, so far observed at z~2, are extremely metal poor (Z<2x10^{-2} Z_{sun}) young galaxies at the 4sigma level. This limit is higher than a simple extrapolation of the observed mass-metallicity relation of z~2 UV-selected galaxies toward lower masses (5x10^{8} M_{sun}), but roughly consistent with a recently proposed fundamental mass-metallicity relation when the LAEs' relatively low SFR is taken into account. The Ha and Lya luminosities of our NB-selected LAEs indicate that the escape fraction of Lya photons is ~12-30 %, much higher than the values derived for other galaxy populations at z~2.

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

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

Red Star Forming Galaxies and Their Environment at z=0.4 Revealed by Panoramic H-alpha Imaging

We present a wide-field H-alpha imaging survey of the rich cluster CL0939+4713 (Abell 851) at z=0.41 with Suprime-Cam on the Subaru Telescope, using the narrow-band filter NB921. The survey is sensitive to active galaxies with star formation rates down to ~0.3Msun/yr throughout the 27'x27' field, and we identified 445 H-alpha emitters along the large-scale structures around the cluster. Using this sample, we find that (1) the fraction of H-alpha emitters is a strong function of environment and shows a clear decline toward the cluster central region. We also find that (2) the color of H-alpha emitters is clearly dependent on environment. The majority of the H-alpha emitters have blue colors with B-I<2, but we find H-alpha emitters with red colors as well. Such red emitters are very rare in the cluster center or its immediate surrounding regions, while they are most frequently found in groups located far away from the cluster center. These groups coincide with the environment where a sharp transition in galaxy color distribution is seen. This may suggest that dusty star formation activity tends to be involved in galaxy truncation processes that are effective in groups, and it is probably related to the "pre-processing" that generates present-day cluster S0 galaxies. Finally, we confirm that (3) the mass-normalized integrated star formation rate in clusters (i.e. the total star formation rate within 0.5xR200 from the cluster center divided by the cluster dynamical mass) rapidly increases with look-back time following approximately ~(1+z)^6, and it is also correlated with the cluster mass.

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.

preprint2011arXiv

The Mass-Dependent Clustering History of K-selected Galaxies at z < 4 in the SXDS/UDS Field

We investigate mass-dependent galaxy evolution based on a large sample of (more than 50,000) K-band selected galaxies in a multi-wavelength catalog of the Subaru/XMM-Newton Deep Survey (SXDS) and the UKIRT Infrared Deep Sky Survey (UKIDSS)/Ultra Deep Survey (UDS). We employ the optical to near-infrared photometry to determine photometric redshifts of these galaxies. Then, we estimate the stellar mass of our sample galaxies using a standard fitting procedure. From the sample galaxies, we obtain the stellar mass function of galaxies and the cosmic stellar mass density up to z<4. Our results are consistent with previous studies and we find a considerable number of low-mass galaxies (M<10^{10.5}) at the redshift range 3<z<4. The stellar-mass dependent correlation functions of our sample galaxies show clear evolution and they connect to that in the local universe consistently. Also, the massive galaxies show strong clustering throughout our studied redshift range. The correlation length of massive galaxies rapidly decreases from z=4 to 2. We also find some high mass density regions of massive galaxies at 1.4<z<2.5 in our sample, which may be candidate progenitors of the present-day clusters of galaxies. At this redshift range, massive star-forming galaxies are the dominant population making up the structures and the passively evolving galaxies show stronger clustering and they may have formed earlier than those star-forming galaxies.

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.

preprint2010arXiv

Stellar Populations of Lyman Alpha Emitters at z~6-7: Constraints on the Escape Fraction of Ionizing Photons from Galaxy Building Blocks

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=5.7 and 6.6 in a 0.65 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field, using deep images taken with Subaru/Suprime-Cam, UKIRT/WFCAM, and Spitzer/IRAC. We produce stacked multiband images at each redshift from 165 (z=5.7) and 91 (z=6.6) IRAC-undetected objects, to derive typical spectral energy distributions (SEDs) of z~6-7 LAEs for the first time. The stacked LAEs have as blue UV continua as the HST/WFC3 z-dropout galaxies of similar Muv, with a spectral slope beta ~ -3, but at the same time they have red UV-to-optical colors with detection in the 3.6um band. Using SED fitting we find that the stacked LAEs have low stellar masses of ~(3-10)*10^7 Msun, very young ages of ~1-3 Myr, negligible dust extinction, and strong nebular emission from the ionized interstellar medium, although the z=6.6 object is fitted similarly well with high-mass models without nebular emission; inclusion of nebular emission reproduces the red UV-to-optical color while keeping the UV color sufficiently blue. We infer that typical LAEs at z~6-7 are building blocks of galaxies seen at lower redshifts. We find a tentative decrease in the Lyman alpha escape fraction from z=5.7 to 6.6, which may imply an increase in the intergalactic medium neutral fraction. From the minimum contribution of nebular emission required to fit the observed SEDs, we place an upper limit on the escape fraction of ionizing photons to be f_esc^ion~0.6 at z=5.7 and ~0.9 at z=6.6. We also compare the stellar populations of our LAEs with that of stacked HST/WFC3 z-dropout galaxies.

preprint2009arXiv

A Search for Molecular Gas toward a BzK-selected Star-forming Galaxy at z = 2.044

We present a search for CO(3-2) emission in SDF-26821, a BzK-selected star-forming galaxy (sBzK) at z = 2.044, using the 45-m telescope of the Nobeyama Radio Observatory and the Nobeyama Millimeter Array. We do not detect significant emission and derive 2 σlimits: the CO luminosity of L'CO < 3.1 x 10^10 K km s^{-1} pc^{-2}, the ratio of far-infrared luminosity to CO luminosity of L_FIR/L'CO > 57 Lsun (K km s^{-1} pc^{-2})^{-1}, and the molecular gas mass of M_H2 < 2.5 x 10^10 Msun, assuming a velocity width of 200 km s^{-1} and a CO-to-H2 conversion factor of alpha_CO=0.8 Msun (K km s^{-1} pc^{-2})^{-1}. The ratio of L_FIR/L'CO, a measure of star formation efficiency (SFE), is comparable to or higher than the two z ~ 1.5 sBzKs detected in CO(2-1) previously, suggesting that sBzKs can have a wide range of SFEs. Comparisons of far-infrared luminosity, gas mass, and stellar mass among the sBzKs suggest that SDF-26821 is at an earlier stage of forming stars with a similar SFE and/or more efficiently forming stars than the two z ~ 1.5 sBzKs. The higher SFEs and specific star formation rates of the sBzKs compared to local spirals are indicative of the difference in star formation modes between these systems, suggesting that sBzKs are not just scaled-up versions of local spirals.

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

High star formation activity in the central region of a distant cluster at z=1.46

We present an unbiased deep [OII] emission survey of a cluster XMMXCS J2215.9-1738 at z=1.46, the most distant cluster to date with a detection of extended X-ray emission. With wide-field optical and near-infrared cameras (Suprime-Cam and MOIRCS, respectively) on Subaru telescope, we performed deep imaging with a narrow-band filter NB912 (lambda_c=9139A, Delta_lambda=134A) as well as broad-band filters (B, z', J and Ks). From the photometric catalogues, we have identified 44 [OII] emitters in the cluster central region of 6'x6' down to a dust-free star formation rate of 2.6 Msun/yr (3 sigma). Interestingly, it is found that there are many [OII] emitters even in the central high density region. In fact, the fraction of [OII] emitters to the cluster members as well as their star formation rates and equivalent widths stay almost constant with decreasing cluster-centric distance up to the cluster core. Unlike clusters at lower redshifts (z<1) where star formation activity is mostly quenched in their central regions, this higher redshift 2215 cluster shows its high star formation activity even at its centre, suggesting that we are beginning to enter the formation epoch of some galaxies in the cluster core eventually. Moreover, we find a deficit of galaxies on the red sequence at magnitudes fainter than ~M*+0.5 on the colour-magnitude diagram. This break magnitude is brighter than that of lower redshift clusters, and it is likely that we are seeing the formation phase of more massive red galaxies in the cluster core at z~1. These results may indicate inside-out and down-sizing propagation of star formation activity in the course of cluster evolution.

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

Panoramic H-alpha and mid-infrared mapping of star formation in a z=0.8 cluster

We present the first wide-field H-alpha imaging survey around the distant cluster RXJ1716.4+6708 at z=0.81 with a narrow-band filter on MOIRCS/Subaru, down to SFR(Ha)~1Msun/yr. Combining with a wide-field MIR imaging survey with AKARI, we compare the unobscured and obscured star formation activities in the cluster. We find that both H-alpha emitters and MIR galaxies avoid the cluster centre and their spatial distribution is quite similar. Most of the H-alpha emitters show blue colours, but we find some H-alpha emitters on the red sequence. The MIR galaxies tend to be systematically redder than the H-alpha emitters probably due to heavy dust extinction. Interestingly, the red H-alpha emitters and the red MIR galaxies (i.e. dusty red galaxies) are most commonly seen in the medium-density environment such as cluster outskirts, groups and filaments, where optical colours of galaxies change. We also find that A(Ha) exceeds ~3 in extreme cases and that such very dusty galaxies are also located in the medium-density environment. These findings suggest that dusty star formation is triggered in the in-fall region of the cluster, implying a probable link between galaxy transition and dusty star formation. We finally calculate the cluster total SFR and find that the cluster total SFR based on H-alpha alone can be underestimated more than factor ~2 even after 1 mag extinction correction. We suggest that the mass-normalized cluster SFR rapidly declines since z~1 following ~(1+z)^6, although the uncertainty is still large.

preprint2009arXiv

Stellar Populations of Lyman Alpha Emitters at z = 3 - 4 Based on Deep Large Area Surveys in the Subaru-SXDS/UKIDSS-UDS Field

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=3.1 and 3.7 in 0.65 deg^2 of the Subaru/XMM-Newton Deep Field, based on rest-frame UV-to-optical photometry obtained from the Subaru/XMM-Newton Deep Survey, the UKIDSS/Ultra Deep Survey, and the Spitzer legacy survey of the UKIDSS/UDS. Among a total of 302 LAEs (224 for z=3.1 and 78 for z=3.7), only 11 are detected in the K band, i.e., brighter than K(3sigma)=24.1 mag. Eight of the 11 K-detected LAEs are spectroscopically confirmed. We find that the K-undetected objects, which should closely represent the LAE population as a whole, have low stellar masses of ~ 10^8 - 10^8.5 Msun, modest SFRs of 1 - 100 Msun yr^-1, and modest dust extinction of E(B-V) < 0.2. The K-detected objects are massive, Mstar ~ 10^9 - 10^10.5 Msun, and have significant dust extinction with a median of E(B-V) ~= 0.3. Four K-detected objects with the reddest spectral energy distributions, two of which are spectroscopically confirmed, are heavily obscured with E(B-V) ~ 0.65, and their continua resemble those of some local ULIRGs. Interestingly, they have large Lyman alpha equivalent widths ~= 70 - 250 A. If these four are excluded, our sample has a weak anti-correlation between EW(Lya) and Mstar. We compare the stellar masses and the specific star formation rates (sSFR) of LAEs with those of Lyman-break galaxies, distant red galaxies, submillimetre galaxies, and I- or K-selected galaxies with z_phot ~ 3. We find that the LAE population is the least massive among all the galaxy populations in question, but with relatively high sSFRs, while NIR-detected LAEs have Mstar and sSFR similar to LBGs. Our reddest four LAEs have very high sSFRs in spite of large Mstar, thus occupying a unique region in the Mstar versus sSFR space. (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

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.

preprint2006arXiv

Systematic Survey of Extended Lyman-alpha Sources over z~3-5

Spatially extended Ly-alpha sources which are faint and compact in coninuum are candidates for extremely young galaxies (age of ~<10e+7 yrs) at high redshifts. We carried out a systematic survey for extended Ly-alpha sources, using deep intermediate-band imaging data taken with the Suprime-Cam on the Subaru Telescope. Our survey covers a field of view of 33' * 25' and a redshift range of 3.24<z<4.95, down to a Ly-alpha flux of ~1e-17 ergs/cm^2/s. We identified 41 extended Ly-alpha sources. The redshift distribution of these sources shows that this kind of objects are common in the early universe over the surveyed redshift range. The identified objects have typical sizes of \~10-15 kpc and luminosities of ~10^42 ergs/s. Follow-up spectroscopy made for 7 of the 41 objects showed that our sample suffers from little contamination. All 7 objects have large equivalent widths of Ly-alpha emission line, all but one exceeding 240A in the rest frame. The large equivalent widths suggest that their extended Ly-alpha emissions are unlikely to be due to normal starbursts, but possibly originated from very young galaxies. All 41 objects in our sample have much smaller Ly-alpha luminosities than the two Ly-alpha Blobs (LABs) found by Steidel et al. (2000) in spite of our much larger survey volume. This suggests that large and luminous extended Ly-alpha objects like the two LABs are very rare and are clustered in overdense regions.

preprint2002arXiv

Subaru Prime Focus Camera -- Suprime-Cam --

We have built an 80-mega pixels (10240 X 8192) mosaic CCD camera, called Suprime-Cam, for the wide-field prime focus of the 8.2 m Subaru telescope. Suprime-Cam covers a field of view 34' X 27', a unique facility among the the 8-10 m class telescopes, with a resolution of 0.''202 per pixel. The focal plane consists of ten high-resistivity 2k X 4k CCDs developed by MIT Lincoln Laboratory, which are cooled by a large stirling-cycle cooler. The CCD readout electronics was designed to be scalable, which allows the multiple read-out of tens of CCDs. It takes 50 seconds to readout entire arrays. We designed a filter-exchange mechanism of the jukebox type that can hold up to ten large filters (205 X 170 X 15 mm^3). The wide-field corrector is basically a three-lens Wynne-type, but has a new type of atmospheric dispersion corrector. The corrector provides a flat focal plane and an un-vignetted field of view of 30' in diameter. The achieved co-planarity of the focal array mosaic is smaller than 30 micron peak-to-peak, which realizes mostly the seeing limited image over the entire field. The median seeing in the Ic-band, measured over one year and a half, is 0.''61. The PSF anisotropy in Suprime-Cam images, estimated by stellar ellipticities, is about 2% under this median seeing condition. At the time of commissioning, Suprime-Cam had the largest survey speed, which is defined as the field of view multiplied by the primary mirror area of the telescope, among those cameras built for sub-arcsecond imaging.

preprint2001arXiv

Galaxy Population in a Cluster of Galaxies around the Radio Galaxy 3C 324 at z=1.2

We discuss the properties of galaxies around the radio galaxy 3C 324 at z=1.2 based on BVRIK' multi-band imaging data. We have applied a photometric-redshift technique to objects in the 3C 324 field, and identified 35 objects as plausible cluster members. We have found that red and luminous members are concentrated in a small region enclosed by a circle of 40" radius (0.33 Mpc at z=1.2 for Ω_0=0.3, λ_0=0.7, H_0=70km/s/Mpc cosmology) from the 3C 324 galaxy. The 3C 324 cluster is probably much more compact in size compared with the local clusters. We constructed a K'-band luminosity function of the cluster members and fit a Schechter function, and found the characteristic magnitude to be K*_{AB}=20.2+-0.6. This value is consistent with the extrapolation of the pure passive evolution seen for z<1 clusters. We have identified eight bright galaxies which form a red color-magnitude sequence. The slope of the sequence is consistent with the passive evolution model down to K'_{AB}<22; we also found that there is no clear age variation in these bright red galaxies. However, seven out of these eight galaxies exhibit a significant excess in the rest UV light with respect to the passive evolution model. This may suggest that the massive early-type galaxies in this high-redshift cluster are still forming stars to some extent. We have confirmed a truncation of the color-magnitude sequence at K'_{AB}\sim22; faint passively-evolving galaxies may not yet be present in this cluster at z\sim1.2. The overall color distribution of the cluster members, selected by the photometric redshift technique, is found to be very broad. We derived the fraction of blue galaxies in this cluster following a

preprint2000arXiv

Evolution of the Blue Luminosity-to-Baryon Mass Ratio of Clusters of Galaxies

We derive the ratio of total blue luminosity to total baryon mass, LB/Mb, for massive (Mgas at the Abell radius is \ge 1 \times 10^{13} h^{-2.5} \Msolar) clusters of galaxies up to z \simeq 1 from the literature. Twenty-two clusters in our sample are at z > 0.1. Assuming that the relative mix of hot gas and galaxies in clusters does not change during cluster evolution, we use LB/Mb to probe the star formation history of the galaxy population as a whole in clusters. We find that LB/Mb of clusters increases with redshift from LB/Mb=0.024 (solar units) at z = 0 to \simeq 0.06 at z=1, indicating a factor of 2-3 brightening (we assume H0=70 km/s/Mpc). This amount of brightening is almost identical to the brightening of the M/LB ratio of early-type galaxies in clusters at 0.02 \le z \le 0.83 reported by van Dokkum et al. (1998). We compare the observed brightening of LB/Mb with luminosity evolution models for the galaxy population as a whole, changing the e-folding time of star formation τby 0.1 \le τ\le 5 Gyr and the formation redshift \zF by 2 \le \zF < \infty. We find that τ=0.1 Gyr 'single burst' models with \zF \ge 3 and τ=5 Gyr 'disk' models with arbitrary \zF are consistent with the observed brightening, while models with τ=1-2 Gyr tend to predict too steep brightening. We also derive the ratio of blue luminosity density to baryon density for field galaxies, adopting Ω_b h^2 = 0.02, and find that blue luminosity per unit baryon is similar in clusters and in fields up to z \simeq 1 within the observational uncertainties.