Source author record

Masahiro Takada

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph gr-qc astro-ph.IM astro-ph.SR hep-ph astro-ph.EP astro-ph.HE physics.data-an

Catalog footprint

What is connected

86works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2024arXiv

The Milky Way Tomography with Subaru Hyper Suprime-Cam. I. Halo substructures

We analyze the photometric data in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) over $\sim 1,200$ deg$^{2}$ to uncover new halo substructures beyond the distance, $D_{\odot}\sim$ 30 kpc, from the Sun. For this purpose, we develop an isochrone filter for an old, metal-poor stellar system to extract the faint main-sequence stars at a range of distances. With this method, we detect, not only the previously discovered substructures such as the Orphan Stream, but also the new overdensity toward Boötes at about $D_{\odot}\sim$ 60 kpc and the new stream-like feature toward Pisces at around $D_{\odot}\sim$ 60 kpc. It has been suggested that a small-scale overdensity exists in this direction of Pisces (the so-called Pisces Overdensity), but our results show that the overdensity is widely spread with a tidally elongated feature. Combining our results with the ongoing Hyper Suprime-Cam narrow-band survey and the near-future spectroscopic survey with Prime Focus Spectrograph (PFS) will allow us to place strong constraints on the origin of these halo substructures.

preprint2022arXiv

Analysis method for 3D power spectrum of projected tensor field with fast estimator and window convolution modelling: an application to intrinsic alignments

Rank-2 tensor fields of large-scale structure, e.g. a tensor field inferred from shapes of galaxies, open up a window to directly access 2-scalar, 2-vector and 2-tensor modes, where the scalar fields can be measured independently from the standard density field that is traced by distribution of galaxies. Here we develop an estimator of the multipole moments of power spectra for the three-dimensional tensor field, taking into account the projection onto plane perpendicular to the line-of-sight direction. To do this, we find that a convenient representation of the power spectrum multipoles can be obtained by the use of the associated Legendre polynomials in the form which allows for the fast Fourier transform estimations under the local plane-parallel (LPP) approximation. The formulation also allows us to obtain the Hankel transforms to connect the two-point statistics in Fourier and configuration space, which are needed to derive theoretical templates of the power spectrum including convolution of a survey window. To validate our estimators, we use the simulation data of the projected tidal field assuming a survey window that mimics the BOSS-like survey footprint. We show that the LPP estimators fairly well recover the multipole moments that are inferred from the global plane-parallel approximation. We find that the survey window causes a more significant change in the multipole moments of projected tensor power spectrum at $k\lesssim 0.1\,h{\rm Mpc}^{-1}$ from the input power spectrum, than in the density power spectrum. Nevertheless, our method to compute the theory template including the survey window effects successfully reproduces the window-convolved multipole moments measured from the simulations. The analysis method presented here paves the way for a cosmological analysis using three-dimensional tensor-type tracers of large-scale structure for current and future surveys.

preprint2022arXiv

Constraining Dark matter annihilation with Dark Energy Survey Y3 LSBG sample

To reveal natures of the dark matter (DM) particles, a gamma-ray signal produced in annihilation processes of DM into the standard model particles has been one of the major probes. The cross-correlation between highly DM dominated structures, such as local dwarf galaxies, and observed photons in the direction of the structures, has been explored and provided stringent constraints on the annihilation rate. In our previous work, we have shown that it is sufficient to know the distance distribution of the galaxy sample and individual distance measurement is not required for constraining the annihilation rate. In this work, we apply the method to low-surface brightness galaxies (LSBGs) with unknown individual redshifts provided from the Dark Energy Survey (DES) Year 3 data. With all DES-LSBGs of about 24,000 objects, we find that the upper limits of the cross section for bb channel with 95% C.L. is 3*10^-25 cm^3/s at DM mass of 100 GeV. To be more conservative, we remove about 7000 LSBGs within 1 degree from resolved gamma-ray point sources and then the constraint becomes 30 % weaker than the one with all samples in all DM mass ranges.

preprint2022arXiv

Cosmological constraints without fingers of God

Non-linear redshift-space distortions ("fingers of God") are challenging to model analytically, a fact that limits the applicability of perturbation theory in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, $Q_0$, which can be easily estimated from the redshift space clustering data and is approximately equal to the real space power spectrum. The new statistic does not suffer from fingers of God and can be accurately described with perturbation theory down to $k_{\rm max}\simeq 0.4~h~\text{Mpc}^{-1}$. It can be straightforwardly included in the likelihood at negligible additional computational cost, and yields noticeable improvements on cosmological parameters compared to standard power spectrum multipole analyses. Using both simulations and observational data from the Baryon Oscillation Spectroscopic Survey, we show that improvements vary from $10\%$ to $100\%$ depending on the cosmological parameter considered, the galaxy sample and the survey volume.

preprint2022arXiv

Dark matter constraints with stacked gamma rays scales with the number of galaxies

Low-surface-brightness galaxies (LSBGs) are interesting targets for searches of dark matter emission due to their low baryonic content. However, predicting their expected dark matter emissivities is difficult because of observational challenges in their distance measurements. Here we present a stacking method that makes use of catalogs of LSBGs and maps of unresolved gamma-ray emission measured by the Fermi Gamma-Ray Space Telescope. We show that, for relatively large number of LSBGs, individual distance measurements to the LSBGs are not necessary, instead the overall distance distribution of the population is sufficient in order to impose dark matter constraints. Further, we demonstrate that the effect of the covariance between two galaxies located closely -- at an angular distance comparable to the size of the Fermi point spread function -- is negligibly small. As a case in point, we apply our pipeline to a sample of 800 faint LSBGs discovered by Hyper Suprime-Cam and find that, the 95 per cent confidence level upper limits on the dark matter annihilation cross-section scales with inverse of the number LSBGs. In light of this linear dependence with the number of objects, we argue this methodology could provide extremely powerful limits if it is applied to the more than 105 LSBGs readily available with the Legacy Survey of Space and Time.

preprint2022arXiv

Full-shape cosmology analysis of SDSS-III BOSS galaxy power spectrum using emulator-based halo model: a $5\%$ determination of $σ_8$

We present the results obtained from the full-shape cosmology analysis of the redshift-space power spectra for 4 galaxy samples of the SDSS-III BOSS DR12 galaxy catalog over $0.2 < z < 0.75$. For the theoretical template, we use an emulator that was built from an ensemble set of $N$-body simulations, which enables fast and accurate computation of the redshift-space power spectrum of halos. Combining with the halo occupation distribution to model the halo-galaxy connection, we can compute the redshift-space power spectrum of BOSS-like galaxies in the flat $Λ$CDM cosmology. In our cosmology inference, we use the power spectrum monopole, quadrupole and hexadecapole and include 7 nuisance parameters to model uncertainties in the halo-galaxy connection for each galaxy sample, but do not use any information on the abundance of galaxies. We demonstrate a validation of our analysis pipeline using the mock catalogs of BOSS-like galaxies, generated using different recipes of the halo-galaxy connection and including the assembly bias effect. Assuming weak priors on cosmological parameters, except for $Ω_{\rm b}h^2$ and $n_{\rm s}$, we show that our model well reproduces the BOSS power spectra. Including the power spectrum information up to $k_{\rm max}=0.25\,h{\rm Mpc}^{-1}$, we find $Ω_\mathrm{m}=0.301^{+0.012}_{-0.011}$, $H_0=68.2 \pm 1.4~\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, and $σ_8=0.786^{+0.036}_{-0.037}$, for the mode and 68\% credible interval, after marginalization over nuisance parameters. We find little improvement in the cosmological parameters beyond a maximum wavelength $k_{\rm max}\simeq 0.2\,h\,{\rm Mpc}^{-1}$ due to the shot noise domination and marginalization of the halo-galaxy connection parameters. Our results are consistent with the Planck CMB results within $1σ$ statistical uncertainties.

preprint2022arXiv

Mock catalogues of emission line galaxies based on the local mass density in dark-matter only simulations

The high-precision measurement of spatial clustering of emission line galaxies (ELGs) is a primary objective for upcoming cosmological spectroscopic surveys. The source of strong emission of ELGs is nebular emission from surrounding ionized gas irradiated by massive short-lived stars in star-forming galaxies. As a result, ELGs are more likely to reside in newly-formed halos and this leads to a nonlinear relation between ELG number density and matter density fields. In order to estimate the covariance matrix of cosmological observables, it is essential to produce many independent realisations to simulate ELG distributions for large survey volumes. To this end, we present a novel and fast scheme to populate ELGs in dark-matter only $N$-body simulations based on local density field. This method enables fast production of mock ELG catalogues suitable for verifying analysis methods and quantifying observational systematics in upcoming spectroscopic surveys and can populate ELGs in moderately high-density regions even though the halo structure cannot be resolved due to low resolution. The power spectrum of simulated ELGs is consistent with results of hydrodynamical simulations up to fairly small scales ($\lesssim 1 h \, \mathrm{Mpc}^{-1}$), and the simulated ELGs are more likely to be found in filamentary structures, which is consistent with results of semi-analytic and hydrodynamical simulations. Furthermore, we address the redshift-space power spectrum of simulated ELGs. The measured multipole moments of simulated ELGs clearly exhibit a weaker Finger-of-God effect than those of matter due to infalling motions towards halo centre, rather than random virial motions inside halos.

preprint2022arXiv

The Subaru HSC weak lensing mass-observable scaling relations of spectroscopic galaxy groups from the GAMA survey

We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) to study the dark matter content of galaxy groups in the Universe using weak lensing. We use galaxy groups from the Galaxy Mass and Assembly galaxy survey in approximately $100$ sq. degrees of the sky that overlap with the HSC survey as lenses. We restrict our analysis to the $1587$ groups with at least five members. We divide these groups into six bins each of group luminosity and group member velocity dispersion and measure the lensing signal with a signal-to-noise ratio of $55$ and $51$ for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our groups binned in the two different observable properties and constrain the power-law scaling relation, and the scatter between mean halo masses and the two group observable properties. We obtain a 5 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\avg{M} = (0.81\pm 0.04)\times10^{14}\hinvMsun$ for $L_{\rm grp}=10^{11.5}\hinvsqLsun$, and a power-law index of $α=1.01\pm 0.07$. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be $\avg{M}=(0.93\pm 0.05)\times10^{14}\hinvMsun$ for $σ=500 \kms$ and a power-law index to be $α=1.52\pm0.10$. However, these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature show that our results are consistent and have significantly reduced errors.

