As cosmic structures form, matter density fluctuations collapse gravitationally and baryonic matter is shock-heated and thermalized. We therefore expect a connection between the mean gravitational potential energy density of collapsed halos, $\Omega_{W}^{\rm halo}$, and the mean thermal energy density of baryons, $\Omega_{\rm th}$. These quantities can be obtained using two fundamentally different estimates: we compute $\Omega_{W}^{\rm halo}$ using the theoretical framework of the halo model which is driven by dark matter statistics, and measure $\Omega_{\rm th}$ using the Sunyaev-Zeldovich (SZ) effect which probes the mean thermal pressure of baryons. First, we derive that, at the present time, about 90% of $\Omega_{W}^{\rm halo}$ originates from massive halos with $M>10^{13}\,M_\odot$. Then, using our measurements of the SZ background, we find that $\Omega_{\rm th}$ accounts for about 80% of the kinetic energy of the baryons available for pressure in halos at $z\lesssim 0.5$. This constrains the amount of non-thermal pressure, e.g., due to bulk and turbulent gas motion sourced by mass accretion, to be about $\Omega_{\rm non-th}\simeq 0.4\times 10^{-8}$ at $z=0$.

The cosmic thermal history, quantified by the evolution of the mean thermal energy density in the universe, is driven by the growth of structures as baryons get shock-heated in collapsing dark matter halos. This process can be probed by redshift-dependent amplitudes of the thermal Sunyaev-Zeldovich (SZ) effect background. To do so, we cross-correlate eight sky intensity maps in the $\it{Planck}$ and Infrared Astronomical Satellite missions with two million spectroscopic-redshift references in the Sloan Digital Sky Surveys. This delivers snapshot spectra for the far-infrared to microwave background light as a function of redshift up to $z\sim3$. We decompose them into the SZ and thermal dust components. Our SZ measurements directly constrain $\langle bP_{\rm e} \rangle$, the halo bias-weighted mean electron pressure, up to $z\sim 1$. This is the highest redshift achieved to date, with uncorrelated redshift bins thanks to the spectroscopic references. We detect a threefold increase in the density-weighted mean electron temperature $\bar{T}_{\rm{e}}$ from $7\times 10^5~{\rm K}$ at $z=1$ to $2\times 10^6~{\rm K}$ today. Over $z=1$-$0$, we witness the build up of nearly $70\%$ of the present-day mean thermal energy density $\rho_{\rm{th}}$, with the corresponding density parameter $\Omega_{\rm th}$ reaching $1.5 \times10^{-8}$. We find the mass bias parameter of $\it{Planck}$'s universal pressure profile of $B=1.27$ (or $1-b=1/B=0.79$), consistent with the magnitude of non-thermal pressure in gas motion and turbulence from mass assembly. We estimate the redshift-integrated mean Compton parameter $y\sim1.2\times10^{-6}$, which will be tested by future spectral distortion experiments. More than half of which originates from large-scale structure at $z<1$, which we detect directly.

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.

Thermal Sunyaev-Zeldovich (tSZ) power spectrum is a powerful probe of the present-day amplitude of matter density fluctuations, and has been measured up to $\ell\approx 10^3$ from the {\it Planck} data. The largest systematic uncertainty in the interpretation of this data is the so-called "mass bias" parameter $B$, which relates the true halo mass to the mass proxy used by the {\it Planck} team as $M_{\rm 500c}^{Planck}=M_{\rm 500c}^{\rm true}/B$. Since the power spectrum of the cosmic weak lensing shear is also sensitive to the amplitude of matter density fluctuations via $S_8\equiv \sigma_8 \Omega_m^{\alpha}$ with $\alpha\sim 0.5$, we can break the degeneracy between the mass bias and the cosmological parameters by combining the tSZ and cosmic shear power spectra. In this paper, we perform a joint likelihood analysis of the tSZ power spectrum from {\it Planck} and the cosmic shear power spectrum from Subaru Hyper Suprime-Cam. Our analysis does not use the primordial cosmic microwave background (CMB) information. We obtain a new constraint on the mass bias as $B = 1.37 ^{+0.15}_{-0.23}$ or $(1-b) = B^{-1}=0.73^{+0.08}_{-0.13}$ (68\%~C.L.), for $\sigma_8 < 0.9$. This value of $B$ is lower than that needed to reconcile the tSZ data with the primordial CMB and CMB lensing data, i.e., $B = 1.64 \pm 0.19$, but is consistent with the mass bias expected from hydrodynamical simulations, $B = 1.28 \pm 0.20$. Our results thus indicate that the mass bias is consistent with the non-thermal pressure support from mass accretion of galaxy clusters via the cosmic structure formation, and that the cosmologies inferred from the tSZ and the cosmic shear are consistent with each other.

