BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Takahiko Matsubara

Takahiko Matsubara contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-phgr-qchep-phhep-thastro-ph.GAmath.STStatistics Theory

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
40works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

40 published item(s)

preprint2023arXiv

Asymptotic expansion of the expected Minkowski functional for isotropic central limit random fields

The Minkowski functionals, including the Euler characteristic statistics, are standard tools for morphological analysis in cosmology. Motivated by cosmic research, we examine the Minkowski functional of the excursion set for an isotropic central limit random field, the $k$-point correlation functions ($k$th order cumulants) of which have the same structure as that assumed in cosmic research. Using 3- and 4-point correlation functions, we derive the asymptotic expansions of the Euler characteristic density, which is the building block of the Minkowski functional. The resulting formula reveals the types of non-Gaussianity that cannot be captured by the Minkowski functionals. As an example, we consider an isotropic chi-square random field and confirm that the asymptotic expansion accurately approximates the true Euler characteristic density.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Clustering of primordial black holes formed in a matter-dominated epoch

In the presence of the local-type primordial non-Gaussianity, it is known that the clustering of primordial black holes (PBHs) emerges even on super-horizon scales at the formation time. This effect has been investigated in the high-peak limit of the PBH formation in the radiation-dominated epoch in the literature. There is another possibility that the PBH formation takes place in the early matter-dominated epoch. In this scenario, the high-peak limit is not applicable because even initially small perturbations grow and can become a PBH. We first derive a general formula to estimate the clustering of PBHs with primordial non-Gaussianity without assuming the high-peak limit, and then apply this formula to a model of PBH formation in a matter-dominated epoch. Clustering is less significant in the case of the PBH formation in the matter-dominated epoch than that in the radiation-dominated epoch. Nevertheless, it is much larger than the Poisson shot noise in many cases. Relations to the constraints of the isocurvature perturbations by the cosmic microwave background radiation are quantitatively discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

de Sitter duality and logarithmic decay of dark energy

We investigate infrared dynamics of four-dimensional Einstein gravity in de Sitter space. We set up a general framework to investigate dynamical scaling relations in quantum/classical gravitational theories. The conformal mode dependence of Einstein gravity is renormalized to the extent that general covariance is not manifest. We point out that the introduction of an inflaton is necessary as a counterterm. We observe and postulate a duality between quantum effects in Einstein gravity and classical evolutions in an inflation (or quintessence) model. The effective action of Einstein gravity can be constructed as an inflation model with manifest general covariance. We show that $g=G_N H^2/π$: the only dimensionless coupling of the Hubble parameter $H^2$ and the Newton's coupling $G_N$ in Einstein gravity is screened by the infrared fluctuations of the conformal mode. We evaluate the one-loop $β$ function of $g$ with respect to the cosmic time $\log Ht$ as $β(g)=-(1/2)g^2$, i.e., $g$ is asymptotically free toward the future. The exact $β$ function with the backreaction of $g$ reveals the existence of the ultraviolet fixed point. It indicates that the de Sitter expansion stared at the Planck scale with a minimal entropy $S=2$. We have identified the de Sitter entropy $1/g$ with the von Neumann entropy of the conformal zero mode. The former evolves according to the screening of $g$ and the Gibbons-Hawking formula. The latter is found to increase by diffusion in the stochastic process at the horizon in a consistent way. Our Universe is located very close to the fixed point $g=0$ with a large entropy. We discuss possible physical implications of our results such as logarithmic decay of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The statistics of peaks of weakly non-Gaussian random fields: Effects of bispectrum in two- and three-dimensions

Analytic expressions for the statistics of peaks of random fields with weak non-Gaussianity are provided. Specifically, the abundance and spatial correlation of peaks are represented by formulas which can be evaluated only by virtually one-dimensional integrals. We assume the non-Gaussianity is weak enough such that it is represented by linear terms of the bispectrum. The formulas are formally given in $N$-dimensional space, and explicitly given in the case of $N=1,2,3$. Some examples of peak statistics in cosmological fields are calculated for the cosmic density field and weak lensing field, assuming the weak non-Gaussianity is induced by gravity. The formulas of this paper would find a fit in many applications to statistical analyses of cosmological fields.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The large-separation expansion of peak clustering in Gaussian random fields

