The observed accelerated expansion of the Universe may be explained by dark energy or the breakdown of general relativity (GR) on cosmological scales. When the latter case, a modified gravity scenario, is considered, it is often assumed that the background evolution is the same as the $\Lambda$CDM model but the density perturbation evolves differently. In this paper, we investigate more general classes of modified gravity, where both the background and perturbation evolutions are deviated from those in the $\Lambda$CDM model. We introduce two phase diagrams, $\alpha{\rm-}f\sigma _8$ and $H{\rm-}f\sigma _8$ diagrams; $H$ is the expansion rate, $f\sigma_8$ is a combination of the growth rate of the Universe and the normalization of the density fluctuation which is directly constrained by redshift-space distortions, and $\alpha$ is a parameter which characterizes the deviation of gravity from GR and can be probed by gravitational lensing. We consider several specific examples of Horndeski's theory, which is a general scalar-tensor theory, and demonstrate how deviations from the $\Lambda$CDM model appears in the $\alpha{\rm-}f\sigma _8$ and $H{\rm-}f\sigma _8$ diagrams. The predicted deviations will be useful for future large-scale structure observations to exclude some of the modified gravity models.

We perform a self-calibration of the richness-to-mass ($N$-$M$) relation of CAMIRA galaxy clusters with richness $N\geq15$ at redshift $0.2\leq z<1.1$ by modeling redshift-space two-point correlation functions, namely, the auto-correlation function $\xi_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $\xi_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $\xi_{\mathrm{cg}}$ between these two samples. We focus on constraining the normalization $A_{\mathrm{N}}$ of the $N$-$M$ relation in a forward-modeling approach, carefully accounting for the redshift-space distortion, the Finger-of-God effect, and the uncertainty in photometric redshifts of CAMIRA clusters. The modeling also takes into account the projection effect on the halo bias of CAMIRA clusters. The parameter constraints are shown to be unbiased according to validation tests using a large set of mock catalogs constructed from N-body simulations. At the pivotal mass $M_{500}=10^{14}h^{-1}M_{\odot}$ and the pivotal redshift $z_{\mathrm{piv}} = 0.6$, the resulting normalization $A_{\mathrm{N}}$ is constrained as $13.8^{+5.8}_{-4.2}$, $13.2^{+3.4}_{-2.7}$, and $11.9^{+3.0}_{-1.9}$ by modeling $\xi_{\mathrm{cc}}$, $\xi_{\mathrm{cc}}+\xi_{\mathrm{cg}}$, and $\xi_{\mathrm{cc}} + \xi_{\mathrm{cg}} + \xi_{\mathrm{gg}}$, with average uncertainties at levels of $36\%$, $23\%$, and $21\%$, respectively. We find that the resulting $A_{\mathrm{N}}$ is statistically consistent with those independently obtained from weak-lensing magnification and from a joint analysis of shear and cluster abundance, with a preference for a lower value at a level of $\lesssim1.9\sigma$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged]

We study an efficient way to enhance the measurability of the galaxy density and/or velocity power spectrum in redshift space. It is based on the angular decomposition with the tripolar spherical harmonic (TripoSH) basis and applicable even to galaxy distributions in wide-angle galaxy surveys. While nontrivial multipole-mode mixings are inevitable in the covariance of the Legendre decomposition coefficient commonly used in the small-angle power spectrum analysis, our analytic computation of the covariance of the TripoSH decomposition coefficient shows that such mixings are absent by virtue of high separability of the TripoSH basis, yielding the minimum variance. Via the simple signal-to-noise ratio assessment, we confirm that the detectability improvement by the TripoSH decomposition approach becomes more significant at higher multipole modes, and, e.g., the hexadecapole of the density power spectrum has two orders of magnitude improvement. The TripoSH decomposition approach is expected to be applied to not only currently available survey data but also forthcoming wide-angle one, and to bring about something new or much more accurate cosmological information.

We show that the spatial correlation of the intrinsic alignments (IAs) of galaxies, measured in galaxy redshift surveys, offers a precision route to improve the geometrical and dynamical constraints on cosmology. The IA has been treated as a contaminant against cosmological probes such as weak gravitational lensing experiments. However, the large-scale correlation of IAs is expected to follow the coherent large-scale matter inhomogeneities. Here, making use of its anisotropic nature, we show that the large-scale IA correlations help to improve the measurements of the geometric distances and growth of structure. In combination with the conventional galaxy clustering statistics, we find that constraints on equation-of-state parameter for dark energy and Hubble parameter can be tighter than those from the clustering statistics alone by a factor of more than 1.5.

Antenna efficiency is one of the most important figures-of-merit of a radio telescope, even for a multibeam radio telescope. To analyze a system where each beam illuminates only a part of the aperture, a lossless antenna consisting of two apertures in series is considered in the frame of the scalar wave approximation. We found that its antenna efficiency is factorized into three factors: efficiencies of beam coupling, transmission spillover, and reception spillover. We numerically confirmed our factorization with an optical simulation. The factorization is applicable to general aperture-type antennas and is useful in design of multibeam radio telescopes.

Investigating intrinsic alignments (IAs) of galaxy shapes is important not only to constrain cosmological parameters unbiasedly from gravitational lensing but also to extract cosmological information complimentary to galaxy clustering analysis. We derive simple and useful formulas for the various IA statistics, including the intrinsic ellipticity--ellipticity correlation, the gravitational shear--intrinsic ellipticity correlation, and the velocity-intrinsic ellipticity correlation functions. The angular dependence of each statistic is explicitly given, namely the angle between the line-of-sight direction and the separation vector of two points. It thus allows us to analyze anisotropies of baryon acoustic oscillations encoded in the IA statistics, and we can extract the maximum cosmological information using the Alcock-Pacyznski and redshift-space distortion effects. We also provide these formulas for the intrinsic ellipticities decomposed into E and B modes.

There is a growing interest of using the intrinsic alignment (IA) of galaxy images as a tool to extract cosmological information complimentary to galaxy clustering analysis. Recently, Okumura & Taruya derived useful formulas for the intrinsic ellipticity--ellipticity correlation, the gravitational shear--intrinsic ellipticity correlation, and the velocity--intrinsic ellipticity correlation functions based on the linear alignment (LA) model. In this paper, using large-volume $N$-body simulations, we measure these alignment statistics in real and redshift space and compare them to the LA and nonlinear alignment model predictions. We find that anisotropic features of baryon acoustic oscillations in the IA statistics can be accurately predicted by our models. The anisotropy due to redshift-space distortions (RSDs) is also well described in the large-scale limit. Our results indicate that one can extract the cosmological information encoded in the IA through the Alcock-Paczynski and RSD effects.

