Search

We have investigated the redshift evolution of the relationship between supermassive black hole (SMBH) mass and host bulge mass using a semi-analytical galaxy formation model $\nu^2$GC. Our model reproduces the relation in the local universe well. We find that, at high redshift ($z \gtrsim 3$), two sequences appear in the SMBH mass--bulge mass plane. The emergence of these two sequences can be attributed to the primary triggers of the growth of the SMBHs and bulges: galaxy mergers and disc instabilities. The growth of SMBHs and bulges as a result of galaxy mergers is responsible for giving rise to the high-mass sequence, in which SMBHs are more massive for a given host bulge mass than in the low-mas sequence. Conversely, disc instabilities are accountable for the emergence of the low-mass sequence. At lower redshifts, galaxy mergers tend to become increasingly deficient in gas, resulting in a preferential increase of bulge mass without a corresponding growth in SMBH mass. This has the effect of causing galaxies in the upper sequence to shift towards the lower one on the SMBH mass-bulge mass plane. The galaxies that undergo dry mergers serve to bridge the gap between the two sequences, eventually leading to convergence into a single relation known in the local universe. Our results suggest that the observations of the SMBH mass-bulge mass relation in high redshifts can provide insight into their growth mechanisms., Comment: 7 pages, 7 figure, accepted for publication in MNRAS

Herein, we present the statistical properties of active galactic nuclei (AGNs) for approximately 1 million member galaxies of galaxy groups and clusters, with 0.1 $<$ cluster redshift ($z_{\rm cl}$) $<$ 1.4, selected using Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction ($f_{\rm AGN}$), which is defined as the proportion of the contribution of AGNs to the total infrared (IR) luminosity, $L_{\rm IR}$ (AGN)/$L_{\rm IR}$, and examined how $f_{\rm AGN}$ depends on (i) $z_{\rm cl}$ and (ii) the distance from the cluster center. We compiled multiwavelength data using the ultraviolet--mid-IR range. Moreover, we performed spectral energy distribution fits to determine $f_{\rm AGN}$ using the CIGALE code with the SKIRTOR AGN model. We found that (i) the value of $f_{\rm AGN}$ in the CAMIRA clusters is positively correlated with $z_{\rm cl}$, with the correlation slope being steeper than that for field galaxies, and (ii) $f_{\rm AGN}$ exhibits a high value at the cluster outskirts. These results indicate that the emergence of AGN population depends on the redshift and environment and that galaxy groups and clusters at high redshifts are important in AGN evolution. Additionally, we demonstrated that cluster--cluster mergers may enhance AGN activity at the outskirts of particularly massive galaxy clusters. Our findings are consistent with a related study on the CAMIRA clusters that was based on the AGN number fraction., Comment: 25 pages, 24 figures, and 3 tables, accepted for publication in ApJ. A value-added CAMIRA member galaxy catalog and the best-fit SED for each member galaxy will be available as FITS or machine-readable tables

We investigate the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of broad-line AGNs at z=4, based on a sample of 52 quasars with i<23.2 at 3.50 < z < 4.25 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A-Wide2 dataset, and 1,462 quasars with i<20.2 in the same redshift range from the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog. Virial BH masses of quasars are estimated using the width of the CIV 1549{\AA} line and the continuum luminosity at 1350{\AA}. To obtain the intrinsic broad-line AGN BHMF and ERDF, we correct for the incompleteness in the low-mass and/or low-Eddington-ratio ranges caused by the flux-limited selection. The resulting BHMF is constrained down to $\log M_{\rm BH}/M_{\odot}\sim7.5$. In comparison with broad-line AGN BHMFs at z=2 in the literature, we find that the number density of massive SMBHs peaks at higher redshifts, consistent with the "down-sizing" evolutionary scenario. Additionally, the resulting ERDF shows a negative dependence on BH mass, suggesting more massive SMBHs tend to accrete at lower Eddington ratios at z=4. With the derived intrinsic broad-line AGN BHMF, we also evaluate the active fraction of broad-line AGNs among the entire SMBH population at z=4. The resulting active fraction may suggest a positive dependence on BH mass. Finally, we examine the time evolution of broad-line AGN BHMF between z=4 and 6 through solving the continuity equation. The results suggest that the broad-line AGN BHMFs at z=4-6 only show evolution in their normalization, but with no significant changes in their shape., Comment: 38 pages, 26 figures. Submitted to ApJ. Data is available in https://github.com/wanqqq31/BHMF_hsc_z4

