Source author record

Ryuichi Takahashi

Ryuichi Takahashi 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

26works
6topics
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

26 published item(s)

preprint2022arXiv

Amplitude and phase fluctuations of gravitational waves magnified by strong gravitational lensing

We discuss how small-scale density perturbations on the Fresnel scale affect amplitudes and phases of gravitational waves that are magnified by gravitational lensing in geometric optics. We derive equations that connect the small-scale density perturbations with the amplitude and phase fluctuations to show that such perturbative wave optics effects are significantly boosted in the presence of macro model magnifications such that amplitude and phase fluctuations can easily be observed for highly magnified gravitational waves. We discuss expected signals due to microlensing by stars, dark matter substructure, fuzzy dark matter, and primordial black holes.

preprint2022arXiv

Probing cosmology and gastrophysics with fast radio bursts: Cross-correlations of dark matter haloes and cosmic dispersion measures

For future surveys of fast radio bursts (FRBs), we clarify information available from cosmic dispersion measures (DMs) through cross-correlation analyses of foreground dark matter haloes (hosting galaxies and galaxy clusters) with their known redshifts. With a halo-model approach, we predict that the cross-correlation with cluster-sized haloes is less affected by the details of gastrophysics, providing robust cosmological information. For less massive haloes, the cross-correlation at angular scales of $<10\, \mathrm{arcmin}$ is sensitive to gas expelled from the halo centre due to galactic feedback. Assuming $20000$ FRBs over $20000 \, {\rm deg}^2$ with a localisation error being 3 arcmin, we expect that the cross-correlation signal at halo masses of $10^{12}-10^{14}\, M_\odot$ can be measured with a level of $\sim 1\%$ precision in a redshift range of $0<z<1$. Such precise measurements enable one to put a 1.5\% level constraint on $σ_8\, (Ω_\mathrm{M}/0.3)^{0.5}$ and a 3\% level constraint on $(Ω_\mathrm{b}/0.049)(h/0.67)(f_\mathrm{e}/0.95)$ ($σ_8$, $Ω_\mathrm{M}$, $Ω_\mathrm{b}$, $h$ and $f_\mathrm{e}$ are the linear mass variance smoothed at $8\, h^{-1}\mathrm{Mpc}$, mean mass density, mean baryon density, the present-day Hubble parameter and fraction of free electrons in cosmic baryons today), whereas the gas-to-halo mass relation in galaxies and clusters can be constrained with a level of $10\%-20\%$. Furthermore the cross-correlation analyses can break the degeneracy among $Ω_\mathrm{b}$, $h$ and $f_\mathrm{e}$, inherent in the DM-redshift relation. Our proposal opens new possibilities for FRB cosmology, while it requires extensive galaxy redshift catalogues and further improvement of the halo model.

preprint2021arXiv

Constraining Primordial Non-Gaussianity with Post-reconstructed Galaxy Bispectrum in Redshift Space

Galaxy bispectrum is a promising probe of inflationary physics in the early universe as a measure of primordial non-Gaussianity (PNG), whereas its signal-to-noise ratio is significantly affected by the mode coupling due to non-linear gravitational growth. In this paper, we examine the standard reconstruction method of linear cosmic mass density fields from non-linear galaxy density fields to de-correlate the covariance in redshift-space galaxy bispectra. In particular, we evaluate the covariance of the bispectrum for massive-galaxy-sized dark matter halos with reconstruction by using 4000 independent $N$-body simulations. Our results show that the bispectrum covariance for the post-reconstructed field approaches the Gaussian prediction at scale of $k<0.2\, h\, {\rm Mpc}^{-1}$. We also verify the leading-order PNG-induced bispectrum is not affected by details of the reconstruction with perturbative theory. We then demonstrate the constraining power of the post-reconstructed bispectrum for PNG at redshift of $\sim0.5$. Further, we perform a Fisher analysis to make a forecast of PNG constraints by galaxy bispectra including anisotropic signals. Assuming a massive galaxy sample in the SDSS Baryon Oscillation Spectroscopic Survey, we find that the post-reconstructed bispectrum can constrain the local-, equilateral- and orthogonal-types of PNG with $Δf_{\rm NL} \sim$13, 90 and 42, respectively, improving the constraints with the pre-reconstructed bispectrum by a factor of $1.3-3.2$. In conclusion, the reconstruction plays an essential role in constraining various types of PNG signatures with a level of $Δf_{\rm NL}<1$ from the galaxy bispectrum based on upcoming galaxy surveys.

