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.

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.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

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.

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<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

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.

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.

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.

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.

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars at high redshift for the first time. In this study, we use the two-point angular cross-correlation function (CCF) for both low-luminosity quasars at $-24<M_{\rm 1450}<-22$ and Lyman-break galaxies (LBGs). Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low-luminosity quasars) in the redshift range $3.1<z<4.5$ and 835 color-selected LBGs with $z^{\prime}_{\rm LBG}<25.0$ at $z\sim4$ in the COSMOS field. We have made our analysis for the following two quasar samples; (1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample (the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity quasars at $z\sim4$ is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. We then obtain the $86\%$ upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the typical dark matter halo masses, log $(M_{\rm DM}/(h^{-1}M_{\odot}))=$$12.7$ and $13.5$, respectively. This result is not inconsistent with the predicted bias for quasars which is estimated by the major merger models.

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.

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.

We present the result of our low-luminosity quasar survey in the redshift range of 4.5 < z < 5.5 in the COSMOS field. Using the COSMOS photometric catalog, we selected 15 quasar candidates with 22 < i' < 24 at z~5, that are ~ 3 mag fainter than the SDSS quasars in the same redshift range. We obtained optical spectra for 14 of the 15 candidates using FOCAS on the Subaru Telescope and did not identify any low-luminosity type-1 quasars at z~5 while a low-luminosity type-2 quasar at z~5.07 was discovered. In order to constrain the faint end of the quasar luminosity function at z~5, we calculated the 1sigma confidence upper limits of the space density of type-1 quasars. As a result, the 1sigma confidence upper limits on the quasar space density are Phi< 1.33*10^{-7} Mpc^{-3} mag^{-1} for -24.52 < M_{1450} < -23.52 and Phi< 2.88*10^{-7} Mpc^{-3} mag^{-1} for -23.52 < M_{1450} < -22.52. The inferred 1sigma confidence upper limits of the space density are then used to provide constrains on the faint-end slope and the break absolute magnitude of the quasar luminosity function at z~5. We find that the quasar space density decreases gradually as a function of redshift at low luminosity (M_{1450} ~ -23), being similar to the trend found for quasars with high luminosity (M_{1450}<-26). This result is consistent with the so-called downsizing evolution of quasars seen at lower redshifts.

We searched for quasars that are ~ 3 mag fainter than the SDSS quasars in the redshift range 3.7 < z < 4.7 in the COSMOS field to constrain the faint end of the quasar luminosity function. Using optical photometric data, we selected 31 quasar candidates with 22 < i' < 24 at z ~ 4. We obtained optical spectra for most of these candidates using FOCAS on the Subaru telescope, and identified 8 low-luminosity quasars at z ~ 4. In order to derive the quasar luminosity function (QLF) based on our spectroscopic follow-up campaign, we estimated the photometric completeness of our quasar survey through detailed Monte Carlo simulations. Our QLF at z ~ 4 has a much shallower faint-end slope beta = -1.67^{+0.11}_{-0.17} than that obtained by other recent surveys in the same redshift. Our result is consistent with the scenario of downsizing evolution of active galactic nuclei inferred by recent optical and X-ray quasar surveys at lower redshifts.

We investigate whether the mean star formation activity of star-forming galaxies from z=0 to z=0.7 in the GOODS-S field can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in galaxies. We selected star-forming galaxies at z=0 and at z=0.7 in IR and in UV to have access to all the recent star formation. We focused on galaxies with a stellar mass ranging between 10^{10} and 10^{11} M_sun for which the results are not biased by the selections. We compared the data to chemical evolution models developed for spiral galaxies and originally built to reproduce the main characteristics of the Milky Way and nearby spirals without fine-tuning them for the present analysis. We find a shallow decrease in the specific star formation rate (SSFR) when the stellar mass increases. The evolution of the SSFR characterizing both UV and IR selected galaxies from z=0 to z=0.7 is consistent with the models built to reproduce the present spiral galaxies. There is no need to strongly modify of the physical conditions in galaxies to explain the average evolution of their star formation from z=0 to z=0.7. We use the models to predict the evolution of the star formation rate and the metallicity on a wider range of redshift and we compare these predictions with the results of semi-analytical models.

A compact binary in an eccentric orbit radiates gravitational waves (GWs) at all integer harmonics of its orbital frequency. In this study, we investigate the effect of orbital eccentricity on the expected gravitational background radiation (GWBR) from supermassive black hole (SMBH) binaries in the nuclei of galaxies. For this purpose, we formulate a power spectrum of the GWBR from cosmological evolving eccentric binaries. Then, we apply this formulation to the case of the GWBR from SMBH binaries. The key to doing this is to correctly estimate the number density of coalescing SMBH binaries. In this study, we use a semi-analytic model of galaxy and SMBH formation. We find that the power spectrum of the GWBR from SMBH binaries on eccentric orbits is suppressed for frequencies $f \lesssim 1~{\rm nHz}$ if the initial eccentricity, $e_0$, satisfies $e_0 > 0.2$ and the initial semi-major axis is 300 times Scwarzschild radius. Our model predicts that while the overall shape and amplitude of the power spectrum depend strongly on the processes of galaxy formation, the eccentricity of binaries can affect the shape of the power spectrum for lower frequencies, i.e., $f \lesssim 1~{\rm nHz}$. Pulsar timing measurements, which can detect GW in this frequency range, could constrain the effect of eccentricity on the power spectrum of the GWBR from SMBH binaries.

We construct the Numerical Galaxy Catalog ($\nu$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $\Lambda$-dominated flat cold dark matter ($\Lambda$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, star formation, heating by supernova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal dust and intervening HI clouds. As the first paper in a series using this model, we focus on basic photometric, structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters are examined, strong and weak supernova feedback models, which are in good agreement with observational luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, HI mass functions, galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isophotal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in much better agreement with near-infrared (K'-band) faint galaxy number counts and redshift distribution than the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter formalism. (Abridged)

We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts rates owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 \mu {\rm Hz})^{-2/3}$ for $f \lesssim 1 \mu {\rm Hz}$. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$.

We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. When galaxies merge, the SMBHs in the host galaxies sink to the center of the new merged galaxy and form a binary system. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 \mu {\rm Hz})^{-2/3}$ for $f \lesssim 1 \mu {\rm Hz}$ just below the detection limit from measurements of the pulsar timing provided that SMBHs coalesce simultaneously when host galaxies merge. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$.

We investigate the environments of quasars such as number distribution of galaxies using a semi-analytic model which includes both galaxy and quasar formations based on the hierarchical clustering scenario. We assume that a supermassive black hole is fueled by accretion of cold gas and that it is a source of quasar activity during a major merger of the quasar host galaxy with another galaxy. This major merger causes spheroid formation of the host galaxy. Our model can reproduce not only general form of the galaxy luminosity functions in the local Universe but also the observed relation between a supermassive black hole mass and a spheroid luminosity, the present black hole mass function and the quasar luminosity functions at different redshifts. Using this model, we predict the mean number of quasars per halo, bias parameter of quasars and the probability distribution of the number of galaxies around quasars. In our model, analysis of the mean number of quasars per halo shows that the spatial distribution of galaxies is different from that of quasars. Furthermore, we found from calculation of the probability distribution of galaxy numbers that at $0.2 \lesssim z \lesssim 0.5$, most quasars are likely to reside in galaxy groups. On the other hand, at $1 \lesssim z \lesssim 2$ most quasars seem to reside in more varied environments than at a lower redshift; quasars reside in environments ranging from small groups of galaxies to clusters of galaxies. Comparing these predictions with observations in future will enable us to constrain our quasar formation model.

We investigate statistical properties of galaxy clusters in the context of hierarchical clustering scenario, taking account of their formation epoch distribution, motivated by the recent finding by Fujita and Takahara that X-ray clusters form a fundamental plane, where the mass and formation epoch are regarded as two independent parameters. Using the formalism which discriminates between major merger and accretion, the epoch of a cluster formation is identified with that of the last major merger. Since tiny mass accretion after the formation does not much affect the core structure of clusters, the properties of X-ray emission from clusters are determined by the total mass and density at their formation time. Under these assumptions, we calculate X-ray luminosity and temperature functions of galaxy clusters. We find that the behavior of luminosity function is different from the model which does not take account of formation epoch distribution, while the behavior of temperature function is not much changed. In our model, luminosity function is shifted to a higher luminosity and shows no significant evolution up to $z \sim 1$, independent of cosmological models. The clusters are populated on the temperature-luminosity plane with a finite dispersion. Since the simple scaling model in which the gas temperature is equal to the virial temperature fails to reproduce the observed luminosity-temperature relation, we also consider a model which takes the effects of preheating into account. The preheating model reproduces the observations much better.

We have investigated the effect of pressure from intracluster medium (ICM) on disk galaxies in merging clusters. The ram-pressure on the galaxies rapidly increases when two clusters collide. This leads to stripping of the interstellar medium (ISM) and decrease of star formation rate (SFR) of the galaxies. On the other hand, the increase of SFR caused by compression of ISM is less significant. Thus, cluster merger does not trigger, but weakens star formation activity of the galaxies. In the central region of the colliding clusters, blue galaxies with high velocity should exist, although most of galaxies become red. Following the decrease of blue galaxy fraction in the clusters, the fraction of post-starburst galaxies increases. After merger, many galaxies in the cluster restart star formation activity and the segregation of blue and red galaxies becomes prominent.