One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction ($f_{\rm AGN}$) on cluster redshift ($z_{\rm cl}$) and distance from the cluster center ($R/R_{\rm 200}$). We focus on more than 27,000 galaxy groups and clusters at $0.1 < z_{\rm cl} < 1.4$ with more than 1 million member galaxies selected from the Subaru Hyper Suprime-Cam. By combining various AGN selection methods based on infrared (IR), radio, and X-ray data, we identify 2,688 AGN. We find that (i) $f_{\rm AGN}$ increases with $z_{\rm cl}$ and (ii) $f_{\rm AGN}$ decreases with $R/R_{\rm 200}$. The main contributors to the rapid increase of $f_{\rm AGN}$ towards high-$z$ and cluster center are IR- and radio-selected AGN, respectively. Those results indicate that the emergence of the AGN population depends on the environment and redshift, and galaxy groups and clusters at high-$z$ play an important role in AGN evolution. We also find that cluster-cluster mergers may not drive AGN activity in at least the cluster center, while we have tentative evidence that cluster-cluster mergers would enhance AGN activity in the outskirts of (particularly massive) galaxy clusters., Comment: 16 pages, 21 figures, and 2 tables, accepted for publication in PASJ

We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $\Phi^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $\beta = -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $\alpha = -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization., Comment: The Astrophysical Journal Letters, in press

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 assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that the lensing mass profile is consistent with that of massive galaxies ($M_{*}\sim 10.8~M_\odot$). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into low and high stellar mass hosts. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tend to suppress AGN activity, probably due to the lack of available gas., Comment: 17 pages,10 figures, accepted for publication in the Astrophysical Journal

We investigate the structural properties, distribution and abundance of LCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify -- for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass -- dependencies of subhalo properties with mass, maximum circular velocity, Vmax, host halo mass and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes in the expected range of values and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. We explore the distribution of subhaloes within their hosts and its evolution over cosmic time for subhaloes located as deep as ~0.1 per cent of the host virial radius. Subhalo structural properties are codified via a concentration parameter, cV, that does not depend on any specific, pre-defined density profile and relies only on Vmax. We derive the cV-Vmax relation in the range 7-1500 km/s and find an important dependence on distance of the subhalo to the host halo centre, as already described in Molin\'e et al. (2017). Interestingly, we also find subhaloes of the same mass to be significantly more concentrated into more massive hosts. Finally, we investigate the redshift evolution of cV, and provide accurate fits that take into account all mentioned dependencies. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Our work enables precision work in any future research involving dark matter halo substructure., Comment: 17 pages, 15 figures, 5 tables. Comments welcome!

We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with particle mass $8.97 \times 10^{5}$ Msun/h. Combining these simulations we can follow the evolution of dark matter halos and subhalos spanning those hosting dwarf galaxies to massive galaxy clusters across an unprecedented volume. In this first paper, we present basic statistics, dark matter power spectra, and the halo and subhalo mass functions, which demonstrate the wide dynamic range and superb statistics of the Uchuu suite. From an analysis of the evolution of the power spectra we conclude that our simulations remain accurate from the Baryon Acoustic Oscillation scale down to the very small. We also provide parameters of a mass-concentration model, which describes the evolution of halo concentration and reproduces our simulation data to within 5 per cent for halos with masses spanning nearly eight orders of magnitude at redshift 00.5, whereas no upturn is detected at z<0.5. We make publicly available various $N$-body products as part of Uchuu Data Release 1 on the Skies & Universes site. Future releases will include gravitational lensing maps and mock galaxy, X-ray cluster, and active galactic nuclei catalogues., Comment: 22 pages, 14 figures, accepted by MNRAS. We release various $N$-body products as data release 1 on http://skiesanduniverses.org/ such as subsets of simulation particles, matter power spectra, halo/subhalo catalogues, and their merger trees