preprint2021arXiv

HSC Year 1 cosmology results with the minimal bias method: HSC$\times$BOSS galaxy-galaxy weak lensing and BOSS galaxy clustering

We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three redshift bins, and use the HSC-Y1 galaxy catalog to define a secure sample of source galaxies at $z_\mathrm{ph}>0.75$ for the $Δ\!Σ$ measurements, selected based on their photometric redshifts. For theoretical template, we use the "minimal bias" model for the cosmological clustering observables for the flat $Λ$CDM cosmological model. We compare the model predictions with the measurements in each redshift bin on large scales, $R>12$ and $8~h^{-1}\mathrm{Mpc}$ for $Δ\!Σ(R)$ and $w_\mathrm{p}(R)$, respectively, where the perturbation theory-inspired model is valid. When we employ weak priors on cosmological parameters, without CMB information, we find $S_8=0.936^{+0.092}_{-0.086}$, $σ_8=0.85^{+0.16}_{-0.11}$, and $Ω_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $Λ$CDM model. Although the central value of $S_8$ appears to be larger than those inferred from other cosmological experiments, we find that the difference is consistent with expected differences due to sample variance, and our results are consistent with the other results to within the statistical uncertainties. (abriged)

preprint2021arXiv

Statistical modelling of the cosmological dispersion measure

We have investigated the basic statistics of the cosmological dispersion measure (DM) -- such as its mean, variance, probability distribution, angular power spectrum and correlation function -- using the state-of-the-art hydrodynamic simulations, IllustrisTNG300, for the fast radio burst (FRB) cosmology. To model the DM statistics, we first measured the free-electron abundance and the power spectrum of its spatial fluctuations. The free-electron power spectrum turns out to be consistent with the dark matter power spectrum at large scales, but it is strongly damped at small scales ($\lesssim 1$Mpc) owing to the stellar and active galactic nucleus feedback. The free-electron power spectrum is well modelled using a scale-dependent bias factor (the ratio of its fluctuation amplitude to that of the dark matter). We provide analytical fitting functions for the free-electron abundance and its bias factor. We next constructed mock sky maps of the DM by performing standard ray-tracing simulations with the TNG300 data. The DM statistics are calculated analytically from the fitting functions of the free-electron distribution, which agree well with the simulation results measured from the mock maps. We have also obtained the probability distribution of source redshift for a given DM, which helps in identifying the host galaxies of FRBs from the measured DMs. The angular two-point correlation function of the DM is described by a simple power law, $ξ(θ) \approx 2400 (θ/{\rm deg})^{-1} \, {\rm pc}^2 \, {\rm cm}^{-6}$, which we anticipate will be confirmed by future observations when thousands of FRBs are available.

preprint2020arXiv

Cosmological information content in redshift-space power spectrum of SDSS-like galaxies in the quasi-nonlinear regime up to $k=0.3\,h\,$Mpc$^{-1}$

Clustering properties and peculiar velocities of halos in large-scale structure carry a wealth of cosmological information over a wide range of scales from linear to nonlinear scales. We use halo catalogs in a suite of high-resolution $N$-body simulations to construct mock catalogs of galaxies that resemble the SDSS-like luminous early-type galaxies at three redshift bins in the range $0.15\le z\le 0.7$. To do this we include 10 nuisance parameters to model variations in halo-galaxy connections for each redshift bin; the halo occupation distribution, and the spatial and velocity distributions of galaxies in the host halos. We evaluate the Fisher information matrix for the redshift-space power spectrum of SDSS-like galaxies using different sets of the mock catalogs that are generated from changes in each of model parameters; cosmological parameters ($σ_8$ and $Ω_{\rm m}$), the halo-galaxy connection parameters, and the cosmological distances ($D_{\rm A}$ and $H$ parameters at each redshift bin) for modeling an apparent geometrical distortion of the redshift-space power spectrum (the Alcock-Paczynski effect). We show that combining the monopole and quadrupole power spectra of galaxies allows for accurate estimations of the cosmological parameters and the cosmological distances, even after marginalization over the halo-galaxy parameters, by lifting the parameter degeneracies that are otherwise inevitable if either of the two spectra alone is used. When including the galaxy power spectrum information up to $k=0.3\,h\,$Mpc$^{-1}$, we find about factor of 6 gain in the cosmological information content of ($σ_8,Ω_{\rm m},D_{\rm A}$'s and $H$'s) compared to $k=0.2\,h\,$Mpc$^{-1}$.

