The long-wavelength coherent overdensity and tidal force, which are not direct observables for a finite-volume survey, affect time evolution of cosmic structure formation and therefore clustering observables through the mode coupling. In this paper we develop an "anisotropic" separate universe (SU) simulation technique to simulate large-scale structure formation taking into account the effect of large-scale tidal force into the anisotropic expansion of local background. We modify the TreePM N-body simulation code to implement the anisotropic SU simulations, and then study the "response" function of matter power spectrum that describes how the matter power spectrum responds to the large-scale tidal effect as a function of wavenumber and redshift for a given global cosmology. We test and validate the SU simulation results from the comparison with the perturbation theory predictions and the results from high-resolution PM simulations. We find that the response function displays characteristic scale dependences over the range of scales down to nonlinear scales, up to k ~ 6 h/Mpc.

We constrain the abundance of primordial black holes (PBH) using 2622 microlensing events obtained from 5-years observations of stars in the Galactic bulge by the Optical Gravitational Lensing Experiment (OGLE). The majority of microlensing events display a single or at least continuous population that has a peak around the light curve timescale $t_{\rm E}\simeq 20~{\rm days}$ and a wide distribution over the range $t_{\rm E}\simeq [1, 300]~{\rm days}$, while the data also indicates a second population of 6 ultrashort-timescale events in $t_{\rm E}\simeq [0.1,0.3]~{\rm days}$, which are advocated to be due to free-floating planets. We confirm that the main population of OGLE events can be well modeled by microlensing due to brown dwarfs, main sequence stars and stellar remnants (white dwarfs and neutron stars) in the standard Galactic bulge and disk models for their spatial and velocity distributions. Using the dark matter (DM) model for the Milky Way (MW) halo relative to the Galactic bulge/disk models, we obtain the tightest upper bound on the PBH abundance in the mass range $M_{\rm PBH}\simeq[10^{-6},10^{-3}]M_\odot$ (Earth-Jupiter mass range), if we employ null hypothesis that the OGLE data does not contain any PBH microlensing event. More interestingly, we also show that Earth-mass PBHs can well reproduce the 6 ultrashort-timescale events, without the need of free-floating planets, if the mass fraction of PBH to DM is at a per cent level, which is consistent with other constraints such as the microlensing search for Andromeda galaxy (M31) and the longer timescale OGLE events. Our result gives a hint of PBH existence, and can be confirmed or falsified by microlensing search for stars in M31, because M31 is towards the MW halo direction and should therefore contain a much less number of free-floating planets, even if exist, than the direction to the MW center.

The peculiar velocity field measured by redshift-space distortions (RSD) in galaxy surveys provides a unique probe of the growth of large-scale structure. However, systematic effects arise when including satellite galaxies in the clustering analysis. Since satellite galaxies tend to reside in massive halos with a greater halo bias, the inclusion boosts the clustering power. In addition, virial motions of the satellite galaxies cause a significant suppression of the clustering power due to nonlinear RSD effects. We develop a novel method to recover the redshift-space power spectrum of halos from the observed galaxy distribution by minimizing the contamination of satellite galaxies. The cylinder grouping method (CGM) we study effectively excludes satellite galaxies from a galaxy sample. However, we find that this technique produces apparent anisotropies in the reconstructed halo distribution over all the scales which mimic RSD. On small scales, the apparent anisotropic clustering is caused by exclusion of halos within the anisotropic cylinder used by the CGM. On large scales, the misidentification of different halos in the large-scale structures, aligned along the line-of-sight, into the same CGM group, causes the apparent anisotropic clustering via their cross-correlation with the CGM halos. We construct an empirical model for the CGM halo power spectrum, which includes correction terms derived using the CGM window function at small scales as well as the linear matter power spectrum multiplied by a simple anisotropic function at large scales. We apply this model to a mock galaxy catalog at z=0.5, designed to resemble SDSS-III BOSS CMASS galaxies, and find that our model can predict both the monopole and quadrupole power spectra of the host halos up to k<0.5 h/Mpc to within 5%.

