The most widely used Galactic extinction map (Schlegel, Finkbeiner, & Davis 1998, SFD) is constructed assuming that the observed FIR fluxes entirely come from the Galactic dust. According to the earlier suggestion by Yahata et al. (2007), we consider how far-infrared (FIR) emission of galaxies affects the SFD map. We first compute the surface number density of SDSS DR7 galaxies as a function of the $r$-band extinction, $A_{r,\rm SFD}$. We confirm that the surface densities of those galaxies positively correlate with $A_{r,\rm SFD}$ for $A_{r,\rm SFD}<0.1$, as first discovered by Yahata et al. (2007) for SDSS DR4 galaxies. Next we construct an analytic model to compute the surface density of galaxies taking account of the contamination of their FIR emission. We adopt a log-normal probability distribution for the ratio of $100\mu {\rm m}$ and $r$-band luminosities of each galaxy, $y \equiv (\nu L)_{100\mu {\rm m}}/(\nu L)_r$. Then we search for the mean and r.m.s values of $y$ that fit the observed anomaly using the analytic model. The required values to reproduce the anomaly are roughly consistent with those measured from the stacking analysis of SDSS galaxies (Kashiwagi, Yahata, & Suto 2013). Due to the limitation of our statistical modeling, we are not yet able to remove the FIR contamination of galaxies from the extinction map. Nevertheless the agreement with the model prediction suggests that the FIR emission of galaxies is mainly responsible for the observed anomaly. While the corresponding systematic error in the Galactic extinction map is 0.1 to 1mmag, it is directly correlated with galaxy clustering, and thus needs to be carefully examined in precision cosmology.

We have performed stacking image analyses of galaxies over the Galactic extinction map constructed by Schlegel, Finkbeiner & Davis (1998). We select ~10^7 galaxies in total from the Sloan Digital Sky Survey (SDSS) DR7 photometric catalog. We detect clear signatures of the enhancement of the extinction in r-band, $\Delta A_r$, around galaxies, indicating that the extinction map is contaminated by their FIR (far infrared) emission. The average amplitude of the contamination per galaxy is well fitted to $\Delta A_r(m_r) = 0.64 \times 10^{0.17(18-m_r)}$ [mmag]. While this value is very small, it is directly associated with galaxies and may have a systematic effect on galaxy statistics. Indeed this correlated contamination leads to a relatively large anomaly of galaxy surface number densities against the SFD extinction A_SFD discovered by Yahata et al. (2007). We model the radial profiles of stacked galaxy images, and find that the FIR signal around each galaxy does not originate from the central galaxy alone, but is dominated by the contributions of nearby galaxies via galaxy angular clustering. The separation of the single galaxy and the clustering terms enables us to infer the statistical relation of the FIR and r-band fluxes of galaxies and also to probe the flux-weighted cross-correlation of galaxies, down to the magnitudes that are difficult to probe directly for individual objects. We repeat the same stacking analysis for SDSS DR6 photometric quasars and discovered the similar signatures but with weaker amplitudes. The implications of the present results for galaxy and quasar statistics and for correction to the Galactic extinction map are briefly discussed.

The peculiar motion of the Earth causes a dipole anisotropy modulation in the distant galaxy distribution due to the aberration effect. However, the amplitude and angular direction of the effect is not necessarily the same as those of the cosmic microwave background (CMB) dipole anisotropy due to the growth of cosmic structures. In other words exploring the aberration effect may give us a clue to the horizon-scale physics perhaps related to the cosmic acceleration. In this paper we develop a method to explore the dipole angular modulation from the pixelized galaxy data on the sky properly taking into account the covariances due to the shot noise and the intrinsic galaxy clustering contamination as well as the partial sky coverage. We applied the method to the galaxy catalogs constructed from the Sloan Digital Sky Survey (SDSS) Data Release 6 data. After constructing the four galaxy catalogs that are different in the ranges of magnitudes and photometric redshifts to study possible systematics, we found that the most robust sample against systematics indicates no dipole anisotropy in the galaxy distribution. This finding is consistent with the expectation from the concordance Lambda-dominated cold dark matter model. Finally we argue that an almost full-sky galaxy survey such as LSST may allow for a significant detection of the aberration effect of the CMB dipole having the precision of constraining the angular direction to ~ 20 degrees in radius. Assuming a hypothetical LSST galaxy survey, we find that this method can confirm or reject the result implied from a stacked analysis of the kinetic Sunyaev-Zel'dovich effect of X-ray luminous clusters in Kashlinsky et al. (2008,2009) if the implied cosmic bulk flow is not extended out to the horizon.