preprint2020arXiv

Cross-correlation of the thermal Sunyaev-Zel'dovich effect and weak gravitational lensing: Planck and Subaru Hyper Suprime-Cam first-year data

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

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

Perturbation theory for the redshift-space matter power spectra after reconstruction

We derive the one-loop perturbative formula of the redshift-space matter power spectrum after density field reconstruction in the Zeldovich approximation. We find that the reconstruction reduces the amplitudes of nonlinear one-loop perturbative terms significantly by partially erasing the nonlinear mode-coupling between density and velocity fields. In comparison with N-body simulations, we find that both the monopole and quadrupole spectra of reconstructed matter density fields agree with the one-loop perturbation theory up to higher wavenumber than those before reconstruction. We also evaluate the impact on cosmic growth rate assuming the survey volume and the number density like the Baryon Oscillation Spectroscopic Survey and find that the total error, including statistical and systematic ones due to one-loop approximation, decreases by half.

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

Constraints on small-scale cosmological fluctuations from SNe lensing dispersion

We provide predictions on small-scale cosmological density power spectrum from supernova lensing dispersion. Parameterizing the primordial power spectrum with running $α$ and running of running $β$ of the spectral index, we exclude large positive $α$ and $β$ parameters which induce too large lensing dispersions over current observational upper bound. We ran cosmological N-body simulations of collisionless dark matter particles to investigate non-linear evolution of the primordial power spectrum with positive running parameters. The initial small-scale enhancement of the power spectrum is largely erased when entering into the non-linear regime. For example, even if the linear power spectrum at $k>10h {\rm Mpc}^{-1}$ is enhanced by $1-2$ orders of magnitude, the enhancement much decreases to a factor of $2-3$ at late time ($z \leq 1.5$). Therefore, the lensing dispersion induced by the dark matter fluctuations weakly constrains the running parameters. When including baryon-cooling effects (which strongly enhance the small-scale clustering), the constraint is comparable or tighter than the PLANCK constraint, depending on the UV cut-off. Further investigations of the non-linear matter spectrum with baryonic processes is needed to reach a firm constraint.

preprint2015arXiv

Constraints on warm dark matter from weak lensing in anomalous quadruple lenses

We investigate the weak lensing effect by line-of-sight structures with a surface mass density of <~10^8 solar mass/arcsec^2 in QSO-galaxy quadruple lens systems. Using high-resolution N-body simulations in warm dark matter (WDM) models and observed four quadruple lenses that show anomalies in the flux ratios, we obtain constraints on the mass of thermal WDM, m_WDM>= 1.3keV(95%CL) assuming that the density of the primary lens is described by a singular isothermal ellipsoid (SIE). The obtained constraint is consistent with those from Lyman-$α$ forests and the number counts of high-redshift galaxies at z>4. Our results show that WDM with a free-streaming comoving wavenumber k_{fs} <= 27 h/Mpc is disfavored as the major component of cosmological density at redshifts 0.5 <~ z <~ 4 provided that the SIE models describe the gravitational potentials of the primary lenses correctly.

preprint2015arXiv

Detection of universality of dark matter profile from Subaru weak lensing measurements of 50 massive clusters

We develop a novel method of measuring the lensing distortion profiles of clusters with stacking the scaled amplitudes of background galaxy ellipticities as a function of the scaled centric radius according to the NFW prediction of each cluster, based on the assumption that the different clusters in a sample follow the universal NFW profile. First we demonstrate the feasibility of this method using both the analytical NFW model and simulated halos in high-resolution $N$-body simulations. We then apply, as a proof of concept, this method to the Subaru weak lensing data and the XMM/Chandra X-ray observables for a sample of 50 massive clusters in the redshift range $0.15\le z\le 0.3$, where their halo masses range over an order of magnitude. To estimate the NFW parameters of each cluster, we use the halo mass proxy relation of X-ray observables, based on either the hydrostatic equilibrium or the gas mass, and then infer the halo concentration from the model $c(M)$ relation. We evaluate a performance of the NFW scaling analysis by measuring the scatters of 50 cluster lensing profiles relative to the NFW predictions over a range of radii, $0.14\le R/[h^{-1}{\rm Mpc}]\le 2.8$. We found a 4 - 6$σ$ level evidence of the universal NFW profile in 50 clusters, for both the X-ray halo mass proxy relations, although the gas mass appears to be a better proxy of the underlying true mass. By comparing the measurements with the simulations of cluster lensing taking into account the statistical errors of intrinsic galaxy shapes in the Subaru data, we argue that additional halo mass errors or intrinsic scatters of $σ_{\ln M_{500c}}\sim 0.2$ - $0.3$ could reconcile a difference between the measurements and the simulations.