In the peaks approach, the formation sites of observable structures in the Universe are identified as peaks in the matter density field. The statistical properties of the clustering of peaks are particularly important in this respect. In this paper, we investigate the large-separation expansion of the correlation function of peaks in Gaussian random fields. The analytic formula up to third order is derived, and the resultant expression can be evaluated by a combination of one-dimensional fast Fourier transforms, which are evaluated very fast. The analytic formula obtained perturbatively in the large-separation limit is compared with a method of Monte-Carlo integrations, and a complementarity between the two methods is demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

CMB probes on the correlated axion isocurvature perturbation

We explore the possible cosmological consequence of the gravitational coupling between the inflaton and axion-like fields. In view of the forthcoming cosmic microwave background (CMB) polarization and lensing data, we study the sensitivity of the CMB data on the cross-correlation between the curvature and axion isocurvature perturbations. Through a concrete example, we illustrate the explicit dependence of the scale dependent cross-correlation power spectrum on the axion parameters.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Impacts of biasing schemes in the one-loop integrated perturbation theory

The impact of biasing schemes on the clustering of tracers of the large-scale structure is analytically studied in the weakly nonlinear regime. For this purpose, we use the one-loop approximation of the integrated perturbation theory together with the renormalized bias functions of various, physically motivated Lagrangian bias schemes. These include the halo, peaks and excursion set peaks model, for which we derive useful formulae for the evaluation of their renormalized bias functions. The shapes of the power spectra and correlation functions are affected by the different bias models at the level of a few percent on weakly nonlinear scales. These effects are studied quantitatively both in real and redshift space. The amplitude of the scale-dependent bias in the presence of primordial non-Gaussianity also depends on the details of the bias models. If left unaccounted for, these theoretical uncertainties could affect the robustness of the cosmological constraints extracted from galaxy clustering data.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Recursive Solutions of Lagrangian Perturbation Theory

In the standard perturbation theory (SPT) of self-gravitating Newtonian fluid in an expanding universe, recurrence relations for higher-order solutions are well known and play an important role both in practical applications and in theoretical investigations. The recurrence relations in Lagrangian perturbation theory (LPT), however, had not been known for a long time. Recently, two different kinds of recurrence relations in LPT have been proposed in limited cases. In this paper, we generalize those methods, and most generally derive the recurrence relations, which are capable to include any initial condition in general models of cosmology. The fastest-growing modes in the general relations are identified, and simplified recurrence relations with accurate approximation for the time-dependence are obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Subaru FMOS Galaxy Redshift Survey (FastSound). I. Overview of the Survey Targeting on H$α$ Emitters at $z \sim 1.4$

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Integrated Perturbation Theory and One-loop Power Spectra of Biased Tracers

General and explicit predictions from the integrated perturbation theory (iPT) for power spectra and correlation functions of biased tracers are derived and presented in the one-loop approximation. The iPT is a general framework of the nonlinear perturbation theory of cosmological density fields in presence of nonlocal bias, redshift-space distortions, and primordial non-Gaussianity. Analytic formulas of auto and cross power spectra of nonlocally biased tracers in both real and redshift spaces are derived and the results are comprehensively summarized. The main difference from previous formulas derived by the present author is to include effects of generally nonlocal Lagrangian bias and primordial non-Gaussianity, and the derivation method of the new formula is fundamentally different from the previous one. Relations to recent work on improved methods of nonlinear perturbation theory in literature are clarified and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Gaussian Minkowski functionals & extrema counts in redshift space

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear stochastic biasing of halos: Analysis of cosmological N-body simulations and perturbation theories