We investigate correlation between the arrival directions of ultra-high-energy cosmic rays (UHECRs) and the large-scale structure (LSS) of the Universe by using statistical quantities which can find the angular scale of the correlation. The Infrared Astronomical Satellite Point Source Redshift Survey (IRAS PSCz) catalog of galaxies is adopted for LSS. We find a positive correlation of the highest energy events detected by the Pierre Auger Observatory (PAO) with the IRAS galaxies inside $z=0.018$ within the angular scale of $\sim 15^{\circ}$. This positive correlation observed in the southern sky implies that a significant fraction of the highest energy events comes from nearby extragalactic objects. We also analyze the data of the Akeno Giant Air Shower Array (AGASA) which observed the northern hemisphere, but the obvious signals of positive correlation with the galaxy distribution are not found. Since the exposure of the AGASA is smaller than the PAO, the cross-correlation in the northern sky should be tested using a larger number of events detected in the future. We also discuss the correlation using the all-sky combined data sets of both the PAO and AGASA, and find a significant correlation within $\sim 8^{\circ}$. These angular scales can constrain several models of intergalactic magnetic field. These cross-correlation signals can be well reproduced by a source model in which the distribution of UHECR sources is related to the IRAS galaxies.

We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a \Lambda Cold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

We critically examine how well the evolution of large-scale density perturbations is followed in cosmological $N$-body simulations. We first run a large volume simulation and perform a mode-by-mode analysis in three-dimensional Fourier space. We show that the growth of large-scale fluctuations significantly deviates from linear theory predictions. The deviations are caused by {\it nonlinear} coupling with a small number of modes at largest scales owing to finiteness of the simulation volume. We then develop an analytic model based on second-order perturbation theory to quantify the effect. Our model accurately reproduces the simulation results. For a single realization, the second-order effect appears typically as ``zig-zag'' patterns around the linear-theory prediction, which imprints artificial ``oscillations'' that lie on the real baryon-acoustic oscillations. Although an ensemble average of a number of realizations approaches the linear theory prediction, the dispersions of the realizations remain large even for a large simulation volume of several hundred megaparsecs on a side. For the standard $\Lambda$CDM model, the deviations from linear growth rate are as large as 10 percent for a simulation volume with $L = 500h^{-1}$Mpc and for a bin width in wavenumber of $\Delta k = 0.005h$Mpc$^{-1}$, which are comparable to the intrinsic variance of Gaussian random realizations. We find that the dispersions scales as $\propto L^{-3/2} \Delta k^{-1/2}$ and that the mean dispersion amplitude can be made smaller than a percent only if we use a very large volume of $L > 2h^{-1}$Gpc. The finite box size effect needs to be appropriately taken into account when interpreting results from large-scale structure simulations for future dark energy surveys using baryon acoustic oscillations.

An acoustic oscillation of the primeval photon-baryon fluid around the decoupling time imprints a characteristic scale in the galaxy distribution today, known as the baryon acoustic oscillation (BAO) scale. Several on-going and/or future galaxy surveys aim at detecting and precisely determining the BAO scale so as to trace the expansion history of the universe. We consider nonlinear and redshift-space distortion effects on the shifts of the BAO scale in $k$-space using perturbation theory. The resulting shifts are indeed sensitive to different choices of the definition of the BAO scale, which needs to be kept in mind in the data analysis. We present a toy model to explain the physical behavior of the shifts. We find that the BAO scale defined as in Percival et al. (2007) indeed shows very small shifts ($\lesssim$ 1%) relative to the prediction in {\it linear theory} in real space. The shifts can be predicted accurately for scales where the perturbation theory is reliable.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We discuss how the baryon acoustic oscillation (BAO) signatures in the galaxy power spectrum can distinguish between modified gravity and the cosmological constant as the source of cosmic acceleration. To this end we consider a model characterized by a parameter $n$, which corresponds to the Dvali-Gabadadze-Porrati (DGP) model if $n=2$ and reduces to the standard spatially flat cosmological constant concordance model for $n$ equal to infinity. We find that the different expansion histories of the modified gravity models systematically shifts the peak positions of BAO. A preliminary analysis using the current SDSS LRG sample indicates that the original DGP model is disfavored unless the matter density parameter exceeds 0.3. The constraints will be strongly tightened with future spectroscopic samples of galaxies at high redshifts. We demonstrate that WFMOS, in collaboration with other surveys such as Planck, will powerfully constrain modified gravity alternatives to dark energy as the explanation of cosmic acceleration.

We compare the most successful and widely used map of Galactic dust extinction, provided by Schlegel, Finkbeiner & Davis (1998; hereafter SFD), to the galaxy number counts in the Sloan Digital Sky Survey (SDSS) photometric/spectroscopic DR4 sample. We divide the SDSS survey area into 69 disjoint subregions according to the dust extinction provided by SFD and compare the surface number density of galaxies in each subregion. As expected, the galaxy surface number density decreases with increasing extinction but only for SFD extinction values above about 0.1 to 0.2 magnitudes (depending on the band). At lower values of the SFD extinction, we find that the sky surface density of galaxies increases with increasing extinction, precisely the opposite of the effect expected from Galactic dust. We suggest that the far infrared (FIR) brightness of the sky in regions of true low dust extinction is significantly ``contaminated'' by the FIR emission from background galaxies. We show that such an explanation is both qualitatively and quantitatively consistent with the available data. Based on this interpretation we conclude that systematic errors in the SFD extinction map due to extragalactic FIR emission are quite small, of order hundredths of a magnitude, but nevertheless statistically detectable. (Abridged)

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)