preprint2022arXiv

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|δm|<9\times 10^{-3}$. The bias residual $δm$ shows no dependence on redshift in the range $0<z\leq 3$. We define the requirements for cosmological weak lensing science for this shear catalog, and quantify potential systematics in the catalog using a series of internal null tests for systematics related to point-spread function modelling and shear estimation. A variety of the null tests are statistically consistent with zero or within requirements, but (i) there is evidence for PSF model shape residual correlations; and (ii) star-galaxy shape correlations reveal additive systematics. Both effects become significant on $>1$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements.

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

Precision analysis of the redshift-space galaxy bispectrum

We study the information content of the angle-averaged (monopole) redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the PT challenge simulations, whose cumulative volume of 566 $h^{-3}$Gpc$^3$ allows for a precise comparison to theoretical predictions. Fitting the power spectrum and bispectrum of our simulated data, and varying all necessary cosmological and nuisance parameters in a consistent Markov chain Monte Carlo analysis, we find that our tree-level bispectrum model is valid up to $k_{\max}=0.08~h{\rm Mpc}^{-1}$ (at $z=0.61$). We also find that inclusion of the bispectrum monopole improves constraints on cosmological parameters by $(5-15)\%$ relative to the power spectrum. The improvement is more significant for the quadratic bias parameters of our simulated galaxies, which we also show to deviate from biases of the host dark matter halos at the $\sim 3σ$ level. Finally, we adjust the covariance and scale cuts to match the volume of the BOSS survey, and estimate that within the minimal $Λ$CDM model the bispectrum data can tighten the constraint on the mass fluctuation amplitude $σ_8$ by roughly $10\%$.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag.

preprint2021arXiv

Super-sample covariance of the thermal Sunyaev-Zel'dovich effect

The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful probe of cosmology. The statistical errors in the tSZ power spectrum measurements are dominated by the presence of massive clusters in a survey volume that are easy to identify on individual cluster basis. First, we study the impact of super sample covariance (SSC) on the tSZ power spectrum measurements, and find that the sample variance is dominated by the connected non-Gaussian (cNG) covariance arising mainly from Poisson number fluctuations of massive clusters in the survey volume. Second, we find that removing such individually-detected, massive clusters from the analysis significantly reduces the cNG contribution, thereby leading the SSC to be a leading source of the sample variance. We then show, based on Fisher analysis, that the power spectrum measured from the remaining diffuse tSZ effects can be used to obtain tight constraints on cosmological parameters as well as the hydrostatic mass bias parameter. Our method offers complementary use of individual tSZ cluster counts and the power spectrum measurements of diffuse tSZ signals for cosmology and intracluster gas physics.

preprint2021arXiv

Testing Stochastic Gravitational Wave Signals from Primordial Black Holes with Optical Telescopes

Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained "sublunar" mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account for the recently reported North American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array data signal, and will be testable in future GW observations by interferometer-type experiments such as Laser Interferometer Space Antenna (LISA). We show that the broad mass function of such PBH DM is already being probed by Subaru Hyper Suprime-Cam (HSC) microlensing data and is consistent with a detected candidate event. Upcoming observations of HSC will be able to provide an independent definitive test of the stochastic GW signals originating from such PBH DM production scenarios.

preprint2021arXiv

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/.

preprint2020arXiv

Anisotropic separate universe simulations

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

preprint2020arXiv

Blinded challenge for precision cosmology with large-scale structure: results from effective field theory for the redshift-space galaxy power spectrum

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

preprint2020arXiv

Constraining dark matter annihilation with HSC Low Surface Brightness Galaxies

Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be less contaminated by extragalactic gamma-ray sources (e.g., blazars) compared to star forming galaxies. We report on an analysis that uses eight LSBGs (detected by Subaru Hyper Suprime-Cam survey data) with known redshifts to conduct a search for gamma-ray emission at the positions of these new objects in Fermi Large Area Telescope data. We found no excesses of gamma-ray emission and set constraints on the dark matter annihilation cross-section. We exclude (at the 95% C.L.) dark matter scenarios predicting a cross-section higher than 10^-23[cm^3/s] for dark matter particles of mass 10 GeV self-annihilating in the b_b channel. Although this constraint is weaker than the ones reported in recent studies using other targets, we note that in the near future, the number of detections of new LSBGs will increase by a few orders of magnitude. We forecast that with the use of the full catalog of soon-to-be-detected LSBGs the constraint will reach cross-section sensitivities of ~ 3*10^-25 [cm^3/s] for dark matter particles with masses less than 10 GeV.

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

Hunting gravitational wave black holes with microlensing

Gravitational microlensing is a powerful tool to search for a population of invisible black holes (BHs) in the Milky Way (MW), including isolated BHs and binary BHs at wide orbits that are complementary to gravitational wave observations. By monitoring highly populated regions of source stars like the MW bulge region, one can pursue microlensing events due to these BHs. We find that if BHs have a Salpeter-like mass function extended beyond $30M_\odot$ and a similar velocity and spatial structure to stars in the Galactic bulge and disk regions, the BH population is a dominant source of the microlensing events at long timescales of the microlensing light curve $\gtrsim 100~$days. This is due to a boosted sensitivity of the microlensing event rate to lens mass, given as $M^2$, for such long-timescale events. A monitoring observation of $2 \times 10^{10}$ stars in the bulge region over 10 years with the Rubin Observatory Legacy Survey of Space and Time (LSST) would enable one to find about $6\times 10^5$ BH microlensing events. We evaluate the efficiency of potential LSST cadences for characterizing the light curves of BH microlensing and find that nearly all events of long timescales can be detected.

preprint2020arXiv

On the wave optics effect on primordial black hole constraints from optical microlensing search

Microlensing of stars, e.g. in the Galactic bulge and Andromeda galaxy (M31), is among the most robust, powerful method to constrain primordial black holes (PBHs) that are a viable candidate of dark matter. If PBHs are in the mass range $M_{\rm PBH} \lower.5ex\hbox{$\; \buildrel < \over \sim \;$} 10^{-10}M_\odot$, its Schwarzschild radius ($r_{\rm Sch}$) becomes comparable with or shorter than optical wavelength ($λ)$ used in a microlensing search, and in this regime the wave optics effect on microlensing needs to be taken into account. For a lensing PBH with mass satisfying $r_{\rm Sch}\sim λ$, it causes a characteristic oscillatory feature in the microlensing light curve, and it will give a smoking gun evidence of PBH if detected, because any astrophysical object cannot have such a tiny Schwarzschild radius. Even in a statistical study, e.g. constraining the abundance of PBHs from a systematic search of microlensing events for a sample of many source stars, the wave effect needs to be taken into account. We examine the impact of wave effect on the PBH constraints obtained from the $r$-band (6210Å) monitoring observation of M31 stars in Niikura et al. (2019), and find that a finite source size effect is dominant over the wave effect for PBHs in the mass range $M_{\rm PBH}\simeq[10^{-11},10^{-10}]M_\odot$. We also discuss that, if a denser-cadence (10~sec), $g$-band monitoring observation for a sample of white dwarfs over a year timescale is available, it would allow one to explore the wave optics effect on microlensing light curve, if it occurs, or improve the PBH constraints in $M_{\rm PBH}\lower.5ex\hbox{$\; \buildrel < \over \sim \;$} 10^{-11}M_\odot$ even from a null detection.

preprint2020arXiv

Proper motion measurements for stars up to $100$ kpc with Subaru HSC and SDSS Stripe 82

We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth and the nearly $12$ years of time baseline between HSC and SDSS enable high-precision measurements of statistical proper motions for stars down to $i\simeq 24$. A validation of our method is demonstrated by the agreement with the $Gaia$ proper motions, to the precision better than $0.1$ mas yr$^{-1}$. To retain the precision, we make a correction of the subtle effects due to the differential chromatic refraction in the SDSS images based on the comparison with the $Gaia$ proper motions against colour of stars, which is validated using the SDSS spectroscopic quasars. Combining with the photometric distance estimates for individual stars based on the precise HSC photometry, we show a significant detection of the net proper motions for stars in each bin of distance out to $100$ kpc. The two-component tangential velocities after subtracting the apparent motions due to our own motion display rich phase-space structures including a clear signature of the Sagittarius stream in the halo region of distance range $[10,\ 35]$ kpc. We also measure the tangential velocity dispersion in the distance range $5-20$ kpc and find that the data are consistent with a constant isotropic dispersion of $80\pm 10\ {\rm km/s}$. More distant stars appear to have random motions with respect to the Galactic centre on average.

preprint2020arXiv

The impact of projection effects on cluster observables: stacked lensing and projected clustering

An optical cluster finder inevitably suffers from projection effects, where it misidentifies a superposition of galaxies in multiple halos along the line-of-sight as a single cluster. Using mock cluster catalogs built from cosmological N-body simulations, we quantify the impact of these projection effects with a particular focus on the observables of interest for cluster cosmology, namely the cluster lensing and the cluster clustering signals. We find that "observed" clusters, i.e. clusters identified by our cluster finder algorithm, exhibit lensing and clustering signals that deviate from expectations based on a statistically isotropic halo model -- while both signals agree with halo model expectations on small scales, they show unexpected boosts on large scales, by up to a factor of 1.2 or 1.4 respectively. We identify the origin of these boosts as the inherent selection bias of optical cluster finders for clusters embedded within filaments aligned with the line-of-sight, and show that a minority ($\sim 30\%$) of such clusters within the entire sample is responsible for this observed boost. We discuss the implications of our results on previous studies of optical cluster, as well as prospects for identifying and mitigating projection effects in future cluster cosmology analyses.

preprint2020arXiv

Towards a non-Gaussian model of redshift space distortions

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments.

preprint2019arXiv

A Search for a Contribution from Axion-Like Particles to the X-Ray Diffuse Background Utilizing the Earth's Magnetic Field

The Axion Like Particle (ALP) is a hypothetical pseudo-scalar particle beyond the Standard Model, with a compelling possible connection to dark matter and early universe physics. ALPs can be converted into photons via interactions with magnetic fields in the universe, i.e., the so-called inverse Primakoff effect. In this paper, we propose a novel method to explore ALP-induced photons from X-ray data obtained from the {\it Suzaku} satellite, arising from a possible interaction of ALPs with the direction-dependent Earth's magnetic field viewed from the satellite. {\it Suzaku} data is suitable for this purpose because its low-altitude Earth orbit result in intrinsically low cosmic-ray background radiation. We study whether the X-ray diffuse background (XDB) spectra estimated from the four deep fields collected over eight years, vary with the integrated Earth's magnetic strength in the direction of each target field at each observation epoch, which amounts to $10^2$ Tm-a value greater than that achieved by terrestrial experiments due to the large coherent length. From the detailed analysis, we did not find evidence of the XDB confidence level spectra having dependence on the Earth's magnetic strength. We obtained 99 % confidence level upper limit on a possible residual contribution to the cosmic X-ray background (CXB) surface brightness to be $1.6\times 10^{-9}~{\rm ergs~s}^{-1}{\rm cm}^{-2}{\rm sr}^{-1}$ normalized at $10^4$ T${}^2$ m${}^2$ in the 2-6 keV range, which corresponds to 6-15 % of the observed CXB brightness, depending on which model of unresolved point sources are used in the interpretation. It is consistent with 80-90 % of the CXB now being resolved into point sources.