It is crucial to understand and model a behavior of galaxy biasing for future ambitious galaxy redshift surveys. Using 40 large cosmological N-body simulations for a standard LambdaCDM cosmology, we study the cross-correlation coefficient between matter and the halo density field, which is an indicator of the stochasticity of bias, over a wide redshift range 0\le z \le 3. The cross-correlation coefficient is important to extract information on the matter density field, e.g., by combining galaxy clustering and galaxy-galaxy lensing measurements. We compare the simulation results with integrated perturbation theory (iPT) proposed by one of the present authors and standard perturbation theory (SPT) combined with a phenomenological model of local bias. The cross-correlation coefficient derived from the iPT agrees with N-body simulation results down to r~15 (10) h^{-1}Mpc within 0.5 (1.0) % for all redshifts and halo masses we consider. The SPT with local bias does not explain complicated behaviors on quasilinear scales at low redshifts, while roughly reproduces the general behavior of the cross-correlation coefficient on fully nonlinear scales. The iPT is powerful to predict the cross-correlation coefficient down to quasilinear regimes with a high precision.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scale-dependent bias due to primordial vector fields

Anisotropic stress perturbations induced by primordial Gaussian vector fields create non-Gaussianity in curvature perturbations. We found that such non-Gaussianity closely resembles the local-type non-Gaussianity parametrized by $f_{\rm NL}$, and generates scale-dependent bias in large-scale structures. We also found a simple relationship between the scale-dependent bias and the power spectrum of the vector fields. When the vector fields are interpreted as primordial magnetic fields, the effective $f_{\rm NL}$ is shown to be always negative. The scale-dependent bias provides a new approach to probing primordial vector fields.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Application of cross correlations between CMB and large scale structure to constraints on the primordial non-Gaussianity

The primordial non-Gaussianity of local type affects the clustering of dark matter halos, and the planned deep and wide photometric surveys are suitable for examining this class of non-Gaussianity. In our previous paper, we investigated the constraint from the cross correlation between CMB lensing potential and galaxy angular distribution on the primordial non-Gaussianity, without taking into account redshift slicing. To improve our previous analysis, in this paper, we add the galaxy lensing shear into our analysis and take into account redshift slicing to follow the redshift evolution of the clustering. By calculating 81 power spectra and using the Fisher matrix method, we find that the constraint on the primordial non-Gaussianity can be improved from ΔfNL \sim 5.4 to 5.1 by including the galaxy-galaxy lensing shear cross correlations expected from the Hyper Suprime-Cam survey (HSC), in comparison with the constraint without any cross correlations. Moreover, the constraint can go down to ΔfNL \sim 4.8 by including the galaxy-CMB lensing cross correlations from the ACTPol and Planck experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Deriving an Accurate Formula of Scale-dependent Bias with Primordial Non-Gaussianity: An Application of the Integrated Perturbation Theory

We apply the integrated perturbation theory (Matsubara 2011, PRD 83, 083518) to evaluate the scale-dependent bias in the presence of primordial non-Gaussianity. The integrated perturbation theory is a general framework of nonlinear perturbation theory, in which a broad class of bias models can be incorporated into perturbative evaluations of biased power spectrum and higher-order polyspectra. Approximations such as the high-peak limit or the peak-background split are not necessary to derive the scale-dependent bias in this framework. Applying the halo approach, previously known formulas are re-derived as limiting cases of a general formula in this work, and it is implied that modifications should be made in general situations. Effects of redshift-space distortions are straightforwardly incorporated. It is found that the slope of the scale-dependent bias on large scales is determined only by the behavior of primordial bispectrum in the squeezed limit, and is not sensitive to bias models in general. It is the amplitude of scale-dependent bias that is sensitive to the bias models. The effects of redshift-space distortions turn out to be quite small for the monopole component of the power spectrum, while the quadrupole component is proportional to the monopole component on large scales, and thus also sensitive to the primordial non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Limits on Second-Order Non-Gaussianity from Minkowski Functionals of WMAP Data