The kinematic Sunyaev-Zel'dovich effect enables us to directly probe the density-weighted velocity field up to very large cosmic scales. We investigate the effects of intrinsic alignments (IAs) of dark-matter halo shapes on cosmic density and velocity fields on such large scales. In the literature IAs have been detected up to $\sim 100~h^{-1}{\rm Mpc}$ using the gravitational shear-intrinsic ellipticity correlation and the alignment correlation function. In this paper we introduce the corresponding various velocity statistics: the (density-weighted) velocity-intrinsic ellipticity correlation as well as the alignment pairwise infall momentum, momentum correlation function, and density-weighted pairwise velocity dispersion. We derive theoretical expressions for these velocity alignment statistics for the first time based on the assumption that the density fluctuation is a Gaussian random field. Using large-volume, high-resolution $N$-body simulations, we measure the alignment statistics of density and velocity fields. The behaviors of IAs in the velocity statistics are similar to those in the density statistics, except that the halo orientations are aligned with the velocity field up to scales larger than those with the density field, $r\gg 100~h^{-1}{\rm Mpc}$, because of a factor of the wave number in the linear relation between the density and velocity fields in Fourier space, $v\propto \delta/k$. We show that the detected IAs of the velocity field can be well predicted by the linear alignment model. We also demonstrate that the baryon acoustic oscillation features can be detected in both the conventional and alignment velocity statistics. Our results indicate that observations of IAs with the velocity field on large scales can provide additional information on cosmological models, complementary to those with the density field.

We show that peculiar velocities of Type Ia supernovae can be used to derive constraints on the sum of neutrino masses, $\Sigma m_{\nu}$, and dark energy equation of state, $w = w_0+w_a(1-a)$, from measurements of the magnitude-redshift relation, complementary to galaxy redshift and weak lensing surveys. Light from a supernova propagates through a perturbed Universe so the luminosity distance is modified from its homogeneous prediction. This modification is proportional to the matter density fluctuation and its time derivative due to gravitational lensing and peculiar velocity respectively. At low redshifts, the peculiar velocity signal dominates while at high redshifts lensing does. We show that using lensing and peculiar velocity of supernovae from the upcoming surveys WFIRST and ZTF, without other observations, we can constrain $\Sigma m_{\nu} \lesssim 0.31$ eV, $\sigma(w_0) \lesssim 0.02$, and ${\sigma(w_a)} \lesssim 0.18$ ($1-\sigma$ CL) in the $\Sigma m_{\nu}$-$w_0$-$w_a$ parameter space, where all the other cosmological parameters are fixed. We find that adding peculiar velocity information from low redshifts shrinks the volume of the parameter ellipsoid in this space by $\sim 33$%. We also allow $\Omega_{\text{CDM}}$ to vary as well as $\Sigma m_{\nu}$, $w_0$ and $w_a$, and demonstrate how these constraints degrade as a consequence.

We review three distance measurement techniques beyond the local universe: (1) gravitational lens time delays, (2) baryon acoustic oscillation (BAO), and (3) HI intensity mapping. We describe the principles and theory behind each method, the ingredients needed for measuring such distances, the current observational results, and future prospects. Time delays from strongly lensed quasars currently provide constraints on $H_0$ with < 4% uncertainty, and with 1% within reach from ongoing surveys and efforts. Recent exciting discoveries of strongly lensed supernovae hold great promise for time-delay cosmography. BAO features have been detected in redshift surveys up to z <~ 0.8 with galaxies and z ~ 2 with Ly-$\alpha$ forest, providing precise distance measurements and $H_0$ with < 2% uncertainty in flat $\Lambda$CDM. Future BAO surveys will probe the distance scale with percent-level precision. HI intensity mapping has great potential to map BAO distances at z ~ 0.8 and beyond with precisions of a few percent. The next years ahead will be exciting as various cosmological probes reach 1% uncertainty in determining $H_0$, to assess the current tension in $H_0$ measurements that could indicate new physics.

We investigate the splashback features of dark-matter halos based on cosmic density and velocity fields. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding velocity statistic, alignment momentum correlation function, to take into account the asphericity of halos. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity. On halo scales, $x\sim R_\mathrm{200m} \sim 1Mpc/h$, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. We also find that the splashback radius determined from the density correlation becomes $\sim 3.5\%$ smaller than that from the momentum correlation, with their correlation coefficient being 0.605. Moreover, the orientation-dependent splashback feature due to halo asphericity is measured when the density profile is determined by dark-matter particles, which can be used as a test of collisional CDM since the halo shape is predicted to be rounder in such a model.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

In this paper, we propose to use the mimetic Horndeski model as a model for the dark universe. Both cold dark matter (CDM) and dark energy (DE) phenomena are described by a single component, the mimetic field. In linear theory, we show that this component effectively behaves like a perfect fluid with zero sound speed and clusters on all scales. For the simpler mimetic cubic Horndeski model, if the background expansion history is chosen to be identical to a perfect fluid DE (PFDE) then the mimetic model predicts the same power spectrum of the Newtonian potential as the PFDE model with zero sound speed. In particular, if the background is chosen to be the same as that of LCDM, then also in this case the power spectrum of the Newtonian potential in the mimetic model becomes indistinguishable from the power spectrum in LCDM on linear scales. A different conclusion may be found in the case of non-adiabatic perturbations. We also discuss the distinguishability, using power spectrum measurements from LCDM N-body simulations as a proxy for future observations, between these mimetic models and other popular models of DE. For instance, we find that if the background has an equation of state equal to -0.95 then we will be able to distinguish the mimetic model from the PFDE model with unity sound speed. On the other hand, it will be hard to do this distinction with respect to the LCDM model.

Intrinsic alignments (IA), the coherent alignment of intrinsic galaxy orientations, can be a source of a systematic error of weak lensing surveys. The redshift evolution of IA also contains information about the physics of galaxy formation and evolution. This paper presents the first measurement of IA at high redshift, $z\sim 1.4$, using the spectroscopic catalog of blue star-forming galaxies of the FastSound redshift survey, with the galaxy shape information from the Canada-Hawaii-France telescope lensing survey. The IA signal is consistent with zero with power-law amplitudes fitted to the projected correlation functions for density-shape and shape-shape correlation components, $A_{\delta+}=-0.0071\pm 0.1340$ and $A_{++}=-0.0505\pm 0.0848$, respectively. These results are consistent with those obtained from blue galaxies at lower redshifts (e.g., $A_{\delta+}=0.0035_{-0.0389}^{+0.0387}$ and $A_{++}=0.0045_{-0.0168}^{+0.0166}$ at $z=0.51$ from the WiggleZ survey). The upper limit of the constrained IA amplitude corresponds to a few percent contamination to the weak-lensing shear power spectrum, resulting in systematic uncertainties on the cosmological parameter estimations by $-0.052<\Delta \sigma_8<0.039$ and $-0.039<\Delta \Omega_m<0.030$.