We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions., Comment: 10 pages, 6 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

The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering., Comment: 20 pages, 13 figures, accepted in MNRAS. Error estimates of the 2PCFs and the fitting to them have been improved. Further discussions and references have been added. Main results remain unchanged

We present the clustering properties of low-$z$ $(z\leq1.4)$ galaxies selected by the Hyper Suprime-Cam Subaru Strategic Program Wide layer over $145$ deg$^{2}$. The wide-field and multi-wavelength observation yields $5,064,770$ galaxies at $0.3\leq z\leq1.4$ with photometric redshifts and physical properties. This enables the accurate measurement of angular correlation functions and subsequent halo occupation distribution (HOD) analysis allows the connection between baryonic properties and dark halo properties. The fraction of less-massive satellite galaxies at $z\lesssim1$ is found to be almost constant at $\sim20\%$, but it gradually decreases beyond $M_{\star} \sim 10^{10.4}h^{-2}M_{\odot}$. However, the abundance of satellite galaxies at $z>1$ is quite small even for less-massive galaxies due to the rarity of massive centrals at high-$z$. This decreasing trend is connected to the small satellite fraction of Lyman break galaxies at $z>3$. The stellar-to-halo mass ratios at $0.3\leq z\leq1.4$ are almost consistent with the predictions obtained using the latest empirical model; however, we identify small excesses from the theoretical model at the massive end. The pivot halo mass is found to be unchanged at $10^{11.9-12.1}h^{-1}M_{\odot}$ at $0.3\leq z\leq1.4$, and we systematically show that $10^{12}h^{-1}M_{\odot}$ is a universal pivot halo mass up to $z\sim5$ that is derived using only the clustering/HOD analyses. Nevertheless, halo masses with peaked instantaneous baryon conversion efficiencies are much smaller than the pivot halo mass regardless of a redshift, and the most efficient stellar-mass assembly is thought to be in progress in $10^{11.0-11.5}h^{-1}M_{\odot}$ dark haloes., Comment: 33 pages, 16 figures, submitted to ApJ

We propose a deep-learning approach based on generative adversarial networks (GANs) to reduce noise in weak lensing mass maps under realistic conditions. We apply image-to-image translation using conditional GANs to the mass map obtained from the first-year data of Subaru Hyper Suprime-Cam (HSC) survey. We train the conditional GANs by using 25000 mock HSC catalogues that directly incorporate a variety of observational effects. We study the non-Gaussian information in denoised maps using one-point probability distribution functions (PDFs) and also perform matching analysis for positive peaks and massive clusters. An ensemble learning technique with our GANs is successfully applied to reproduce the PDFs of the lensing convergence. About $60\%$ of the peaks in the denoised maps with height greater than $5\sigma$ have counterparts of massive clusters within a separation of 6 arcmin. We show that PDFs in the denoised maps are not compromised by details of multiplicative biases and photometric redshift distributions, nor by shape measurement errors, and that the PDFs show stronger cosmological dependence compared to the noisy counterpart. We apply our denoising method to a part of the first-year HSC data to show that the observed mass distribution is statistically consistent with the prediction from the standard $\Lambda$CDM model., Comment: 16 pages, 17 figures, 1 table. Accepted for publication in MNRAS

We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the "slowing down" of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios., Comment: 11 pages, 13 figures. Accepted for publication in MNRAS

The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between red-spiral galaxies and galaxy clusters has a dip at $v_{\rm{los}} = 0$, which is a very odd feature, at 93\% confidence level. To understand its origin, we construct a model of the phase space distribution of galaxies surrounding galaxy clusters in three-dimensional space by using cosmological $N$-body simulations. We adopt a two component model that consists of the infall component, which corresponds to galaxies that are now falling into galaxy clusters, and the splashback component, which corresponds to galaxies that are on their first (or more) orbit after falling into galaxy clusters. We find that we can reproduce the distribution of the line-of-sight velocity of red-spiral galaxies with the dip with a very simple assumption that red-spiral galaxies reside predominantly in the infall component, regardless of the choice of the functional form of their spatial distribution. Our results constrain the quenching timescale of red-spiral galaxies to a few Gyrs, and the radius where the morphological transformation is effective as $r \sim 0.2 h^{-1} \rm{Mpc}$., Comment: 10 pages, 9 figures, accepted for publication in MNRAS

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds., Comment: 46 pages, 47 figures; version accepted for publication in ApJ

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