preprint2019arXiv

Mitigating the impact of fiber assignment on clustering measurements from deep galaxy redshift surveys

We examine the impact of fiber assignment on clustering measurements from fiber-fed spectroscopic galaxy surveys. We identify new effects which were absent in previous, relatively shallow galaxy surveys such as Baryon Oscillation Spectroscopic Survey . Specifically, we consider deep surveys covering a wide redshift range from z=0.6 to z=2.4, as in the Subaru Prime Focus Spectrograph survey. Such surveys will have more target galaxies than we can place fibers on. This leads to two effects. First, it eliminates fluctuations with wavelengths longer than the size of the field of view, as the number of observed galaxies per field is nearly fixed to the number of available fibers. We find that we can recover the long-wavelength fluctuation by weighting galaxies in each field by the number of target galaxies. Second, it makes the preferential selection of galaxies in under-dense regions. We mitigate this effect by weighting galaxies using the so-called individual inverse probability. Correcting these two effects, we recover the underlying correlation function at better than 1 percent accuracy on scales greater than 10 Mpc/h.

preprint2016arXiv

A New Milky Way Satellite Discovered In The Subaru/Hyper Suprime-Cam Survey

We report the discovery of a new ultra-faint dwarf satellite companion of the Milky Way based on the early survey data from the Hyper Suprime-Cam Subaru Strategic Program. This new satellite, Virgo I, which is located in the constellation of Virgo, has been identified as a statistically significant (5.5 sigma) spatial overdensity of star-like objects with a well-defined main sequence and red giant branch in their color-magnitude diagram. The significance of this overdensity increases to 10.8 sigma when the relevant isochrone filter is adopted for the search. Based on the distribution of the stars around the likely main sequence turn-off at r ~ 24 mag, the distance to Virgo I is estimated as 87 kpc, and its most likely absolute magnitude calculated from a Monte Carlo analysis is M_V = -0.8 +/- 0.9 mag. This stellar system has an extended spatial distribution with a half-light radius of 38 +12/-11 pc, which clearly distinguishes it from a globular cluster with comparable luminosity. Thus, Virgo I is one of the faintest dwarf satellites known and is located beyond the reach of the Sloan Digital Sky Survey. This demonstrates the power of this survey program to identify very faint dwarf satellites. This discovery of VirgoI is based only on about 100 square degrees of data, thus a large number of faint dwarf satellites are likely to exist in the outer halo of the Milky Way.

preprint2016arXiv

Detection of Stacked Filament Lensing Between SDSS Luminous Red Galaxies

We search for the lensing signal of massive filaments between 135,000 pairs of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey. We develop a new estimator that cleanly removes the much larger shear signal of the neighboring LRG halos, relying only on the assumption of spherical symmetry. We consider two models: a thick filament model constructed from ray-tracing simulations for $Λ$CDM model, and a thin filament model which models the filament by a string of halos along the line connecting the two LRGs. We show that the filament lensing signal is in nice agreement with the thick simulation filament, while strongly disfavoring the thin model. The magnitude of the lensing shear due to the filament is below $10^{-4}$. Employing the likelihood ratio test, we find a 4.5$σ$ significance for the detection of the filament lensing signal, corresponding to a null hypothesis fluctuation probability of $3 \times 10^{-6}$. We also carried out several null tests to verify that the residual shear signal from neighboring LRGs and other shear systematics are minimized.

preprint2016arXiv

Detection of the Splashback Radius and Halo Assembly bias of Massive Galaxy Clusters

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$σ$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $Λ$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

preprint2016arXiv

Evidence of Halo Assembly Bias in Massive Clusters

We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$σ$ in the weak lensing and 4.4$σ$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

preprint2016arXiv

Prime Focus Spectrograph (PFS) for the Subaru Telescope: Overview, recent progress, and future perspectives

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

preprint2016arXiv

Separate Universe Consistency Relation and Calibration of Halo Bias

Linear halo bias is the response of dark matter halo number density to a long wavelength fluctuation in the dark matter density. Using abundance matching between separate universe simulations which absorb the latter into a change in the background, we test the consistency relation between the change in a one point function, the halo mass function, and a two point function, the halo-matter cross correlation in the long wavelength limit. We find excellent agreement between the two at the $1-2\%$ level for average halo biases between $1 \lesssim \bar b_1 \lesssim 4$ and no statistically significant deviations at the $4-5\%$ level out to $\bar b_1 \approx 8$. Halo bias inferred assuming instead a universal mass function is significantly different and inaccurate at the 10\% level or more. The separate universe technique provides a way of calibrating linear halo bias efficiently for even highly biased rare halos in the $Λ$CDM model. Observational violation of the consistency relation would indicate new physics, e.g.~in the dark matter, dark energy or primordial non-Gaussianity sectors.

preprint2016arXiv

Subaru high-z exploration of low-luminosity quasars (SHELLQs). I. Discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

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.

preprint2015arXiv

Properties of Weak Lensing Clusters Detected on Hyper Suprime-Cam 2.3 Square Degree Field

We present properties of moderately massive clusters of galaxies detected by the newly developed Hyper Suprime-Cam on the Subaru telescope using weak gravitational lensing. Eight peaks exceeding a S/N ratio of 4.5 are identified on the convergence S/N map of a 2.3 square degree field observed during the early commissioning phase of the camera. Multi-color photometric data is used to generate optically selected clusters using the CAMIRA algorithm. The optical cluster positions were correlated with the peak positions from the convergence map. All eight significant peaks have optical counterparts. The velocity dispersion of clusters are evaluated by adopting the Singular Isothemal Sphere (SIS) fit to the tangential shear profiles, yielding virial mass estimates, M500c, of the clusters which range from 2.7x10^13 to 4.4x10^14 solar mass. The number of peaks is considerably larger than the average number expected from LambdaCDM cosmology but this is not extremely unlikely if one takes the large sample variance in the small field into account. We could, however, safely argue that the peak count strongly favours the recent Planck result suggesting high sigma8$value of 0.83. The ratio of stellar mass to the dark matter halo mass shows a clear decline as the halo mass increases. If the gas mass fraction, fg, in halos is universal, as has been suggested in the literature, the observed baryon mass in stars and gas shows a possible deficit compared with the total baryon density estimated from the baryon oscillation peaks in anisotropy of the cosmic microwave background.

preprint2015arXiv

The Weak Lensing Signal and the Clustering of BOSS Galaxies I: Measurements

A joint analysis of the clustering of galaxies and their weak gravitational lensing signal is well-suited to simultaneously constrain the galaxy-halo connection as well as the cosmological parameters by breaking the degeneracy between galaxy bias and the amplitude of clustering signal. In a series of two papers, we perform such an analysis at the highest redshift ($z\sim0.53$) in the literature using CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Eleventh Data Release (SDSS-III/BOSS DR11) catalog spanning 8300~deg$^2$. In this paper, we present details of the clustering and weak lensing measurements of these galaxies. We define a subsample of 400,916 CMASS galaxies based on their redshifts and stellar mass estimates so that the galaxies constitute an approximately volume-limited and similar population over the redshift range $0.47\le z\le 0.59$. We obtain a signal-to-noise ratio $S/N\simeq 56$ for the galaxy clustering measurement. We also explore the redshift and stellar mass dependence of the clustering signal. For the weak lensing measurement, we use existing deeper imaging data from the CFHTLS with publicly available shape and photometric redshift catalogs from CFHTLenS, but only in a 105~deg$^2$ area which overlaps with BOSS. This restricts the lensing measurement to only 5,084 CMASS galaxies. After careful systematic tests, we find a highly significant detection of the CMASS weak lensing signal, with total $S/N\simeq 26$. These measurements form the basis of the halo occupation distribution and cosmology analysis presented in More et al. (Paper II).

preprint2015arXiv

The Weak Lensing Signal and the Clustering of BOSS Galaxies II: Astrophysical and Cosmological Constraints

We perform a joint analysis of the abundance, the clustering and the galaxy-galaxy lensing signal of galaxies measured from Data Release 11 of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS III-BOSS) in our companion paper, Miyatake et al. (2014). The lensing signal was obtained by using the shape catalog of background galaxies from the Canada France Hawaii Telescope Legacy Survey, which was made publicly available by the CFHTLenS collaboration, with an area overlap of about 105 deg$^2$. We analyse the data in the framework of the halo model in order to fit halo occupation parameters and cosmological parameters ($Ω_{\rm m}$ and $σ_8$) to these observables simultaneously, and thus break the degeneracy between galaxy bias and cosmology. Adopting a flat $Λ$CDM cosmology with priors on $Ω_b h^2$, $n_{\rm s}$ and $h$ from the analysis of WMAP 9-year data, we obtain constraints on the stellar mass-halo mass relation of galaxies in our sample. Marginalizing over the halo occupation distribution parameters and a number of other nuisance parameters in our model, we obtain $Ω_{\rm m}=0.310^{+0.019}_{-0.020}$ and $σ_8=0.785^{+0.044}_{-0.044}$ (68% confidence). We demonstrate the robustness of our results with respect to sample selection and a variety of systematics such as the halo off-centering effect and possible incompleteness in our sample. Our constraints are consistent, complementary and competitive with those obtained using other independent probes of these cosmological parameters. The cosmological analysis is the first of its kind to be performed at a redshift as high as $0.53$.

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

Joint analysis of cluster number counts and weak lensing power spectrum to correct for the super-sample covariance