We study the dependence of surface mass density profiles, which can be directly measured by weak gravitational lensing, on the orientation of haloes with respect to the line-of-sight direction, using a suite of $N$-body simulations. We find that, when major axes of haloes are aligned with the line-of-sight direction, surface mass density profiles have higher amplitudes than those averaged over all halo orientations, over all scales from $0.1$ to $100\,\mathrm{Mpc}/h$ we studied. While the orientation dependence at small scales is ascribed to the halo triaxiality, our results indicate even stronger orientation dependence in the so-called two-halo regime, up to $100\,\mathrm{Mpc}/h$. The orientation dependence for the two-halo term is well approximated by a multiplicative shift of the amplitude and therefore a shift in the halo bias parameter value. The halo bias from the two-halo term can be overestimated or underestimated by up to $\sim 30 \%$ depending on the viewing angle, which translates into the bias in estimated halo masses by up to a factor of two from halo bias measurements. The orientation dependence at large scales originates from the anisotropic halo-matter correlation function, which has an elliptical shape with the axis ratio of $\sim 0.55$ up to $100\, \mathrm{Mpc}/h$. We discuss potential impacts of halo orientation bias on other observables such as optically selected cluster samples and a clustering analysis of large-scale structure tracers such as quasars.

We present the measurement of the kinematic Sunyaev-Zel'dovich (kSZ) effect in Fourier space, rather than in real space. We measure the density-weighted pairwise kSZ power spectrum, the first use of this promising approach, by cross-correlating a cleaned Cosmic Microwave Background (CMB) temperature map, which jointly uses both Planck Release 2 and Wilkinson Microwave Anisotropy Probe nine-year data, with the two galaxy samples, CMASS and LOWZ, derived fr om the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12. With the current data, we constrain the average optical depth $\tau$ multiplied by the ratio of the Hubble parameter at redshift $z$ and the present day, $E=H/H_0$; we find $\tau E = (3.95\pm1.62)\times10^{-5}$ for LOWZ and $\tau E = ( 1.25\pm 1.06)\times10^{-5}$ for CMASS, with the optimal angular radius of an aperture photometry filter to estimate the CMB temperature distortion associ ated with each galaxy. By repeating the pairwise kSZ power analysis for various aperture radii, we measure the optical depth as a function of aperture ra dii. While this analysis results in the kSZ signals with only evidence for a detection, ${\rm S/N}=2.54$ for LOWZ and $1.24$ for CMASS, the combination of future CMB and spectroscopic galaxy surveys should enable precision measurements. We estimate that the combination of CMB-S4 and data from DESI shoul d yield detections of the kSZ signal with ${\rm S/N}=70-100$, depending on the resolution of CMB-S4.

We measure statistically anisotropic signatures imprinted in three-dimensional galaxy clustering using bipolar spherical harmonics (BipoSHs) in both Fourier space and configuration space. We then constrain a well-known quadrupolar anisotropy parameter $g_{2M}$ in the primordial power spectrum, parametrized by $P(\vec{k}) = \bar{P}(k) [ 1 + \sum_{M} g_{2M} Y_{2M}(\hat{k}) ]$, with $M$ determining the direction of the anisotropy. Such an anisotropic signal is easily contaminated by artificial asymmetries due to specific survey geometry. We precisely estimate the contaminated signal and finally subtract it from the data. Using the galaxy samples obtained by the Baryon Oscillation Spectroscopic Survey Data Release 12, we find no evidence for violation of statistical isotropy, $g_{2M}$ for all $M$ to be of zero within the $2\sigma$ level. The $g_{2M}$-type anisotropy can originate from the primordial curvature power spectrum involving a directional-dependent modulation $g_* (\hat{k} \cdot \hat{p})^2$. The bound on $g_{2M}$ is translated into $g_*$ as $-0.09 < g_* < 0.08$ with a $95\%$ confidence level when $\hat{p}$ is marginalized over.

We investigate the effects of intrinsic alignments (IA) of dark-matter halo shapes on cosmic density and velocity fields from cluster to cosmic scales beyond 100 Mpc/h. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding two velocity statistics, the angle-binned pairwise infall momentum and momentum correlation function. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity, both in real and redshift space. We find that the alignment signal is not amplified by redshift-space distortions at linear scales. Behaviors of IA in the velocity statistics are similar to those in the density statistics, except that the halo orientations are aligned with the velocity field up to a scale larger than those with the density field, x>100 Mpc/h. On halo scales, x~ R_{200m} ~ 1 Mpc/h, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. Our results indicate that observations of IA with the velocity field can provide additional information on cosmological models from large scales and on physical sizes of halos from small scales.