We present morphological properties of dusty star-forming galaxies at z=1-3 determined with high-resolution (FWHM~0"19) Atacama Large Milllimeter/submilimeter Array (ALMA) 1-mm band maps of our ASAGAO survey covering a 26-arcmin^2 area in GOODS-S. In conjunction with the ALMA archival data, the present sample consists of 42 ALMA sources with a wide rest-frame far-infrared (FIR) luminosity L_FIR range of ~10^11-10^13 Lo. To obtain an average rest-frame FIR profile, we perform individual measurements and careful stacking of the ALMA sources using the uv-visibility method that includes positional-uncertainty and smoothing-effect evaluations through Monte-Carlo simulations. We find that the dusty star-forming galaxies have the average FIR-wavelength Sersic index and effective radius of n_FIR=1.2+/-0.2 and R_e,FIR=1.0-1.3 kpc, respectively, additionally with a point source at the center, indicative of the existence of AGN. The average FIR profile agrees with a morphology of an exponential-disk clearly distinguished from a spheroidal profile (Sersic index of 4). We also examine the rest-frame optical Sersic index n_opt and effective radius R_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain n_opt=0.9+/-0.3 (~n_FIR) and R_e,opt=3.2+/-0.6 kpc (>R_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs (n_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation., Comment: 13 pages, 6 figures, ApJ in press

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location., Comment: 18pages, 10 figures, accepted as PASJ special issue

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 examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73

We present measurements of the clustering properties of a sample of infrared (IR) bright dust-obscured galaxies (DOGs). Combining 125 deg$^2$ of wide and deep optical images obtained with the Hyper Suprime-Cam on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer, we have discovered 4,367 IR-bright DOGs with $(i - [22])_{\rm AB}$ $>$ 7.0 and flux density at 22 $\mu$m $>$ 1.0 mJy. We calculate the angular autocorrelation function (ACF) for a uniform subsample of 1411 DOGs with 3.0 mJy $<$ flux (22 $mu$m) $<$ 5.0 mJy and $i_{\rm AB}$ $<$ 24.0. The ACF of our DOG subsample is well-fit with a single power-law, $\omega (\theta)$ = (0.010 $\pm$ 0.003) $\theta^{-0.9}$, where $\theta$ in degrees. The correlation amplitude of IR-bright DOGs is larger than that of IR-faint DOGs, which reflects a flux-dependence of the DOG clustering, as suggested by Brodwin et al. (2008). We assume that the redshift distribution for our DOG sample is Gaussian, and consider 2 cases: (1) the redshift distribution is the same as IR-faint DOGs with flux at 22 $\mu$m $<$ 1.0 mJy, mean and sigma $z$ = 1.99 $\pm$ 0.45, and (2) $z$ = 1.19 $\pm$ 0.30, as inferred from their photometric redshifts. The inferred correlation length of IR-bright DOGs is $r_0$ = 12.0 $\pm$ 2.0 and 10.3 $\pm$ 1.7 $h^{-1}$ Mpc, respectively. IR-bright DOGs reside in massive dark matter halos with a mass of $\log [\langle M_{\mathrm{h}} \rangle / (h^{-1} M_{\odot})]$ = 13.57$_{-0.55}^{+0.50}$ and 13.65$_{-0.52}^{+0.45}$ in the two cases, respectively., Comment: 15 pages, 8 figures, and 1 table, accepted for publication in ApJ

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 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 perform dry merger simulations to investigate the role of dry mergers in the size growth of early-type galaxies in high density environments. We replace the virialized dark matter haloes obtained by a large cosmological $N$-body simulation with $N$-body galaxy models consisting of two components, a stellar bulge and a dark matter halo, which have higher mass resolution than the cosmological simulation. We then re-simulate nine cluster forming regions, whose masses range from 1e+14 Msun to 5e+14 Msun. Masses and sizes of stellar bulges are also assumed to satisfy the stellar mass--size relation of high-z compact massive early-type galaxies. We find that dry major mergers considerably contribute to the mass and size growth of central massive galaxies. One or two dry major mergers double the average stellar mass and quadruple the average size between $z=2$ and $z=0$. These growths favorably agree with observations. Moreover, the density distributions of our simulated central massive galaxies grow from the inside-out, which is consistent with recent observations. The mass--size evolution is approximated as R propto M_{*}^{alpha}, with alpha ~ 2.24. Most of our simulated galaxies are efficiently grown by dry mergers, and their stellar mass--size relations match the ones observed in the local Universe. Our results show that the central galaxies in the cluster haloes are potential descendants of high-z (z ~ 2-3) compact massive early-type galaxies. This conclusion is consistent with previous numerical studies which investigate the formation and evolution of compact massive early-type galaxies., Comment: 15 pages, 10 figures, accepted for MNRAS; adjusted to match the published version

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 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 report the discovery of a new ultra-bright submillimeter galaxy (SMG) behind the Small Magellanic Cloud (SMC). This SMG is detected as a 43.3+-8.4 mJy point source (MM J01071-7302, hereafter MMJ0107) in the 1.1 mm continuum survey of the SMC by AzTEC on the ASTE telescope. MMJ0107 is also detected in the radio (843 MHz), Herschel/SPIRE, Spitzer MIPS 24 {\mu}m, all IRAC bands, Wide-field Infrared Survey Explorer, and near-infrared (J, H, KS). We find an optical (U, B, V) source, which might be the lensing object, at a distance of 1.4 arcsec from near-infrared and IRAC sources. Photometric redshift estimates for the SMG using representative spectral energy distribution templates show the redshifts of 1.4-3.9. We estimate total far-infrared luminosity of (0.3-2.2)x10^14 {\mu}^-1 L_sun and a star formation rate of 5600-39, 000 {\mu}^-1 M_sun yr^-1, where {\mu} is the gravitational magnification factor. This apparent extreme star formation activity is likely explained by a highly magnified gravitational lens system., Comment: 14 pages, 2 figures, published in ApJ Letters

Recent observations show evidence that high-z (z\sim 2 - 3) early-type galaxies (ETGs) are more compact than those with comparable mass at z\sim 0. Such a size evolution is most likely explained by the `Dry Merger Sceanario'. However, previous studies based on this scenario are not able to consistantly explain both the properties of the high-z compact massive ETGs and the local ETGs. We investigate the effect of multiple sequential dry minor mergers on the size evolution of the compact massive ETGs. From an analysis of the Millennium Simulation Database, we show that such minor (stellar mass ratio $M_{2}/M_{1} < 1/4$) mergers are extremely common during hierarchical structure formation. We perform N-body simulations of sequential minor mergers with parabolic and head-on orbits, including a dark matter component and a stellar component. Typical mass ratios of the minor mergers are $1/20 < M_{2}/M_{1} < 1/10$. We show that sequential minor mergers of compact satellite galaxies are the most efficient at promoting size growth and decreasing the velocity dispersion of the compact massive ETGs in our simulations. The change of stellar size and density of the merger remnants is consistent with recent observations. Furthermore, we construct the merger histories of candidates for the high-z compact massive ETGs using the Millennium Simulation Database, and estimate the size growth of the galaxies by the dry minor merger scenario. We can reproduce the mean size growth factor between $z=2$ and $z=0$, assuming the most efficient size growth obtained during sequential minor mergers in our simulations. However, we note that our numerical result is only valid for merger histories with typical mass ratios between 1/20 and 1/10 with parabolic and head-on orbits, and that our most efficient size growth efficiency is likely to an upper limit., Comment: 23 pages, 15 figures, accepted for MNRAS, Fig. 5 is fixed