Source author record

Atsushi Taruya

Atsushi Taruya appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

69works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

69 published item(s)

preprint2022arXiv

Detectability of the gravitational redshift effect from the asymmetric galaxy clustering

It has been recently recognized that the observational relativistic effects, mainly arising from the light propagation in an inhomogeneous universe, induce the dipole asymmetry in the cross-correlation function of galaxies. In particular, the dipole asymmetry at small scales is shown to be dominated by the gravitational redshift effects. In this paper, we exploit a simple analytical description for the dipole asymmetry in the cross-correlation function valid at quasi-linear regime. In contrast to the previous model, a new prescription involves only one dimensional integrals, providing a faster way to reproduce the results obtained by Saga et al. (2020). Using the analytical model, we discuss the detectability of the dipole signal induced by the gravitational redshift effect from upcoming galaxy surveys. The gravitational redshift effect at small scales enhances the signal-to-noise ratio (S/N) of the dipole, and in most of the cases considered, the S/N is found to reach a maximum at $z\approx0.5$. We show that current and future surveys such as DESI and SKA provide an idealistic data set, giving a large S/N of $10\sim 20$. Two potential systematics arising from off-centered galaxies are also discussed (transverse Doppler effect and diminution of the gravitational redshift effect), and their impacts are found to be mitigated by a partial cancellation between two competitive effects. Thus, the detection of the dipole signal at small scales is directly linked to the gravitational redshift effect, and should provide an alternative route to test gravity.

preprint2022arXiv

Tightening geometric and dynamical constraints on dark energy and gravity: galaxy clustering, intrinsic alignment and kinetic Sunyaev-Zel'dovich effect

Conventionally, in galaxy surveys, cosmological constraints on the growth and expansion history of the universe have been obtained from the measurements of redshift-space distortions and baryon acoustic oscillations embedded in the large-scale galaxy density field. In this paper, we study how well one can improve the cosmological constraints from the combination of the galaxy density field with velocity and tidal fields, which are observed via the kinetic Sunyaev-Zel'dovich (kSZ) and galaxy intrinsic alignment (IA) effects, respectively. For illustration, we consider the deep galaxy survey by Subaru Prime Focus Spectrograph, whose survey footprint perfectly overlaps with the imaging survey of the Hyper Suprime-Cam and the CMB-S4 experiment. We find that adding the kSZ and IA effects significantly improves cosmological constraints, particularly when we adopt the non-flat cold dark matter model which allows both time variation of the dark energy equation-of-state and deviation of the gravity law from general relativity. Under this model, we achieve $31\%$ improvement for the growth index $γ$ and $>35\%$ improvement for other parameters except for the curvature parameter, compared to the case of the conventional galaxy-clustering-only analysis. As another example, we also consider the wide galaxy survey by the {\it Euclid} satellite, in which shapes of galaxies are noisier but the survey volume is much larger. We demonstrate that when the above model is adopted, the clustering analysis combined with kSZ and IA from the deep survey can achieve tighter cosmological constraints than the clustering-only analysis from the wide survey.

preprint2021arXiv

Covariance of the matter power spectrum including the survey window function effect: N-body simulations vs. fifth-order perturbation theory on grid

We present a Next-to-next-to-leading (fifth or NNLO) order calculation for the covariance matrix of the matter power spectrum, taking into account the effect of survey window functions. Using the grid-based calculation scheme for the standard perturbation theory, GridSPT, we quickly generate multiple realizations of the nonlinear density fields to fifth order in perturbation theory, then estimate the power spectrum and the covariance matrix from the sample. To the end, we have obtained the non-Gaussian covariance originated from the one-loop trispectrum without explicitly computing the trispectrum. By comparing the GridSPT calculations with the N-body results, we show that NNLO GridSPT result reproduces the N-body results on quasi-linear scales where SPT accurately models nonlinear matter power spectrum. Incorporating the survey window function effect to GridSPT is rather straightforward, and the resulting NNLO covariance matrix also matches well with the N-body results.

preprint2021arXiv

Grid-based calculations of redshift-space matter fluctuations from perturbation theory: UV sensitivity and convergence at the field level

Perturbation theory (PT) has been used to interpret the observed nonlinear large-scale structure statistics at the quasi-linear regime. To facilitate the PT-based analysis, we have presented the GridSPT algorithm, a grid-based method to compute the nonlinear density and velocity fields in standard perturbation theory (SPT) from a given linear power spectrum. Here, we further put forward the approach by taking the redshift-space distortions into account. With the new implementation, we have, for the first time, generated the redshift-space density field to the fifth order and computed the next-to-next-to-leading order (2 loop) power spectrum and the next-to-leading order (1 loop) bispectrum of matter clustering in redshift space. By comparing the result with corresponding analytical SPT calculation and $N$-body simulations, we find that the SPT calculation (A) suffers much more from the UV sensitivity due to the higher-derivative operators and (B) deviates from the $N$-body results from the Fourier wavenumber smaller than real space $k_{\rm max}$. Finally, we have shown that while Padé approximation removes spurious features in morphology, it does not improve the modeling of power spectrum and bispectrum.

preprint2021arXiv

Wide-angle effects on galaxy ellipticity correlations

We show an efficient way to compute wide-angle or all-sky statistics of galaxy intrinsic alignment in three-dimensional configuration space. For this purpose, we expand the two-point correlation function using a newly introduced spin-dependent tripolar spherical harmonic basis. Therefore, the angular dependences on the two line-of-sight (LOS) directions pointing to each pair of objects, which are degenerate with each other in the conventional analysis under the small-angle or plane-parallel (PP) approximation, are unambiguously decomposed. By means of this, we, for the first time, compute the wide-angle auto and cross correlations between intrinsic ellipticities, number densities and velocities of galaxies, and compare them with the PP-limit results. For the ellipticity-ellipticity and density-ellipticity correlations, we find more than $10\%$ deviation from the PP-limit results if the opening angle between two LOS directions exceeds $30^\circ - 50^\circ$. It is also shown that even if the PP-limit result is strictly zero, the non-vanishing correlation is obtained over the various scales, arising purely from the curved-sky effects. Our results indicate the importance of the data analysis not relying on the PP approximation in order to determine the cosmological parameters more precisely and/or find new physics via ongoing and forthcoming wide-angle galaxy surveys.

preprint2020arXiv

Anisotropies of galaxy ellipticity correlations in real and redshift space: angular dependence in linear tidal alignment model

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.

preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

preprint2020arXiv

Fitting the nonlinear matter bispectrum by the Halofit approach

We provide a new fitting formula of the matter bispectrum in the nonlinear regime calibrated by high-resolution cosmological $N$-body simulations of $41$ cold dark matter ($w$CDM, $w=$ constant) models around the Planck 2015 best-fit parameters. As the parameterization in our fitting function is similar to that in Halofit, our fitting is named BiHalofit. The simulation volume is sufficiently large ($> 10 \, {\rm Gpc}^3$) to cover almost all measurable triangle bispectrum configurations in the universe. The function is also calibrated using one-loop perturbation theory at large scales ($k<0.3 \, h \, {\rm Mpc}^{-1}$). Our formula reproduced the matter bispectrum to within $10 \, (15) \, \%$ accuracy in the Planck 2015 model at wavenumber $k< 3 \, (10) \, h \, {\rm Mpc}^{-1}$ and redshifts $z=0-3$. The other $40$ $w$CDM models obtained poorer fits, with accuracy approximating $20 \, \%$ at $k<3 \, h \, {\rm Mpc}^{-1}$ and $z=0-1.5$ (the deviation includes the $10 \, \%$-level sample variance of the simulations). We also provide a fitting formula that corrects the baryonic effects such as radiative cooling and active galactic nucleus feedback, using the latest hydrodynamical simulation IllustrisTNG. We demonstrate that our new formula more accurately predicts the weak-lensing bispectrum than the existing fitting formulas. This formula will assist current and future weak-lensing surveys and cosmic microwave background lensing experiments. Numerical codes of the formula are available, written in Python, C and Fortran.