We examine the impact of fiber assignment on clustering measurements from fiber-fed spectroscopic galaxy surveys. We identify new effects which were absent in previous, relatively shallow galaxy surveys such as Baryon Oscillation Spectroscopic Survey . Specifically, we consider deep surveys covering a wide redshift range from z=0.6 to z=2.4, as in the Subaru Prime Focus Spectrograph survey. Such surveys will have more target galaxies than we can place fibers on. This leads to two effects. First, it eliminates fluctuations with wavelengths longer than the size of the field of view, as the number of observed galaxies per field is nearly fixed to the number of available fibers. We find that we can recover the long-wavelength fluctuation by weighting galaxies in each field by the number of target galaxies. Second, it makes the preferential selection of galaxies in under-dense regions. We mitigate this effect by weighting galaxies using the so-called individual inverse probability. Correcting these two effects, we recover the underlying correlation function at better than 1 percent accuracy on scales greater than 10 Mpc/h.

Galaxies evolve from a blue star-forming phase into a red quiescent one by quenching their star formation activity. In high density environments, this galaxy evolution proceeds earlier and more efficiently. Therefore, local galaxy clusters are dominated by well-evolved red, elliptical galaxies. The fraction of blue galaxies in clusters monotonically declines with decreasing redshift, i.e., the Butcher-Oemler effect. In the local Universe, observed blue fractions of massive clusters are as small as $\lesssim$ 0.2. Here we report a discovery of a \lq \lq blue cluster\rq \rq, that is a local galaxy cluster with an unprecedentedly high fraction of blue star-forming galaxies yet hosted by a massive dark matter halo. The blue fraction is 0.57, which is 4.0 $\sigma$ higher than those of the other comparison clusters under the same selection and identification criteria. The velocity dispersion of the member galaxies is 510 km s$^{-1}$, which corresponds to a dark matter halo mass of 2.0$^{+1.9}_{-1.0}\times 10^{14}$ M$_{\odot}$. The blue fraction of the cluster is more than 4.7 $\sigma$ beyond the standard theoretical predictions including semi-analytic models of galaxy formation. The probability to find such a high blue fraction in an individual cluster is only 0.003\%, which challenges the current standard frameworks of the galaxy formation and evolution in the $\Lambda$CDM Universe. The spatial distribution of galaxies around the blue cluster suggests that filamentary cold gas streams can exist in massive halos even in the local Universe. However these cold streams have already disappeared in the theoretically simulated local universes.

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.

We present a first joint analysis of the power spectra of the thermal Sunyaev-Zeldovich (tSZ) effect measured by the Planck and the number density fluctuations of galaxies in the 2MASS redshift survey (2MRS) catalog, including their cross-correlation. Combining these measurements with the cosmic microwave background (CMB) data and CMB lensing of Planck assuming a flat $\Lambda$CDM model, we constrain the mass bias parameter as $B = 1.54 \pm 0.098\;(1\sigma)$ [$(1-b) = 0.649 \pm 0.041$, where $(1-b) \equiv B^{-1}$], i.e., the Planck cluster mass should be $35\%$ lower than the true mass. The mass bias determined by the 2MRS-tSZ cross-power spectrum alone is consistent with that determined by the tSZ auto-power spectrum alone, suggesting that this large mass bias is not due to obvious systematics in the tSZ data. We find that the 2MRS-tSZ cross-power spectrum is more sensitive to less massive halos than the tSZ auto-power spectrum and it significantly improves a constraint on the mass dependence of the mass bias. The redshift dependence is not strongly constrained since the multipole range in which high redshift clusters mainly contribute to the tSZ auto is dominated by the contaminating sources. We conclude that no strong mass or redshift evolution of the mass bias is needed to explain the data.