A coherent over- or under-density contrast across a finite survey volume causes an upward- or downward-fluctuation in the observed number of halos. This fluctuation in halo number adds a significant co-variant scatter in the observed amplitudes of weak lensing power spectrum at nonlinear, small scales -- the so-called super-sample variance or the halo sample variance. In this paper, we show that by measuring both the number counts of clusters and the power spectrum in the same survey region, we can mitigate this loss of information and significantly enhance the scientific return from the upcoming surveys. First, using the halo model approach, we derive the cross-correlation between the halo number counts and the weak lensing power spectrum, taking into account the super-sample covariance effect, which well matches the distributions measured from 1000 realizations for a Lambda-dominated cold dark matter model. Then we show that adding the observed number counts of massive halos with M> 10^14 Ms/h can significantly improve the information content of weak lensing power spectrum, almost recovering the Gaussian information up to lmax~1000, if the average mass profiles of the massive halos are known, which can be estimated from stacked lensing. When combined with the halo number counts for M > 3 or 1 x 10^14 Ms/h, the improvement is up to a factor of 1.4 or 2 at lmax~1000 --2000, equivalent to a factor 2 or 4 times larger survey volume, compared to the power spectrum measurement alone.

preprint2014arXiv

Joint likelihood function of cluster counts and $n$-point correlation functions: Improving their power through including halo sample variance

Naive estimates of the statistics of large scale structure and weak lensing power spectrum measurements that include only Gaussian errors exaggerate their scientific impact. Non-linear evolution and finite volume effects are both significant sources of non-Gaussian covariance that reduce the ability of power spectrum measurements to constrain cosmological parameters. Using a halo model formalism, we derive an intuitive understanding of the various contributions to the covariance and show that our analytical treatment agrees with simulations. This approach enables an approximate derivation of a joint likelihood for the cluster number counts, the weak lensing power spectrum and the bispectrum. We show that this likelihood is a good description of the ray-tracing simulation. Since all of these observables are sensitive to the same finite volume effects and contain information about the non-linear evolution, a combined analysis recovers much of the "lost" information and obviates the non-Gaussian covariance. For upcoming weak lensing surveys, we estimate that a joint analysis of power spectrum, number counts and bispectrum will produce an improvement of about $30-40\%$ in determinations of the matter density and the scalar amplitude. This improvement is equivalent to doubling the survey area.

preprint2014arXiv

Power Spectrum Super-Sample Covariance

We provide a simple, unified approach to describing the impact of super-sample covariance, or beat coupling, on power spectrum estimation in a finite-volume survey. For a wide range of survey volumes, the sample variance that arises from modes that are larger than the survey dominates the covariance of power spectrum estimators for modes much smaller than the survey. The deeply nonlinear version of this effect is known as halo sample variance. We show that all variants are unified by the matter trispectrum of squeezed configurations and that such configurations obey a consistency relation which relates them to the response of the power spectrum to a change in the background density. Our method also applies to statistics that are based on radial projections of the density field such as weak lensing shear. While we use the halo model for an analytic description to expose the nature of the effect, the consistency description enables an accurate calibration of the full effect directly from simulations. It also suggests that super-sample covariance may be viewed as an additional interesting signal rather than excess noise.

preprint2014arXiv

Statistical challenges in weak lensing cosmology

Cosmological weak lensing is the powerful probe of cosmology. Here we address one of the most fundamental, statistical questions inherent in weak lensing cosmology: whether or not we can recover the initial Gaussian information content of large-scale structure by combining the weak lensing observables, here focused on the weak lensing power spectrum and bispectrum. To address this question we fully take into account correlations between the power spectra of different multipoles and the bispectra of different triangle configurations, measured from a finite area survey. In particular we show that super-survey modes whose length scale is larger than or comparable with the survey size cause significant sample variance in the weak lensing correlations via the mode-coupling with sub-survey modes due to nonlinear gravitational clustering -- the so-called super-sample variance. In this paper we discuss the origin of the super-sample variance and then study the information content inherent in the weak lensing correlation functions up to three-point level.

preprint2014arXiv

Super-Sample Covariance in Simulations

Using separate universe simulations, we accurately quantify super-sample covariance (SSC), the typically dominant sampling error for matter power spectrum estimators in a finite volume, which arises from the presence of super survey modes. By quantifying the power spectrum response to a background mode, this approach automatically captures the separate effects of beat coupling in the quasilinear regime, halo sample variance in the nonlinear regime and a new dilation effect which changes scales in the power spectrum coherently across the survey volume, including the baryon acoustic oscillation scale. It models these effects at typically the few percent level or better with a handful of small volume simulations for any survey geometry compared with directly using many thousands of survey volumes in a suite of large volume simulations. The stochasticity of the response is sufficiently small that in the quasilinear regime, SSC can be alternately included by fitting the mean density in the volume with these fixed templates in parameter estimation. We also test the halo model prescription and find agreement typically at better than the 10% level for the response.

preprint2014arXiv

Super-Sample Signal

When extracting cosmological information from power spectrum measurements, we must consider the impact of super-sample density fluctuations whose wavelengths are larger than the survey scale. These modes contribute to the mean density fluctuation $δ_b$ in the survey and change the power spectrum in the same way as a change in the cosmological background. They can be simply included in cosmological parameter estimation and forecasts by treating $δ_b$ as an additional cosmological parameter enabling efficient exploration of its impact. To test this approach, we consider here an idealized measurement of the matter power spectrum itself in the $Λ$CDM cosmology though our techniques can readily be extended to more observationally relevant statistics or other parameter spaces. Using sub-volumes of large-volume $N$-body simulations for power spectra measured with respect to either the global or local mean density, we verify that the minimum variance estimator of $δ_b$ is both unbiased and has the predicted variance. Parameter degeneracies arise since the response of the matter power spectrum to $δ_b$ and cosmological parameters share similar properties in changing the growth of structure and dilating the scale of features especially in the local case. For matter power spectrum measurements, these degeneracies can lead in certain cases to substantial error degradation and motivates future studies of specific cosmological observables such as galaxy clustering and weak lensing statistics with these techniques.

preprint2013arXiv

Baryon acoustic oscillations with the cross-correlation of spectroscopic and photometric samples

The baryon acoustic oscillation (BAO) experiment requires a sufficiently dense sampling of large-scale structure tracers with spectroscopic redshift, which is observationally expensive especially at high redshifts $z\simgt 1$. Here we present an alternative route of the BAO analysis that uses the cross-correlation of sparse spectroscopic tracers with a much denser photometric sample, where the spectroscopic tracers can be quasars or bright, rare galaxies that are easier to access spectroscopically. We show that measurements of the cross-correlation as a function of the transverse comoving separation rather than the angular separation avoid a smearing of the BAO feature without mixing the different scales at different redshifts in the projection, even for a wide redshift slice $Δz\simeq 1$. The bias, scatter, and catastrophic redshift errors of the photometric sample affect only the overall normalization of the cross-correlation which can be marginalized over when constraining the angular diameter distance. As a specific example, we forecast an expected accuracy of the BAO geometrical test via the cross-correlation of the SDSS and BOSS spectroscopic quasar sample with a dense photometric galaxy sample that is assumed to have a full overlap with the SDSS/BOSS survey region. We show that this cross-correlation BAO analysis allows us to measure the angular diameter distances to a fractional accuracy of about 10% at each redshift bin over $1\simlt z\simlt 3$, if the photometric redshift errors of the galaxies, $σ_z/(1+z)$, are better than 10-20% level.

preprint2013arXiv

Cosmological parameters from weak lensing power spectrum and bispectrum tomography: including the non-Gaussian errors

We re-examine a genuine power of weak lensing bispectrum tomography for constraining cosmological parameters, when combined with the power spectrum tomography, based on the Fisher information matrix formalism. To account for the full information at two- and three-point levels, we include all the power spectrum and bispectrum information built from all-available combinations of tomographic redshift bins, multipole bins and different triangle configurations over a range of angular scales (up to lmax=2000 as our fiducial choice). For the parameter forecast, we use the halo model approach in Kayo, Takada & Jain (2013) to model the non-Gaussian error covariances as well as the cross-covariance between the power spectrum and the bispectrum, including the halo sample variance or the nonlinear version of beat-coupling. We find that adding the bispectrum information leads to about 60% improvement in the dark energy figure-of-merit compared to the lensing power spectrum tomography alone, for three redshift-bin tomography and a Subaru-type survey probing galaxies at typical redshift of zs~1. The improvement is equivalent to a 1.6 larger survey area. Thus our results show that the bispectrum or more generally any three-point correlation based statistics carries complementary information on cosmological parameters to the power spectrum. However, the improvement is modest compared to the previous claim derived using the Gaussian error assumption, and therefore our results imply less additional information in even higher-order moments such as the four-point correlation function.

preprint2013arXiv

Extragalactic Science, Cosmology and Galactic Archaeology with the Subaru Prime Focus Spectrograph (PFS)

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

preprint2013arXiv

Information content of weak lensing power spectrum and bispectrum: including the non-Gaussian error covariance matrix