preprint2020arXiv

Improving geometric and dynamical constraints on cosmology with intrinsic alignments of galaxies

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.

preprint2020arXiv

Modelling the asymmetry of the halo cross-correlation function with relativistic effects at quasi-linear scales

The observed galaxy distribution via galaxy redshift surveys appears distorted due to redshift-space distortions (RSD). While one dominant contribution to RSD comes from the Doppler effect induced by the peculiar velocity of galaxies, the relativistic effects, including the gravitational redshift effect, are recently recognized to give small but important contributions. Such contributions lead to an asymmetric galaxy clustering along the line of sight, and produce non-vanishing odd multipoles when cross-correlating between different biased objects. However, non-zero odd multipoles are also generated by the Doppler effect beyond the distant-observer approximation, known as the wide-angle effect, and at quasi-linear scales, the interplay between wide-angle and relativistic effects becomes significant. In this paper, based on the formalism developed by Taruya et al., we present a quasi-linear model of the cross-correlation function taking a proper account of both the wide-angle and gravitational redshift effects, as one of the major relativistic effects. Our quasi-linear predictions of the dipole agree well with simulations even at the scales below $20\,h^{-1}\,$Mpc, where non-perturbative contributions from the halo potential play an important role, flipping the sign of the dipole amplitude. When increasing the bias difference and redshift, the scale where the sign flip happens is shifted to a larger scale. We derive a simple approximate formula to quantitatively account for the behaviors of the sign flip.

preprint2020arXiv

Observing Baryonic Acoustic Oscillations in tomographic cosmic shear surveys

We show that it is possible to build effective matter density power spectra in tomographic cosmic shear observations that exhibit the Baryonic Acoustic Oscillations (BAO) features once a nulling transformation has been applied to the data. The precision with which the amplitude and position of these features can be reconstructed is quantified in terms of sky coverage, intrinsic shape noise, median source redshift and number density of sources. BAO detection in Euclid or LSST like wide surveys will be possible with a modest signal-to-noise ratio. It would improve dramatically for slightly deeper surveys.

preprint2020arXiv

Phase-space structure of cold dark matter halos inside splashback: multi-stream flows and self-similar solution

Using the motion of accreting particles onto halos in cosmological $N$-body simulations, we study the radial phase-space structures of cold dark matter (CDM) halos. In CDM cosmology, formation of virialized halos generically produces radial caustics, followed by multi-stream flows of accreted dark matter inside the halos. In particular, the radius of the outermost caustic called the splashback radius exhibits a sharp drop in the slope of the density profile. Here, we focus on the multi-stream structure of CDM halos inside the splashback radius. To analyze this, we use and extend the SPARTA algorithm developed by Diemer. By tracking the particle trajectories accreting onto the halos, we count their number of apocenter passages, which is then used to reveal the multi-stream flows of the dark matter particles. The resultant multi-stream structure in radial phase space is compared with the prediction of the self-similar solution by Fillmore & Goldreich for each halo. We find that $\sim30\%$ of the simulated halos satisfy our criteria to be regarded as being well fitted to the self-similar solution. The fitting parameters in the self-similar solution characterizes physical properties of the halos, including the mass accretion rate and the size of the outermost caustic (i.e., the splashback radius). We discuss in detail the correlation of these fitting parameters and other measures directly extracted from the $N$-body simulation.

preprint2020arXiv

Testing tidal alignment models for anisotropic correlations of halo ellipticities with N-body simulations

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.

preprint2019arXiv

Wide-angle redshift-space distortions at quasi-linear scales: cross-correlation functions from Zel'dovich approximation

Redshift-space distortions (RSD) in galaxy redshift surveys generally break both the isotropy and homogeneity of galaxy distribution. While the former aspect is particularly highlighted as a probe of growth of structure induced by gravity, the latter aspect, often quoted as wide-angle RSD but ignored in most of the cases, will become important and critical to account for as increasing the statistical precision in next-generation surveys. However, the impact of wide-angle RSD has been mostly studied using linear perturbation theory. In this paper, employing the Zel'dovich approximation, i.e., first-order Lagrangian perturbation theory for gravitational evolution of matter fluctuations, we present a quasi-linear treatment of wide-angle RSD, and compute the cross-correlation function. The present formalism consistently reproduces linear theory results, and can be easily extended to incorporate relativistic corrections (e.g., gravitational redshift).

preprint2016arXiv

Anisotropies of gravitational-wave standard sirens as a new cosmological probe without redshift information

Gravitational waves (GWs) from compact binary stars at cosmological distances are promising and powerful cosmological probes, referred to as the GW standard sirens. With future GW detectors, we will be able to precisely measure source luminosity distances out to a redshift $z\sim5$. To extract cosmological information, previously proposed cosmological studies using the GW standard sirens rely on source redshift information obtained through an extensive electromagnetic follow-up campaign. However, the redshift identification is typically time consuming and rather challenging. Here, we propose a novel method for cosmology with the GW standard sirens free from the redshift measurements. Utilizing the anisotropies of the number density and luminosity distances of compact binaries originated from the large-scale structure, we show that, once GW observations will be well established in the future, (i) these anisotropies can be measured even at very high redshifts ($z\geq 2$), where the identification of the electromagnetic counterpart is difficult, (ii) the expected constraints on the primordial non-Gaussianity with the Einstein Telescope would be comparable to or even better than the other large-scale structure probes at the same epoch, and (iii) the cross-correlation with other cosmological observations is found to have high-statistical significance, providing additional cosmological information at very high redshifts.

preprint2016arXiv

Constraining higher-order parameters for primordial non-Gaussianities from power spectra and bispectra of imaging survey

We investigate the statistical power of higher-order statistics and cross-correlation statistics to constrain the primordial non-Gaussianity from the imaging surveys. In particular, we consider the local-type primordial non- Gaussianity and discuss how well one can tightly constrain the higher-order non-Gaussian parameters ($g_{\rm NL}$ and $τ_{\rm NL}$) as well as the leading order parameter $f_{\rm NL}$ from the halo/galaxy clustering and weak gravitational lensing measurements. Making use of a strong scale-dependent behavior in the galaxy/halo clustering, Fisher matrix analysis reveals that the bispectra can break the degeneracy between non-Gaussian parameters ($f_{\rm NL}$, $g_{\rm NL}$ and $τ_{\rm NL}$) and this will give simultaneous constraints on those three parameters. The combination of cross-correlation statistics further improves the constraints by factor of 2. As a result, upcoming imaging surveys like the Large Synoptic Survey Telescope have the potential to improve the constraints on the primordial non-Gaussianity much tighter than those obtained from the CMB measurement by Planck, giving us an opportunity to test the single-sourced consistency relation, $τ_{\rm NL} \ge (36/25) f_{\rm NL}^2$.

preprint2016arXiv

Constructing perturbation theory kernels for large-scale structure in generalized cosmologies

We present a simple numerical scheme for perturbation theory (PT) calculations of large-scale structure. Solving the evolution equations for perturbations numerically, we construct the PT kernels as building blocks of statistical calculations, from which the power spectrum and/or correlation function can be systematically computed. The scheme is especially applicable to the generalized structure formation including modified gravity, in which the analytic construction of PT kernels is intractable. As an illustration, we show several examples for power spectrum calculations in $f(R)$ gravity and $Λ$CDM models.

preprint2016arXiv

Detecting Black-Hole Binary Clustering via the Second-Generation Gravitational-Wave Detectors

The first discovery of the gravitational wave (GW) event, GW150914, suggests a higher merger rate of black-hole (BH) binaries. If this is true, a number of BH binaries will be observed via the second-generation GW detectors, and the statistical properties of the observed BH binaries can be scrutinized. A naive but important question to ask is whether the spatial distribution of BH binaries faithfully traces the matter inhomogeneities in the Universe or not. Although the BH binaries are thought to be formed inside the galaxies in most of the scenarios, there is no observational evidence to confirm such a hypothesis. Here, we estimate how well the second-generation GW detectors can statistically confirm the BH binaries to be a tracer of the large-scale structure by looking at the auto- and cross-correlation of BH binaries with photometric galaxies and weak lensing measurements, finding that, with a 3 year observation, the $>3σ$ detection of non-zero signal is possible if the BH merger rate today is $\dot{n}_0\gtrsim100$ Gpc$^{-3}$yr$^{-1}$ and the clustering bias of BH binaries is $b_{\rm BH,0}\gtrsim1.5$.