We present our ALMA Cycle 4 measurements of the [CII] emission line and the underlying far-infrared (FIR) continuum emission from four optically low-luminosity ($M_{\rm 1450} > -25$) quasars at $z \gtrsim 6$ discovered by the Subaru Hyper Suprime Cam (HSC) survey. The [CII] line and FIR continuum luminosities lie in the ranges $L_{\rm [CII]} = (3.8-10.2) \times 10^8~L_\odot$ and $L_{\rm FIR} = (1.2-2.0) \times 10^{11}~L_\odot$, which are at least one order of magnitude smaller than those of optically-luminous quasars at $z \gtrsim 6$. We estimate the star formation rates (SFR) of our targets as $\simeq 23-40~M_\odot ~{\rm yr}^{-1}$. Their line and continuum-emitting regions are marginally resolved, and found to be comparable in size to those of optically luminous quasars, indicating that their SFR or likely gas mass surface densities (key controlling parameter of mass accretion) are accordingly different. The $L_{\rm [CII]}/L_{\rm FIR}$ ratios of the hosts, $\simeq (2.2-8.7) \times 10^{-3}$, are fully consistent with local star-forming galaxies. Using the [CII] dynamics, we derived their dynamical masses within a radius of 1.5-2.5 kpc as $\simeq (1.4-8.2) \times 10^{10}~M_\odot$. By interpreting these masses as stellar ones, we suggest that these faint quasar hosts are on or even below the star-forming main sequence at $z \sim 6$, i.e., they appear to be transforming into quiescent galaxies. This is in contrast to the optically luminous quasars at those redshifts, which show starburst-like properties. Finally, we find that the ratios of black hole mass to host galaxy dynamical mass of the most of low-luminosity quasars including the HSC ones are consistent with the local value. The mass ratios of the HSC quasars can be reproduced by a semi-analytical model that assumes merger-induced black hole-host galaxy evolution.

We present a public code to generate a mock galaxy catalog in redshift space assuming a log-normal probability density function (PDF) of galaxy and matter density fields. We draw galaxies by Poisson-sampling the log-normal field, and calculate the velocity field from the linearised continuity equation of matter fields, assuming zero vorticity. This procedure yields a PDF of the pairwise velocity fields that is qualitatively similar to that of N-body simulations. We check fidelity of the catalog, showing that the measured two-point correlation function and power spectrum in real space agree with the input precisely. We find that a linear bias relation in the power spectrum does not guarantee a linear bias relation in the density contrasts, leading to a cross-correlation coefficient of matter and galaxies deviating from unity on small scales. We also find that linearising the Jacobian of the real-to-redshift space mapping provides a poor model for the two-point statistics in redshift space. That is, non-linear redshift-space distortion is dominated by non-linearity in the Jacobian. The power spectrum in redshift space shows a damping on small scales that is qualitatively similar to that of the well-known Fingers-of-God (FoG) effect due to random velocities, except that the log-normal mock does not include random velocities. This damping is a consequence of non-linearity in the Jacobian, and thus attributing the damping of the power spectrum solely to FoG, as commonly done in the literature, is misleading.

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 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 present stacking analyses on our ALMA deep 1.1 mm imaging in the SXDF using 1.6 {\mu}m and 3.6 {\mu}m selected galaxies in the CANDELS WFC3 catalog. We detect a stacked flux of ~0.03-0.05 mJy, corresponding to LIR < 10^11 Lsun and a star formation rate (SFR) of ~ 15 Msun/yr at z = 2. We find that galaxies brighter in the rest-frame near-infrared tend to be also brighter at 1.1 mm, and galaxies fainter than m[3.6um] = 23 do not produce detectable 1.1 mm emission. This suggests a correlation between stellar mass and SFR, but outliers to this correlation are also observed, suggesting strongly boosted star formation or extremely large extinction. We also find tendencies that redder galaxies and galaxies at higher redshifts are brighter at 1.1 mm. Our field contains z ~ 2.5 H-alpha emitters and a bright single-dish source. However, we do not find evidence of bias in our results caused by the bright source. By combining the fluxes of sources detected by ALMA and fluxes of faint sources detected with stacking, we recover a 1.1 mm surface brightness of up to 20.3 +/- 1.2 Jy/deg, comparable to the extragalactic background light measured by COBE. Based on the fractions of optically faint sources in our and previous ALMA studies and the COBE measurements, we find that approximately half of the cosmic star formation may be obscured by dust and missed by deep optical surveys, Much deeper and wider ALMA imaging is therefore needed to better constrain the obscured cosmic star formation history.