We present a first application of the subhalo abundance matching (SHAM) method to describe the redshift-space clustering of galaxies including the non-linear redshift-space distortion, i.e., the Fingers-of-God. We find that the standard SHAM connecting the luminosity of galaxies to the maximum circular velocity of subhalos well reproduces the luminosity dependence of redshift-space clustering of galaxies in the Sloan Digital Sky Survey in a wide range of scales from 0.3 to 40 Mpc/h. The result indicates that the SHAM approach is very promising for establishing a theoretical model of redshift-space galaxy clustering without additional parameters. We also test color abundance matching using two different proxies for colors: subhalo age and local dark matter density following the method by Masaki et al. (2013b). Observed clustering of red galaxies exhibits much stronger Fingers-of-God effect than blue galaxies. We find that the subhalo age model describes the observed color-dependent redshift-space clustering much better than the local dark matter density model. The result infers that the age of subhalos is a key ingredient to determine the color of galaxies.

We extend the subhalo abundance matching method to assign galaxy color to subhalos. We separate a luminosity-binned subhalo sample into two groups by a secondary subhalo property which is presumed to be correlated with galaxy color. The two subsamples then represent red and blue galaxy populations. We explore two models for the secondary property, namely subhalo assembly time and local dark matter density around each subhalo. The model predictions for the galaxy two-point correlation functions are compared with the recent results from the Sloan Digital Sky Survey. We show that the observed color dependence of galaxy clustering can be reproduced well by our method applied to cosmological N-body simulations without baryonic processes. We then compare the model predictions for the color-dependent galaxy-mass cross correlation functions with the results from gravitational lensing observations. The comparison allows us to distinguish the models, and also to discuss what subhalo property should be used to assign color to subhalos accurately. We show that the extended abundance matching method using the local dark matter density as a color proxy provides an accurate description of the galaxy populations in the local universe. We also study impacts of scatter in the local dark matter density-color relations. Introducing scatter improves agreements of our model predictions with the observed red and blue galaxy clustering and is needed to explain observed correlation functions in finer color bins. Finally, we study red galaxy fraction profiles in galaxy group- and cluster-sized halos and find the red fraction profiles have a relatively strong dependence on our model parameters. We argue that the red fraction profiles can be an important observational clue, in addition to galaxy clustering and galaxy-galaxy lensing, to explore the galaxy-(sub)halo connections.

We develop a novel abundance matching method to construct a mock catalog of luminous red galaxies (LRGs) in SDSS, using catalogs of halos and subhalos in N-body simulations for a LCDM model. Motivated by observations suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age >5Gyr, we assume that simulated halos at z=2 (z2-halo) are progenitors for LRG-host subhalos observed today, and we label the most tightly bound particles in each progenitor z2-halo as LRG ``stars''. We then identify the subhalos containing these stars to z=0.3 (SDSS redshift) in descending order of the masses of z2-halos until the comoving number density of the matched subhalos becomes comparable to the measured number density of SDSS LRGs, n=10^{-4} (h/Mpc)^3. Once the above prescription is determined, our only free parameter is the number density of halos identified at z=2 and this parameter is fixed to match the observed number density at z = 0.3. By tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly compute, from the mock catalog, the distributions of central and satellite LRGs and their internal motions in each host halo at z=0.3. While the SDSS LRGs are galaxies selected by the magnitude and color cuts from the SDSS images and are not necessarily a stellar-mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo occupation distribution for central and satellite LRGs, the projected auto-correlation function of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, in the angle-averaged redshift-space power spectrum.