preprint2016arXiv

Response function of the large-scale structure of the universe to the small scale inhomogeneities

In order to infer the impact of the small-scale physics to the large-scale properties of the universe, we use a series of cosmological $N$-body simulations of self-gravitating matter inhomogeneities to measure, for the first time, the response function of such a system defined as a functional derivative of the nonlinear power spectrum with respect to its linear counterpart. Its measured shape and amplitude are found to be in good agreement with perturbation theory predictions except for the coupling from small to large-scale perturbations. The latter is found to be significantly damped, following a Lorentzian form. These results shed light on validity regime of perturbation theory calculations giving a useful guideline for regularization of small scale effects in analytical modeling. Most importantly our result indicates that the statistical properties of the large-scale structure of the universe are remarkably insensitive to the details of the small-scale physics, astrophysical or gravitational, paving the way for the derivation of robust estimates of theoretical uncertainties on the determination of cosmological parameters from large-scale survey observations.

preprint2016arXiv

The Subaru FMOS galaxy redshift survey (FastSound). IV. New constraint on gravity theory from redshift space distortions at $z\sim 1.4$

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 $σ_{\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)σ_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)σ_8(z)$. This corresponds to $4.2σ$ detection of RSD. Our constraint is consistent with the prediction of general relativity $fσ_8\sim 0.392$ within the $1-σ$ confidence level. When we allow $σ_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $fσ_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $fσ_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

preprint2015arXiv

Consistent Modified Gravity Analysis of Anisotropic Galaxy Clustering Using BOSS DR11

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.

preprint2015arXiv

Cosmology with anisotropic galaxy clustering from the combination of power spectrum and bispectrum

The apparent anisotropies of the galaxy clustering in observable redshift space provide a unique opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales via the Alcock--Paczynski effect and redshift-space distortions. While the improved theoretical models have been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear scales, the applicability of these models is still limited in the presence of the non--perturbative smearing effect caused by the randomness of the relative velocities. Although the cosmological constraint from the anisotropic clustering will be certainly improved with a more elaborate theoretical model, we here consider an alternative way by using the statistical power of both the power spectrum and bispectrum at large scales. Based on the Fisher matrix analysis, we estimate the benefit of combining the power spectra and bispectra, finding that the constraints on the cosmic expansion and growth of structure will be improved by a factor of two. This compensates for the loss of constraining power using the power spectrum alone due to the randomness of the relative velocities.

preprint2015arXiv

Full-sky lensing reconstruction of gradient and curl modes from CMB maps

We present a method of lensing reconstruction on the full sky, by extending the optimal quadratic estimator proposed by Okamoto & Hu (2003) to the case including the curl mode of deflection angle. The curl mode is induced by the vector and tensor metric perturbations, and the reconstruction of the curl mode would be a powerful tool to not only check systematics in the estimated gradient mode but also probe any vector and tensor sources. We find that the gradient and curl modes can be reconstructed separately, thanks to the distinctive feature in the parity symmetry between the gradient and curl modes. We compare our estimator with the flat-sky estimator proposed by Cooray et al (2005). Based on the new formalism, the expected signal-to-noise ratio of the curl mode produced by the primordial gravitational-waves and a specific model of cosmic strings are estimated, and prospects for future observations are discussed.

preprint2015arXiv

Future detectability of gravitational-wave induced lensing from high-sensitivity CMB experiments

We discuss the future detectability of gravitational-wave induced lensing from high-sensitivity cosmic microwave background (CMB) experiments. Gravitational waves can induce a rotational component of the weak-lensing deflection angle, usually referred to as the curl mode, which would be imprinted on the CMB maps. Using the technique of reconstructing lensing signals involved in CMB maps, this curl mode can be measured in an unbiased manner, offering an independent confirmation of the gravitational waves complementary to the B-mode polarization experiments. Based on the Fisher matrix analysis, we first show that with the noise levels necessary to confirm the consistency relation for the primordial gravitational waves, the future CMB experiments will be able to detect the gravitational-wave induced lensing signals. For a tensor-to-scalar ratio of $r < 0.1$, even if the consistency relation is difficult to confirm with a high significance, the gravitational-wave induced lensing would be detected at more than $3\,σ$ significance level. Further, we point out that high-sensitivity experiments will be also powerful to constrain the gravitational waves generated after the recombination epoch. Compared to the B-mode polarization, the curl mode is particularly sensitive to gravitational waves generated at low redshifts ($z < 10$) with a low frequency ($k < 10^{-3}$ Mpc$^{-1}$), and it could give a much tighter constraint on their energy density $Ω_{\rm GW}$ by more than three orders of magnitude.

preprint2014arXiv

Beyond consistency test of gravity with redshift-space distortions at quasi-linear scales

Redshift-space distortions (RSD) offer an attractive method to measure the growth of cosmic structure on large scales, and combining with the measurement of the cosmic expansion history, it can be used as cosmological tests of gravity. With the advent of future galaxy redshift surveys aiming at precisely measuring the RSD, an accurate modeling of RSD going beyond linear theory is a critical issue in order to detect or disprove small deviations from general relativity (GR). While several improved models of RSD have been recently proposed based on the perturbation theory (PT), the framework of these models heavily relies on GR. Here, we put forward a new PT prescription for RSD in general modified gravity models. As a specific application, we present theoretical predictions of the redshift-space power spectra in f(R) gravity model, and compare them with N-body simulations. Using the PT template that takes into account the effects of both modifications of gravity and RSD properly, we successfully recover the fiducial model parameter in N-body simulations in an unbiased way. On the other hand, we found it difficult to detect the scale dependence of the growth rate in a model-independent way based on GR templates.

preprint2014arXiv

Broadband Alcock-Paczynski test exploiting redshift 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, $ξ(σ,π)$, can form a great circle (hereafter 2D BAO circle) in the $σ$ and $π$ plane, where $σ$ and $π$ 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 $ξ(σ,π)$ 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.

preprint2014arXiv

Cosmic shear full nulling: sorting out dynamics, geometry and systematics

An explicit full nulling scheme for cosmic shear observations is presented. It makes possible the construction of shear maps from extended source distributions for which the lens distance distribution is restricted to a definite interval. Such a construction allows to build totally independent shear maps, at all scales, that can be taken advantage of to constrain background cosmological parameters and systematics using the full statistical power of cosmic shear observations. Another advantage of such construction is that, as the lens redshift distribution can be made arbitrarily narrow, scale mixing due to projection effects can be limited allowing controlled predictions on the large scale shear power spectrum from perturbation theory calculations.

preprint2014arXiv

Halo/Galaxy Bispectrum with Primordial non-Gaussianity from integrated Perturbation Theory (iPT)

We derive a formula for the halo/galaxy bispectrum on the basis of the integrated Perturbation Theory (iPT). In addition to the gravity-induced non-Gaussianity, we consider the non-Gaussianity of the primordial curvature perturbations, and investigate in detail the effect of such primordial non-Gaussianity on the large-scale halo/galaxy bispectrum. In iPT, the effects of primordial non-Gaussianity are wholly encapsulated in the linear (primordial) polyspectra, and we systematically calculate the contributions to the large-scale behaviors arising from the three types of primordial bispectrum (local-, equilateral-, and orthogonal-types), and primordial trispectrum of the local-type non-Gaussianity. We find that the equilateral- and orthogonal-type non-Gaussianities show distinct scale-dependent behaviors which can dominate the gravity-induced non-Gaussianity at very large scales. For the local-type non-Gaussianity, higher-order loop corrections are found to give a significantly large contribution to the halo/galaxy bispectrum of the squeezed shape, and eventually dominate over the other contributions on large scales. A diagrammatic approach based on the iPT helps us to systematically investigate an impact of such higher-order contributions to the large-scale halo/galaxy bispectrum.

