Search

Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$\nu^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$\nu^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST)., Comment: 18 pages, 16 figures, accepted for publication in MNRAS, Figs. 3, 4, 11, 13, and 14 revised, Figs. 14 and 16 in the original manuscript moved to the Appendix, All Uchuu-$\nu^2$GC data used in this paper are publicly available on our Skies & Universes website https://www.skiesanduniverses.org/Simulations/Uchuu/GalaxyCatalogues/

We present a new catalog of $9318$ Ly$\alpha$ emitter (LAE) candidates at $z = 2.2$, $3.3$, $4.9$, $5.7$, $6.6$, and $7.0$ that are photometrically selected by the SILVERRUSH program with a machine learning technique from large area (up to $25.0$ deg$^2$) imaging data with six narrowband filters taken by the Subaru Strategic Program with Hyper Suprime-Cam (HSC SSP) and a Subaru intensive program, Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). We construct a convolutional neural network that distinguishes between real LAEs and contaminants with a completeness of $94$% and a contamination rate of $1$%, enabling us to efficiently remove contaminants from the photometrically selected LAE candidates. We confirm that our LAE catalogs include $177$ LAEs that have been spectroscopically identified in our SILVERRUSH programs and previous studies, ensuring the validity of our machine learning selection. In addition, we find that the object-matching rates between our LAE catalogs and our previous results are $\simeq 80$-$100$% at bright NB magnitudes of $\lesssim 24$ mag. We also confirm that the surface number densities of our LAE candidates are consistent with previous results. Our LAE catalogs will be made public on our project webpage., Comment: 19 pages, 10 figures, accepted for publication in ApJ. Our LAE catalogs will become available at http://cos.icrr.u-tokyo.ac.jp/rush.html

We present the quasar luminosity function at $z \sim 5$ derived from the optical wide-field survey data obtained as a part of the Subaru strategic program (SSP) with Hyper Suprime-Cam (HSC). From $\sim$81.8 deg$^2$ area in the Wide layer of the HSC-SSP survey, we selected 224 candidates of low-luminosity quasars at $z \sim 5$ by adopting the Lyman-break method down to $i = 24.1$ mag. Based on our candidates and spectroscopically-confirmed quasars from the Sloan Digital Sky Survey (SDSS), we derived the quasar luminosity function at $z \sim 5$ covering a wide luminosity range of $-28.76 < M_{\rm 1450} < -22.32$ mag. We found that the quasar luminosity function is fitted by a double power-law model with a break magnitude of $M^{*}_{1450} = -25.05^{+0.10}_{-0.24}$ mag. The inferred number density of low-luminosity quasars is lower, and the derived faint-end slope, $-1.22^{+0.03}_{-0.10}$, is flatter than those of previous studies at $z \sim 5$. A compilation of the quasar luminosity function at $4 \leq z \leq 6$ from the HSC-SSP suggests that there is little redshift evolution in the break magnitude and in the faint-end slope within this redshift range, although previous studies suggest that the faint-end slope becomes steeper at higher redshifts. The number density of low-luminosity quasars decreases more rapidly from $z \sim 5$ to $z \sim 6$ than from $z \sim 4$ to $z \sim 5$., Comment: 18 pages, 11 figures, accepted for publication in ApJ

We construct a model of H$\alpha$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($\nu^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$\alpha$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$\alpha$ luminosity function (LF) at $z=0.4$. The large volume of the $\nu^2$GC makes it possible to examine the spatial distribution of HAEs over a region of (411.8 Mpc)$^3$ in the comoving scale. The surface number density of $z=0.4$ HAEs with $L_{\rm H\alpha} \geq 10^{40}$ erg s$^{-1}$ is 308.9 deg$^{-2}$. We have confirmed that the HAE is a useful tracer for the large-scale structure of the Universe because of their significant overdensity ($>$ 5$\sigma$) at clusters and the filamentary structures. The H$\alpha$ LFs within a survey area of $\sim$2 deg$^2$ (typical for previous observational studies) show a significant field variance up to $\sim$1 dex. Based on our model, one can estimate the variance on the H$\alpha$ LFs within given survey areas., Comment: 13 pages, 11 figures, accepted for publication in the ApJ

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

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ 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$\alpha$ 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}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$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$\alpha$ emission of the LABs. The Ly$\alpha$ 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$\alpha$ luminosity relation of the Ly$\alpha$ 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., Comment: Submitted to ApJ