We analyze non-Gaussianity (NG) due to the primordial bispectrum and trispectrum using CMB temperature maps of WMAP 7-year data. We first apply the perturbative formulae of Minkowski functionals up to second-order NG derived by Matsubara (2010), which enable us to give limits on cubic NG parametrized with tau_NL and g_NL as well as various types of quadratic NG parametrized with f_NL. We find no signature of primordial NG in WMAP 7-year data, but give constraints on the local-type, equilateral-type, orthogonal-type f_NL: f_NL(loc)=20+-42, f_NL(eq)=-121+-208, f_NL(ort)=-129+-171, respectively, and tau_NL/10^4=-7.6+-8.7, and g_NL/10^5=-1.9+-6.4. We also find that these constraints are consistent with the limits from skewness and kurtosis parameters which characterize the perturbative corrections of MFs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scale-dependent bias with higher order primordial non-Gaussianity: Use of the Integrated Perturbation Theory

We analytically derive a more accurate formula for the power spectrum of the biased objects with the primordial non-Gaussianity parameterized not only by the non-linearity parameter fNL, but also by gNL and tauNL which characterize the trispectrum of the primordial curvature perturbations. We adopt the integrated perturbation theory which was constructed in Matsubara (2011). We discuss an inequality between fNL and tauNL in the context of the scale-dependent bias, by introducing a stochasticity parameter. We also mention higher order loop corrections into the scale-dependency of the bias parameter.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonlinear Biasing and Redshift-Space Distortions in Lagrangian Resummation Theory and N-body Simulations

Understanding a behavior of galaxy biasing is crucial for future galaxy redshift surveys. One aim is to measure the baryon acoustic oscillations (BAOs) within the precision of a few percent level. Using 30 large cosmological N-body simulations for a standard LCDM cosmology, we study the halo biasing over a wide redshift range. We compare the simulation results with theoretical predictions proposed by Matsubara (2008) which naturally incorporate the halo bias and redshift-space distortions into their formalism of perturbation theory with a resummation technique via the Lagrangian picture. The power spectrum and correlation function of halos obtained from Lagrangian resummation theory (LRT) well agree with N-body simulation results on scales of BAOs. Especially nonlinear effects on the baryon acoustic peak of the halo correlation function are accurately explained both in real and redshift space. We find that nonlinearity and scale dependence of bias are fairly well reproduced by 1-loop LRT up to k=0.35hMpc^{-1} (z=2 and 3) within a few percent level in real space and up to k=0.1hMpc^{-1} (z=2) and 0.15hMpc^{-1} (z=3) in redshift space. Thus, the LRT is very powerful for accurately extracting cosmological information in upcoming high redshift BAO surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonlinear Perturbation Theory Integrated with Nonlocal Bias, Redshift-space Distortions, and Primordial Non-Gaussianity

The standard nonlinear perturbation theory of the gravitational instability is extended to incorporate the nonlocal bias, redshift-space distortions, and primordial non-Gaussianity. We show that local Eulerian bias is not generally compatible to local Lagrangian bias in nonlinear regime. The Eulerian and Lagrangian biases are nonlocally related order by order in the general perturbation theory. The relation between Eulerian and Lagrangian kernels of density perturbations with biasing are derived. The effects of primordial non-Gaussianity and redshift-space distortions are also incorporated in our general formalism, and diagrammatic methods are introduced. Vertex resummations of higher-order perturbations in the presence of bias are considered. Resummations of Lagrangian bias are shown to be essential to handle biasing schemes in a general framework.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulations of Wide-Field Weak Lensing Surveys II: Covariance Matrix of Real Space Correlation Functions