The anthropic principle is one of the possible explanations for the cosmological constant ($\Lambda$) problem. In previous studies, a dark halo mass threshold comparable with our Galaxy must be assumed in galaxy formation to get a reasonably large probability of finding the observed small value, $P(<$$\Lambda_{\rm obs})$, though stars are found in much smaller galaxies as well. Here we examine the anthropic argument by using a semi-analytic model of cosmological galaxy formation, which can reproduce many observations such as galaxy luminosity functions. We calculate the probability distribution of $\Lambda$ by running the model code for a wide range of $\Lambda$, while other cosmological parameters and model parameters for baryonic processes of galaxy formation are kept constant. Assuming that the prior probability distribution is flat per unit $\Lambda$, and that the number of observers is proportional to stellar mass, we find $P(<$$\Lambda_{\rm obs}) = 6.7 \%$ without introducing any galaxy mass threshold. We also investigate the effect of metallicity; we find $P(<$$\Lambda_{\rm obs}) = 9.0 \%$ if observers exist only in galaxies whose metallicity is higher than the solar abundance. If the number of observers is proportional to metallicity, we find $P(<$$\Lambda_{\rm obs}) = 9.7 \%$. Since these probabilities are not extremely small, we conclude that the anthropic argument is a viable explanation, if the value of $\Lambda$ observed in our universe is determined by a probability distribution.