Recently a new method to discover obscured active galactic nuclei (AGNs) by utilizing X-ray and Infrared data has been developed. We carried out a survey of H2O maser emission toward ten obscured AGNs with the Nobeyama 45-m telescope. We newly detected the maser emission with the signal-noise-ratio (SNR) of above 4 from two AGNs, NGC 1402 and NGC 7738. We also found a tentative detection with its SNR > 3 in NGC 5037. The detection rate of 20% is higher than those of previous surveys (usually several percents).

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 propose an efficient way to test rotational invariance in the cosmological perturbations by use of galaxy correlation functions. In symmetry-breaking cases, the galaxy power spectrum can have extra angular dependence in addition to the usual one due to the redshift-space distortion, $ \hat{k} \cdot \hat{n}$. We confirm that, via the decomposition into not the usual Legendre basis ${\cal L}_\ell(\hat{k} \cdot \hat{n})$ but the bipolar spherical harmonic one $\{Y_{\ell}(\hat{k}) \otimes Y_{\ell'}(\hat{n})\}_{LM}$, the symmetry-breaking signal can be completely distinguished from the usual isotropic one since the former yields nonvanishing $L \geq 1$ modes but the latter is confined to the $L = 0$ one. As a demonstration, we analyze the signatures due to primordial-origin symmetry breakings such as the well-known quadrupolar-type and dipolar-type power asymmetries and find nonzero $L = 2$ and $1$ modes, respectively. Fisher matrix forecasts of their constraints indicate that the $Planck$-level sensitivity could be achieved by the SDSS or BOSS-CMASS data, and an order-of-magnitude improvement is expected in a near future survey as PFS or Euclid by virtue of an increase in accessible Fourier mode. Our methodology is model-independent and hence applicable to the searches for various types of statistically anisotropic fluctuations.

Yes. Future CMB experiments such as Advanced ACTPol and CMB-S4 should achieve measurements with S/N of $> 0.1$ for the typical galaxies in redshift surveys. These measurements will provide complementary measurements of the growth rate of large scale structure $f$ and the expansion rate of the Universe $H$ to galaxy clustering measurements. This paper emphasizes that there is significant information in the anisotropy of the relative pairwise kSZ measurements. We expand the relative pairwise kSZ power spectrum in Legendre polynomials and consider up to its octopole. Assuming that the noise in the filtered maps is uncorrelated between the positions of galaxies in the survey, we derive a simple analytic form for the power spectrum covariance of the relative pairwise kSZ temperature in redshift space. While many previous studies have assumed optimistically that the optical depth of the galaxies $\tau_{\rm T}$ in the survey is known, we marginalize over $\tau_{\rm T}$, to compute constraints on the growth rate $f$ and the expansion rate $H$. For realistic sure parameters, we find that combining kSZ and galaxy redshift survey data reduces the marginalized $1$-$\sigma$ errors on $H$ and $f$ by $\sim50$-$70\%$ compared to the galaxy-only analysis.

We present basic properties of $\sim$3,300 emission line galaxies detected by the FastSound survey, which are mostly H$\alpha$ emitters at $z \sim$ 1.2-1.5 in the total area of about 20 deg$^2$, with the H$\alpha$ flux sensitivity limit of $\sim 1.6 \times 10^{-16} \rm erg \ cm^{-2} s^{-1}$ at 4.5 sigma. This paper presents the catalogs of the FastSound emission lines and galaxies, which will be open to the public in the near future. We also present basic properties of typical FastSound H$\alpha$ emitters, which have H$\alpha$ luminosities of $10^{41.8}$-$10^{43.3}$ erg/s, SFRs of 20--500 $M_\odot$/yr, and stellar masses of $10^{10.0}$--$10^{11.3}$ $M_\odot$. The 3D distribution maps for the four fields of CFHTLS W1--4 are presented, clearly showing large scale clustering of galaxies at the scale of $\sim$ 100--600 comoving Mpc. Based on 1,105 galaxies with detections of multiple emission lines, we estimate that contamination of non-H$\alpha$ lines is about 4% in the single-line emission galaxies, which are mostly [OIII]$\lambda$5007. This contamination fraction is also confirmed by the stacked spectrum of all the FastSound spectra, in which H$\alpha$, [NII]$\lambda \lambda$6548,6583, [SII]$\lambda \lambda$6717, 6731, and [OI]$\lambda \lambda$6300,6364 are seen.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

Observations of the kinetic Sunyaev-Zel'dovich (kSZ) effect measure the density-weighted velocity field, a potentially powerful cosmological probe. This paper presents an analytical method to predict the power spectrum and two-point correlation function of the density-weighted velocity in redshift space, the direct observables in kSZ surveys. We show a simple relation between the density power spectrum and the density-weighted velocity power spectrum that holds for both dark matter and halos. Using this relation, we can then extend familiar perturbation expansion techniques to the kSZ power spectrum. One of the most important features of density-weighted velocity statistics in redshift space is the change in sign of the cross-correlation between the density and density-weighted velocity at mildly small scales due to nonlinear redshift space distortions. Our model can explain this characteristic feature without any free parameters. As a result, our results can precisely predict the non-linear behavior of the density-weighted velocity field in redshift space up to $\sim30\ h^{-1} {\rm Mpc}$ for dark matter particles at the redshifts of $z=0.0$, $0.5$, and $1.0$.

Theoretical modeling of the redshift-space power spectrum of galaxies is crucially important to correctly extract cosmological information from redshift surveys. The task is complicated by the nonlinear biasing and redshift space distortion effects, which change with halo mass, and by the wide distribution of halo masses and their occupations by galaxies. One of the main modeling challenges is the existence of satellite galaxies that have both radial distribution and large virial velocities inside halos, a phenomenon known as the Finger-of-God effect. We present a model for the galaxy power spectrum of in which we decompose a given galaxy sample into central and satellite galaxies and relate different contributions to 1- and 2-halo terms in a halo model. Our primary goal is to ensure that any parameters that we introduce have physically meaningful values, and are not just fitting parameters. For the 2-halo terms we use the previously developed RSD modeling of halos in the context of distribution function and perturbation theory approach. This term needs to be multiplied by the effect of radial distances and velocities of satellites inside the halo. To this one needs to add the 1-halo terms, which are non-perturbative. We show that the real space 1-halo terms can be modeled as almost constant, with the finite extent of the satellites inside the halo inducing a small k^2R^2 term, where R is related to the size of the halo. We adopt a similar model for FoG in redshift space, ensuring that FoG velocity dispersion is related to the halo mass. For FoG k^2 type expansions do not work and FoG resummation must be used instead. We test several damping functions to model the velocity dispersion FoG effect. Applying the formalism to mock galaxies modeled after the "CMASS" sample of the BOSS survey, we find that our predictions for the redshift-space power spectra are accurate up to k~0.4Mpc/h.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

Understanding the relation between underlying matter distribution and biased tracers such as galaxy or dark matter halo is essential to extract cosmological information from ongoing or future galaxy redshift surveys. At sufficiently large scales such as the BAO scale, a standard approach for the bias problem on the basis of the perturbation theory (PT) is to assume the `local bias' model in which the density field of biased tracers is deterministically expanded in terms of matter density field at the same position. The higher-order bias parameters are then determined by combining the power spectrum with higher-order statistics such as the bispectrum. As is pointed out by recent studies, however, nonlinear gravitational evolution naturally induces nonlocal bias terms even if initially starting only with purely local bias. As a matter of fact, previous works showed that the second-order nonlocal bias term, which corresponds to the gravitational tidal field, is important to explain the characteristic scale-dependence of the bispectrum. In this paper we extend the nonlocal bias term up to third order, and investigate whether the PT-based model including nonlocal bias terms can simultaneously explain the power spectrum and the bispectrum of simulated halos in $N$-body simulations. We show that the power spectrum, including density and momentum, and the bispectrum between halo and matter in $N$-body simulations can be simultaneously well explained by the model including up to third-order nonlocal bias terms up to k~0.1h/Mpc. Also, the results seem in a good agreement with theoretical predictions of a simple coevolution picture, although the agreement is not perfect. These demonstration clearly shows a failure of the local bias model even at such large scales, and we conclude that nonlocal bias terms should be consistently included in order to model statistics of halos. [abridged]

We analyze the clustering of large scale structure in the Universe in a model independent method, accounting for anisotropic effects along and transverse to the line of sight. The Baryon Oscillation Spectroscopy Survey Data Release 11 provides a large sample of 690,000 galaxies, allowing determination of the Hubble expansion H, angular distance D_A, and growth rate G_T at an effective redshift of z=0.57. After careful bias and convergence studies of the effects from small scale clustering, we find that cutting transverse separations below 40 Mpc/h delivers robust results while smaller scale data leads to a bias due to unmodelled nonlinear and velocity effects. The converged results are in agreement with concordance LCDM cosmology, general relativity, and minimal neutrino mass, all within the 68% confidence level. We also present results separately for the northern and southern hemisphere sky, finding a slight tension in the growth rate -- potentially a signature of anisotropic stress, or just covariance with small scale velocities -- but within 68% CL.

Direct measurements of peculiar velocities of galaxies and clusters of galaxies can in principle provide explicit information on the three dimensional mass distribution, but this information is modulated by the fact that velocity field is sampled at galaxy positions, and is thus probing galaxy momentum. We derive expressions for the cross power spectrum between the density and momentum field and the auto spectrum of the momentum field in redshift space, by extending the distribution function method to these statistics. The resulting momentum cross and auto power spectra in redshift space are expressed as infinite sums over velocity moment correlators in real space, as is the case for the density power spectrum in redshift space. We compute each correlator using Eulerian perturbation theory (PT) and halo biasing model and compare the resulting redshift-space velocity statistics to those measured from N-body simulations for both dark matter and halos. We find that in redshift space linear theory predictions for the density-momentum cross power spectrum as well as for the momentum auto spectrum fail to predict the N-body results at very large scales. On the other hand, our nonlinear PT prediction for these velocity statistics, together with real-space power spectrum for dark matter from simulations, improves the accuracy for both dark matter and halos. We also present the same analysis in configuration space, computing the redshift-space pairwise mean infall velocities and velocity correlation function and compare to nonlinear PT.

We examine the signature of dynamic (redshift-space) distortions and geometric distortions (including the Alcock-Paczynski effect) in the context of the galaxy power spectrum measured in upcoming galaxy redshift surveys. Information comes from both the baryon acoustic oscillation (BAO) feature and the broadband power spectrum shape. Accurate modeling is required to extract this information without introducing systematic biases in the result. We consider an analytic model for the power spectrum of dark matter halos in redshift space, based on the distribution function expansion, and compare with halo clustering measured in N-body simulations. We forecast that the distribution function model is sufficiently accurate to allow the inclusion of broadband information on scales down to k~0.2 h/Mpc, with somewhat better accuracy for higher bias halos. Compared with a BAO-only analysis with reconstruction, including broadband shape information can improve unbiased constraints on distance measures H(z) and D_A(z) by ~30% and 20%, respectively, for a galaxy sample similar to the DESI luminous red galaxies. The gains in precision are larger in the absence of BAO reconstruction. Furthermore, including broadband shape information allows the measurement of structure growth, through redshift-space distortions. For the same galaxy sample, the distribution function model is able to constrain f*sigma_8 to ~2%, when simultaneously fitting for H(z) and D_A(z). We discuss techniques to optimize the analysis of the power spectrum, including removing modes near the line-of-sight that are particularly challenging to model, and whether these approaches can improve parameter constraints. We find that such techniques are unlikely to significantly improve constraints on geometry, although they may allow higher precision measurements of redshift-space distortions.

Baryon acoustic oscillations (BAO), known as one of the largest cosmological objects, is now recognized as standard cosmological tool to measure geometric distances via the Alcock-Paczynski effect, by which the observed BAO exhibits characteristic anisotropies in addition to the redshift distortions. This implies that once we know the correct distances to the observed BAO, the tip points of baryon acoustic peaks in the anisotropic correlation function of galaxies, $\xi(\sigma,\pi)$, can form a great circle (hereafter 2D BAO circle) in the $\sigma$ and $\pi$ plane, where $\sigma$ and $\pi$ are the separation of galaxy pair parallel and perpendicular to the line-of-sight, respectively. This 2D BAO circle remains unchanged under the variations of the unknown galaxy bias and/or coherent motion, while it varies transversely and radially with respect to the variations of $D_A$ and $H^{-1}$, respectively. Hereby the ratio between transverse distance $D_A$ and the radial distance $H^{-1}$ reproduces the intrinsic shape of 2D BAO circle, which is {\it a priori} given by the known broadband shape of spectra. All BAO peaks of $\xi(\sigma,\pi)$ are precisely calculated with the improved theoretical model of redshift distortion. We test this broadband Alcock--Paczynski method using BOSS--like mock catalogues. The transverse and radial distances are probed in precision of several percentage fractional errors, and the coherent motion is observed to match with the fiducial values accurately.

Our observations of the Universe are fundamentally anisotropic, with data from galaxies separated transverse to the line of sight coming from the same epoch while that from galaxies separated parallel to the line of sight coming from different times. Moreover, galaxy velocities along the line of sight change their redshift, giving redshift space distortions. We perform a full two-dimensional anisotropy analysis of galaxy clustering data, fitting in a substantially model independent manner the angular diameter distance D_A, Hubble parameter H, and growth rate ddelta/dln a without assuming a dark energy model. The results demonstrate consistency with LCDM expansion and growth, hence also testing general relativity. We also point out the interpretation dependence of the effective redshift z_eff, and its cosmological impact for next generation surveys.

Numerical simulations show that redshift space distortions (RSD) introduce strong scale dependence in the power spectra of halos, with ten percent deviations relative to linear theory predictions even on relatively large scales (k<0.1h/Mpc) and even in the absence of satellites (which induce Fingers-of-God, FoG, effects). If unmodeled these effects prevent one from extracting cosmological information from RSD surveys. In this paper we use perturbation theory (PT) and halo biasing model and apply it to the distribution function approach to RSD, in which RSD is decomposed into several correlators of density weighted velocity moments. We model each of these correlators using PT and compare the results to simulations over a wide range of halo masses and redshifts. We find that with an introduction of a physically motivated halo biasing, and using dark matter power spectra from simulations, we can reproduce the simulation results at a percent level on scales up to k~0.15h/Mpc at z=0, without the need to have free FoG parameters in the model.

We develop a perturbative approach to redshift space distortions (RSD) using the phase space distribution function approach and apply it to the dark matter redshift space power spectrum and its moments. RSD can be written as a sum over density weighted velocity moments correlators, with the lowest order being density, momentum density and stress energy density. We use standard and extended perturbation theory (PT) to determine their auto and cross correlators, comparing them to N-body simulations. We show which of the terms can be modeled well with the standard PT and which need additional terms that include higher order corrections which cannot be modeled in PT. Most of these additional terms are related to the small scale velocity dispersion effects, the so called finger of god (FoG) effects, which affect some, but not all, of the terms in this expansion, and which can be approximately modeled using a simple physically motivated ansatz such as the halo model. We point out that there are several velocity dispersions that enter into the detailed RSD analysis with very different amplitudes, which can be approximately predicted by the halo model. In contrast to previous models our approach systematically includes all of the terms at a given order in PT and provides a physical interpretation for the small scale dispersion values. We investigate RSD power spectrum as a function of \mu, the cosine of the angle between the Fourier mode and line of sight, focusing on the lowest order powers of \mu and multipole moments which dominate the observable RSD power spectrum. Overall we find considerable success in modeling many, but not all, of the terms in this expansion.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

It was recently shown that the power spectrum in redshift space can be written as a sum of cross-power spectra between number weighted velocity moments. We investigate the properties of these power spectra for simulated galaxies and dark matter halos and compare them to the dark matter power spectra, generalizing the concept of the bias. Because all of the quantities are number weighted this approach is well defined even for sparse systems such as massive halos, in contrasts to the previous approaches to RSD where velocity correlations have been explored. We find that the number density weighting leads to a strong scale dependence of the bias terms for momentum density auto-correlation and cross-correlation with density. This trend becomes more significant for the more biased halos and leads to an enhancement of RSD power relative to the linear theory. Fingers-of-god effects, which in this formalism come from the correlations of the higher order moments beyond the momentum density, lead to smoothing of the power spectrum and can reduce this enhancement of power, but are relatively small for halos with no small-scale velocity dispersion. In comparison, for a more realistic galaxy sample with satellites the velocity dispersion generated by satellite motions inside the halos leads to a larger power suppression on small scales, but this depends on the satellite fraction. We investigate several statistics such as the two-dimensional power spectrum, its multipole moments, its powers of mu^2, and configuration space statistics. Overall we find that the nonlinear effects in realistic galaxy samples such as luminous red galaxies affect the redshift space clustering on very large scales: for example, the quadrupole moment is affected by 10% for k<0.1h/mpc, which means that these effects need to be understood if we want to extract cosmological information from the redshift space distortions.

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-1101/MOA-2011-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness profile of source stars by measuring the deviation of the light curves near the peak affected by the finite-source effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of $\theta_{\rm E}\sim 0.08$ mas combined with the short time scale of $t_{\rm E}\sim 2.7$ days suggests the possibility that the lens is a free-floating planet. For MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of the lens. We also find that the measured lens mass of $\sim 0.84\ M_\odot$ is consistent with that of a star blended with the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a function of the separation and mass ratio.

Measurement of redshift-space distortions (RSD) offers an attractive method to directly probe the cosmic growth history of density perturbations. A distribution function approach where RSD can be written as a sum over density weighted velocity moment correlators has recently been developed. We use Nbody simulations to investigate the individual contributions and convergence of this expansion for dark matter. If the series is expanded as a function of powers of mu, cosine of the angle between the Fourier mode and line of sight, there are a finite number of terms contributing at each order. We present these terms and investigate their contribution to the total as a function of wavevector k. For mu^2 the correlation between density and momentum dominates on large scales. Higher order corrections, which act as a Finger-of-God (FoG) term, contribute 1% at k~0.015h/Mpc, 10% at k~0.05h/Mpc at z=0, while for k>0.15h/Mpc they dominate and make the total negative. These higher order terms are dominated by density-energy density correlations which contribute negatively to the power, while the contribution from vorticity part of momentum density auto-correlation is an order of magnitude lower. For mu^4 term the dominant term on large scales is the scalar part of momentum density auto-correlation, while higher order terms dominate for k>0.15h/Mpc. For mu^6 and mu^8 we find it has very little power for k<0.15h/Mpc. We also compare the expansion to the full 2D P^ss(k,mu) as well as to their multipoles. For these statistics an infinite number of terms contribute and we find that the expansion achieves percent level accuracy for kmu<0.15h/Mpc at 6th order, but breaks down on smaller scales because the series is no longer perturbative. We explore resummation of the terms into FoG kernels, which extend the convergence up to a factor of 2 in scale. We find that the FoG kernels are approximately Lorentzian.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

We develop a new method which measures the projected density distribution w_p(r_p)n of photometric galaxies surrounding a set of spectroscopically-identified galaxies, and simultaneously the projected correlation function w_p(r_p) between the two populations. In this method we are able to divide the photometric galaxies into subsamples in luminosity intervals when redshift information is unavailable, enabling us to measure w_p(r_p)n and w_p(r_p) as a function of not only the luminosity of the spectroscopic galaxy, but also that of the photometric galaxy. Extensive tests show that our method can measure w_p(r_p) in a statistically unbiased way. The accuracy of the measurement depends on the validity of the assumption in the method that the foreground/background galaxies are randomly distributed and thus uncorrelated with those galaxies of interest. Therefore, our method can be applied to the cases where foreground/background galaxies are distributed in large volumes, which is usually valid in real observations. We applied our method to data from SDSS including a sample of 10^5 LRGs at z~0.4 and a sample of about half a million galaxies at z~0.1, both of which are cross-correlated with a deep photometric sample drawn from the SDSS. On large scales, the relative bias factor of galaxies measured from w_p(r_p) at z~0.4 depends on luminosity in a manner similar to what is found at z~0.1, which are usually probed by autocorrelations of spectroscopic samples. On scales smaller than a few Mpc and at both z~0.4 and z~0.1, the photometric galaxies of different luminosities exhibit similar density profiles around spectroscopic galaxies at fixed luminosity and redshift. This provides clear support for the assumption commonly-adopted in HOD models that satellite galaxies of different luminosities are distributed in a similar way, following the dark matter distribution within their host halos.

Since 1995, more than 500 exoplanets have been detected using different techniques, of which 11 were detected with gravitational microlensing. Most of these are gravitationally bound to their host stars. There is some evidence of free-floating planetary mass objects in young star-forming regions, but these objects are limited to massive objects of 3 to 15 Jupiter masses with large uncertainties in photometric mass estimates and their abundance. Here, we report the discovery of a population of unbound or distant Jupiter-mass objects, which are almost twice (1.8_{-0.8}^{+1.7}) as common as main-sequence stars, based on two years of gravitational microlensing survey observations toward the Galactic Bulge. These planetary-mass objects have no host stars that can be detected within about ten astronomical units by gravitational microlensing. However a comparison with constraints from direct imaging suggests that most of these planetary-mass objects are not bound to any host star. An abrupt change in the mass function at about a Jupiter mass favours the idea that their formation process is different from that of stars and brown dwarfs. They may have formed in proto-planetary disks and subsequently scattered into unbound or very distant orbits.

The linear growth factor of density perturbations is believed to be a powerful observable of future redshift surveys to probe physical properties of dark energy and to distinguish among gravity theories. We investigate systematic effects on determination of the growth factor f from a measurement of redshift-space distortions. Using N-body simulations we identify dark matter halos over a broad mass range. We compute the power spectra and correlation functions for the halos and then examine how well the redshift distortion parameter beta=f/b can be reconstructed as a function of halo mass. We find that beta measured for a fixed halo mass is generally a function of scale even on large scales, in contrast with the common expectation that beta approaches a constant described by Kaiser's formula on such scales. The scale dependence depends on the halo mass, being stronger for smaller halos. It also cannot be easily explained with the well-known distribution function of the halo peculiar velocities. We demonstrate that the biasing for smaller halos has larger nonlinearity and stochasticity, thus the linear bias assumption becomes worse for smaller halos. Only for massive halos with b>1.5, beta approaches the linear theory prediction on scales of r or pi/k>30Mpc/h. Luminous red galaxies (LRG), targeted by the SDSS-III's BOSS survey, tend to reside in very massive halos. Our results indicate that if the LRG is used for the measurement of redshift distortions, f can be measured unbiasedly. On the other hand, if one considers to use emission line galaxies, which are targeted by the BigBOSS survey and inhabited in halos of a broad mass range, the scale dependence of beta must be taken into account carefully; otherwise one might give incorrect constraints on dark energy or modified gravity theories. We also find that beta reconstructed in Fourier space behaves better than that in configuration space.

We present a new analysis of the Jupiter+Saturn analog system, OGLE-2006-BLG-109Lb,c, which was the first double planet system discovered with the gravitational microlensing method. This is the only multi-planet system discovered by any method with measured masses for the star and both planets. In addition to the signatures of two planets, this event also exhibits a microlensing parallax signature and finite source effects that provide a direct measure of the masses of the star and planets, and the expected brightness of the host star is confirmed by Keck AO imaging, yielding masses of M_* = 0.51(+0.05-0.04) M_sun, M_b = 231+-19 M_earth, M_c = 86+-7 M_earth. The Saturn-analog planet in this system had a planetary light curve deviation that lasted for 11 days, and as a result, the effects of the orbital motion are visible in the microlensing light curve. We find that four of the six orbital parameters are tightly constrained and that a fifth parameter, the orbital acceleration, is weakly constrained. No orbital information is available for the Jupiter-analog planet, but its presence helps to constrain the orbital motion of the Saturn-analog planet. Assuming co-planar orbits, we find an orbital eccentricity of eccentricity = 0.15 (+0.17-0.10) and an orbital inclination of i = 64 (+4-7) deg. The 95% confidence level lower limit on the inclination of i > 49 deg. implies that this planetary system can be detected and studied via radial velocity measurements using a telescope of >30m aperture.

Gravitational microlensing is not only a successful tool for discovering distant exoplanets, but it also enables characterization of the lens and source stars involved in the lensing event. In high magnification events, the lens caustic may cross over the source disk, which allows a determination of the angular size of the source and additionally a measurement of its limb darkening. When such extended-source effects appear close to maximum magnification, the resulting light curve differs from the characteristic Paczynski point-source curve. The exact shape of the light curve close to the peak depends on the limb darkening of the source. Dense photometric coverage permits measurement of the respective limb-darkening coefficients. In the case of microlensing event OGLE 2008-BLG-290, the K giant source star reached a peak magnification of about 100. Thirteen different telescopes have covered this event in eight different photometric bands. Subsequent light-curve analysis yielded measurements of linear limb-darkening coefficients of the source in six photometric bands. The best-measured coefficients lead to an estimate of the source effective temperature of about 4700 +100-200 K. However, the photometric estimate from colour-magnitude diagrams favours a cooler temperature of 4200 +-100 K. As the limb-darkening measurements, at least in the CTIO/SMARTS2 V and I bands, are among the most accurate obtained, the above disagreement needs to be understood. A solution is proposed, which may apply to previous events where such a discrepancy also appeared.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

We analyze OGLE-2007-BLG-050, a high magnification microlensing event (A ~ 432) whose peak occurred on 2 May, 2007, with pronounced finite-source and parallax effects. We compute planet detection efficiencies for this event in order to determine its sensitivity to the presence of planets around the lens star. Both finite-source and parallax effects permit a measurement of the angular Einstein radius \theta_E = 0.48 +/- 0.01 mas and the parallax \pi_E = 0.12 +/- 0.03, leading to an estimate of the lens mass M = 0.50 +/- 0.14 M_Sun and its distance to the observer D_L = 5.5 +/- 0.4 kpc. This is only the second determination of a reasonably precise (<30%) mass estimate for an isolated unseen object, using any method. This allows us to calculate the planetary detection efficiency in physical units (r_\perp, m_p), where r_\perp is the projected planet-star separation and m_p is the planet mass. When computing planet detection efficiency, we did not find any planetary signature and our detection efficiency results reveal significant sensitivity to Neptune-mass planets, and to a lesser extent Earth-mass planets in some configurations. Indeed, Jupiter and Neptune-mass planets are excluded with a high confidence for a large projected separation range between the planet and the lens star, respectively [0.6 - 10] and [1.4 - 4] AU, and Earth-mass planets are excluded with a 10% confidence in the lensing zone, i.e. [1.8 - 3.1] AU.

We investigate the orientation correlation of giant elliptical galaxies by measuring the intrinsic ellipticity correlation function of 83,773 luminous red galaxies (LRGs) at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We have accurately determined the correlation up to 30 $h^{-1}$Mpc. Luminosity dependence of the ellipticity correlation is also detected although the error bars are large, while no evidence is found for its redshift evolution between $z=0.2$ and $z=0.4$. Then we use a cosmological $N$-body simulation to examine misalignment between the central LRGs and their parent dark matter halos. Central and satellite galaxies are assigned to simulated halos by employing a halo occupation distribution model for the LRGs. The ellipticity correlation is predicted to have the same shape as but an amplitude about 4 times higher than our observation if the central LRGs are perfectly aligned with their host halos. This indicates that the central LRG galaxies are preferentially but not perfectly aligned with their host halos. With the assumption that there is a misalignment angle between a central LRG and its host halo which follows a Gaussian distribution with a zero mean and a width $\sigma_\theta$, we obtain a tight constraint on the misalignment parameter, $\sigma_\theta={35.4}^{+4.0}_{-3.3}$ deg. This type of intrinsic ellipticity correlation, if not corrected, can lead to contamination at 5% level to the shear power spectrum in weak lensing surveys of limiting magnitude $R_{AB}=24.5$ if the source central galaxies follow the same misalignment distribution as the LRGs.

We examine whether the gravitational shear--intrinsic ellipticity (GI) correlation function of the luminous red galaxies (LRGs) can be modeled with the distribution function of a misalignment angle advocated recently by Okumura et al. For this purpose, we have accurately measured the GI correlation for the LRGs in the Data Release 6 (DR6) of the Sloan Digital Sky Survey (SDSS), which confirms the results of Hirata et al. who used the DR4 data. By comparing the GI correlation functions in the simulation and in the observation, we find that the GI correlation can be modeled in the current $\Lambda$CDM model if the misalignment follows a Gaussian distribution with a zero mean and a typical misalignment angle $\sigma_\theta=34.9^{+1.9}_{-2.1}$ degrees. We also find a correlation between the axis ratios and intrinsic alignments of LRGs. This effect should be taken into account in theoretical modeling of the GI and intrinsic ellipticity--ellipticity correlations for weak lensing surveys.

We report the detection of the cool, Jovian-mass planet MOA-2007-BLG-400Lb. The planet was detected in a high-magnification microlensing event (with peak magnification A_max = 628) in which the primary lens transited the source, resulting in a dramatic smoothing of the peak of the event. The angular extent of the region of perturbation due to the planet is significantly smaller than the angular size of the source, and as a result the planetary signature is also smoothed out by the finite source size. Thus the deviation from a single-lens fit is broad and relatively weak (~ few percent). Nevertheless, we demonstrate that the planetary nature of the deviation can be unambiguously ascertained from the gross features of the residuals, and detailed analysis yields a fairly precise planet/star mass ratio of q = 0.0026+/-0.0004, in accord with the large significance (\Delta\chi^2=1070) of the detection. The planet/star projected separation is subject to a strong close/wide degeneracy, leading to two indistinguishable solutions that differ in separation by a factor of ~8.5. Upper limits on flux from the lens constrain its mass to be M < 0.75 M_Sun (assuming it is a main-sequence star). A Bayesian analysis that includes all available observational constraints indicates a primary in the Galactic bulge with a mass of ~0.2-0.5 M_Sun and thus a planet mass of ~ 0.5-1.3 M_Jupiter. The separation and equilibrium temperature are ~0.6-1.1AU (~5.3-9.7AU) and ~103K (~34K) for the close (wide) solution. If the primary is a main-sequence star, follow-up observations would enable the detection of its light and so a measurement of its mass and distance.

We report the detection of an extrasolar planet of mass ratio q ~ 2 x 10^(-4) in microlensing event MOA-2007-BLG-192. The best fit microlensing model shows both the microlensing parallax and finite source effects, and these can be combined to obtain the lens masses of M = 0.060 (+0.028 -0.021) M_sun for the primary and m = 3.3 (+4.9 -1.6) M_earth for the planet. However, the observational coverage of the planetary deviation is sparse and incomplete, and the radius of the source was estimated without the benefit of a source star color measurement. As a result, the 2-sigma limits on the mass ratio and finite source measurements are weak. Nevertheless, the microlensing parallax signal clearly favors a sub-stellar mass planetary host, and the measurement of finite source effects in the light curve supports this conclusion. Adaptive optics images taken with the Very Large Telescope (VLT) NACO instrument are consistent with a lens star that is either a brown dwarf or a star at the bottom of the main sequence. Follow-up VLT and/or Hubble Space Telescope (HST) observations will either confirm that the primary is a brown dwarf or detect the low-mass lens star and enable a precise determination of its mass. In either case, the lens star, MOA-2007-BLG-192L, is the lowest mass primary known to have a companion with a planetary mass ratio, and the planet, MOA-2007-BLG-192Lb, is probably the lowest mass exoplanet found to date, aside from the lowest mass pulsar planet.

We study the large-scale anisotropic two-point correlation function using 46,760 luminous red galaxies at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We measure the correlation function as a function of separations parallel and perpendicular to the line-of-sight in order to take account of anisotropy of the large-scale structure in redshift space. We find a slight signal of baryonic features in the anisotropic correlation function, i.e., a ``baryon ridge'' which corresponds to a baryon acoustic peak in the spherically averaged correlation function which has already been reported using the same sample. The baryon ridge has primarily a spherical structure with a known radius in comoving coordinates. It enables us to divide the redshift distortion effects into dynamical and geometrical components and provides further constraints on cosmological parameters, including the dark energy equation-of-state. With an assumption of a flat $\Lambda$ cosmology, we find the best-fit values of $\Omega_{\rm m} = 0.218^{+0.047}_{-0.037}$ and $\Omega_{\rm b} = 0.047^{+0.016}_{-0.016}$ (68% C.L.) when we use the overall shape of the anisotropic correlation function of $40<s<200\himpc$ including a scale of baryon acoustic oscillations. When an additional assumption $\Omega_{\rm b}h^2=0.024$ is adopted, we obtain $\Omega_{\rm DE}=0.770^{+0.051}_{-0.040}$ and $w=-0.93^{+0.45}_{-0.35}$. These constraints are estimated only from our data of the anisotropic correlation function, and they agree quite well with values both from the cosmic microwave background (CMB) anisotropies and from other complementary statistics using the LRG sample. With the CMB prior from the 3 year WMAP results, we give stronger constraints on those parameters.

Searches for extrasolar planets have uncovered an astonishing diversity of planetary systems, yet the frequency of solar system analogs remains unknown. The gravitational microlensing planet search method is potentially sensitive to multiple-planet systems containing analogs of all the solar system planets except Mercury. We report the detection of a multiple-planet system with microlensing. We identify two planets with masses of ~0.71 and ~0.27 times the mass of Jupiter and orbital separations of ~2.3 and ~4.6 astronomical units orbiting a primary star of mass ~0.50 solar masses at a distance of ~1.5 kiloparsecs. This system resembles a scaled version of our solar system in that the mass ratio, separation ratio, and equilibrium temperatures of the planets are similar to those of Jupiter and Saturn. These planets could not have been detected with other techniques; their discovery from only six confirmed microlensing planet detections suggests that solar system analogs may be common.