We address the amount of information in the non-Gaussian regime of weak lensing surveys by modelling all relevant covariances of the power spectra and bispectra, using 1000 ray-tracing simulation realizations for a Lambda-CDM model and an analytical halo model. We develop a formalism to describe the covariance matrices of power spectra and bispectra of all triangle configurations. In addition to the known contributions which extend up to six-point correlation functions, we propose a new contribution `the halo sample variance (HSV)' arising from the coupling of the lensing Fourier modes with large-scale mass fluctuations on scales comparable with the survey region via halo bias theory. We show that the model predictions are in good agreement with the simulation once we take the HSV into account. The HSV gives a dominant contribution to the covariance matrices at multipoles l > 10^3, which arises from massive haloes with a mass of > 10^14 solar mass and at relatively low redshifts z < 0.4. Since such haloes are easily identified from a multi-colour imaging survey, the effect can be estimated from the data. By adding the bispectrum to the power spectrum, the total information content or the cumulative signal-to-noise ratio up to a certain maximum multipole of a few 10^3 is improved by 20--50 per cent, which is equivalent to a factor of 1.4--2.3 larger survey area for the power spectrum measurement alone. However, it is still smaller than the case of a Gaussian field by a factor of 3 mostly due to the HSV. Thus bispectrum measurements are useful for cosmology, but using information from upcoming surveys requires that non-Gaussian covariances are carefully estimated.

preprint2013arXiv

Modeling the Phase-Space Distribution around Massive Halos

The comparison between dynamical mass and lensing mass provides a targeted test for a wide range of modified gravity models. In our previous paper we showed, through numerical simulations, that the measurement of the line-of-sight velocity dispersion around stacked massive clusters whose lensing masses are known allows for stringent constraints on modified gravity on scales of 2 - 15 Mpc/h. In this work we develop a semi-analytical approach based on the halo model to describe the phase-space distribution and the line-of-sight velocity dispersion for different tracers. The model distinguishes contributions from the halo pairwise velocity and the virial velocity within halos. We also discuss observational complications, in particular the contribution from Hubble flow, and show how our model can incorporate these complications. We then incorporate the effects of modified gravity (specifically, f(R) and braneworld models), and show that the model predictions are in excellent agreement with modified gravity simulations. More broadly, the phase-space distribution provides a sensitive test of our understanding of hierarchical structure formation when confronted with observations via this model.

preprint2013arXiv

Novel Probes of Gravity and Dark Energy

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

preprint2013arXiv

Perturbation theory for nonlinear halo power spectrum: the renormalized bias and halo bias

We revisit an analytical model to describe the halo-matter cross-power spectrum and the halo auto-power spectrum in the weakly nonlinear regime, by combining the perturbation theory (PT) for matter clustering, the local bias model, and the halo bias. Nonlinearities in the power spectra arise from the nonlinear clustering of matter as well as the nonlinear relation between the matter and halo density fields. By using the "renormalization" approach, we express the nonlinear power spectra by a sum of the two contributions: the nonlinear matter power spectrum with the effective linear bias parameter, and the higher-order PT spectra having the halo bias parameters as the coefficients. The halo auto-power spectrum includes the residual shot noise contamination that needs to be treated as additional free parameter. The term(s) of the higher-order PT spectra and the residual shot noise cause a scale-dependent bias function relative to the nonlinear matter power spectrum in the weakly nonlinear regime. We show that the model predictions are in good agreement with the spectra measured from a suit of high-resolution $N$-body simulations up to $k\simeq 0.2 h$/Mpc at $z=0.35$, for different halo mass bins.

preprint2013arXiv

The Synergy between Weak Lensing and Galaxy Redshift Surveys

We study the complementarity of weak lensing (WL) and spectroscopic galaxy clustering (GC) surveys, by forecasting dark energy and modified gravity constraints for three upcoming survey combinations: SuMIRe (Subaru Measurement of Images and Redshifts, the combination of the Hyper Suprime-Cam lensing survey and the Prime Focus Spectrograph redshift survey), EUCLID and WFIRST. From the WL surveys, we take into account both the shear and clustering of the source galaxies and from the GC surveys, we use the three-dimensional clustering of spectroscopic galaxies, including redshift space distortions. A CMB prior is included in all cases. Focusing on the large-scale, two-point function information, we find strong synergy between the two probes. The dark energy figure of merit from WL+GC is up to a factor ~2.5 larger than from either probe alone. Considering modified gravity, if the growth factor f(z) is treated as a free function, it is very poorly constrained by WL or GC alone, but can be measured at the few percent level by the combination of the two. On the other hand, for cosmological constraints derived from (angular) power spectra and considering statistical errors only, it hardly matters whether the surveys overlap on the sky or not. For instance, the dark energy figure of merit for overlapping surveys is at most ~20 % better than in the disjoint case. This modest gain can be traced to the fundamental fact that only a small fraction of the total number of modes sampled by the GC survey (or by the WL survey) contributes to the cross-correlations between WL and GC.

preprint2013arXiv

Where are the Luminous Red Galaxies (LRGs)? Using correlation measurements and lensing to relate LRGs to dark matter halos

Nonlinear redshift-space distortions, the Finger-of-God (FoG) effect, can complicate the interpretation of the galaxy power spectrum. Here, we demonstrate the method proposed by Hikage et al. (2012) to use complimentary observations to directly constrain this effect on the data. We use catalogs of Luminous Red Galaxies (LRGs) and photometric galaxies from the SDSS DR7 to measure the redshift-space power spectrum of LRGs, the cross-correlation of LRGs with the shapes of background photometric galaxies (galaxy-galaxy weak lensing), and the projected cross-correlation of LRGs with photometric galaxies having similar photometric redshifts to the LRG spectroscopic redshift. All of these measurements use a reconstructed halo field. While we use the position of each LRG for single LRG systems, we compare the measurements using different halo-center proxies for multiple-LRG systems (4.5 per cent of all the halos): the brightest LRG position (BLRG), the faintest LRG position (FLRG) and their arithmetical mean position (Mean), respectively, in each system. We find significant differences in the measured correlations of different centers, showing consistent off-centering effects in the three observables. By comparing the measurements with a halo model that treats the satellite photometric galaxies as being distributed according to a generalized NFW profile, we find that about 40 (70) per cent of BLRGs (FLRGs) are off-centered satellite galaxies in the multiple-LRG systems. The satellite LRGs have typical off-centering radius of about 400 kpc/h, and velocity dispersion of about 500 km/s in host halos with a mean mass of 1.6x10^14 Ms/h. We show that, if LRGs in the single LRG systems have similar offsets, the residual FoG contamination in the LRG power spectrum can be significant at k>0.1 h/Mpc, which may cause a bias in cosmological parameters such as the neutrino mass.

preprint2012arXiv

Cross-correlating the Thermal Sunyaev-Zel'dovich Effect and the Distribution of Galaxy Clusters

We present the analytical formulas, derived based on the halo model, to compute the cross-correlation between the thermal Sunyaev-Zel'dovich (SZ) effect and the distribution of galaxy clusters. By binning the clusters according to their redshifts and masses, this cross-correlation, the so-called stacked SZ signal, reveals the average SZ profile around the clusters. The stacked SZ signal is obtainable from a joint analysis of an arcminute-resolution cosmic microwave background (CMB) experiment and an overlapping optical survey, which allows for detection of the SZ signals for clusters whose masses are below the individual cluster detection threshold. We derive the error covariance matrix for measuring the stacked SZ signal, and then forecast for its detection from ongoing and forthcoming combined CMB-optical surveys. We find that, over a wide range of mass and redshift, the stacked SZ signal can be detected with a significant signal to noise ratio (total S/N \gsim 10), whose value peaks for the clusters with intermediate masses and redshifts. Our calculation also shows that the stacking method allows for probing the clusters' SZ profiles over a wide range of scales, even out to projected radii as large as the virial radius, thereby providing a promising way to study gas physics at the outskirts of galaxy clusters.

preprint2012arXiv

LoCuSS: The Sunyaev-Zel'dovich Effect and Weak Lensing Mass Scaling Relation

We present the first weak-lensing-based scaling relation between galaxy cluster mass, M_wl, and integrated Compton parameter Y_sph. Observations of 18 galaxy clusters at z~0.2 were obtained with the Subaru 8.2-m telescope and the Sunyaev-Zel'dovich Array. The M_wl-Y_sph scaling relations, measured at Delta=500, 1000, and 2500 rho_c, are consistent in slope and normalization with previous results derived under the assumption of hydrostatic equilibrium (HSE). We find an intrinsic scatter in M_wl at fixed Y_sph of 20%, larger than both previous measurements of M_HSE-Y_sph scatter as well as the scatter in true mass at fixed Y_sph found in simulations. Moreover, the scatter in our lensing-based scaling relations is morphology dependent, with 30-40% larger M_wl for undisturbed compared to disturbed clusters at the same Y_sph at r_500. Further examination suggests that the segregation may be explained by the inability of our spherical lens models to faithfully describe the three-dimensional structure of the clusters, in particular, the structure along the line-of-sight. We find that the ellipticity of the brightest cluster galaxy, a proxy for halo orientation, correlates well with the offset in mass from the mean scaling relation, which supports this picture. This provides empirical evidence that line-of-sight projection effects are an important systematic uncertainty in lensing-based scaling relations.

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

Subaru weak-lensing measurement of a z = 0.81 cluster discovered by the Atacama Cosmology Telescope Survey

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\barρ_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the Λ-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

preprint2012arXiv

Testing Gravity with the Stacked Phase Space around Galaxy Clusters

In General Relativity, the average velocity field of dark matter around galaxy clusters is uniquely determined by the mass profile. The latter can be measured through weak lensing. We propose a new method of measuring the velocity field (phase space density) by stacking redshifts of surrounding galaxies from a spectroscopic sample. In combination with lensing, this yields a direct test of gravity on scales of 1-30 Mpc. Using N-body simulations, we show that this method can improve upon current constraints on f(R) and DGP model parameters by several orders of magnitude when applied to upcoming imaging and redshift surveys.

preprint2012arXiv

The impact of massive neutrinos on the abundance of massive clusters

We study the spherical, top-hat collapse model for a mixed dark matter model including cold dark matter (CDM) and massive neutrinos of mass scales ranging from m_nu= 0.05 to a few 0.1eV, the range of lower- and upper-bounds implied from the neutrino oscillation experiments and the cosmological constraints. To develop this model, we properly take into account relative differences between the density perturbation amplitudes of different components (radiation, baryon, CDM and neutrinos) around the top-hat CDM overdensity region assuming the adiabatic initial conditions. Furthermore, we solve the linearized Boltzmann hierarchy equations to obtain time evolution of the lineariezed neutrino perturbations, yet including the effect of nonlinear gravitational potential due to the nonlinear CDM and baryon overdensities in the late stage. We find that the presence of massive neutrinos slows down the collapse of CDM (plus baryon) overdensity, however, that the neutrinos cannot fully catch up with the the nonlinear CDM perturbation due to its large free-streaming velocity for the ranges of neutrino masses and halo masses we consider. We find that, just like CDM models, the collapse time of CDM overdensity is well monitored by the linear-theory extrapolated overdensity of CDM plus baryon perturbation, smoothed with a given halo mass scale, if taking into account the suppression effect of the massive neutrinos on the linear growth rate. Using these findings, we argue that the presence of massive neutrinos of mass scales 0.05 or 0.1eV may cause a significant decrease in the abundance of massive halos compared to the model without the massive neutrinos; e.g., by 25% or factor 2, respectively, for halos with 10^15Ms and at z=1.

preprint2012arXiv

Using galaxy-galaxy weak lensing measurements to correct the Finger-of-God

For decades, cosmologists have been using galaxies to trace the large-scale distribution of matter. At present, the largest source of systematic uncertainty in this analysis is the challenge of modeling the complex relationship between galaxy redshift and the distribution of dark matter. If all galaxies sat in the centers of halos, there would be minimal Finger-of-God (FoG) effects and a simple relationship between the galaxy and matter distributions. However, many galaxies, even some of the luminous red galaxies (LRGs), do not lie in the centers of halos. Because the galaxy-galaxy lensing is also sensitive to the off-centered galaxies, we show that we can use the lensing measurements to determine the amplitude of this effect and to determine the expected amplitude of FoG effects. We develop an approach for using the lensing data to model how the FoG suppresses the power spectrum amplitudes and show that the current data implies a 30% suppression at wavenumber k=0.2h/Mpc. Our analysis implies that it is important to complement a spectroscopic survey with an imaging survey with sufficient depth and wide field coverage. Joint imaging and spectroscopic surveys allow a robust, unbiased use of the power spectrum amplitude information: it improves the marginalized error of growth rate fg=dln D/dln a by up to a factor of 2 over a wide range of redshifts z<1.4. We also find that the dark energy equation-of-state parameter, w0, and the neutrino mass, fnu, can be unbiasedly constrained by combining the lensing information, with an improvement of 10--25% compared to a spectroscopic survey without lensing calibration.

preprint2011arXiv

Dark Energy from the log-transformed convergence field

A logarithmic transform of the convergence field improves `the information content', ie., the overall precision associated with the measurement of the amplitude of the convergence power spectrum by improving the covariance matrix properties. The translation of this improvement in the information content to that in cosmological parameters, such as those associated with dark energy, requires knowing the sensitivity of the log-transformed field to those cosmological parameters. In this paper we use N-body simulations with ray tracing to generate convergence fields at multiple source redshifts as a function of cosmology. The gain in information associated with the log-transformed field does lead to tighter constraints on dark energy parameters, but only if shape noise is neglected. The presence of shape noise quickly diminishes the advantage of the log mapping, more quickly than we would expect based on the information content. With or without shape noise, using a larger pixel size allows for a more efficient log-transformation.