The observed value $\Lambda_{\rm obs}$ of the cosmological constant $\Lambda$ is extremely smaller than theoretical expectations, and the anthropic argument has been proposed as a solution to this problem because galaxies do not form when $\Lambda \gg \Lambda_{\rm obs}$. However, the contemporary galaxy formation theory predicts that stars form even with a high value of $\Lambda / \Lambda_{\rm obs} \sim$ 50, which makes the anthropic argument less persuasive. Here we calculate the probability distribution of $\Lambda$ using a model of cosmological galaxy formation, considering extinction of observers caused by radiation from nearby supernovae. The life survival probability decreases in a large $\Lambda$ universe because of higher stellar density. Using a reasonable rate of lethal supernovae, we find that the mean expectation value of $\Lambda$ can be close to $\Lambda_{\rm obs}$, and hence this effect may be essential to understand the small but nonzero value of $\Lambda$. It is predicted that we are located on the edge of habitable regions about stellar density in the Galaxy, which may be tested by future exoplanet studies., Comment: 13 pages, 2 figures. the version accepted to Astrobiology

We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray luminosity 10^{ 44 }< erg/s is underestimated compared with recent observational estimates when we assume the exponentially decreasing accretion rate and the accretion timescale which is proportional to the dynamical time of the host halo or the bulge, as is often assumed in semi-analytic models. We show that to solve this discrepancy, the accretion timescale of such less luminous AGNs instead should be a function of the black hole mass and the accreted gas mass. This timescale can be obtained from a phenomenological modelling of the gas angular momentum loss in the circumnuclear torus and/or the accretion disc. Such models predict a longer accretion timescale for less luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be 10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can explain the observed luminosity functions of AGNs at z < 6.0., Comment: 31 pages and 29 figures. Accepted for publication in MNRAS

Super-Eddington mass accretion has been suggested as an efficient mechanism to grow supermassive black holes (SMBHs). We investigate the imprint left by the radiative efficiency of the super-Eddington accretion process on the clustering of quasars using a new semi-analytic model of galaxy and quasar formation based on large-volume cosmological $N$-body simulations. Our model includes a simple model for the radiative efficiency of a quasar, which imitates the effect of photon trapping for a high mass accretion rate. We find that the model of radiative efficiency affects the relation between the quasar luminosity and the quasar host halo mass. The quasar host halo mass has only weak dependence on quasar luminosity when there is no upper limit for quasar luminosity. On the other hand, it has significant dependence on quasar luminosity when the quasar luminosity is limited by its Eddington luminosity. In the latter case, the quasar bias also depends on the quasar luminosity, and the quasar bias of bright quasars is in agreement with observations. Our results suggest that the quasar clustering studies can provide a constraint on the accretion disc model., Comment: 6 pages, 4 figures, to appear in MNRAS Letters

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016., Comment: 21 pages, 14 figures, accepted for publication in PASJ special issue

We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3., Comment: 26 pages, 12 figures. Accepted for publication in PASJ special issue

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs., Comment: 19 pages, 11 figures. Accepted for publication in PASJ special issue

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe., Comment: 16 pages, 10 figures. Accepted for publication in PASJ special issue