preprint2014arXiv

An optimal survey geometry of weak lensing survey: minimizing super-sample covariance

Upcoming wide-area weak lensing surveys are expensive both in time and cost and require an optimal survey design in order to attain maximum scientific returns from a fixed amount of available telescope time. The super-sample covariance (SSC), which arises from unobservable modes that are larger than the survey size, significantly degrades the statistical precision of weak lensing power spectrum measurement even for a wide-area survey. Using the 1000 mock realizations of the log-normal model, which approximates the weak lensing field for a $Λ$-dominated cold dark matter model, we study an optimal survey geometry to minimize the impact of SSC contamination. For a continuous survey geometry with a fixed survey area, a more elongated geometry such as a rectangular shape of 1:400 side-length ratio reduces the SSC effect and allows for a factor 2 improvement in the cumulative signal-to-noise ratio ($S/N$) of power spectrum measurement up to $\ell_{\rm max}\simeq $ a few $10^3$, compared to compact geometries such as squares or circles. When we allow the survey geometry to be disconnected but with a fixed total area, assuming $1\times 1$ sq. degrees patches as the fundamental building blocks of survey footprints, the best strategy is to locate the patches with $\sim 15$ degrees separation. This separation angle corresponds to the scale at which the two-point correlation function has a negative minimum. The best configuration allows for a factor 100 gain in the effective area coverage as well as a factor 2.5 improvement in the $S/N$ at high multipoles, yielding a much wider coverage of multipoles than in the compact geometry.

preprint2014arXiv

On the Systematic Errors of Cosmological-Scale Gravity Tests using Redshift Space Distortion: Non-linear Effects and the Halo Bias

Redshift space distortion (RSD) observed in galaxy redshift surveys is a powerful tool to test gravity theories on cosmological scales, but the systematic uncertainties must carefully be examined for future surveys with large statistics. Here we employ various analytic models of RSD and estimate the systematic errors on measurements of the structure growth-rate parameter, $fσ_8$, induced by non-linear effects and the halo bias with respect to the dark matter distribution, by using halo catalogues from 40 realisations of $3.4 \times 10^8$ comoving $h^{-3}$Mpc$^3$ cosmological N-body simulations. We consider hypothetical redshift surveys at redshifts z=0.5, 1.35 and 2, and different minimum halo mass thresholds in the range of $5.0 \times 10^{11}$ -- $2.0 \times 10^{13} h^{-1} M_\odot$. We find that the systematic error of $fσ_8$ is greatly reduced to ~5 per cent level, when a recently proposed analytical formula of RSD that takes into account the higher-order coupling between the density and velocity fields is adopted, with a scale-dependent parametric bias model. Dependence of the systematic error on the halo mass, the redshift, and the maximum wavenumber used in the analysis is discussed. We also find that the Wilson-Hilferty transformation is useful to improve the accuracy of likelihood analysis when only a small number of modes are available in power spectrum measurements.

preprint2014arXiv

Weak Lensing by Intergalactic Mini-Structures in Quadruple Lens Systems: Simulation and Detection