preprint2014arXiv

Modeling the anomaly of surface number densities of galaxies on the Galactic extinction map due to their FIR emission contamination

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μ{\rm m}$ and $r$-band luminosities of each galaxy, $y \equiv (νL)_{100μ{\rm m}}/(ν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.

preprint2014arXiv

Regularized cosmological power spectrum and correlation function in modified gravity models

Based on the multi-point propagator expansion, we present resummed perturbative calculations for cosmological power spectra and correlation functions in the context of modified gravity. In a wide class of modified gravity models that have a screening mechanism to recover general relativity (GR) on small scales, we apply the eikonal approximation to derive the governing equation for resummed propagator that partly includes the non-perturbative effect in the high-k limit. The resultant propagator in the high-k limit contains the new corrections arising from the screening mechanism as well as the standard exponential damping. We explicitly derive the expression for new high-k contributions in specific modified gravity models, and find that in the case of f(R) gravity for a currently constrained model parameter, the corrections are basically of the sub-leading order and can be neglected. Thus, in f(R) gravity, similarly to the GR case, we can analytically construct the regularized propagator that reproduces both the resummed high-k behavior and the low-k results computed with standard perturbation theory, consistently taking account of the nonlinear modification of gravity valid at large scales. With the regularized multi-point propagators, we give predictions for power spectrum and correlation function at one-loop order, and compare those with N-body simulations in f(R) gravity model. As an important application, we also discuss the redshift-space distortions and compute the anisotropic power spectra and correlation functions.

preprint2014arXiv

Simultaneous constraints on the growth of structure and cosmic expansion from the multipole power spectra of the SDSS DR7 LRG sample

The anisotropic galaxy clustering on large scales provides us with a unique opportunity to probe into the gravity theory through the redshift-space distortions (RSDs) and the Alcock-Paczynski effect. Using the multipole power spectra up to hexadecapole (ell=4), of the Luminous Red Galaxy (LRG) sample in the data release 7 (DR7) of the Sloan Digital Sky Survey II (SDSS-II), we obtain simultaneous constraints on the linear growth rate f, angular diameter distance D_A, and Hubble parameter H at redshift z = 0.3. For this purpose, we first extensively examine the validity of a theoretical model for the non-linear RSDs using mock subhalo catalogues from N-body simulations, which are constructed to match with the observed multipole power spectra. We show that the input cosmological parameters of the simulations can be recovered well within the error bars by comparing the multipole power spectra of our theoretical model and those of the mock subhalo catalogues. We also carefully examine systematic uncertainties in our analysis by testing the dependence on prior assumption of the theoretical model and the range of wavenumbers to be used in the fitting. These investigations validate that the theoretical model can be safely applied to the real data. Thus, our results from the SDSS DR7 LRG sample are robust including systematics of theoretical modeling; f(z = 0.3) sigma_8(z = 0.3) =0.49+-0.08(stat.)+-0.04(sys.), D_A (z = 0.3) =968+-42(stat.)+-17(sys.)[Mpc], H (z = 0.3) =81.7+-5.0(stat.)+-3.7(sys.)[km/s/Mpc]. We believe that our method to constrain the cosmological parameters using subhaloes catalogues will be useful for more refined samples like CMASS and LOWZ catalogues in the Baryon Oscillation Spectroscopic Survey in SDSS-III.

preprint2014arXiv

Tidal evolution of the spin-orbit angle in exoplanetary systems

The angle between the stellar spin and the planetary orbit axes (spin-orbit angle) is supposed to carry valuable information on the initial condition of the planet formation and the subsequent migration history. Indeed current observations of the Rossiter- McLaughlin effect have revealed a wide range of spin-orbit misalignments for transiting exoplanets. We examine in detail the tidal evolution of a simple system comprising a Sun-like star and a hot Jupiter adopting the equilibrium tide and the inertial wave dissipation effects simultaneously. We find that the combined tidal model works as a very efficient realignment mechanism; it predicts three distinct states of the spin-orbit angle (i.e., parallel, polar, and anti-parallel orbits) for a while, but the latter two states eventually approach the parallel spin-orbit configuration. The intermediate spin-orbit angles as measured in recent observations are difficult to be achieved. Therefore the current model cannot reproduce the observed broad distribution of the spin-orbit angles, at least in its simple form. This indicates that the observed diversity of the spin-orbit angles may emerge from more complicated interactions with outer planets and/or may be the consequence of the primordial misalignment between the proto-planetary disk and the stellar spin, which requires future detailed studies.

preprint2013arXiv

Characterization of the KOI-94 System with Transit Timing Variation Analysis: Implication for the Planet-Planet Eclipse

The KOI-94 system is a closely-packed, multi-transiting planetary system discovered by the Kepler space telescope. It is known as the first system that exhibited a rare event called a "planet-planet eclipse (PPE)," in which two planets partially overlap with each other in their double-transit phase. In this paper, we constrain the parameters of the KOI-94 system with an analysis of the transit timing variations (TTVs). Such constraints are independent of the radial velocity (RV) analysis recently performed by Weiss and coworkers, and valuable in examining the reliability of the parameter estimate using TTVs. We numerically fit the observed TTVs of KOI-94c, KOI-94d, and KOI-94e for their masses, eccentricities, and longitudes of periastrons, and obtain the best-fit parameters including $m_{\rm c} = 9.4_{-2.1}^{+2.4} M_{\oplus}$, $m_{\rm d} = 52.1_{-7.1}^{+6.9} M_{\oplus}$, $m_{\rm e} = 13.0_{-2.1}^{+2.5} M_{\oplus}$, and $e \lesssim 0.1$ for all the three planets. While these values are mostly in agreement with the RV result, the mass of KOI-94d estimated from the TTV is significantly smaller than the RV value $m_{\rm d} = 106 \pm 11 M_{\oplus}$. In addition, we find that the TTV of the outermost planet KOI-94e is not well reproduced in the current modeling. We also present analytic modeling of the PPE and derive a simple formula to reconstruct the mutual inclination of the two planets from the observed height, central time, and duration of the brightening caused by the PPE. Based on this model, the implication of the results of TTV analysis for the time of the next PPE is discussed.

preprint2013arXiv

Chasing Unbiased Spectra of the Universe

The cosmological power spectrum of the coherent matter flow is measured exploiting an improved prescription for the apparent anisotropic clustering pattern in redshift space. New statistical analysis is presented to provide an optimal observational platform to link the improved redshift distortion theoretical model to future real datasets. The statistical power as well as robustness of our method are tested against 60 realizations of 8 Gpc/h^3 dark matter simulation maps mocking the precision level of upcoming wide--deep surveys. We showed that we can accurately extract the velocity power spectrum up to quasi linear scales of k~0.1 h/Mpc at z = 0.35 and up to k~0.15 h/Mpc at higher redshifts within a couple of percentage precision level. Our understanding of redshift space distortion is proved to be appropriate for precision cosmology, and our statistical method will guide us to righteous path to meet the real world.

preprint2013arXiv

Constraining cosmic string parameters with curl mode of CMB lensing

We present constraints on a cosmic string network with a measurement of weak gravitational lensing from CMB temperature map. The cosmic string network between observer and last scattering surface of CMB photons generates vector and/or tensor metric perturbations, and the deflection of CMB photons by these gravitational fields has curl mode which is not produced by the scalar metric perturbations. In this paper, we use the power spectrum of curl mode obtained from Planck to constrain the string tension, Gμ, and the reconnection probability, P. In demonstrating the parameter constraints with Planck curl mode, we also measure the lensing power spectrum from the Atacama Cosmology Telescope (ACT) 2008 season data, which have better angular resolution with lower instrumental noise on a much smaller chunk of the sky. Assuming P=1, the upper bound on tension is Gμ= 6.6\times 10^{-5} with 2σ(95% C.L.), using curl mode from Planck, which is weaker than that from the small-scale temperature power spectrum. For small values of P, however, the constraint from curl mode becomes tighter compared to that from temperature power spectrum. For P\lsim 10^{-2}, we obtain the constraint on the combination of the string parameters as GμP^{-1}\leq 3.4\times 10^{-5} at more than 2σ(95% C.L.).

preprint2013arXiv

Full-sky formulae for weak lensing power spectra from total angular momentum method

We systematically derive full-sky formulae for the weak lensing power spectra generated by scalar, vector and tensor perturbations from the total angular momentum (TAM) method. Based on both the geodesic and geodesic deviation equations, we first give the gauge-invariant expressions for the deflection angle and Jacobi map as observables of the CMB lensing and cosmic shear experiments. We then apply the TAM method, originally developed in the theoretical studies of CMB, to a systematic derivation of the angular power spectra. The TAM representation, which characterizes the total angular dependence of the spatial modes projected along a line-of-sight, can carry all the information of the lensing modes generated by scalar, vector, and tensor metric perturbations. This greatly simplifies the calculation, and we present a complete set of the full-sky formulae for angular power spectra in both the E-/B-mode cosmic shear and gradient-/curl-mode lensing potential of deflection angle. Based on the formulae, we give illustrative examples of non-vanishing B-mode cosmic shear and curl-mode of deflection angle in the presence of the vector and tensor perturbations, and explicitly compute the power spectra.

preprint2013arXiv

Matter power spectrum from a Lagrangian-space regularization of perturbation theory

We present a new approach to computing the matter density power spectrum, from large linear scales to small highly nonlinear scales. Instead of explicitly computing a partial series of high-order diagrams, as in perturbative resummation schemes, we embed the standard perturbation theory within a realistic nonlinear Lagrangian-space ansatz. We also point out that an "adhesion-like" regularization of the shell-crossing regime is more realistic than a "Zel'dovich-like" behavior, where particles freely escape to infinity. This provides a "cosmic web" power spectrum with good small-scale properties that provide a good matching with a halo model on mildly nonlinear scales. We obtain a good agreement with numerical simulations on large scales, better than 3% for $k\leq 1 h$Mpc$^{-1}$, and on small scales, better than 10% for $k \leq 10 h$Mpc$^{-1}$, at $z \geq 0.35$, which improves over previous methods.

preprint2013arXiv

Precision modeling of redshift-space distortions from multi-point propagator expansion

Using a full implementation of resummed perturbation theory (PT) from a multi-point propagator expansion, we put forward new theoretical predictions for the two-point statistics of matter fluctuations in redshift space. The predictions consistently include PT corrections up to the two-loop order and are based on an improved prescription of the redshift-space distortions (RSD) that properly takes into account their non-Gaussian impact from a systematic low-k expansion. In contrast to the previous studies that partly used standard PT calculations, the present treatment is able to provide a consistent prediction for both power spectra and correlation functions. These results are compared with $N$-body simulations with which a very good agreement is found up to the quadrupole moment. The theoretical predictions for the hexadecapole moment of the power spectra are however found to significantly departs from the numerical results at low redshift. We examine this issue and found it to be likely related to an improper modeling of the RSD damping effects on which this moment shows large dependence.

preprint2012arXiv

Cosmic propagators at two-loop order

We explore the properties of two-point cosmic propagators when Perturbation Theory (PT) loop corrections are consistently taken into account. We show in particular how the interpolation scheme proposed in arXiv:1112.3895 can be explicitly used up to two-loop order introducing the notion of regular parts for the contributing terms. Extending the one-loop results, we then derive and give semi analytical forms of the two-loop contributions for both the cosmic density and velocity propagators. These results are tested against numerical simulations and shown to significantly improve upon one-loop results for redshifts above 0.5. We found however that at lower redshift two-loop order corrections are too large partly due to a strong sensitivity of those terms to the small scale modes. We show that this dependence is expected to be even larger for higher order loop corrections both from theoretical investigations and numerical tests, the latter obtained with Monte Carlo evaluations of the three-loop contributions. This makes small-scale regularization schemes necessary for the use of such higher order corrections.

preprint2012arXiv

Cosmology with space-based gravitational-wave detectors --- dark energy and primordial gravitational waves ---

Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect ~10^6 neutron-star (NS) binaries and determine the luminosity distances to the binaries with high precision. Combining the luminosity distances with cosmologically-induced phase corrections on the GWs, cosmological expansion out to high redshift can be measured without the redshift determinations of host galaxies by electromagnetic observation and be a unique probe for dark energy. On the other hand, such a NS-binary foreground should be subtracted to detect primordial GWs produced during inflation. Thus, the constraining power on dark energy and the detectability of the primordial gravitational waves strongly depend on the detector sensitivity and are in close relation with one another. In this paper, we investigate the constraints on the equation of state of dark energy with future space-based GW detectors with/without identifying the redshifts of host galaxies. We also study the sensitivity to the primordial GWs, properly dealing with the residual of the NS binary foreground. Based on the results, we discuss the detector sensitivity required to achieve the forementioned targeted study of cosmology.

preprint2012arXiv

Gravitational-wave standard siren without redshift identification

Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect $\sim10^6$ neutron-star (NS) binaries and determine the luminosity distances to the binaries with high precision. Combining the luminosity distances with cosmologically-induced phase corrections on the GWs, cosmological expansion out to high redshift can be measured without the redshift determinations of host galaxies by electromagnetic observation and can be a unique probe for dark energy. This article is based on the results obtained in [1] where we investigated constraining power of the GW standard siren without redshift information on the equation of state of dark energy with future space-based GW detectors. We also compare the results with those obtained with other instruments and methods.

preprint2012arXiv

Measurements of Stellar Inclinations for Kepler Planet Candidates

We present an investigation of spin-orbit angles for planetary system candidates reported by Kepler. By combining the rotational period $P_s$ inferred from the flux variation due to starspots and the projected rotational velocity $V\sin I_s$ and stellar radius obtained by a high resolution spectroscopy, we attempt to estimate the inclination $I_s$ of the stellar spin axis with respect to the line-of-sight. For transiting planetary systems, in which planetary orbits are edge-on seen from us, the stellar inclination $I_s$ can be a useful indicator of a spin-orbit alignment/misalignment. We newly conducted spectroscopic observations with Subaru/HDS for 15 KOI systems, whose lightcurves show periodic flux variations. After detailed analyses of their lightcurves and spectra, it turned out that some of them are binaries, or the flux variations are too coherent to be caused by starspots, probably representing ellipsoidal variations, and consequently we could constrain stellar inclinations $I_s$ for eight systems. Among them, KOI-262 and 280 are in good agreement with $I_s=90^\circ$ suggesting a spin-orbit alignment, while at least one system, KOI-261, shows a possible spin-orbit misalignment. We also obtain a small $I_s$ for KOI-1463, but the transiting companion seems to be a star rather than a planet. The results for KOI-257, 269, 367, and 974 are ambiguous, and can be explained with either misalignments or moderate differential rotation. Since our method can be applied to any system having starspots regardless of the planet size, future observations will allow for the expansion of the parameter space in which the spin-orbit relations are investigated.

preprint2012arXiv

RegPT: Direct and fast calculation of regularized cosmological power spectrum at two-loop order

We present a specific prescription for the calculation of cosmological power spectra, exploited here at two-loop order in perturbation theory (PT), based on the multi-point propagator expansion. In this approach power spectra are constructed from the regularized expressions of the propagators that reproduce both the resummed behavior in the high-k limit and the standard PT results at low-k. With the help of N-body simulations, we show that such a construction gives robust and accurate predictions for both the density power spectrum and the correlation function at percent-level in the weakly non-linear regime. We then present an algorithm that allows accelerated evaluations of all the required diagrams by reducing the computational tasks to one-dimensional integrals. This is achieved by means of pre-computed kernel sets defined for appropriately chosen fiducial models. The computational time for two-loop results is then reduced from a few minutes, with the direct method, to a few seconds with the fast one. The robustness and applicability of this method are tested against the power spectrum cosmic emulator from which a wide variety of cosmological models can be explored. The fortran program with which direct and fast calculations of power spectra can be done, RegPT, is publicly released as part of this paper.

preprint2012arXiv

Revising the Halofit Model for the Nonlinear Matter Power Spectrum

Based on a suite of state-of-the-art high-resolution $N$-body simulations, we revisit the so-called halofit model (Smith et al. 2003) as an accurate fitting formula for the nonlinear matter power spectrum. While the halofit model has been frequently used as a standard cosmological tool to predict the nonlinear matter power spectrum in a universe dominated by cold dark matter, its precision has been limited by the low-resolution of $N$-body simulations used to determine the fitting parameters, suggesting the necessity of improved fitting formula at small scales for future cosmological studies. We run high-resolution $N$-body simulations for 16 cosmological models around the Wilkinson Microwave Anisotropy Probe (WMAP) best-fit cosmological parameters (1, 3, 5, and 7 year results), including dark energy models with a constant equation of state. The simulation results are used to re-calibrate the fitting parameters of the halofit model so as to reproduce small-scale power spectra of the $N$-body simulations, while keeping the precision at large scales. The revised fitting formula provides an accurate prediction of the nonlinear matter power spectrum in a wide range of wavenumber ($k \leq 30h$\,Mpc$^{-1}$) at redshifts $0 \leq z \leq 10$, with 5% precision for $k\leq1 h$ Mpc$^{-1}$ at $0 \leq z \leq 10$ and 10% for $1 \leq k\leq 10 h$ Mpc$^{-1} $ at $0 \leq z \leq 3$. We discuss the impact of the improved halofit model on weak lensing power spectra and correlation functions, and show that the improved model better reproduces ray-tracing simulation results.

preprint2012arXiv

Weak lensing generated by vector perturbations and detectability of cosmic strings

We study the observational signature of vector metric perturbations through the effect of weak gravitational lensing. In the presence of vector perturbations, the non-vanishing signals for B-mode cosmic shear and curl-mode deflection angle, which have never appeared in the case of scalar metric perturbations, naturally arise. Solving the geodesic and geodesic deviation equations, we drive the full-sky formulas for angular power spectra of weak lensing signals, and give the explicit expressions for E-/B-mode cosmic shear and gradient-/curl-mode deflection angle. As a possible source for seeding vector perturbations, we then consider a cosmic string network, and discuss its detectability from upcoming weak lensing and CMB measurements. Based on the formulas and a simple model for cosmic string network, we calculate the angular power spectra and expected signal-to-noise ratios for the B-mode cosmic shear and curl-mode deflection angle. We find that the weak lensing signals are enhanced for a smaller intercommuting probability of the string network, $P$, and they are potentially detectable from the upcoming cosmic shear and CMB lensing observations. For $P\sim 10^{-1}$, the minimum detectable tension of the cosmic string will be down to $Gμ\sim 5\times 10^{-8}$. With a theoretically inferred smallest value $P\sim 10^{-3}$, we could even detect the string with $Gμ\sim 5\times 10^{-10}$.

preprint2011arXiv

Baryon Acoustic Oscillations in 2D II: Redshift-space halo clustering in N-body simulations

We measure the halo power spectrum in redshift space from cosmological N-body simulations, and test the analytical models of redshift distortions particularly focusing on the scales of baryon acoustic oscillations (BAOs). Remarkably, the measured halo power spectrum in redshift space exhibits a large-scale enhancement in amplitude relative to the real-space clustering, and the effect becomes significant for the massive or highly biased halo samples. These findings cannot be simply explained by the so-called streaming model frequently used in the literature. By contrast, a physically-motivated perturbation theory model developed in the previous paper reproduces the halo power spectrum very well, and the model combining a simple linear scale-dependent bias can accurately characterize the clustering anisotropies of halos in two dimensions, i.e., line-of-sight and its perpendicular directions. The results highlight the significance of non-linear coupling between density and velocity fields associated with two competing effects of redshift distortions, i.e., Kaiser and Finger-of-God effects, and a proper account of this effect would be important in accurately characterizing the BAOs in two dimensions.

preprint2011arXiv

Forecasting the Cosmological Constraints with Anisotropic Baryon Acoustic Oscillations from Multipole Expansion

Baryon acoustic oscillations (BAOs) imprinted in the galaxy power spectrum can be used as a standard ruler to determine angular diameter distance and Hubble parameter at high redshift galaxies. Combining redshift distortion effect which apparently distorts the galaxy clustering pattern, we can also constrain the growth rate of large-scale structure formation. Usually, future forecast for constraining these parameters from galaxy redshift surveys has been made with a full 2D power spectrum characterized as function of wavenumber $k$ and directional cosine $μ$ between line-of-sight direction and wave vector, i.e., $P(k,μ)$. Here, we apply the multipole expansion to the full 2D power spectrum, and discuss how much cosmological information can be extracted from the lower-multipole spectra, taking a proper account of the non-linear effects on gravitational clustering and redshift distortion. The Fisher matrix analysis reveals that compared to the analysis with full 2D spectrum, a partial information from the monopole and quadrupole spectra generally degrades the constraints by a factor of $\sim1.3$ for each parameter. The additional information from the hexadecapole spectrum helps to improve the constraints, which lead to an almost comparable result expected from the full 2D spectrum.

preprint2011arXiv

Further Observations of the Tilted Planet XO-3: A New Determination of Spin-Orbit Misalignment, and Limits on Differential Rotation

We report on observations of the Rossiter-McLaughlin (RM) effect for the XO-3 exoplanetary system. The RM effect for the system was previously measured by two different groups, but their results were statistically inconsistent. To obtain a decisive result we observed two full transits of XO-3b with the Subaru 8.2-m telescope. By modeling these data with a new and more accurate analytic formula for the RM effect, we find the projected spin-orbit angle to be λ=37.3 deg \pm 3.0 deg, in good agreement with the previous finding by Winn et al. (2009). In addition, an offset of ~22 m/s was observed between the two transit datasets. This offset could be a signal of a third body in the XO-3 system, a possibility that should be checked with future observations. We also attempt to search for a possible signature of the stellar differential rotation in the RM data for the first time, and put weak upper limits on the differential rotation parameters.

preprint2011arXiv

Improved Modeling of the Rossiter-McLaughlin Effect for Transiting Exoplanets

We present an improved formula for the anomalous radial velocity of the star during planetary transits due to the Rossiter-McLaughlin (RM) effect. The improvement comes from a more realistic description of the stellar absorption line profiles, taking into account stellar rotation, macroturbulence, thermal broadening, pressure broadening, and instrumental broadening. Although the formula is derived for the case in which radial velocities are measured by cross-correlation, we show through numerical simulations that the formula accurately describes the cases where the radial velocities are measured with the iodine absorption-cell technique. The formula relies on prior knowledge of the parameters describing macroturbulence, instrumental broadening and other broadening mechanisms, but even 30% errors in those parameters do not significantly change the results in typical circumstances. We show that the new analytic formula agrees with previous ones that had been computed on a case-by-case basis via numerical simulations. Finally, as one application of the new formula, we reassess the impact of the differential rotation on the RM velocity anomaly. We show that differential rotation of a rapidly rotating star may have a significant impact on future RM observations.

preprint2011arXiv

Magnification effect on the detection of primordial non-Gaussianity from photometric surveys

We present forecast results for constraining the primordial non-Gaussianity from photometric surveys through a large-scale enhancement of the galaxy clustering amplitude. In photometric surveys, the distribution of observed galaxies at high redshifts suffers from the gravitational-lensing magnification, which systematically alters the number density for magnitude-limited galaxy samples. We estimate size of the systematic bias in the best-fit cosmological parameters caused by the magnification effect, particularly focusing on the primordial non-Gaussianity. For upcoming deep and/or wide photometric surveys like HSC, DES and LSST, the best-fit value of the non-Gaussian parameter, fNL, obtained from the galaxy count data is highly biased, and the true values of fNL would typically go outside the 3-sigma error of the biased confidence region, if we ignore the magnification effect in the theoretical template of angular power spectrum. The additional information from cosmic shear data helps not only to improve the constraint, but also to reduce the systematic bias. As a result, the size of systematic bias on fNL would become small enough compared to the expected 1-sigma error for HSC and DES, but it would be still serious for deep surveys with z_m > 1.5, like LSST. Tomographic technique improves the constraint on fNL by a factor of 2-3 compared to the one without tomography, but the systematic bias would increase.

preprint2011arXiv

Neutrino mass constraint with the Sloan Digital Sky Survey power spectrum of luminous red galaxies and perturbation theory

We compare the model power spectrum, computed based on perturbation theory (PT) with the power spectrum of luminous red galaxies (LRG) measured from the SDSSDR7 catalog, assuming a flat, CDM-dominated cosmology. The model includes the effects of massive neutrinos, nonlinear matter clustering and nonlinear, scale-dependent galaxy bias in a self-consistent manner. We first test the accuracy of PT-model by comparing the model predictions with the halo power spectrum in real- and redshift-space measured from simulations without massive neutrinos. We show that the PT-model with bias parameters being properly adjusted can fairly well reproduce the simulation results. As a result the best-fit parameters obtained from the hypothetical parameter fitting recover, within statistical uncertainties, the input cosmological parameters in simulations, including an upper bound on neutrino mass, if the power spectrum information up to k~0.15h/Mpc is used. However, for the redshift-space power spectrum, the best-fit cosmological parameters show a sizable bias from the input values if using the information up to k~0.2h/Mpc, probably due to nonlinear redshift distortion effect. Given these tests, we decided, as a conservative choice, to use the LRG power spectrum up to k=0.1h/Mpc in order to minimize possible unknown nonlinearity effects. In combination with the recent results from Wilkinson Microwave Background Anisotropy Probe (WMAP), we derive a robust upper-bound on the sum of neutrino masses, given as m_nu,tot < 0.81eV (95% C.L.), marginalized over other parameters including nonlinear bias parameters and dark energy equation of state parameter. The neutrino mass limit is improved by a factor of 1.85 compared to the limit from the WMAP5 alone, m_nu,tot < 1.5eV.

preprint2011arXiv

Next-to-leading resummation of cosmological perturbations via the Lagrangian picture: 2-loop correction in real and redshift spaces

We present an improved prediction of the nonlinear perturbation theory (PT) via the Lagrangian picture, which was originally proposed by Matsubara (2008). Based on the relations between the power spectrum in standard PT and that in Lagrangian PT, we derive analytic expressions for the power spectrum in Lagrangian PT up to 2-loop order in both real and redshift spaces. Comparing the improved prediction of Lagrangian PT with $N$-body simulations in real space, we find that the 2-loop corrections can extend the valid range of wave numbers where we can predict the power spectrum within 1% accuracy by a factor of 1.0 ($z=0.5$), 1.3 (1), 1.6 (2) and 1.8 (3) vied with 1-loop Lagrangian PT results. On the other hand, in all redshift ranges, the higher-order corrections are shown to be less significant on the two-point correlation functions around the baryon acoustic peak, because the 1-loop Lagrangian PT is already accurate enough to explain the nonlinearity on those scales in $N$-body simulations.

preprint2011arXiv

Tracing the redshift evolution of Hubble parameter with gravitational-wave standard sirens

Proposed space-based gravitational-wave detectors such as BBO and DECIGO can detect ~10^6 neutron-star binaries and determine luminosity distance to the binaries with a high precision. Combining the luminosity distance and electromagnetically-derived redshift, one would be able to probe cosmological expansion out to high redshift. In this paper, we show that the Hubble parameter as a function of redshift can be directly measured with monopole and dipole components of the luminosity distance on the sky. As a result, the measurement accuracies of the Hubble parameter in each redshift bin up to z=1 are 3-14 %, 1.5-8 %, and 0.8-4% for the observation time 1 yr, 3 yr, and 10 yr, respectively.

preprint2010arXiv

Baryon Acoustic Oscillations in 2D: Modeling Redshift-space Power Spectrum from Perturbation Theory

We present an improved prescription for matter power spectrum in redshift space taking a proper account of both the non-linear gravitational clustering and redshift distortion, which are of particular importance for accurately modeling baryon acoustic oscillations (BAOs). Contrary to the models of redshift distortion phenomenologically introduced but frequently used in the literature, the new model includes the corrections arising from the non-linear coupling between the density and velocity fields associated with two competitive effects of redshift distortion, i.e., Kaiser and Finger-of-God effects. Based on the improved treatment of perturbation theory for gravitational clustering, we compare our model predictions with monopole and quadrupole power spectra of N-body simulations, and an excellent agreement is achieved over the scales of BAOs. Potential impacts on constraining dark energy and modified gravity from the redshift-space power spectrum are also investigated based on the Fisher-matrix formalism. We find that the existing phenomenological models of redshift distortion produce a systematic error on measurements of the angular diameter distance and Hubble parameter by 1~2%, and the growth rate parameter by ~5%, which would become non-negligible for future galaxy surveys. Correctly modeling redshift distortion is thus essential, and the new prescription of redshift-space power spectrum including the non-linear corrections can be used as an accurate theoretical template for anisotropic BAOs.

preprint2010arXiv

Cosmological test of gravity with polarizations of stochastic gravitational waves around 0.1-1 Hz

In general relativity, a gravitational wave has two polarization modes (tensor mode), but it could have additional polarizations (scalar and vector modes) in the early stage of the universe, where the general relativity may not strictly hold and/or the effect of higher-dimensional gravity may become significant. In this paper, we discuss how to detect extra-polarization modes of stochastic gravitational wave background (GWB), and study the separability of each polarization using future space-based detectors such as BBO and DECIGO. We specifically consider two plausible setups of the spacecraft constellations consisting of two and four clusters, and estimate the sensitivity to each polarization mode of GWBs. We find that a separate detection of each polarization mode is rather sensitive to the geometric configuration and distance between clusters and that the clusters should be, in general, separated by an appropriate distance. This seriously degrades the signal sensitivity, however, for suitable conditions, space-based detector can separately detect scalar, vector and tensor modes of GWBs with energy density as low as ~10^-15.

preprint2010arXiv

Non-Gaussian Error Contribution to Likelihood Analysis of the Matter Power Spectrum

We study the sample variance of the matter power spectrum for the standard Lambda Cold Dark Matter universe. We use a total of 5000 cosmological N-body simulations to study in detail the distribution of best-fit cosmological parameters and the baryon acoustic peak positions. The obtained distribution is compared with the results from the Fisher matrix analysis with and without including non-Gaussian errors. For the Fisher matrix analysis, we compute the derivatives of the matter power spectrum with respect to cosmological parameters using directly full nonlinear simulations. We show that the non-Gaussian errors increase the unmarginalized errors by up to a factor 5 for k_{max}=0.4h/Mpc if there is only one free parameter provided other parameters are well determined by external information. On the other hand, for multi-parameter fitting, the impact of the non-Gaussian errors is significantly mitigated due to severe parameter degeneracies in the power spectrum. The distribution of the acoustic peak positions is well described by a Gaussian distribution, with its width being consistent with the statistical interval predicted from the Fisher matrix. We also examine systematic bias in the best-fit parameter due to the non-Gaussian errors. The bias is found to be smaller than the 1 sigma statistical error for both the cosmological parameters and the acoustic scale positions.

preprint2010arXiv

Scale Dependence of Halo Bispectrum from Non-Gaussian Initial Conditions in Cosmological N-body Simulations

We study the halo bispectrum from non-Gaussian initial conditions. Based on a set of large $N$-body simulations starting from initial density fields with local type non-Gaussianity, we find that the halo bispectrum exhibits a strong dependence on the shape and scale of Fourier space triangles near squeezed configurations at large scales. The amplitude of the halo bispectrum roughly scales as $f_nl^2$. The resultant scaling on the triangular shape is consistent with that predicted by Jeong & Komatsu based on perturbation theory. We systematically investigate this dependence with varying redshifts and halo mass thresholds. It is shown that the $f_nl$ dependence of the halo bispectrum is stronger for more massive haloes at higher redshifts. This feature can be a useful discriminator of inflation scenarios in future deep and wide galaxy redshift surveys.

preprint2009arXiv

Analytic Description of the Rossiter-McLaughlin Effect for Transiting Exoplanets: Cross-Correlation Method and Comparison with Simulated Data

We obtain analytical expressions for the velocity anomaly due to the Rossiter- McLaughlin effect, for the case when the anomalous radial velocity is obtained by cross-correlation with a stellar template spectrum. In the limit of vanishing width of the stellar absorption lines, our result reduces to the formula derived by Ohta et al. (2005), which is based on the first moment of distorted stellar lines. Our new formula contains a term dependent on the stellar linewidth, which becomes important when rotational line broadening is appreciable. We generate mock transit spectra for four existing exoplanetary systems (HD17156, TrES-2, TrES- 4, and HD209458) following the procedure of Winn et al. (2005), and find that the new formula is in better agreement with the velocity anomaly extracted from the mock data. Thus, our result provides a more reliable analytical description of the velocity anomaly due to the Rossiter-McLaughlin effect, and explains the previously observed dependence of the velocity anomaly on the stellar rotation velocity.

preprint2009arXiv

Non-linear Evolution of Baryon Acoustic Oscillations from Improved Perturbation Theory in Real and Redshift Spaces

We study the non-linear evolution of baryon acoustic oscillations in the matter power spectrum and correlation function from the improved perturbation theory (PT). Based on the framework of renormalized PT, we apply the {\it closure approximation} that truncates the infinite series of loop contributions at one-loop order, and obtain a closed set of integral equations for power spectrum and non-linear propagator. The resultant integral expressions keep important non-perturbative properties which can dramatically improve the prediction of non-linear power spectrum. Employing the Born approximation, we then derive the analytic expressions for non-linear power spectrum and the predictions are made for non-linear evolution of baryon acoustic oscillations in power spectrum and correlation function. A detailed comparison between improved PT results and N-body simulations shows that a percent-level agreement is achieved in a certain range in power spectrum and in a rather wider range in correlation function. Combining a model of non-linear redshift-space distortion, we also evaluate the power spectrum and correlation function in correlation function. In contrast to the results in real space, the agreement between N-body simulations and improved PT predictions tends to be worse, and a more elaborate modeling for redshift-space distortion needs to be developed. Nevertheless, with currently existing model, we find that the prediction of correlation function has a sufficient accuracy compared with the cosmic-variance errors for future galaxy surveys with volume of a few (Gpc/h)^3 at z>=0.5.

preprint2008arXiv

Modeling Nonlinear Evolution of Baryon Acoustic Oscillations: Convergence Regime of N-body Simulations and Analytic Models

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 Λ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.

preprint2007arXiv

Characteristic Scales of Baryon Acoustic Oscillations from Perturbation Theory: Non-linearity and Redshift-Space Distortion Effects

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.

preprint2007arXiv

Measurement of the Rossiter--McLaughlin Effect in the Transiting Exoplanetary System TrES-1

We report a measurement of the Rossiter--McLaughlin effect in the transiting extrasolar planetary system TrES-1, via simultaneous spectroscopic and photometric observations with the Subaru and MAGNUM telescopes. By modeling the radial velocity anomaly that was observed during a transit, we determine the sky-projected angle between the stellar spin axis and the planetary orbital axis to be $λ= 30 \pm 21$ [deg]. This is the third case for which $λ$ has been measured in a transiting exoplanetary system, and the first demonstration that such measurements are possible for relatively faint host stars ($V \sim 12$, as compared to $V \sim 8$ for the other systems). We also derive a time of mid-transit, constraints on the eccentricity of the TrES-1b orbit ($e = 0.048 \pm 0.025$), and upper limits on the mass of the Trojan companions ($\lesssim$14 $M_{\oplus}$) at the 3$σ$ level.

preprint2006arXiv

Bispectrum and Nonlinear Biasing of Galaxies: Perturbation Analysis, Numerical Simulation and SDSS Galaxy Clustering

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.

preprint2005arXiv

Measurement of Spin-Orbit Alignment in an Extrasolar Planetary System

We determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms.

preprint2005arXiv

Minkowski Functionals of SDSS galaxies I : Analysis of Excursion Sets

We present a first morphometric investigation of a preliminary sample from the SDSS of 154287 galaxies with apparent magnitude 14.5<m_r<17.5 and redshift 0.001<z<0.4. We measure the Minkowski Functionals, which are a complete set of morphological descriptors. To account for the complicated wedge--like geometry of the present survey data, we construct isodensity contour surfaces from the galaxy positions in redshift space and employ two complementary methods of computing the Minkowski Functionals. We find that the observed Minkowski Functionals for SDSS galaxies are consistent with the prediction of a Lambda--dominated spatially--flat Cold Dark Matter model with random--Gaussian initial conditions, within the cosmic variance estimated from the corresponding mock catalogue. We expect that future releases of the SDSS survey will allow us to distinguish morphological differences in the galaxy distribution with regard to different morphological type and luminosity ranges.

preprint2003arXiv

Biasing and Genus Statistics of Dark Matter Halos in the Hubble Volume Simulation

We present a numerical analysis of genus statistics for dark matter halo catalogs from the Hubble volume simulation. The huge box-size of the Hubble volume simulation enables us to carry out a reliable statistical analysis of the biasing properties of halos at a Gaussian smoothing scale of R_G>30Mpc/h with a cluster-mass scale of between 7*10^{13}Msolar/h and 6*10^{15}Msolar/h. A detailed comparison of the genus for dark matter halos with that for the mass distribution shows that the non-Gaussianity induced by the halo biasing is comparable to that by nonlinear gravitational evolution, and both the shape and the amplitude of the genus are almost insensitive to the halo mass at R_G>30Mpc/h. In order to characterize the biasing effect on the genus, we apply a perturbative formula developed by Matsubara (1994). We find that the perturbative formula well describes the simulated halo genus at R_G>50Mpc/h. The result indicates that the biasing effect on the halo genus is well approximated by nonlinear deterministic biasing up to the second-order term in the mass density fluctuation. The two parameters describing the linear and quadratic terms in the nonlinear biasing accurately specify the genus for galaxy clusters.

preprint2003arXiv

Evolution of gravitational waves from inflationary brane-world : numerical study of high-energy effects

We study the evolution of gravitational waves(GWs) after inflation in a brane-world cosmology embedded in five-dimensional anti-de Sitter spacetime. Contrary to the standard four-dimensional results, the GWs at the high-energy regime in brane-world model suffer from the effects of the non-standard cosmological expansion and the excitation of the Kaluza-Klein modes(KK-modes), which can affect the amplitude of stochastic gravitational wave background significantly. To investigate these two high-energy effects quantitatively, we numerically solve the wave equation of the GWs in the radiation dominated epoch at relatively low-energy scales. We show that the resultant GWs is suppressed by the excitation of the KK modes. The created KK modes are rather soft and escape away from the brane to the bulk gravitational field. The results are also compared to the semi-analytic prediction from the low-energy approximation and the evolved amplitude of GWs on the brane reasonably matches the numerical simulations.

preprint2002arXiv

Modeling Pairwise Peculiar Velocity Distribution Function of Dark Matter from Halo Density Profiles

We derive the pairwise peculiar velocity distribution function of dark matter particles applying the dark matter halo approach. Unlike the previous work, we do not assume a Gaussian velocity distribution function of dark matter in a single halo, but compute it self-consistently with the assumed density profile for dark matter halo. The resulting distribution function is well approximated by an exponential distribution which is consistent with the previous observational, numerical and theoretical results. We also compute the pairwise peculiar velocity dispersion for different density profiles, and provide a practical fitting formula. We apply an empirical biasing scheme into our model and present prediction for pairwise peculiar velocity dispersion of {\it galaxies}, and reproduce the previous results of simulations using our semi-analytical method.