In order to investigate the physical relationship between strong Ly-alpha absorbers (logN_HI>20.0 cm^-2) such as damped Ly-alpha absorption systems (DLAs) and young star-forming galaxies at high redshift, we have conducted narrow-band observations of Ly-alpha emitters (LAEs) in a concentrated region of strong Ly-alpha absorbers at z= 2.3, the J1230+34 field. Using a catalog of Ly-alpha absorbers with logN_HI>20.0 cm^-2 based on the baryon oscillation spectroscopic survey (BOSS), we found 6 fields where 3 or more absorbers are concentrated within a (50 Mpc)^3 cubic box in the comoving scale. Among them, we focus on the J1230+34 field, where 2 DLAs and 2 sub-DLAs present. Our narrow-band imaging observations with Subaru/Suprime-Cam using a custom-made filter, NB400 (lambda_c=4003 A and FWHM=92 A) yield a sample of 149 LAEs in this field. In the large scale (~50 Mpc), we have found no differences between the obtained Ly-alpha luminosity function and those in the blank fields at similar redshifts. We also compare the frequency distribution of the Ly-alpha rest-frame equivalent width (EW_0) in the target field and other fields including both overdensity region and blank field, but find no differences. On the other hand, in the small scale (~10 Mpc), we have found a possible overdensity of LAEs around a DLA with the highest HI column density (N_HI=21.08 cm^-2) in the target field while there are no density excess around the other absorbers with a lower N_HI., Comment: 15 pages, 15 figures, Accepted for publication to PASJ

We detect 20 $z=7.0$ Ly$\alpha$ emitter (LAE) candidates to $L({\rm Ly}\alpha) \geq 2 \times 10^{42}$ erg s$^{-1}$ or $0.3$ $L^*_{z=7}$ and in $6.1\times 10^5$ Mpc$^3$ volume in the Subaru Deep Field and the Subaru/XMM-Newton Deep Survey field by 82 and 37 hours of Subaru Suprime-Cam narrowband NB973 and reddest optical $y$-band imaging. We compare their Ly$\alpha$ and UV luminosity functions (LFs) and densities and Ly$\alpha$ equivalent widths (EWs) to those of $z=5.7$, 6.6 and 7.3 LAEs from previous Suprime-Cam surveys. The Ly$\alpha$ LF (density) rapidly declines by a factor of $\times$1.5 (1.9) in $L({\rm Ly}\alpha)$ at $z=5.7-6.6$ (160 Myr), $\times$1.5 (1.6) at $z=6.6-7.0$ (60 Myr) at the faint end and $\times$2.0 (3.8) at $z=7.0-7.3$ (40 Myr). Also, in addition to the systematic decrease in EW at $z=5.7-6.6$ previously found, 2/3 of the $z=7.0$ LAEs detected in the UV continuum exhibit lower EWs than the $z=6.6$ ones. Moreover, while the UV LF and density do not evolve at $z=5.7-6.6$, they modestly decline at $z=6.6-7.0$, implying galaxy evolution contributing to the decline of the Ly$\alpha$ LF. Comparison of the $z=7.0$ Ly$\alpha$ LF to the one predicted by an LAE evolution model further reveals that galaxy evolution alone cannot explain all the decline of Ly$\alpha$ LF. If we attribute the discrepancy to Ly$\alpha$ attenuation by neutral hydrogen, the intergalactic medium transmission of Ly$\alpha$ photons at $z=7.0$ would be $T_{{\rm Ly}\alpha}^{\rm IGM} \leq 0.6-0.7$. It is lower (higher) than $T_{{\rm Ly}\alpha}^{\rm IGM}$ at $z=6.6$ (7.3) derived by previous studies, suggesting rapid increase in neutral fraction at $z > 6$., Comment: 37 Pages, 22 Figures, 10 Tables, Accepted for publication in ApJ

Recent studies have derived quasar luminosity functions (QLFs) at various redshifts. However, the faint side of the QLF at high redshifts is still too uncertain. An accurate estimate of the survey completeness is essential to derive an accurate QLF for use in studying the luminosity-dependent density evolution of the quasar population. Here we investigate how the luminosity dependence of quasar spectra (the Baldwin effect) and the attenuation model for the inter-galactic medium (IGM) affect the completeness estimates. For this purpose, we revisit the completeness of quasar surveys specifically at $z\sim4-5$, using the COSMOS images observed with Subaru/Suprime-Cam. As the result, we find that the completeness estimates are sensitive to the luminosity dependence of the quasar spectrum and difference in the IGM attenuation models. At $z\sim4$, the number density of quasars when we adopt the latest IGM model and take the luminosity dependence of spectra into account are $(3.49\pm1.62)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-24.09