preprint2011arXiv

Likelihood reconstruction method of real-space density and velocity power spectra from a redshift galaxy survey

We develop a maximum likelihood based method of reconstructing band powers of the density and velocity power spectra at each wavenumber bins from the measured clustering features of galaxies in redshift space, including marginalization over uncertainties inherent in the Fingers-of-God (FoG) effect. The reconstruction can be done assuming that the density and velocity power spectra depend on the redshift-space power spectrum having different angular modulations of mu with mu^{2n} (n=0,1,2) and that the model FoG effect is given as a multiplicative function in the redshift-space spectrum. By using N-body simulations and the halo catalogs, we test our method by comparing the reconstructed power spectra with the simulations. For the spectrum of mu^0 or equivalently the density power spectrum P_dd(k), our method recovers the amplitudes to a few percent accuracies up to k=0.3 h/Mpc for both dark matter and halos. For the power spectrum of mu^2, which is equivalent to the density-velocity spectrum P_dv(k) in the linear regime, our method can recover the input power spectrum for dark matter up to k=0.2 h/Mpc and at both z=0 and 1, if using the adequate FoG model. However, for the halo spectrum, the reconstructed spectrum shows greater amplitudes than the simulation P_dv(k). We argue that the disagreement is ascribed to nonlinearity effect that arises from the cross-bispectra of density and velocity perturbations. Using the perturbation theory, we derive the nonlinear correction term, and find that the leading-order correction term is proportional to mu^2 and increases the mu^2-power spectrum amplitudes at larger k, at lower redshifts and for more massive halos. We find that adding the nonlinearity correction term to the simulation P_dv(k) can fairly well reproduce the reconstructed P_dv(k) for halos up to k~0.2 h/Mpc.

preprint2011arXiv

LoCuSS: Subaru Weak Lensing Study of 30 Galaxy Clusters

(Abridged) We use Subaru data to conduct a detailed weak-lensing study of the dark matter distribution in a sample of 30 X-ray luminous galaxy clusters at 0.15<z<0.3. A weak-lensing signal is detected at high statistical significance in each cluster, the total detection S/N ranging from 5 to 13. In this paper we concentrate on fitting spherical models to the tangential distortion profiles of the clusters. When the models are fitted to the clusters individually, we are unable to discriminate statistically between SIS and NFW models. However when the tangential distortion profiles of the individual clusters are combined, and models fitted to the stacked profile, the SIS model is rejected at 6- and 11-sigma, respectively, for low- and high-mass bins. We also use the individual cluster NFW model fits to investigate the relationship between cluster mass (M_vir) and concentration (c_vir), finding an anti-correlation of c_vir and M_vir. The best-fit c_vir-M_vir relation is: c_vir(M_vir) propto M_vir^{-alpha} with alpha=0.41+/-0.19 -- i.e. a non-zero slope is detected at 2sigma significance. We then investigate the optimal radius within which to measure cluster mass, finding that the typical fractional errors are improved to sigma(M_Delta)/M_Delta ~ 0.1-0.2 for cluster masses at higher over-densities Delta=500-2000, from 0.2-0.3 for the virial over-density (~110). Further comparisons between mass measurements based on spherical model fitting and the model-independent aperture mass method reveal that the 2D aperture mass enclosed within a cylinder of a given aperture radius is systematically greater than the 3D spherical mass obtained from NFW model fitting: M_2D/M_3D= 1.34 and 1.40 for Delta=500 and 110, respectively. The amplitude of this effect agrees well with that predicted by integrating the NFW model along the line-of-sight.

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

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.

preprint2010arXiv

A Clipping Method to Mitigate the Impact of Catastrophic Photometric Redshift Errors on Weak Lensing Tomography

We use the mock catalog of galaxies, constructed based on the COSMOS galaxy catalog including information on photometric redshifts (photo-z) and SED types of galaxies, in order to study how to define a galaxy subsample suitable for weak lensing tomography feasible with optical (and NIR) multi-band data. Since most of useful cosmological information arises from the sample variance limited regime for upcoming lensing surveys, a suitable subsample can be obtained by discarding a large fraction of galaxies that have less reliable photo-z estimations. We develop a method to efficiently identify photo-z outliers by monitoring the width of posterior likelihood function of redshift estimation for each galaxies. This clipping method may allow to obtain clean tomographic redshift bins (here three bins considered) that have almost no overlaps between different bins, by discarding more than 70% galaxies of ill-defined photo-z's corresponding to the number densities of remaining galaxies less than 20 per square arcminutes for a Subaru-type deep survey. Restricting the ranges of magnitudes and redshifts and/or adding near infrared data help obtain a cleaner redshift binning. By using the Fisher information matrix formalism, we propagate photo-z errors into biases in the dark energy equation of state parameter w. We found that, by discarding most of ill-defined photo-z galaxies, the bias in w can be reduced to the level comparable to the marginalized statistical error, however, the residual, small systematic bias remains due to asymmetric scatters around the relation between photometric and true redshifts. We also use the mock catalog to estimate the cumulative signal-to-noise (S/N) ratios for measuring the angular cross-correlations of galaxies between finner photo-z bins, finding the higher S/N values for photo-z bins including photo-z outliers.

preprint2010arXiv

A dipole anisotropy of galaxy distribution: Does the CMB rest-frame exist in the local universe?

The peculiar motion of the Earth causes a dipole anisotropy modulation in the distant galaxy distribution due to the aberration effect. However, the amplitude and angular direction of the effect is not necessarily the same as those of the cosmic microwave background (CMB) dipole anisotropy due to the growth of cosmic structures. In other words exploring the aberration effect may give us a clue to the horizon-scale physics perhaps related to the cosmic acceleration. In this paper we develop a method to explore the dipole angular modulation from the pixelized galaxy data on the sky properly taking into account the covariances due to the shot noise and the intrinsic galaxy clustering contamination as well as the partial sky coverage. We applied the method to the galaxy catalogs constructed from the Sloan Digital Sky Survey (SDSS) Data Release 6 data. After constructing the four galaxy catalogs that are different in the ranges of magnitudes and photometric redshifts to study possible systematics, we found that the most robust sample against systematics indicates no dipole anisotropy in the galaxy distribution. This finding is consistent with the expectation from the concordance Lambda-dominated cold dark matter model. Finally we argue that an almost full-sky galaxy survey such as LSST may allow for a significant detection of the aberration effect of the CMB dipole having the precision of constraining the angular direction to ~ 20 degrees in radius. Assuming a hypothetical LSST galaxy survey, we find that this method can confirm or reject the result implied from a stacked analysis of the kinetic Sunyaev-Zel'dovich effect of X-ray luminous clusters in Kashlinsky et al. (2008,2009) if the implied cosmic bulk flow is not extended out to the horizon.

preprint2010arXiv

Combining cluster observables and stacked weak lensing to probe dark energy: Self-calibration of systematic uncertainties

We develop a new method of combining cluster observables (number counts and cluster-cluster correlation functions) and stacked weak lensing signals of background galaxy shapes, both of which are available in a wide-field optical imaging survey. Assuming that the clusters have secure redshift estimates, we show that the joint experiment enables a self-calibration of important systematic errors including the source redshift uncertainty and the cluster mass-observable relation, by adopting a single population of background source galaxies for the lensing analysis. It allows us to use the relative strengths of stacked lensing signals at different cluster redshifts for calibrating the source redshift uncertainty, which in turn leads to accurate measurements of the mean cluster mass in each bin. In addition, our formulation of stacked lensing signals in Fourier space simplifies the Fisher matrix calculations, as well as the marginalization over the cluster off-centering effect, the most significant uncertainty in stacked lensing. We show that upcoming wide-field surveys yield stringent constraints on cosmological parameters including dark energy parameters, without any priors on nuisance parameters that model systematic uncertainties. Specifically, the stacked lensing information improves the dark energy FoM by a factor of 4, compared to that from the cluster observables alone. The primordial non-Gaussianity parameter can also be constrained with a level of f_NL~10. In this method, the mean source redshift is well calibrated to an accuracy of 0.1 in redshift, and the mean cluster mass in each bin to 5-10% accuracies, which demonstrates the success of the self-calibration of systematic uncertainties from the joint experiment. (Abridged)

preprint2010arXiv

Direct measurement of dark matter halo ellipticity from two-dimensional lensing shear maps of 25 massive clusters

We present new measurements of dark matter distributions in 25 X-ray luminous clusters by making a full use of the two-dimensional (2D) weak lensing signals obtained from high-quality Subaru/Suprime-Cam imaging data. Our approach to directly compare the measured lensing shear pattern with elliptical model predictions allows us to extract new information on the mass distributions of individual clusters, such as the halo ellipticity and mass centroid. We find that these parameters on the cluster shape are little degenerate with cluster mass and concentration parameters. By combining the 2D fitting results for a subsample of 18 clusters, the elliptical shape of dark matter haloes is detected at 7σsignificance level. The mean ellipticity is found to be e = 0.46 \pm 0.04 (1σ), which is in excellent agreement with the standard collisionless CDM model prediction. The mass centroid can be constrained with a typical accuracy of ~20" (~50 kpc/h) in radius for each cluster with some significant outliers, enabling to assess one of the most important systematic errors inherent in the stacked cluster weak lensing technique, the mass centroid uncertainty. In addition, the shape of the dark mass distribution is found to be only weakly correlated with that of the member galaxy distribution. We carefully examine possible sources of systematic errors in our measurements, finding none of them to be significant. Our results demonstrate the power of high-quality imaging data for exploring the detailed spatial distribution of dark matter (Abridged).