Using 1000 ray-tracing simulations for a Λ-dominated cold dark model in Sato et al. (2009), we study the covariance matrix of cosmic shear correlation functions, which is the standard statistics used in the previous measurements. The shear correlation function of a particular separation angle is affected by Fourier modes over a wide range of multipoles, even beyond a survey area, which complicates the analysis of the covariance matrix. To overcome such obstacles we first construct Gaussian shear simulations from the 1000 realizations, and then use the Gaussian simulations to disentangle the Gaussian covariance contribution to the covariance matrix we measured from the original simulations. We found that an analytical formula of Gaussian covariance overestimates the covariance amplitudes due to an effect of finite survey area. Furthermore, the clean separation of the Gaussian covariance allows to examine the non-Gaussian covariance contributions as a function of separation angles and source redshifts. For upcoming surveys with typical source redshifts of z_s=0.6 and 1.0, the non-Gaussian contribution to the diagonal covariance components at 1 arcminute scales is greater than the Gaussian contribution by a factor of 20 and 10, respectively. Predictions based on the halo model qualitatively well reproduce the simulation results, however show a sizable disagreement in the covariance amplitudes. By combining these simulation results we develop a fitting formula to the covariance matrix for a survey with arbitrary area coverage, taking into account effects of the finiteness of survey area on the Gaussian covariance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Analytic Minkowski Functionals of the Cosmic Microwave Background: Second-order Non-Gaussianity with Bispectrum and Trispectrum

Analytic formulas of Minkowski functionals in two-dimensional random fields are derived, including effects of second-order non-Gaussianity in the presence of both the bispectrum and trispectrum. The set of formulas provides a promising method to constrain the primordial non-Gaussianity of the universe by temperature fluctuations in the cosmic microwave background radiation. In a case of local-type non-Gaussianity, the Minkowski functionals are analytically given by powers of quadratic and cubic parameters, $f_{\rm NL}$ and $g_{\rm NL}$. Our formulas are not restricted to this particular model, and applicable to a wide class of non-Gaussian models. The analytic formulas are compared to numerical evaluations from non-Gaussian realizations of temperature maps, showing very good agreements.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraints on primordial non-Gaussianity from Galaxy-CMB lensing cross-correlation

Recent studies have shown that the primordial non-Gaussianity affects clustering of dark matter halos through a scale-dependent bias and various constraints on the non-Gaussianity through this scale-dependent bias have been placed. Here we introduce the cross-correlation between the CMB lensing potential and the galaxy angular distribution to effectively extract information about the bias from the galaxy distribution. Then, we estimate the error of non-linear parameter, f_NL, for the on-going CMB experiments and galaxy surveys, such as Planck and Hyper Suprime-Cam (HSC). We found that for the constraint on f_NL with Planck and HSC, the wide field galaxy survey is preferable to the deep one, and the expected error on f_NL can be as small as: Δf_NL ~ 20 for b_0 = 2 and Δf_NL ~ 10 for b_0 = 4, where b_0 is the linear bias parameter. It is also found that future wide field galaxy survey could achieve Δf_NL ~ 5 with CMB prior from Planck if one could observe highly biased objects at higher redshift (z ~ 2).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Nonlinear perturbation theory with halo bias and redshift-space distortions via the Lagrangian picture

The nonlinear perturbation theory of gravitational instability is extended to include effects of both biasing and redshift-space distortions, which are inevitable in predicting observable quantities in galaxy surveys. The precise determination of scales of baryon acoustic oscillations is crucial to investigate the nature of dark energy by galaxy surveys. We find that a local Lagrangian bias and redshift-space distortions are naturally incorporated in our formalism of perturbation theory with a resummation technique via the Lagrangian picture. Our formalism is applicable to any biasing scheme which is local in Lagrangian space, including the halo bias as a special case. Weakly nonlinear effects on halo clustering in redshift space are analytically given.We assume only a fundamental idea of the halo model: haloes form according to the extended Press-Schechter theory, and the spatial distributions are locally biased in Lagrangian space. There is no need for assuming the spherical collapse model to follow the dynamical evolution, which is additionally assumed in standard halo prescriptions. One-loop corrections to the power spectrum and correlation function of haloes in redshift space are explicitly derived and presented. Instead of relying on expensive numerical simulations, our approach provides an analytic way of investigating the weakly nonlinear effects, simultaneously including the nonlinear biasing and nonlinear redshift-space distortions. Nonlinearity introduces a weak scale dependence in the halo bias. The scale dependence is a smooth function in Fourier space, and the bias does not critically change the feature of baryon acoustic oscillations in the power spectrum. The same feature in the correlation function is less affected by nonlinear effects of biasing.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Effects of Smoothing Functions in Cosmological Counts-in-Cells Analysis