We explore the effect of varying the mass of the seed black hole on the resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model by requiring the observed properties of galaxies at z ~ 0 are reproduced. In keeping with previous semi-analytic models, we place a seed black hole immediately after a galaxy forms. When the mass of the seed is set at 10^5 M_sun, we find that the model results become inconsistent with recent observational results of the black hole mass - bulge mass relation for dwarf galaxies. In particular, the model predicts that bulges with ~ 10^9 M_sun harbour larger black holes than observed. On the other hand, when we employ seed black holes with 10^3 M_sun, or randomly select their mass within a 10^(3-5) M_sun range, the resulting relation is consistent with observation estimates, including the observed dispersion. We find that to obtain stronger constraints on the mass of seed black holes, observations of less massive bulges at z ~ 0 are a more powerful comparison than the relations at higher redshifts., Comment: 6 pages, 2 tables, 4 figures. Accepted for publication by MNRAS Letter

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies., Comment: 20 pages, 7 figures, accepted for publication in ApJ

We investigate morphological properties of 61 Lyman-alpha emitters (LAEs) at z = 4.86 identified in the COSMOS field, based on Hubble Space Telescope Advanced Camera for Surveys (ACS) imaging data in the F814W-band. Out of the 61 LAEs, we find the ACS counterparts for the 54 LAEs. Eight LAEs show double-component structures with a mean projected separation of 0."63 (~ 4.0 kpc at z = 4.86). Considering the faintness of these ACS sources, we carefully evaluate their morphological properties, that is, size and ellipticity. While some of them are compact and indistinguishable from the PSF half-light radius of 0."07 (~ 0.45 kpc), the others are clearly larger than the PSF size and spatially extended up to 0."3 (~ 1.9 kpc). We find that the ACS sources show a positive correlation between ellipticity and size and that the ACS sources with large size and round shape are absent. Our Monte Carlo simulation suggests that the correlation can be explained by (1) the deformation effects via PSF broadening and shot noise or (2) the source blending in which two or more sources with small separation are blended in our ACS image and detected as a single elongated source. Therefore, the 46 single-component LAEs could contain the sources which consist of double (or multiple) components with small spatial separation (i.e., < 0."3 or 1.9 kpc). Further observation with high angular resolution at longer wavelengths (e.g., rest-frame wavelengths of > 4000 A) is inevitable to decipher which interpretation is adequate for our LAE sample., Comment: 20 pages, 18 figures, 4 tables, accepted for publication in ApJ on 19 January 2016

We investigate clustering properties of quasars using a new version of our semi-analytic model of galaxy and quasar formation with state-of-the-art cosmological N-body simulations. In this study, we assume that a major merger of galaxies triggers cold gas accretion on to a supermassive black hole and quasar activity. Our model can reproduce the downsizing trend of the evolution of quasars. We find that the median mass of quasar host dark matter haloes increases with cosmic time by an order of magnitude from z=4 (a few 1e+11 Msun) to z=1 (a few 1e+12 Msun), and depends only weakly on the quasar luminosity. Deriving the quasar bias through the quasar--galaxy cross-correlation function in the model, we find that the quasar bias does not depend on the quasar luminosity, similar to observed trends. This result reflects the fact that quasars with a fixed luminosity have various Eddington ratios and thus have various host halo masses that primarily determine the quasar bias. We also show that the quasar bias increases with redshift, which is in qualitative agreement with observations. Our bias value is lower than the observed values at high redshifts, implying that we need some mechanisms that make quasars inactive in low-mass haloes and/or that make them more active in high-mass haloes., Comment: 6 pages, 3 figures, to appear in MNRAS Letters