We investigate the weak lensing effects of line-of-sight structures on quadruple images in quasar-galaxy strong lens systems based on N-body and ray-tracing simulations that can resolve halos with a mass of 10^5 solar mass. The intervening halos and voids disturb the magnification ratios of lensed images as well as their relative positions due to lensing. The magnification ratios typically change by O(10%) when the shifts of relative angular positions of lensed images are constrained to <0.004 arcsec. The constrained amplitudes of projected density perturbations due to line-of-sight structures are O(10^8) solar mass per arcsec^2. These results are consistent with our new analytical estimate based on the two-point correlation of density fluctuations. The observed mid-infrared (MIR) flux ratios for 6 quasar-galaxy lens systems with quadruple images agree well with the numerically estimated values without taking into account of subhalos residing in the lensing galaxies. We find that the constrained mean amplitudes of projected density perturbations in the line-of-sight are negative, which suggests that the fluxes of lensed images are perturbed mainly by minivoids and minihalos in underdense regions. We derive a new fitting formula for estimating the probability distribution function of magnification perturbation. We also find that the mean amplitude of magnification perturbation roughly equals the standard deviation regardless of the model parameters.

preprint2013arXiv

Bias-Hardened CMB Lensing

We present new methods for lensing reconstruction from CMB temperature fluctuations which have smaller mean-field and reconstruction noise bias corrections than current lensing estimators, with minimal loss of signal-to-noise. These biases are usually corrected using Monte Carlo simulations, and to the extent that these simulations do not perfectly mimic the underlying sky there are uncertainties in the bias corrections. The bias-hardened estimators which we present can have reduced sensitivity to such uncertainties, and provide a desirable cross-check on standard results. To test our approach, we also show the results of lensing reconstruction from simulated temperature maps given on 100 deg^2, and confirm that our approach works well to reduce biases for a typical masked map in which 70 square masks each having 10 arcminute on a side exist, covering 2% of the simulated map, which is similar to the masks used in the current SPT lensing analysis.

preprint2013arXiv

Observational Upper Bound on the Cosmic Abundances of Negative-mass Compact Objects and Ellis Wormholes from the Sloan Digital Sky Survey Quasar Lens Search

The latest result in the Sloan Digital Sky Survey Quasar Lens Search (SQLS) has set the first cosmological constraints on negative-mass compact objects and Ellis wormholes. There are no multiple images lensed by the above two exotic objects for $\sim 50000$ distant quasars in the SQLS data. Therefore, an upper bound is put on the cosmic abundances of these lenses. The number density of negative mass compact objects is $n<10^{-8} (10^{-4}) h^3 {\rm Mpc}^{-3}$ at the mass scale $|M| > 10^{15} (10^{12}) M_\odot$, which corresponds to the cosmological density parameter $|Ω| < 10^{-4}$ at the galaxy-scale mass range $|M|=10^{12-15}M_\odot$. The number density of the Ellis wormhole is $n<10^{-4} h^3 {\rm Mpc}^{-3}$ for a range of the throat radius $a = 10^{1-4}$pc, which is much smaller than the Einstein ring radius.

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

Simulations of Wide-Field Weak Lensing Surveys I: Basic Statistics and Non-Gaussian Effects

We study the lensing convergence power spectrum and its covariance for a standard LCDM cosmology. We run 400 cosmological N-body simulations and use the outputs to perform a total of 1000 independent ray-tracing simulations. We compare the simulation results with analytic model predictions. The semi-analytic model based on Smith et al.(2003) fitting formula underestimates the convergence power by ~30% at arc-minute angular scales. For the convergence power spectrum covariance, the halo model reproduces the simulation results remarkably well over a wide range of angular scales and source redshifts. The dominant contribution at small angular scales comes from the sample variance due to the number fluctuations of halos in a finite survey volume. The signal-to-noise ratio for the convergence power spectrum is degraded by the non-Gaussian covariances by up to a factor 5 for a weak lensing survey to z_s ~1. The probability distribution of the convergence power spectrum estimators, among the realizations, is well approximated by a chi-square distribution with broadened variance given by the non-Gaussian covariance, but has a larger positive tail. The skewness and kurtosis have non-negligible values especially for a shallow survey. We argue that a prior knowledge on the full distribution may be needed to obtain an unbiased estimate on the ensemble averaged band power at each angular scale from a finite volume survey.

preprint2012arXiv

Weak Lensing by Line-of-sight Halos as the Origin of Flux-ratio Anomalies in Quadruply Lensed QSOs