A method of counts-in-cells analysis of galaxy distribution is investigated with arbitrary smoothing functions in obtaining the galaxy counts. We explore the possiblity of optimizing the smoothing function, considering a series of $m$-weight Epanechnikov kernels. The popular top-hat and Gaussian smoothing functions are two special cases in this series. In this paper, we mainly consider the second moments of counts-in-cells as a first step. We analytically derive the covariance matrix among different smoothing scales of cells, taking into account possible overlaps between cells. We find that the Epanechnikov kernel of $m=1$ is better than top-hat and Gaussian smoothing functions in estimating cosmological parameters. As an example, we estimate expected parameter bounds which comes only from the analysis of second moments of galaxy distributions in a survey which is similar to the Sloan Digital Sky Survey.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Primordial Non-Gaussianity and Analytical Formula for Minkowski Functionals of the Cosmic Microwave Background and Large-scale Structure

We derive analytical formulae for the Minkowski Functions of the cosmic microwave background (CMB) and large-scale structure (LSS) from primordial non-Gaussianity. These formulae enable us to estimate a non-linear coupling parameter, f_NL, directly from the CMB and LSS data without relying on numerical simulations of non-Gaussian primordial fluctuations. One can use these formulae to estimate statistical errors on f_NL from Gaussian realizations, which are much faster to generate than non-Gaussian ones, fully taking into account the cosmic/sampling variance, beam smearing, survey mask, etc. We show that the CMB data from the Wilkinson Microwave Anisotropy Probe should be sensitive to |f_NL|\simeq 40 at the 68% confidence level. The Planck data should be sensitive to |f_NL|\simeq 20. As for the LSS data, the late-time non-Gaussianity arising from gravitational instability and galaxy biasing makes it more challenging to detect primordial non-Gaussianity at low redshifts. The late-time effects obscure the primordial signals at small spatial scales. High-redshift galaxy surveys at z>2 covering \sim 10Gpc^3 volume would be required for the LSS data to detect |f_NL|\simeq 100. Minkowski Functionals are nicely complementary to the bispectrum because the Minkowski Functionals are defined in real space and the bispectrum is defined in Fourier space. This property makes the Minksowski Functionals a useful tool in the presence of real-world issues such as anisotropic noise, foreground and survey masks. Our formalism can be extended to scale-dependent f_NL easily.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Fourier Phase Analysis of SDSS Galaxies

We present a first analysis of the clustering of SDSS galaxies using the distribution function of the sum of Fourier phases. This statistic was recently proposed by one of authors as a new method to probe phase correlations of cosmological density fields. Since the Fourier phases are statistically independent of Fourier amplitudes, the phase statistic plays a complementary role to the conventional two-point statistics of galaxy clustering. In particular, we focus on the distribution functions of phase sum over three closed wavevectors as a function of the configurations of triangle wavevectors. We find that the observed distribution functions of phase sum are in good agreement with the lowest-order approximation from perturbation theory. For direct comparison with observations, we construct mock catalogs from N-body simulations taking account of several observational effects such as the survey geometry, the redshift distortion, and the discreteness due to the limited number of data. Indeed the observed phase correlations for the galaxies in the range of absolute magnitude -22<Mr<-18 agree well with those for Lambda CDM predictions with sigma8=0.9 evolved from Gaussian initial condition. This agreement implies that the galaxy biasing is approximately linear in redshift space. Instead, assuming that the galaxy biasing is described by a quadratic deterministic function at scale larger than 30Mpc/h, we can constrain the ratio of the quadratic biasing parameter b2 to the linear biasing parameter b1. We find that resulting b2/b1 is mainly dependent on sigma8 and that it is well fitted by b2/b1=0.54(+-0.06)-0.62(+-0.08)sigma8 in CDM models. Indeed, b2/b1 is nearly zero when sigma8=0.9.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Large Scale Clustering of Sloan Digital Sky Survey Quasars: Impact of the Baryon Density and the Cosmological Constant

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

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galaxy Structures and External Perturbations in Gravitational Lenses