We present a new cosmological galaxy formation model, $\nu^2$GC, as an updated version of our previous model $\nu$GC. We adopt the so-called "semi-analytic" approach, in which the formation history of dark matter halos is computed by ${\it N}$-body simulations, while the baryon physics such as gas cooling, star formation and supernova feedback are simply modeled by phenomenological equations. Major updates of the model are as follows: (1) the merger trees of dark matter halos are constructed in state-of-the-art ${\it N}$-body simulations, (2) we introduce the formation and evolution process of supermassive black holes and the suppression of gas cooling due to active galactic nucleus (AGN) activity, (3) we include heating of the intergalactic gas by the cosmic UV background, and (4) we tune some free parameters related to the astrophysical processes using a Markov chain Monte Carlo method. Our ${\it N}$-body simulations of dark matter halos have unprecedented box size and mass resolution (the largest simulation contains 550 billion particles in a 1.12 Gpc/h box), enabling the study of much smaller and rarer objects. The model was tuned to fit the luminosity functions of local galaxies and mass function of neutral hydrogen. Local observations, such as the Tully-Fisher relation, size-magnitude relation of spiral galaxies and scaling relation between the bulge mass and black hole mass were well reproduced by the model. Moreover, the model also well reproduced the cosmic star formation history and the redshift evolution of rest-frame ${\it K}$-band luminosity functions. The numerical catalog of the simulated galaxies and AGNs is publicly available on the web., Comment: 25 pages, 21 figures, 2 tables, accepted to PASJ. The numerical catalog of galaxies and AGNs are available at http://cdsarc.u-strasbg.fr/cgi-bin/VizieR?-source=J/PASJ/68/25

We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$\alpha$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history., Comment: Accepted for publication in ApJ Letters on 15th July, 2015. 6 pages including 3 figures and 2 tables

We present the evolution of dark matter halos in six large cosmological N-body simulations, called the $\nu^2$GC (New Numerical Galaxy Catalog) simulations on the basis of the LCDM cosmology consistent with observational results obtained by the Planck satellite. The largest simulation consists of $8192^3$ (550 billion) dark matter particles in a box of $1.12 \, h^{-1} \rm Gpc$ (a mass resolution of $2.20 \times 10^{8} \, h^{-1} M_{\odot}$). Among simulations utilizing boxes larger than $1 \, h^{-1} \rm Gpc$, our simulation yields the highest resolution simulation that has ever been achieved. A $\nu^2$GC simulation with the smallest box consists of eight billions particles in a box of $70 \, h^{-1} \rm Mpc$ (a mass resolution of $3.44 \times 10^{6} \, h^{-1} M_{\odot}$). These simulations can follow the evolution of halos over masses of eight orders of magnitude, from small dwarf galaxies to massive clusters. Using the unprecedentedly high resolution and powerful statistics of the $\nu^2$GC simulations, we provide statistical results of the halo mass function, mass accretion rate, formation redshift, and merger statistics, and present accurate fitting functions for the Planck cosmology. By combining the $\nu^2$GC simulations with our new semi-analytic galaxy formation model, we are able to prepare mock catalogs of galaxies and active galactic nuclei, which will be made publicly available in the near future., Comment: 18 pages, 16 figures, accepted, simulation data can be downloaded on http://hpc.imit.chiba-u.jp/~nngc/

We present a new chemical evolution model for dwarf spheroidal galaxies (dSphs) in the Local Universe. Our main aim is to explain both their observed star formation histories and metallicity distribution functions simultaneously. Applying our new model for the four local dSphs, that is, Fornax, Sculptor, Leo II, and Sextans, we find that our new model reproduces the observed chemical properties of the dSphs consistently. Our results show that the dSphs have evolved with both a low star formation efficiency and a large gas outflow efficiency compared with the Milky Way, as suggested by previous works. Comparing the observed [alpha/Fe]--[Fe/H] relation of the dSphs with the model predictions, we find that our model favors a longer onset time of Type Ia supernovae (i.e., 0.5 Gyr) than that suggested in previous studies (i.e., 0.1 Gyr). We discuss the origin of this discrepancy in detail., Comment: 18 pages, 13 figures, Accepted for publication in ApJ

We have investigated effects of dust attenuation on quasar luminosity functions using a semi-analytic galaxy formation model combined with a large cosmological N-body simulation. We estimate the dust attenuation of quasars self-consistently with that of galaxies by considering the dust in their host bulges. We find that the luminosity of the bright quasars is strongly dimmed by the dust attenuation, about 2 mag in the B-band. Assuming the empirical bolometric corrections for active galactic nuclei (AGNs) by Marconi et al., we find that this dust attenuation is too strong to explain the B-band and X-ray quasar luminosity functions simultaneously. We consider two possible mechanisms that weaken the dust attenuation. As such a mechanism, we introduce a time delay for AGN activity, that is, gas fueling to a central black hole starts some time after the beginning of the starburst induced by a major merger. The other is the anisotropy in the dust distribution. We find that in order to make the dust attenuation of the quasars negligible, either the gas accretion into the black holes has to be delayed at least three times the dynamical timescale of their host bulges or the dust covering factor is as small as 0.1., Comment: 5 pages, 4 figures, 1 tables. Accepted for publication by MNRAS Letters