We explore the weak lensing effect by line-of-sight halos and sub-halos with a mass of M < 10^7 solar mass in QSO-galaxy strong lens systems with quadruple images in a concordant LCDM universe. Using a polynomially fitted non-linear power spectrum P(k) obtained from N-body simulations that can resolve halos with a mass of M ~ 10^5 solar mass, or structures with a comoving wavenumber of k ~ 3*10^2 h/Mpc, we find that the ratio of magnification perturbation due to intervening halos to that of a primary lens is typically ~10 per cent and the predicted values agree well with the estimated values for 6 observed QSO-galaxy lens systems with quadruple images in the mid-infrared band without considering the effects of substructures inside a primary lens. We also find that the estimated amplitudes of convergence perturbation for the 6 lenses increase with the source redshift as predicted by theoretical models. Using an extrapolated matter power spectrum, we demonstrate that small halos or sub-halos in the line-of-sight with a mass of M=10^3-10^7 solar mass, or structures with a comoving wavenumber of k=3*10^2-10^4 h/Mpc can significantly affect the magnification ratios of the lensed images. Flux ratio anomalies in QSO-galaxy strong lens systems offer us a unique probe into clustering property of mini-halos with a mass of M < 10^6 solar mass.

preprint2011arXiv

Probability Distribution Functions of Cosmological Lensing: Convergence, Shear, and Magnification

We perform high resolution ray-tracing simulations to investigate probability distribution functions (PDFs) of lensing convergence, shear, and magnification on distant sources up to the redshift of $z_s=20$. We pay particular attention to the shot noise effect in $N$-body simulations by explicitly showing how it affects the variance of the convergence. We show that the convergence and magnification PDFs are closely related with each other via the approximate relation $μ=(1-κ)^{-2}$, which can reproduce the behavior of PDFs surprisingly well up to the high magnification tail. The mean convergence measured in the source plane is found to be systematically negative, rather than zero as often assumed, and is correlated with the convergence variance. We provide simple analytical formulae for the PDFs, which reproduce simulated PDFs reasonably well for a wide range of redshifts and smoothing sizes. As explicit applications of our ray-tracing simulations, we examine the strong lensing probability and the magnification effects on the luminosity functions of distant galaxies and quasars.

preprint2010arXiv

High-precision predictions for the acoustic scale in the non-linear regime

We measure shifts of the acoustic scale due to nonlinear growth and redshift distortions to a high precision using a very large volume of high-force-resolution simulations. We compare results from various sets of simulations that differ in their force, volume, and mass resolution. We find a consistency within 1.5-sigma for shift values from different simulations and derive shift alpha(z) -1 = (0.300\pm 0.015)% [D(z)/D(0)]^{2} using our fiducial set. We find a strong correlation with a non-unity slope between shifts in real space and in redshift space and a weak correlation between the initial redshift and low redshift. Density-field reconstruction not only removes the mean shifts and reduces errors on the mean, but also tightens the correlations: after reconstruction, we recover a slope of near unity for the correlation between the real and redshift space and restore a strong correlation between the low and the initial redshifts. We derive propagators and mode-coupling terms from our N-body simulations and compared with Zeldovich approximation and the shifts measured from the chi^2 fitting, respectively. We interpret the propagator and the mode-coupling term of a nonlinear density field in the context of an average and a dispersion of its complex Fourier coefficients relative to those of the linear density field; from these two terms, we derive a signal-to-noise ratio of the acoustic peak measurement. We attempt to improve our reconstruction method by implementing 2LPT and iterative operations: we obtain little improvement. The Fisher matrix estimates of uncertainty in the acoustic scale is tested using 5000 (Gpc/h)^3 of cosmological PM simulations from Takahashi et al. (2009). (abridged)

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.

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

Weak Lensing of Galaxy Clusters in MOND

We study weak gravitational lensing of galaxy clusters in terms of the MOND (MOdified Newtonian Dynamics) theory. We calculate shears and convergences of background galaxies for three clusters (A1689, CL0024+1654, CL1358+6245) and the mean profile of 42 SDSS (Sloan Digital Sky Survey) clusters and compare them with observational data. The mass profile is modeled as a sum of X-ray gas, galaxies and dark halo. For the shear as a function of the angular radius, MOND predicts a shallower slope than the data irrespective of the critical acceleration parameter $g_0$. The dark halo is necessary to explain the data for any $g_0$ and for three interpolation functions. If the dark halo is composed of massive neutrinos, its mass should be heavier than 2 eV. However the constraint still depends on the dark halo model and there are systematic uncertainties, and hence the more careful study is necessary to put a stringent constraint.