We present a multi-wavelength analysis of five submillimeter sources (S_1.1mm = 0.54-2.02 mJy) that were detected during our 1.1-mm-deep continuum survey in the SXDF-UDS-CANDELS field (2 arcmin^2, 1sigma = 0.055 mJy beam^-1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be ~ 4.1^{+5.4}_{-3.0} Jy deg^{-2}, which corresponds to ~ 16^{+22}_{-12}% of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multi-wavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity (L_IR ~ 1e12 L_sun or SFR ~ 100 M_sun yr^{-1}). By fitting the spectral energy distributions (SEDs) at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (z_photo) and stellar masses (M_*): z_photo ~ 1.3-2.5, M_* ~ (3.5-9.5)e10 M_sun. We also derived their star formation rates (SFRs) and specific SFRs (sSFRs) as ~ 30-200 M_sun yr^{-1} and ~ 0.8-2 Gyr^{-1}, respectively. These values imply that they are main-sequence star-forming galaxies.

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 report 1.1 mm number counts revealed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey Field (SXDF). The advent of ALMA enables us to reveal millimeter-wavelength number counts down to the faint end without source confusion. However, previous studies are based on the ensemble of serendipitously-detected sources in fields originally targeting different sources and could be biased due to the clustering of sources around the targets. We derive number counts in the flux range of 0.2-2 mJy by using 23 (>=4sigma) sources detected in a continuous 2.0 arcmin$^2$ area of the SXDF. The number counts are consistent with previous results within errors, suggesting that the counts derived from serendipitously-detected sources are not significantly biased, although there could be field-to-field variation due to the small survey area. By using the best-fit function of the number counts, we find that ~40% of the extragalactic background light at 1.1 mm is resolved at S(1.1mm) > 0.2 mJy.

We have conducted 1.1 mm ALMA observations of a contiguous $105'' \times 50''$ or 1.5 arcmin$^2$ window in the SXDF-UDS-CANDELS. We achieved a 5$\sigma$ sensitivity of 0.28 mJy, providing a flat sensus of dusty star-forming galaxies with $L_{\rm IR} \sim6\times10^{11}$ $L_\odot$ (for $T_{\rm dust}$ =40K) up to $z\sim10$ thanks to the negative K-correction at this wavelength. We detected 5 brightest sources (S/N$>$6) and 18 low-significance sources (5$>$S/N$>$4; these may contain spurious detections, though). One of the 5 brightest ALMA sources ($S_{\rm 1.1mm} = 0.84 \pm 0.09$ mJy) is extremely faint in the WFC3 and VLT/HAWK-I images, demonstrating that a contiguous ALMA imaging survey is able to uncover a faint dust-obscured population that is invisible in deep optical/near-infrared surveys. We found a possible [CII]-line emitter at $z=5.955$ or a low-$z$ CO emitting galaxy within the field, which may allow us to constrain the [CII] and/or the CO luminosity functions across the history of the universe.

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.

We report the results of $1^{\prime}.5 \times3^{\prime}$ mapping at 1.1~mm with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the central region of the $z=3.09$ SSA22 protocluster. By combining our source catalog with archival spectroscopic redshifts, we find that eight submillimeter galaxies (SMGs) with flux densities, $S_{\rm 1.1~mm}=0.7-6.4$~mJy ($L_{\rm IR}\sim10^{12.1}-10^{13.1}L_\odot$) are at $z=3.08-3.10$. Not only are these SMGs members of the protocluster but they in fact reside within the node at the junction of the 50 Mpc-scale filamentary three-dimensional structure traced by Lyman-$\alpha$ emitters (LAEs) in this field. The eight SMGs account for a star formation rate density (SFRD) $\sim$10 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ in the node, which is two orders of magnitudes higher than the global SFRD at this redshift. We find that four of the eight SMGs host a X-ray luminous active galactic nuclei (AGN). Our results suggest that the vigorous star formation activity and the growth of super massive black holes (SMBHs) occurred simultaneously in the densest regions at $z\sim3$, which may correspond to the most active historical phase of the massive galaxy population found in the core of the clusters in the present universe. Two SMGs are associated with Lyman-$\alpha$ blobs (LABs), implying that the two populations coexist in high density environments for a few cases.

We present first results from the SXDF-ALMA 1.5 arcmin^2 deep survey at 1.1 mm using Atacama Large Millimeter Array (ALMA). The map reaches a 1sigma depth of 55 uJy/beam and covers 12 Halpha-selected star-forming galaxies at z = 2.19 or z=2.53. We have detected continuum emission from three of our Halpha-selected sample, including one compact star-forming galaxy with high stellar surface density, NB2315-07. They are all red in the rest-frame optical and have stellar masses of log (M*/Msun)>10.9 whereas the other blue, main-sequence galaxies with log(M*/Msun)=10.0-10.8 are exceedingly faint, <290 uJy (2sigma upper limit). We also find the 1.1 mm-brightest galaxy, NB2315-02, to be associated with a compact (R_e=0.7+-0.1 kpc), dusty star-forming component. Given high gas fraction (44^{+20}_{-8}% or 37^{+25}_{-3}%) and high star formation rate surface density (126^{+27}_{-30} Msun yr^{-1}kpc^{-2}), the concentrated starburst can within less than 50^{+12}_{-11} Myr build up a stellar surface density matching that of massive compact galaxies at z~2, provided at least 19+-3% of the total gas is converted into stars in the galaxy centre. On the other hand, NB2315-07, which already has such a high stellar surface density core, shows a gas fraction (23+-8%) and is located in the lower envelope of the star formation main-sequence. This compact less star-forming galaxy is likely to be in an intermediate phase between compact dusty star-forming and quiescent galaxies.

We report optical-infrared (IR) properties of faint 1.3 mm sources (S_1.3mm = 0.2-1.0 mJy) detected with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey (SXDS) field. We searched for optical/IR counterparts of 8 ALMA-detected sources (>=4.0 sigma, the sum of the probability of spurious source contamination is ~1) in a K-band source catalog. Four ALMA sources have K-band counterpart candidates within a 0.4" radius. Comparison between ALMA-detected and undetected K-band sources in the same observing fields shows that ALMA-detected sources tend to be brighter, more massive, and more actively forming stars. While many of the ALMA-identified submillimeter-bright galaxies (SMGs) in previous studies lie above the sequence of star-forming galaxies in stellar mass--star-formation rate plane, our ALMA sources are located in the sequence, suggesting that the ALMA-detected faint sources are more like `normal' star-forming galaxies rather than `classical' SMGs. We found a region where multiple ALMA sources and K-band sources reside in a narrow photometric redshift range (z ~ 1.3-1.6) within a radius of 5" (42 kpc if we assume z = 1.45). This is possibly a pre-merging system and we may be witnessing the early phase of formation of a massive elliptical galaxy.

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.

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.

We present a new model of the extragalactic background light (EBL) and corresponding gamma-gamma opacity for intergalactic gamma-ray absorption from z = 0 up to z = 10, based on a semi-analytical model of hierarchical galaxy formation that reproduces key observed properties of galaxies at various redshifts. Including the potential contribution from Population III stars and following the cosmic reionization history in a simplified way, the model is also broadly consistent with available data concerning reionization, particularly the Thomson scattering optical depth constraints from WMAP. In comparison with previous EBL studies up to z ~ 3-5, our predicted gamma-gamma opacity is in general agreement for observed gamma-ray energy below 400/(1 + z) GeV, whereas it is a factor of ~ 2 lower above this energy because of a correspondingly lower cosmic star formation rate, even though the observed UV luminosity is well reproduced by virtue of our improved treatment of dust obscuration and direct estimation of star formation rate. The horizon energy at which the gamma-ray opacity is unity does not evolve strongly beyond z ~ 4 and approaches ~ 20 GeV. The contribution of Population III stars is a minor fraction of the EBL at z = 0, and is also difficult to distinguish through gamma-ray absorption in high-z objects, even at the highest levels allowed by the WMAP constraints. Nevertheless, the attenuation due to Population II stars should be observable in high-z gamma-ray sources by telescopes such as Fermi or CTA and provide a valuable probe of the evolving EBL in the rest-frame UV. The detailed results of our model are publicly available in numerical form at the URL http://www.slac.stanford.edu/~yinoue/Download.html.

We present deep optical and X-ray follow-up observations of the bright unassociated Fermi-LAT gammaray source 1FGL J1311.7-3429. The source was already known as an unidentified EGRET source (3EG J1314-3431, EGR J1314-3417), hence its nature has remained uncertain for the past two decades. For the putative counterpart, we detected a quasi-sinusoidal optical modulation of delta_m\sim2 mag with a period of ~1.5 hr in the Rc, r' and g' bands. Moreover, we found that the amplitude of the modulation and peak intensity changed by > 1 mag and 0.5 mag respectively, over our total six nights of observations from 2012 March and May. Combined with Swif t UVOT data, the optical-UV spectrum is consistent with a blackbody temperature, kT \sim1 eV, and the emission volume radius Rbb\sim 1.5x10^4 km. In contrast, deep Suzaku observations conducted in 2009 and 2011 revealed strong X-ray flares with a lightcurve characterized with a power spectrum density of P(f) propto f^(-2) but the folded X-ray light curves suggest an orbital modulation also in X-rays. Together with the non-detection of a radio counterpart, and significant curved spectrum and non-detection of variability in gamma-rays, the source may be the second radio-quiet gamma-ray emitting milli-second pulsar candidate after 1FGL J2339.7-0531, although the origin of flaring X-ray and optical variability remains an open question.

We report X-ray observations of two galaxy clusters Abell 1555 and Abell 1558 with Suzaku, which are included in a large scale filamentary structure and a supercluster, to search for non-thermal emission driven by shocks produced in structure formation. These two clusters are detected by Suzaku/XIS for the first time in the X-ray band of 0.5-7 keV. No significant flux is detected by HXD in the energy band of 13-40 keV, and upper limits are reported. From the analysis of the XIS data, we find that the spectrum of A1555 is fit by a thermal plus power-law model, significantly better than a single-temperature pure thermal spectrum. If this power-law component is due to inverse-Compton scattering, the fraction of total baryon energy imparted to non-thermal electrons is consistent with the typical value inferred from the observation of other clusters. However, other scenarios (e.g., under lying AGNs, multi-temperature thermal models) cannot be excluded and further investigation of this system is desired. Basic physical properties of A1555 (e.g., total mass) are also reported.

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.

We report on our second-year campaign of X-ray follow-up observations of unidentified Fermi-LAT \gamma-ray sources at high Galactic latitudes (|b|>10 degree) using the X-ray Imaging Spectrometer onboard the Suzaku X-ray Observatory. In this second year of the project, seven new targets were selected from the First Fermi-LAT Catalog, and studied with 20-40 ks effective Suzaku exposures. We detected an X-ray point source coincident with the position of the recently discovered millisecond pulsar PSR J2302+4442 within the 95% confidence error circle of 1FGL J2302.8+4443. The X-ray spectrum of the detected counterpart was well fit by a blackbody model with temperature of kT ~0.3 keV, consistent with an origin of the observed X-ray photons from the surface of a rotating magnetized neutron star. For four other targets which were also recently identified with a normal pulsar (1FGL J0106.7+4853) and millisecond pulsars (1FGL J1312.6+0048, J1902.0-5110, and J2043.2+1709), only upper limits in the 0.5-10 keV band were obtained at the flux levels of ~10^{-14} erg cm^{-2} s^{-1}. A weak X-ray source was found in the field of 1FGL J1739.4+8717, but its association with the variable \gamma-ray emitter could not be confirmed with the available Suzaku data alone. For the remaining Fermi-LAT object 1FGL J1743.8-7620 no X-ray source was detected within the LAT 95% error ellipse. We briefly discuss the general properties of the observed high Galactic-latitude Fermi-LAT objects by comparing their multiwavelength properties with those of known blazars and millisecond pulsars.

We report on the results of deep X-ray follow-up observations of four unidentified Fermi/LAT gamma-ray sources at high Galactic latitudes using Suzaku. The studied objects were detected with high significance during the first 3 months of Fermi/LAT operation, and subsequently better localized in the Fermi/LAT 1 year catalog (1FGL). Possible associations with pulsars and active galaxies have subsequently been discussed, and our observations provide an important contribution to this debate. In particular, an X-ray point source was found within the 95% confidence error circle of 1FGL J1231.1-1410. X-ray spectrum is well-fitted by a blackbody with an additional power-law. This supports the recently claimed identification of this source with a millisecond pulsar (MSP) PSR J1231-1411. Concerning 1FGL J1311.7-3429, two X-ray sources were found within the LAT error circle. Even though the X-ray spectral and variability properties were accessed, their nature and relationship with the gamma-ray source remain uncertain. We found several weak X-ray sources in the field of 1FGL J1333.2+5056, one coinciding with CLASS J1333+5057. We argue the available data are consistent with the association between these two objects. Finally, we have detected an X-ray source in the vicinity of 1FGL J2017.3+0603. This object was recently suggested to be associated with a newly discovered MSP PSR J2017+0603, because of the spatial-coincidence and the gamma-ray pulse detection. We have only detected the X-ray counterpart of the CLASS J2017+0603, while we determined an X-ray flux upper limit at the pulsar position. All in all, our studies indicate while a significant fraction of unidentified high Galactic latitude gamma-ray sources is related to the pulsar and blazar phenomena, associations with other classes of astrophysical objects are still valid options.

We present a new theoretical calculation of the contribution to the extragalactic gamma-ray background radiation (EGRB) from star-forming galaxies, based on a state-of-the-art model of hierarchical galaxy formation that is in quantitative agreement with a variety of observations of local and high-redshift galaxies. Gamma-ray luminosity ($L_\gamma$) and spectrum of galaxies are related to star formation rate ($\psi$), gas mass ($M_gas$), and star formation mode (quiescent or starburst) of model galaxies using latest observed data of nearby galaxies. We try the two limiting cases about gamma-ray production: the escape limit ($L_\gamma \propto \psi M_gas$) and the calorimetric limit ($L_\gamma \propto \psi$), and our standard model predicts 7 and 4% contribution from star-forming galaxies to the total EGRB flux (including bright resolved sources) recently reported by the Fermi Gamma-Ray Space Telescope. Systematic uncertainties do not allow us to determine the EGRB flux better than by a factor of ~2. The predicted number of nearby galaxies detectable by Fermi is consistent with the observation. Intergalactic absorption by pair-production attenuates the EGRB flux only by a modest factor of ~1.3 at the highest Fermi energy band, and the reprocessed cascade emission does not significantly alter EGRB at lower photon energies. The sum of the known contributions from AGNs and star-forming galaxies can explain a large part of EGRB, with a remarkable agreement between the predicted model spectrum and observation.