We develop an analytic halo model for the distribution of dust around galaxies. The model results are compared with the observed surface dust density profile measured through reddening of background quasars in the Sloan Digital Sky Survey (SDSS) reported by Menard et al.(2010). We assume that the dust distribution around a galaxy is described by a simple power law, similarly to the mass distribution, but with a sharp cut-off at $\alpha R_{\rm vir}$ where $R_{\rm vir}$ is the galaxy's virial radius and $\alpha$ is a model parameter. Our model reproduces the observed dust distribution profile very well over a wide range of radial distance of $10 - 10^{4} h^{-1}$kpc. For the characteristic galaxy halo mass of $2\times 10^{12} h^{-1}M_{\odot}$ estimated for the SDSS galaxies, the best fit model is obtained if $\alpha$ is greater than unity, which suggests that dust is distributed to over a few hundred kilo-parsecs from the galaxies. The observed large-scale dust distribution profile is reproduced if we assume the total amount of dust is equal to that estimated from the integrated stellar evolution over the cosmic time.

We explore the mass distribution of material associated with galaxies from the observation of gravitational weak lensing for the galaxy mass correlation function with the aid of $N$-body simulations of dark matter. The latter is employed to unfold various contributions that contribute to the integrated line of sight mass density. We conclude that galaxies have no definite edges of the matter distribution, extending to the middle to neighbouring galaxies with the density profile roughly $r^{-2.4}$ beyond the virial radius. The mass distributed beyond the virial radius (gravitationally bound radius) explains the gap seen in the mass density estimates, the global value $\Omega_m\sim 0.27$ and typically $\Omega_{\rm gal} \sim 0.15$ from the luminosity density multiplied by the mass to light ratio. We suggest to use a physical method of gravitational lensing to characterise galaxy samples rather than characterise them with photometric means.

A systematic study of gravitational waves from galaxy mergers, through N-body simulations, was performed. In particular, we investigated the relative importance of galaxy components (disk, bulge and halo) and effects of initial relative velocity, relative angular momentum and mass ratio of the galaxies. We found that the features of light curve of gravitational waves, such as peak width and luminosity, are reliably simulated with particle numbers larger than ~10^4. Dominant contribution to gravitational wave emission came from the halo component, while peak luminosity amounted to 10^31 erg/sec for the collision of two halos with mass 3.8 x10^12Msun/h. We also found that the initial relative velocity in the direction of the initial separation did not significantly affect gravitational wave emission, while the initial relative angular momentum broadened the peak width and suppressed the luminosity. Mass dependence of the peak luminosity was also investigated, and we obtained evidence that the luminosity is proportional to the cubic mass when the scaling relation is satisfied. This behavior was considered by a simple analysis.

We present a clustering analysis of ~60,000 massive (stellar mass Mstar > 10^{11} Msun) galaxies out to z = 1 drawn from 55.2 deg2 of the UKIRT Infrared Deep Sky Survey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II Supernova Survey. Strong clustering is detected for all the subsamples of massive galaxies characterized by different stellar masses (Mstar = 10^{11.0-11.5} Msun, 10^{11.5-12.0} Msun) or rest-frame colors (blue: U - V < 1.0, red: U - V > 1.0). We find that more mature (more massive or redder) galaxies are more clustered, which implies that more mature galaxies have started stellar-mass assembly earlier within the highly-biased region where the structure formation has also started earlier. By means of halo occupation distribution (HOD) models fitted to the observed angular correlation function, we infer the properties of the underlying host dark halos. We find that the estimated bias factors and host halo masses are systematically larger for galaxies with larger stellar masses, which is consistent with the general agreement that the capability of hosting massive galaxies depends strongly on halo mass. The estimated effective halo masses are ~10^{14} Msun, which gives the stellar-mass to halo-mass ratios of ~0.003. The observed evolution of bias factors indicates rapid evolution of spatial distributions of cold dark matter relative to those traced by the massive galaxies, while the transition of host halo masses might imply that the fractional mass growth rate of halos is less than those of stellar systems. The inferred halo masses and high fractions of central galaxies indicate that the massive galaxies in the current sample are possibly equivalent to central galaxies of galaxy clusters.