preprint2010arXiv

Halo-Galaxy Lensing: A Full Sky Approach

The halo-galaxy lensing correlation function or the average tangential shear profile over sampled halos is a very powerful means of measuring the halo masses, the mass profile, and the halo-mass correlation function of very large separations in the linear regime. We reformulate the halo-galaxy lensing correlation in harmonic space. We find that, counter-intuitively, errors in the conventionally used flat-sky approximation remain at a % level even at very small angles. The errors increase at larger angles and for lensing halos at lower redshifts: the effect is at a few % level at the baryonic acoustic oscillation scales for lensing halos of $z\sim 0.2$, and comparable with the effect of primordial non-Gaussianity with $f_{\rm NL}\sim 10$ at large separations. Our results allow to readily estimate/correct for the full-sky effect on a high-precision measurement of the average shear profile available from upcoming wide-area lensing surveys.

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

Re-capturing cosmic information

Gravitational lensing of distant galaxies can be exploited to infer the convergence field as a function of angular position on the sky. The statistics of this field, much like that of the cosmic microwave background (CMB), can be studied to extract information about fundamental parameters in cosmology, most notably the dark energy in the Universe. Unlike the CMB, the distribution of matter in the Universe which determines the convergence field is highly non-Gaussian, reflecting the nonlinear processes which accompanied structure formation. Much of the cosmic information contained in the initial field is therefore unavailable to the standard power spectrum measurements. Here we propose a method for re-capturing cosmic information by using the power spectrum of a simple function of the observed (nonlinear) convergence field. We adapt the approach of Neyrinck et al. (2009) to lensing by using a modified logarithmic transform of the convergence field. The Fourier transform of the log-transformed field has modes that are nearly uncorrelated, which allows for additional cosmological information to be extracted from small-scale modes.

preprint2010arXiv

Shear Power Spectrum Reconstruction using Pseudo-Spectrum Method

We develop a pseudo power spectrum technique for measuring the lensing power spectrum from weak lensing surveys in both the full sky and flat sky limits. The power spectrum approaches have a number of advantages over the traditional correlation function approach. We test the pseudo spectrum method by using numerical simulations with square-shape boundary that include masked regions with complex configuration due to bright stars and saturated spikes. Even when 25% of total area of the survey is masked, the method recovers the E-mode power spectrum at a sub-percent precision over a wide range of multipoles 100<l<10000. The systematic error is smaller than the statistical errors expected for a 2000 square degree survey. The residual B-mode spectrum is well suppressed in the amplitudes at less than a percent level relative to the E-mode. We also find that the correlated errors of binned power spectra caused by the survey geometry effects are not significant. Our method is applicable to the current and upcoming wide-field lensing surveys.

preprint2009arXiv

Constraints on Neutrino Masses from Weak Lensing

The weak lensing (WL) distortions of distant galaxy images are sensitive to neutrino masses by probing the suppression effect on clustering strengths of total matter in large-scale structure. We use the latest measurement of WL correlations, the CFHTLS data, to explore constraints on neutrino masses. We find that, while the WL data alone cannot place a stringent limit on neutrino masses due to parameter degeneracies, the constraint can be significantly improved when combined with other cosmological probes, the WMAP 5-year (WMAP5) data and the distance measurements of type-Ia supernovae (SNe) and baryon acoustic oscillations (BAO). The upper bounds on the sum of neutrino masses are m_tot = 1.1, 0.76 and 0.54 eV (95% CL) for WL+WMAP5, WMAP5+SNe+BAO, and WL+WMAP5+SNe+BAO, respectively, assuming a flat LCDM model with finite-mass neutrinos. In deriving these constraints, our analysis includes the non-Gaussian covariances of the WL correlation functions to properly take into account significant correlations between different angles.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

Tests of Gravity from Imaging and Spectroscopic Surveys

Tests of gravity on large-scales in the universe can be made using both imaging and spectroscopic surveys. The former allow for measurements of weak lensing, galaxy clustering and cross-correlations such as the ISW effect. The latter probe galaxy dynamics through redshift space distortions. We use a set of basic observables, namely lensing power spectra, galaxy-lensing and galaxy-velocity cross-spectra in multiple redshift bins (including their covariances), to estimate the ability of upcoming surveys to test gravity theories. We use a two-parameter description of gravity that allows for the Poisson equation and the ratio of metric potentials to depart from general relativity. We find that the combination of imaging and spectroscopic observables is essential in making robust tests of gravity theories. The range of scales and redshifts best probed by upcoming surveys is discussed. We also compare our parametrization to others used in the literature, in particular the gamma parameter modification of the growth factor.

preprint2009arXiv

The Impact of Non-Gaussian Errors on Weak Lensing Surveys

The weak lensing power spectrum carries cosmological information via its dependence on the growth of structure and on geometric factors. Since much of the cosmological information comes from scales affected by nonlinear clustering, measurements of the lensing power spectrum can be degraded by non-Gaussian covariances. Recently there have been conflicting studies about the level of this degradation. We use the halo model to estimate it and include new contributions related to the finite size of lensing surveys, following Rimes and Hamilton's study of 3D simulations. We find that non-Gaussian correlations between different multipoles can degrade the cumulative signal-to-noise for the power spectrum amplitude by up to a factor of 2 (or 5 for a worst-case model that exceeds current N-body simulation predictions). However, using an eight-parameter Fisher analysis we find that the marginalized errors on individual parameters are degraded by less than 10% (or 20% for the worst-case model). The smaller degradation in parameter accuracy is primarily because: individual parameters in a high-dimensional parameter space are degraded much less than the volume of the full Fisher ellipsoid; lensing involves projections along the line of sight, which reduce the non-Gaussian effect; some of the cosmological information comes from geometric factors which are not degraded at all. We contrast our findings with those of Lee & Pen (2008) who suggested a much larger degradation in information content. Finally, our results give a useful guide for exploring survey design by giving the cosmological information returns for varying survey area, depth and the level of some systematic errors.

preprint2007arXiv

Color Tomography

Lensing tomography with multi-color imaging surveys can probe dark energy and the cosmological power spectrum. However accurate photometric redshifts for tomography out to high redshift require imaging in five or more bands, which is expensive to carry out over thousands of square degrees. Since lensing makes coarse, statistical use of redshift information, we explore the prospects for tomography using limited color information from two or three band imaging. With an appropriate calibration sample, we find that it is feasible to create up to four redshift bins using imaging data in just the g, r and i bands. We construct such redshift sub-samples from mock catalogs by clustering galaxies in color space and discarding regions with poorly-defined redshift distributions. The loss of galaxy number density decreases the accuracy of lensing measurements, but even losing half or more of the galaxies is not a severe loss for large area surveys. We estimate the errors on lensing power spectra and dark energy parameters with color tomography and discuss trade-offs in survey area and filter choice. We discuss the systematic errors that may change our conclusions, especially the information needed to tackle intrinsic alignments.

preprint2005arXiv

Calibrating the Nonlinear Matter Power Spectrum: Requirements for Future Weak Lensing Surveys

Uncertainties in predicting the nonlinear clustering of matter are among the most serious theoretical systematics facing the upcoming wide-field weak gravitational lensing surveys. We estimate the accuracy with which the matter power spectrum will need to be calibrated in order not to contribute appreciably to the error budget for future weak lensing surveys. We consider the random statistical errors and the systematic biases in P(k), as well as some estimates based on current N-body simulations. While the power spectrum on relevant scales (0.1 < k/h Mpc^{-1} < 10) is currently calibrated with N-body simulations to about 5-10%, in the future it will have to be calibrated to about 1-2% accuracy, depending on the specifications of the survey. Encouragingly, we find that even the worst-case error that mimics the effect of cosmological parameters needs to be calibrated to no better than about 0.5-1%. These goals require a suite of high resolution N-body simulations on a relatively fine grid in cosmological parameter space, and should be achievable in the near future.

preprint2005arXiv

Cosmological parameters from lensing power spectrum and bispectrum tomography

We examine how lensing tomography with the bispectrum and power spectrum can constrain cosmological parameters and the equation of state of dark energy. Our analysis uses the full information at the two- and three-point level from angular scales of a few degrees to 5 arcminutes (50 < l < 3000), which will be probed by lensing surveys. We use all triangle configurations, cross-power spectra and bispectra constructed from up to three redshift bins with photometric redshifts, and relevant covariances in our analysis. We find that the parameter constraints from bispectrum tomography are comparable to those from power spectrum tomography. Combining the two improves parameter accuracies by a factor of three due to their complementarity. For the dark energy parameterization w(a) = w0 + wa(1-a), the marginalized errors from lensing alone are sigma(w0) = 0.03 fsky^{-1/2} and sigma(wa) = 0.1 fsky^{-1/2}. We show that these constraints can be further improved when combined with measurements of the cosmic microwave background or Type Ia supernovae. The amplitude and shape of the mass power spectrum are also shown to be precisely constrained. We use hyper-extended perturbation theory to compute the nonlinear lensing bispectrum for dark energy models. Accurate model predictions of the bispectrum in the moderately nonlinear regime, calibrated with numerical simulations, will be needed to realize the parameter accuracy we have estimated. Finally, we estimate how well the lensing bispectrum can constrain a model with primordial non-Gaussianity.

preprint2005arXiv

Systematic Errors in Future Weak Lensing Surveys: Requirements and Prospects for Self-Calibration

We study the impact of systematic errors on planned weak lensing surveys and compute the requirements on their contributions so that they are not a dominant source of the cosmological parameter error budget. The generic types of error we consider are multiplicative and additive errors in measurements of shear, as well as photometric redshift errors. In general, more powerful surveys have stronger systematic requirements. For example, for a SNAP-type survey the multiplicative error in shear needs to be smaller than 1%(fsky/0.025)^{-1/2} of the mean shear in any given redshift bin, while the centroids of photometric redshift bins need to be known to better than 0.003(fsky/0.025)^{-1/2}. With about a factor of two degradation in cosmological parameter errors, future surveys can enter a self-calibration regime, where the mean systematic biases are self-consistently determined from the survey and only higher-order moments of the systematics contribute. Interestingly, once the power spectrum measurements are combined with the bispectrum, the self-calibration regime in the variation of the equation of state of dark energy w_a is attained with only a 20-30% error degradation.

preprint2004arXiv