Recent observations show that the space density of luminous active galactic nuclei (AGNs) peaks at higher redshifts than that of faint AGNs. This downsizing trend in the AGN evolution seems to be contradictory to the hierarchical structure formation scenario. In this study, we present the AGN space density evolution predicted by a semi-analytic model of galaxy and AGN formation based on the hierarchical structure formation scenario. We demonstrate that our model can reproduce the downsizing trend of the AGN space density evolution. The reason for the downsizing trend in our model is a combination of the cold gas depletion as a consequence of star formation, the gas cooling suppression in massive halos and the AGN lifetime scaling with the dynamical timescale. We assume that a major merger of galaxies causes a starburst, spheroid formation, and cold gas accretion onto a supermassive black hole (SMBH). We also assume that this cold gas accretion triggers AGN activity. Since the cold gas is mainly depleted by star formation and gas cooling is suppressed in massive dark halos, the amount of cold gas accreted onto SMBHs decreases with cosmic time. Moreover, AGN lifetime increases with cosmic time. Thus, at low redshifts, major mergers do not always lead to luminous AGNs. Because the luminosity of AGNs is correlated with the mass of accreted gas onto SMBHs, the space density of luminous AGNs decreases more quickly than that of faint AGNs. We conclude that the anti-hierarchical evolution of the AGN space density is not contradictory to the hierarchical structure formation scenario., Comment: 10 pages, 5 figures. Accepted for publication in ApJ. Minor editorial changes and typos corrected

We present a new version of a semi-analytic model of cosmological galaxy formation, incorporating a star formation law with a feedback depending on the galaxy-scale mean dust opacity and metallicity, motivated by recent observations of star formation in nearby galaxies and theoretical considerations. This new model is used to investigate the effect of such a feedback on shaping the galaxy luminosity function and its evolution. Star formation activity is significantly suppressed in dwarf galaxies by the new feedback effect, and the faint-end slope of local luminosity functions can be reproduced with a reasonable strength of supernova feedback, which is in contrast to the previous models that require a rather extreme strength of supernova feedback. Our model can also reproduce the early appearance of massive galaxies manifested in the bright-end of high redshift K-band luminosity functions. Though some of the previous models also succeeded in reproducing this, they assumed a star formation law depending on the galaxy-scale dynamical time, which is not supported by observations. We argue that the feedback depending on dust opacity (or metal column density) is essential, rather than that simply depending on gas column density, to get these results., Comment: 10 pages, 7 figures, Accepted for publication in MNRAS. v2: minor revision, new references added

We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) is tightly correlated by a power-law to specific SFR (SSFR), and L_TIR is a good SFR indicator only for galaxies with the largest SSFR. We discovered a tight linear correlation for normal galaxies between the radiation field strength of dust heating, estimated by infrared SED fits (U_h), and that of galactic-scale infrared emission (U_TIR ~ L_TIR/R^2), where R is the optical size of a galaxy. The dispersion of U_h along this relation is 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the U_h/U_TIR ratio can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the dust column density, supporting this interpretation. In the large U_TIR limit, the data of circumnuclear starbursts indicate the existence of an upper limit on U_h, corresponding to the maximum SFR per gas mass of ~ 10 Gyr^{-1}. We find that the number of galaxies sharply drops when they become optically thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. Implications are discussed for the M_HI-size relation, the Kennicutt-Schmidt relation, and galaxy formation in the cosmological context., Comment: 29 pages including 28 figures. matches the published version (PASJ 2011 Dec. 25 issue). The E-open option was chosen for this article, i.e., the official version available from PASJ site (http://pasj.asj.or.jp/v63/n6/630613/630613-frame.html) without restriction