In modeling strong gravitational lens systems, one often adopts simple models, such as singular isothermal elliptical density plus lowest-order external perturbation. However, such simple models may mislead us if the real mass distribution is more complicated than that in the assumed models. In particular, assumptions on mass models are crucial in studying flux ratio anomalies that have been suggested as evidence for a cold dark-matter substructure. We reinvestigated four quadruple lens systems using power-law Fourier models, which have advantages of clear physical meanings and the applicability of a linear method, as well as a simple singular isothermal elliptical density model. We also investigated the effect of external perturbations, including a singular isothermal sphere, lowest-order expansion, and next-order expansion. We have found that the $\cos3θ$ terms of the primary galaxy and/or of external perturbation significantly reduce the χ^2 in PG 1115+080 and B1422+231. In particular, we could reproduce the flux ratios of B1422+231 with next-order external perturbation assuming 5% flux uncertainties, suggesting that external perturbation cannot be described by a simple singular isothermal sphere approximation. On the other hand, we could not fit B0712+472 and B2045+265 very well even with our models, although the χ^2 were reduced compared with the case of using simple models. Our results clearly demonstrate that both the primary lens galaxy and the external perturbation are often more complicated than we usually assume.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Three-point Correlation Functions of SDSS Galaxies in Redshift Space: Morphology, Color, and Luminosity Dependence

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Diagrammatic Methods in Statistics and Biasing in the Large-scale Structure of the Universe

We develop technique concerning statistical analyses of the large-scale structure of the universe in the presence of biasing in the structure formation. We formulate the diagrammatic method to calculate the correlation functions of the nonlocally biased field from a generally non-Gaussian density field. Our method is based on generalized Wiener-Hermite expansion of the density field. The present formalism has not only general applicability, but also practical significance, too. To show the effectiveness of this method, we revisit the problems on biasing that have been considered previously, i.e., various approximations for peak statistics and hierarchical structure of correlation functions of locally biased field. Further analyses which have not been possible so far on these problems can be performed by our formalism. The gravitational evolution of primordial fluctuation or other structure-forming processes can be described by nonlocal biasing, so can be treated by our new formalism in principle.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Statistics of Isodensity Contours in Redshift Space

The peculiar velocities of galaxies distort the clustering pattern anisotropically in redshift space. This effect on the statistics of isodensity contours is examined by linear theory. The statistics considered in this paper are the three- and two-dimensional genus of isodensity contours, the area of isodensity contours, the length of isodensity contours in the 2-dimensional slice and the level crossing statistics on the line. We find that all these statistics in redshift space as functions of density threshold of contours have the same shape as in real space. The redshift space distortion affects only amplitudes of these statistics. The three-dimensional genus and the area hardly suffer from the redshift space distortion for $0 \leq Ωb^{-5/3} \leq 1$, where $b$ is a linear bias parameter. The other statistics are defined in one- or two-dimensional slices of the sample volume and depend on the direction of these slices relative to the line of sight. This dependence on direction of these statistics provides a way to determine the density parameter of the universe.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Weakly Nonlinear Evolution of Topology of Large-scale Structure

The gravitational evolution of the genus and other statistics of isodensity contours of the density field is derived analytically in a weakly nonlinear regime using second-order perturbation theory. The effect of final smoothing in perturbation theory on the statistics of isodensity contours is also evaluated. The resulting analytic expression for the genus is compared with $N$-body numerical simulations and exhibits a good agreement.

0 citations0 reviewsNo code linkedOpen access