The Surprisingly Steep Mass Profile of Abell 1689, from a Lensing Analysis of Subaru Images

Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner profile, r<200 kpc/h, derived from deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NFW profile with virial mass, (1.93 \pm 0.20)10^15 M_sun, and surprisingly high concentration, c_vir=13.7^{+1.4}_{-1.1}, significantly larger than theoretically expected (c_vir~4), corresponding to a relatively steep overall profile. A slightly better fit is achieved with a steep power-law model that has its 2D logarithmic slope -3 and core radius theta_c~1.7' (r_c~210 kpc/h), whereas an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies selected to be redder than the cluster E/S0 sequence. By including the faint blue galaxy population a much smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal for r~400 kpc/h. This contamination is likely to affect most weak lensing results to date.

preprint2002arXiv

Probing Intracluster Magnetic Fields with Cosmic Microwave Background Polarization

Intracluster magnetic fields with $ \sim μG $ strength induce Faraday rotation on the cosmic microwave background (CMB) polarization. Measurements of this effect can potentially probe the detailed structure of intracluster magnetic fields across clusters, since the CMB polarization is a continuously varying field on the sky, in contrast to the conventional method restricted by the limited number of radio sources behind or inside a cluster. We here construct a method for extracting information on magnetic fields from measurements of the effect, combined with possible observations of the Sunyaev-Zel'dovich effect and X-ray emission for the same cluster which are needed to reconstruct the electron density fields. Employing the high-resolution magneto-hydrodynamic simulations performed by Dolag, Bartelmann & Lesch (1999) as a realistic model of magnetized intracluster gas distribution, we demonstrate how our reconstruction technique can well reproduce the magnetic fields, i.e., the spherically average radial profiles of the field strength and the coherence length.

preprint2002arXiv

The Kurtosis of the Cosmic Shear Field

We study the fourth-order moment of the cosmic shear field using the dark matter halo approach to describe the nonlinear gravitational evolution of structure in the universe. Since the third-order moment of the shear field vanishes because of symmetry, non-Gaussian signatures in its one-point statistics emerge at the fourth-order level. We argue that the shear kurtosis parameter S_4 = <gamma^4>/<gamma^2>^3 may be more directly applicable to realistic data than the well-studied higher-order statistics of the convergence field, since obtaining the convergence requires a non-local reconstruction from the measured shear field. We compare our halo model predictions for the variance, skewness and kurtosis of lensing fields with ray-tracing simulations of cold dark matter models and find good agreement. The shear kurtosis calculation is made tractable by developing approximations for fast and accurate evaluations of the 8-dimensional integrals needed to obtain the kurtosis. We show that on small scales it is dominated by correlations within halos more massive than 10^14 solar masses. The shear kurtosis is sensitive to the mass density parameter of the universe, Omega, and has relatively weak dependences on other parameters. The approximations we develop for the third- and fourth-order moments allow for accurate halo model predictions for the 3-dimensional mass distribution as well. We demonstrate their accuracy in the small scale regime, below 2 Mpc, where analytical approaches used in the literature so far cease to be accurate.

preprint2000arXiv

Gravitational Lensing Effect on the Two-point Correlation of Hotspots in the Cosmic Microwave Background

We investigate the weak gravitational lensing effect due to the large-scale structure of the universe on two-point correlations of local maxima ({\em hotspots}) in the 2D sky map of the cosmic microwave background (CMB) anisotropy. According to the Gaussian random statistics as most inflationary scenarios predict, the hotspots are discretely distributed with some {\em characteristic} angular separations on the last scattering surface owing to oscillations of the CMB angular power spectrum. The weak lensing then causes pairs of hotspots which are separated with the characteristic scale to be observed with various separations. We found that the lensing fairly smoothes the oscillatory features of the two-point correlation function of hotspots. This indicates that the hotspots correlations can be a new statistical tool for measuring shape and normalization of the power spectrum of matter fluctuations from the lensing signatures.

preprint1999arXiv

How is the local-scale gravitational instability influenced by the surrounding large-scale structure formation?

We develop the formalism to investigate the relation between the evolution of the large-scale (quasi) linear structure and that of the small-scale nonlinear structure in Newtonian cosmology within the Lagrangian framework. In doing so, we first derive the standard Friedmann expansion law using the averaging procedure over the present horizon scale. Then the large-scale (quasi) linear flow is defined by averaging the full trajectory field over a large-scale domain, but much smaller than the horizon scale. The rest of the full trajectory field is supposed to describe small-scale nonlinear dynamics. We obtain the evolution equations for the large-scale and small-scale parts of the trajectory field. These are coupled to each other in most general situations. It is shown that if the shear deformation of fluid elements is ignored in the averaged large-scale dynamics, the small-scale dynamics is described by Newtonian dynamics in an effective Friedmann-Robertson-Walker (FRW) background with a local scale factor. The local scale factor is defined by the sum of the global scale factor and the expansion deformation of the averaged large-scale displacement field. This means that the evolution of small-scale fluctuations is influenced by the surrounding large-scale structure through the modification of FRW scale factor. The effect might play an important role in the structure formation scenario. Furthermore, it is argued that the so-called {\it optimized} or {\it truncated} Lagrangian perturbation theory is a good approximation in investigating the large-scale structure formation up to the quasi nonlinear regime, even when the small-scale fluctuations are in the non-linear regime.

Masahiro Takada

What is connected

Connect this record

See the researcher in context

Building this map preview

86 published item(s)

The Milky Way Tomography with Subaru Hyper Suprime-Cam. I. Halo substructures

Analysis method for 3D power spectrum of projected tensor field with fast estimator and window convolution modelling: an application to intrinsic alignments

Constraining Dark matter annihilation with Dark Energy Survey Y3 LSBG sample

Cosmological constraints without fingers of God

Dark matter constraints with stacked gamma rays scales with the number of galaxies

Full-shape cosmology analysis of SDSS-III BOSS galaxy power spectrum using emulator-based halo model: a $5\%$ determination of $σ_8$

Mock catalogues of emission line galaxies based on the local mass density in dark-matter only simulations

The Subaru HSC weak lensing mass-observable scaling relations of spectroscopic galaxy groups from the GAMA survey

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

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

Precision analysis of the redshift-space galaxy bispectrum

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

Super-sample covariance of the thermal Sunyaev-Zel'dovich effect

Testing Stochastic Gravitational Wave Signals from Primordial Black Holes with Optical Telescopes

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Anisotropic separate universe simulations

Blinded challenge for precision cosmology with large-scale structure: results from effective field theory for the redshift-space galaxy power spectrum

Constraining dark matter annihilation with HSC Low Surface Brightness Galaxies

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}$

Hunting gravitational wave black holes with microlensing

On the wave optics effect on primordial black hole constraints from optical microlensing search

Proper motion measurements for stars up to $100$ kpc with Subaru HSC and SDSS Stripe 82

The impact of projection effects on cluster observables: stacked lensing and projected clustering

Towards a non-Gaussian model of redshift space distortions

A Search for a Contribution from Axion-Like Particles to the X-Ray Diffuse Background Utilizing the Earth's Magnetic Field

Mitigating the impact of fiber assignment on clustering measurements from deep galaxy redshift surveys

A New Milky Way Satellite Discovered In The Subaru/Hyper Suprime-Cam Survey

Detection of Stacked Filament Lensing Between SDSS Luminous Red Galaxies

Detection of the Splashback Radius and Halo Assembly bias of Massive Galaxy Clusters

Evidence of Halo Assembly Bias in Massive Clusters

Prime Focus Spectrograph (PFS) for the Subaru Telescope: Overview, recent progress, and future perspectives

Separate Universe Consistency Relation and Calibration of Halo Bias

Subaru high-z exploration of low-luminosity quasars (SHELLQs). I. Discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9

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

Properties of Weak Lensing Clusters Detected on Hyper Suprime-Cam 2.3 Square Degree Field

The Weak Lensing Signal and the Clustering of BOSS Galaxies I: Measurements

The Weak Lensing Signal and the Clustering of BOSS Galaxies II: Astrophysical and Cosmological Constraints

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

Joint analysis of cluster number counts and weak lensing power spectrum to correct for the super-sample covariance

Joint likelihood function of cluster counts and $n$-point correlation functions: Improving their power through including halo sample variance

Power Spectrum Super-Sample Covariance

Statistical challenges in weak lensing cosmology

Super-Sample Covariance in Simulations

Super-Sample Signal

Baryon acoustic oscillations with the cross-correlation of spectroscopic and photometric samples

Cosmological parameters from weak lensing power spectrum and bispectrum tomography: including the non-Gaussian errors

Extragalactic Science, Cosmology and Galactic Archaeology with the Subaru Prime Focus Spectrograph (PFS)

Information content of weak lensing power spectrum and bispectrum: including the non-Gaussian error covariance matrix

Modeling the Phase-Space Distribution around Massive Halos

Novel Probes of Gravity and Dark Energy

Perturbation theory for nonlinear halo power spectrum: the renormalized bias and halo bias

The Synergy between Weak Lensing and Galaxy Redshift Surveys

Where are the Luminous Red Galaxies (LRGs)? Using correlation measurements and lensing to relate LRGs to dark matter halos

Cross-correlating the Thermal Sunyaev-Zel'dovich Effect and the Distribution of Galaxy Clusters

LoCuSS: The Sunyaev-Zel'dovich Effect and Weak Lensing Mass Scaling Relation

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

Subaru weak-lensing measurement of a z = 0.81 cluster discovered by the Atacama Cosmology Telescope Survey

Testing Gravity with the Stacked Phase Space around Galaxy Clusters

The impact of massive neutrinos on the abundance of massive clusters

Using galaxy-galaxy weak lensing measurements to correct the Finger-of-God

Dark Energy from the log-transformed convergence field

Likelihood reconstruction method of real-space density and velocity power spectra from a redshift galaxy survey

LoCuSS: Subaru Weak Lensing Study of 30 Galaxy Clusters

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

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

A Clipping Method to Mitigate the Impact of Catastrophic Photometric Redshift Errors on Weak Lensing Tomography

A dipole anisotropy of galaxy distribution: Does the CMB rest-frame exist in the local universe?

Combining cluster observables and stacked weak lensing to probe dark energy: Self-calibration of systematic uncertainties

Direct measurement of dark matter halo ellipticity from two-dimensional lensing shear maps of 25 massive clusters

Halo-Galaxy Lensing: A Full Sky Approach

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

Re-capturing cosmic information

Shear Power Spectrum Reconstruction using Pseudo-Spectrum Method

Constraints on Neutrino Masses from Weak Lensing