Using $N$-body simulations for the first time, we show that the anisotropic primordial non-Gaussianity (PNG) causes a scale-dependent modification, given by $1/k^2$ at small $k$ limit, in the three-dimensional power spectra of halo shapes (intrinsic alignments), whilst the conventional power spectrum of halo number density field remains unaffected. We discuss that wide-area imaging and spectroscopic surveys observing the same region of the sky allow us to constrain the quadrupole PNG coefficient $f_{\rm NL}^{s=2}$ at a precision comparable with that of the cosmic microwave background.

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.

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 10$^9$ stars in the bulge region over 10 years with the Rubin Observatory Legacy Survey of Space and Time (LSST) would enable one to find more than $3\times 10^4$ 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.

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.

An accurate theoretical template of the redshift-space galaxy power spectrum, if applicable out to nonlinear scales, enables us to extract more stringent and robust constraints on cosmological parameters from the measured galaxy clustering. In this work we develop a simulation-based template, so-called emulator, for the redshift-space power spectrum of dark matter halos. Using the redshift-space halo power spectra measured from the Dark Quest $N$-body simulation suite that covers 101 flat-geometry $w$CDM cosmologies around the Planck $\Lambda$CDM model, we feed these data into a feed-forward neural network to build the fast and accurate emulation of the power spectrum from the linear to nonlinear scales up to $k \simeq 0.6 \, h \, {\rm Mpc}^{-1}$. Our emulator achieves about 1 and 5% fractional accuracies in predicting the monopole and quadrupole moments of the power spectrum, respectively, for halos of $\sim 10^{13}h^{-1}M_\odot$ that correspond to host halos of the SDSS LOWZ- and CMASS-like galaxies, where the achieved accuracies are sufficient compared to the statistical errors of SDSS volume. The validation and performance of the emulator are given by the comparison of the emulator predictions with the power spectra directly measured from the simulations for validation sets that are not used in the training. We demonstrate that the emulator outputs can be used to make model predictions for the redshift-space power spectrum of galaxies by employing user-fed models for the halo-galaxy connection, such as the halo occupation distribution. The emulator allows us to easily incorporate the Finger-of-God effect due to the virial motions of galaxies and the Alcock--Paczy\'{n}ski distortions. Our code can compute the redshift-space galaxy power spectrum in a CPU subseconds, and is ready to perform the emulator-based cosmological analysis for the exiting and upcoming galaxy redshift surveys.

We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $\Lambda$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the $E$ mode, $P^{(\ell)}_{EE}(k)$, based on the $E$/$B$-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to nonlinear scales, and the large-scale amplitudes at $k\lesssim 0.1~h~{\rm Mpc}^{-1}$ are related to the matter power spectrum via a constant coefficient ($A_{\rm IA}$), similar to the linear bias parameter. We find that the cross- and auto-power spectra $P_{\delta E}$ and $P_{EE}$ at nonlinear scales, $k\gtrsim 0.1~h~{\rm Mpc}^{-1}$, show different $k$-dependences relative to the matter power spectrum, suggesting a violation of the nonlinear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts and cosmological parameters ($\Omega_{\rm m}, S_8$). The cumulative signal-to-noise ratio for the IA power spectra is about 60\% of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be powerful tools to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.

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.

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 ($\lambda)$ 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 \lambda$, 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\AA) 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.

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.

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.

Primordial black holes (PBHs) are a viable candidate for dark matter if the PBH masses are in the currently unconstrained "sublunar" mass range. We revisit the possibility that PBHs were produced by nucleation of false vacuum bubbles during inflation. We show that this scenario can produce a population of PBHs that simultaneously accounts for all dark matter, explains the candidate event in Subaru Hyper Suprime-Cam (HSC) data, and contains both heavy black holes as observed by LIGO and very heavy seeds of super-massive black holes. We demonstrate with numerical studies that future observations of HSC, as well as other optical surveys, such as LSST, will be able to provide a definitive test for this generic PBH formation mechanism if it is the dominant source of dark matter.

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.

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 ($\sigma_8$ and $\Omega_{\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 ($\sigma_8,\Omega_{\rm m},D_{\rm A}$'s and $H$'s) compared to $k=0.2\,h\,$Mpc$^{-1}$.

We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $\Delta z =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and 2.4 arcmin$^{-2}$ from lower to higher redshifts, respectively. We adopt the standard TPCF estimators, $\xi_\pm$, for our cosmological analysis, given that we find no evidence of the significant B-mode shear. The TPCFs are detected at high significance for all ten combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7'<\theta<56'$ for $\xi_+$ and $28'<\theta<178'$ for $\xi_-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}=0.804_{-0.029}^{+0.032}$, and $\Omega_m=0.346_{-0.100}^{+0.052}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints.

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.

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.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

Constraining the relation between the richness $N$ and the halo mass $M$ over a wide redshift range for optically-selected clusters is a key ingredient for cluster-related science in optical surveys, including the Subaru Hyper Suprime-Cam (HSC) survey. We measure stacked weak lensing profiles around 1747 HSC CAMIRA clusters over a redshift range of $0.1\leq z_{\rm cl}\leq 1.0$ with $N\geq 15$ using the HSC first-year shear catalog covering $\sim$$140$ ${\rm deg^2}$. The exquisite depth and image quality of the HSC survey allow us to measure lensing signals around the high-redshift clusters at $0.7\leq z_{\rm cl}\leq 1.0$ with a signal-to-noise ratio of 19 in the comoving radius range $0.5\lesssim R\lesssim 15 h^{-1}{\rm Mpc}$. We constrain richness-mass relations $P(\ln N|M,z)$ of the HSC CAMIRA clusters assuming a log-normal distribution without informative priors on model parameters, by jointly fitting to the lensing profiles and abundance measurements under both Planck and WMAP cosmological models. We show that our model gives acceptable $p$-values when we add redshift dependent terms which are proportional to $\ln (1+z)$ and $[\ln (1+z)]^{2}$ into the mean and scatter relations of $P(\ln N|M,z)$. Such terms presumably originate from the variation of photometric redshift errors as a function of the redshift. We show that the constraints on the mean relation $\langle M|N \rangle$ are consistent between the Planck and WMAP models, whereas the scatter values $\sigma_{\ln M|N}$ for the Planck model are systematically larger than those for the WMAP model. We also show that the scatter values for the Planck model increase toward lower richness values, whereas those for the WMAP model are consistent with constant values as a function of richness. This result highlights the importance of the scatter in the mass-richness relation for cluster cosmology.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360<lambda<1330 nm with options for fiber-fed spectrographs at high resolution and a panoramic IFU at a separate focus. This would enable transformational progress via its ability to access a larger fraction of objects from Gaia, LSST, Euclid, and WFIRST than any currently funded or planned spectroscopic facility. An ESO-sponsored study (arXiv:1701.01976) discussed the scientific potential in ambitious new spectroscopic surveys in Galactic astronomy, extragalactic astronomy, and cosmology. The US community should establish links with European and other international communities to plan for such a powerful facility and maximize the potential of large aperture multi-object spectroscopy given the considerable investment in deep imaging surveys.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

We constrain the abundance of primordial black holes (PBH) using 2622 microlensing events obtained from 5-years observations of stars in the Galactic bulge by the Optical Gravitational Lensing Experiment (OGLE). The majority of microlensing events display a single or at least continuous population that has a peak around the light curve timescale $t_{\rm E}\simeq 20~{\rm days}$ and a wide distribution over the range $t_{\rm E}\simeq [1, 300]~{\rm days}$, while the data also indicates a second population of 6 ultrashort-timescale events in $t_{\rm E}\simeq [0.1,0.3]~{\rm days}$, which are advocated to be due to free-floating planets. We confirm that the main population of OGLE events can be well modeled by microlensing due to brown dwarfs, main sequence stars and stellar remnants (white dwarfs and neutron stars) in the standard Galactic bulge and disk models for their spatial and velocity distributions. Using the dark matter (DM) model for the Milky Way (MW) halo relative to the Galactic bulge/disk models, we obtain the tightest upper bound on the PBH abundance in the mass range $M_{\rm PBH}\simeq[10^{-6},10^{-3}]M_\odot$ (Earth-Jupiter mass range), if we employ null hypothesis that the OGLE data does not contain any PBH microlensing event. More interestingly, we also show that Earth-mass PBHs can well reproduce the 6 ultrashort-timescale events, without the need of free-floating planets, if the mass fraction of PBH to DM is at a per cent level, which is consistent with other constraints such as the microlensing search for Andromeda galaxy (M31) and the longer timescale OGLE events. Our result gives a hint of PBH existence, and can be confirmed or falsified by microlensing search for stars in M31, because M31 is towards the MW halo direction and should therefore contain a much less number of free-floating planets, even if exist, than the direction to the MW center.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

We use the full-sky ray-tracing weak lensing simulations to generate 2268 mock catalogues for the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalogue. Our mock catalogues take into account various effects as in the real data: the survey footprints, inhomogeneous angular distribution of source galaxies, statistical uncertainties in photometric redshift (photo-$z$) estimate, variations in the lensing weight, and the statistical noise in galaxy shape measurements including both intrinsic shapes and the measurement errors. We then utilize our mock catalogues to evaluate statistical uncertainties expected in measurements of cosmic shear two-point correlations $\xi_{\pm}$ with tomographic redshift information for the HSC survey. We develop a quasi-analytical formula for the Gaussian sample variance properly taking into account the number of source pairs in the survey footprints. The standard Gaussian formula significantly overestimates or underestimates the mock results by $50\%$ level. We also show that different photo-$z$ catalogues or the six disconnected fields, rather than a consecutive geometry, cause variations in the covariance by $\sim 5\%$. The mock catalogues enable us to study the chi-square distribution for $\xi_{\pm}$. We find the wider distribution than that naively expected for the distribution with the degrees-of-freedom of data vector used. Finally, we propose a method to include non-zero multiplicative bias in mock shape catalogue and show the non-zero multiplicative bias can change the effective shape noise in cosmic shear analyses. Our results suggest an importance of estimating an accurate form of the likelihood function (and therefore the covariance) for robust cosmological parameter inference from the precise measurements.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

The galaxy bispectrum, measuring excess clustering of galaxy triplets, offers a probe of dark energy via baryon acoustic oscillations (BAOs). However up to now it has been severely underused due to the combinatorically explosive number of triangles. Here we exploit interference in the bispectrum to identify triangles that amplify BAOs. This approach reduces the computational cost of estimating covariance matrices, offers an improvement in BAO constraints equivalent to lengthening BOSS by 30%, and simplifies adding bispectrum BAO information to future large-scale redshift survey analyses.

We present a picture of the matter bispectrum in a novel "interferometric" basis designed to highlight interference of the baryon acoustic oscillations (BAO) in the power spectra composing it. Triangles where constructive interference amplifies BAO provide stronger cosmic distance constraints than triangles with destructive interference. We show that the amplitude of the BAO feature in the full cyclically summed bispectrum can be decomposed into simpler contributions from single terms or pairs of terms in the perturbation theory bispectrum, and that across large swathes of our parameter space the full BAO amplitude is described well by the amplitude of BAO in a single term. The dominant term is determined largely by the $F^{(2)}$ kernel of Eulerian standard perturbation theory. We present a simple physical picture of the BAO amplitude in each term; the BAO signal is strongest in triangle configurations where two wavenumbers differ by a multiple of the BAO fundamental wavelength.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

In this study, we measure the response of matter and halo projected power spectra $P^{\rm 2D}_{\rm XY}(k)$ (X, Y are matter and/or halos), to a large-scale density contrast, $\delta_{\rm b}$, using separate universe simulations. We show that the fractional response functions, i.e., $\mathrm{d}\ln P^{\rm 2D}_{\rm XY}(k)/\mathrm{d}\delta_{\rm b}$, are identical to their respective three-dimensional power spectra within simulation measurement errors. Then, using various $N$-body simulation combinations (small-box simulations with periodic boundary conditions and sub-volumes of large-box simulations) to construct {mock observations of projected fields}, we study how super-survey modes, in both parallel and perpendicular directions to the projection direction, affect the covariance matrix of $P^{\rm 2D}_{\rm XY}(k)$, known as super-sample covariance (SSC). Our results indicate that the SSC term provides dominant contributions to the covariances of matter-matter and matter-halo spectra at small scales but does not provide significant contributions in the halo-halo spectrum. We observe that the large-scale density contrast in each redshift shell causes most of the SSC effect, and we did not observe a SSC signature arising from large-scale tidal field within the levels of measurement accuracy. We also develop a response approach to calibrate the SSC term for cosmic shear correlation function and galaxy--galaxy weak lensing, and validate the method by comparison with the light-cone ray-tracing simulations. Our method provides a reasonably accurate, albeit computationally inexpensive, way to calibrate the covariance matrix for clustering observables available from wide-area galaxy surveys.

Primordial black holes (PBHs) have long been suggested as a viable candidate for the elusive dark matter (DM). The abundance of such PBHs has been constrained using a number of astrophysical observations, except for a hitherto unexplored mass window of $M_{\rm PBH}=[10^{-14},10^{-9}]M_\odot$. Here we carry out a dense-cadence (2~min sampling rate), 7 hour-long observation of the Andromeda galaxy (M31) with the Subaru Hyper Suprime-Cam to search for microlensing of stars in M31 by PBHs lying in the halo regions of the Milky Way (MW) and M31. Given our simultaneous monitoring of tens of millions of stars in M31, if such light PBHs make up a significant fraction of DM, we expect to find many microlensing events for the PBH DM scenario. However, we identify only a single candidate event, which translates into the most stringent upper bounds on the abundance of PBHs in the mass range $M_{\rm PBH}\simeq [10^{-11}, 10^{-6}]M_\odot$.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

Stacked lensing is a powerful means of measuring the average mass distribution around large-scale structure tracers. There are two stacked lensing estimators used in the literature, denoted as $\Delta\Sigma$ and $\gamma_+$, which are related as $\Delta\Sigma=\Sigma_{\rm cr}\gamma_+$, where $\Sigma_{\rm cr}(z_l,z_s)$ is the critical surface mass density for each lens-source pair ($z_l$ and $z_s$ are lens and source redshifts, respectively). In this paper we derive a formula for the covariance matrix of $\Delta\Sigma$-estimator focusing on `weight' function to improve the signal-to-noise ($S/N$). We assume that the lensing fields and the distribution of lensing objects obey the Gaussian statistics. With this formula, we show that, if background galaxy shapes are weighted by an amount of $\Sigma_{\rm cr}^{-2}(z_l,z_s)$, the $\Delta\Sigma$-estimator maximizes the $S/N$ in the shot noise limited regime. We also show that the $\Delta\Sigma$-estimator with the weight $\Sigma_{\rm cr}^{-2}$ gives a greater $(S/N)^2$ than that of the $\gamma_+$-estimator by about 5--25\% for lensing objects at redshifts comparable with or higher than the median of source galaxy redshifts for hypothetical Subaru HSC and DES surveys. However, for low-redshift lenses such as $z_l<0.3$, the $\gamma_+$-estimator has higher $(S/N)^2$ than $\Delta\Sigma$. We also discuss that the $(S/N)^2$ for $\Delta\Sigma$ at large separations in the sample variance limited regime can be boosted, by up to a factor of 1.5, if one adopts a weight of $\Sigma_{\rm cr}^{-\alpha}$ with $\alpha>2$. Our formula allows one to explore how the combination of the different estimators can approach an optimal estimator in all regimes of redshifts and separation scales.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

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

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

The spatial curvature ($K$ or $\Omega_K$) is one of the most fundamental parameters of an isotropic and homogeneous universe and has a close link to the physics of the early Universe. Combining the radial and angular diameter distances measured via the baryon acoustic oscillation (BAO) experiments allows us to unambiguously constrain the curvature. The method is primarily based on the metric theory, but is less sensitive to the theory of structure formation other than the existence of the BAO scale and is free of any model of dark energy. In this paper, we estimate a best achievable accuracy of constraining the curvature with the BAO experiments. We show that an all-sky, cosmic-variance-limited galaxy survey covering the Universe up to $z> 4$ enables a precise determination of the curvature to an accuracy of $\sigma(\Omega_K)\simeq 10^{-3}$. When we assume a model of dark energy - either the cosmological constant or the $(w_0,w_a)$ model - it can achieve a precision of $\sigma(\Omega_K)\simeq \mbox{a few}\times 10^{-4}$. These forecasts require a high sampling density of galaxies, and are degraded by up to a factor of a few for a survey with a finite number density of $\sim 10^{-3}(h/{\rm Mpc})^3$.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We constrain the scaling relation between optical richness ($\lambda$) and halo mass ($M$) for a sample of SDSS redMaPPer galaxy clusters within the context of the {\it Planck} cosmological model. We use a forward modeling approach where we model the probability distribution of optical richness for a given mass, $P(\ln \lambda| M)$. To model the abundance and the stacked lensing profiles, we use an emulator specifically built to interpolate the halo mass function and the stacked lensing profile for an arbitrary set of halo mass and redshift, which is calibrated based on a suite of high-resolution $N$-body simulations. We apply our method to 8,312 SDSS redMaPPer clusters with $20\le \lambda \le 100$ and $0.10\le z_{\lambda}\le0.33$, and show that the log-normal distribution model for $P(\lambda|M)$, with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, $\left\langle \ln{\lambda}\right\rangle(M)=A+B\ln(M/M_{\rm pivot})$, with $A=3.207^{+0.044}_{-0.046}$ and $B=0.993^{+0.041}_{-0.055}$ (68\%~CL), where the pivot mass scale $M_{\rm pivot}=3\times 10^{14} h^{-1}M_\odot$, and the scatter $\sigma_{\mathrm{\ln\lambda}|M}=\sigma_0+q\ln(M/M_{\rm pivot})$ with $\sigma_0=0.456^{+0.047}_{-0.039}$ and $q=-0.169^{+0.035}_{-0.026}$. We find that a large scatter in halo masses is required at the lowest richness bins ($20\le \lambda \lesssim 30$) in order to reproduce the measurements. Without such a large scatter, the model prediction for the lensing profiles tends to overestimate the measured amplitudes. This might imply a possible contamination of intrinsically low-richness clusters due to the projection effects. Such a low-mass halo contribution is significantly reduced when applying our method to the sample of $30\le \lambda \le 100$.

We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.

Although the large-scale perturbations beyond a finite-volume survey region are not direct observables, these affect measurements of clustering statistics of small-scale (sub-survey) perturbations in the large-scale structure, compared with the ensemble average, via the mode-coupling effect. In this paper we show that the large-scale tides induced by scalar perturbations cause apparent anisotropic distortions in the redshift-space power spectrum of galaxies in a way depending on an alignment between the tides, the wavevector of small-scale modes and the line-of-sight direction. Using the perturbation theory of structure formation, we derive the $\textit{response}$ function of the redshift-space power spectrum to the large-scale tides. We then investigate the impact of the large-scale tides on estimation of cosmological distances and the redshift-space distortion parameter via the measured redshift-space power spectrum for a hypothetical large-volume survey, based on the Fisher matrix formalism. To do this, we treat the large-scale tides as a signal, rather than an additional source of the statistical errors, and show that a degradation in the parameter is restored if we can employ the prior on the rms amplitude expected for the standard cold dark matter (CDM) model. We also discuss whether the large-scale tides can be constrained at an accuracy better than the CDM prediction, if the effects up to a larger wavenumber in the nonlinear regime can be included.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We present wide-field (167 deg$^2$) weak lensing mass maps from the Hyper Supreme-Cam Subaru Strategic Program (HSC-SSP). We compare these weak lensing based dark matter maps with maps of the distribution of the stellar mass associated with luminous red galaxies. We find a strong correlation between these two maps with a correlation coefficient of $\rho=0.54\pm0.03$ (for a smoothing size of $8'$). This correlation is detected even with a smaller smoothing scale of $2'$ ($\rho=0.34\pm 0.01$). This detection is made uniquely possible because of the high source density of the HSC-SSP weak lensing survey ($\bar{n}\sim 25$ arcmin$^{-2}$). We also present a variety of tests to demonstrate that our maps are not significantly affected by systematic effects. By using the photometric redshift information associated with source galaxies, we reconstruct a three-dimensional mass map. This three-dimensional mass map is also found to correlate with the three-dimensional galaxy mass map. Cross-correlation tests presented in this paper demonstrate that the HSC-SSP weak lensing mass maps are ready for further science analyses.

We develop a method to simulate galaxy-galaxy weak lensing by utilizing all-sky, light-cone simulations and their inherent halo catalogs. Using the mock catalog to study the error covariance matrix of galaxy-galaxy weak lensing, we compare the full covariance with the "jackknife" (JK) covariance, the method often used in the literature that estimates the covariance from the resamples of the data itself. We show that there exists the variation of JK covariance over realizations of mock lensing measurements, while the average JK covariance over mocks can give a reasonably accurate estimation of the true covariance up to separations comparable with the size of JK subregion. The scatter in JK covariances is found to be $\sim10\%$ after we subtract the lensing measurement around random points. However, the JK method tends to underestimate the covariance at the larger separations, more increasingly for a survey with a higher number density of source galaxies. We apply our method to the the Sloan Digital Sky Survey (SDSS) data, and show that the 48 mock SDSS catalogs nicely reproduce the signals and the JK covariance measured from the real data. We then argue that the use of the accurate covariance, compared to the JK covariance, allows us to use the lensing signals at large scales beyond a size of the JK subregion, which contains cleaner cosmological information in the linear regime.

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.1<z<1.1$ and richness $\hat{N}_{\rm mem}>15$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.

We utilize the HSC CAMIRA cluster catalog and the photo-$z$ galaxy catalog constructed in the HSC wide field (S16A), covering $\sim$ 174 deg$^{2}$, to study the star formation activity of galaxies in different environments over 0.2 $<$ $z$ $<$ 1.1. We probe galaxies down to $i \sim$ 26, corresponding to a stellar mass limit of log$_{10}$(M$_*$/M$_{\odot}$) $\sim$ 8.2 and $\sim$ 8.6 for star-forming and quiescent populations, respectively, at $z$ $\sim$ 0.2. The existence of the red sequence for low stellar mass galaxies in clusters suggests that the environmental quenching persists to halt the star formation in the low-mass regime. In addition, star-forming galaxies in groups or clusters are systematically biased toward lower values of specific star formation rate by 0.1 -- 0.3 dex with respect to those in the field and the offsets shows no strong redshift evolution over our redshift range, implying a universal slow quenching mechanism acting in the dense environments since $z$ $\sim$ 1.1. Moreover, the environmental quenching dominates the mass quenching in low mass galaxies, and the quenching dominance reverses in high mass ones. The transition mass is greater in clusters than in groups, indicating that the environmental quenching is more effective for massive galaxies in clusters compared to groups.

Long-wavelength matter inhomogeneities contain cleaner information on the nature of primordial perturbations as well as the physics of the early universe. The large-scale coherent overdensity and tidal force, not directly observable for a finite-volume galaxy survey, are both related to the Hessian of large-scale gravitational potential and therefore of equal importance. We show that the coherent tidal force causes a homogeneous anisotropic distortion of the observed distribution of galaxies in all three directions, perpendicular and parallel to the line-of-sight direction. This effect mimics the redshift-space distortion signal of galaxy peculiar velocities, as well as a distortion by the Alcock-Paczynski effect. We quantify its impact on the redshift-space power spectrum to the leading order, and discuss its importance for the ongoing and upcoming galaxy surveys.

The peculiar velocity field measured by redshift-space distortions (RSD) in galaxy surveys provides a unique probe of the growth of large-scale structure. However, systematic effects arise when including satellite galaxies in the clustering analysis. Since satellite galaxies tend to reside in massive halos with a greater halo bias, the inclusion boosts the clustering power. In addition, virial motions of the satellite galaxies cause a significant suppression of the clustering power due to nonlinear RSD effects. We develop a novel method to recover the redshift-space power spectrum of halos from the observed galaxy distribution by minimizing the contamination of satellite galaxies. The cylinder grouping method (CGM) we study effectively excludes satellite galaxies from a galaxy sample. However, we find that this technique produces apparent anisotropies in the reconstructed halo distribution over all the scales which mimic RSD. On small scales, the apparent anisotropic clustering is caused by exclusion of halos within the anisotropic cylinder used by the CGM. On large scales, the misidentification of different halos in the large-scale structures, aligned along the line-of-sight, into the same CGM group, causes the apparent anisotropic clustering via their cross-correlation with the CGM halos. We construct an empirical model for the CGM halo power spectrum, which includes correction terms derived using the CGM window function at small scales as well as the linear matter power spectrum multiplied by a simple anisotropic function at large scales. We apply this model to a mock galaxy catalog at z=0.5, designed to resemble SDSS-III BOSS CMASS galaxies, and find that our model can predict both the monopole and quadrupole power spectra of the host halos up to k<0.5 h/Mpc to within 5%.

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

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.

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.

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 $\Lambda$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.

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 $\Lambda$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$\sigma$ 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.

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$\sigma$ in the weak lensing and 4.4$\sigma$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

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-$\sigma$ 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 $\Lambda$-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.

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$\sigma$ 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 $\sigma_{\ln M_{500c}}\sim 0.2$ - $0.3$ could reconcile a difference between the measurements and the simulations.

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.

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 ($\Omega_{\rm m}$ and $\sigma_8$) to these observables simultaneously, and thus break the degeneracy between galaxy bias and cosmology. Adopting a flat $\Lambda$CDM cosmology with priors on $\Omega_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 $\Omega_{\rm m}=0.310^{+0.019}_{-0.020}$ and $\sigma_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$.

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

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 $\delta_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 $\delta_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 $\Lambda$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 $\delta_b$ is both unbiased and has the predicted variance. Parameter degeneracies arise since the response of the matter power spectrum to $\delta_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.

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 $\Lambda$-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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We develop a novel abundance matching method to construct a mock catalog of luminous red galaxies (LRGs) in SDSS, using catalogs of halos and subhalos in N-body simulations for a LCDM model. Motivated by observations suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age >5Gyr, we assume that simulated halos at z=2 (z2-halo) are progenitors for LRG-host subhalos observed today, and we label the most tightly bound particles in each progenitor z2-halo as LRG ``stars''. We then identify the subhalos containing these stars to z=0.3 (SDSS redshift) in descending order of the masses of z2-halos until the comoving number density of the matched subhalos becomes comparable to the measured number density of SDSS LRGs, n=10^{-4} (h/Mpc)^3. Once the above prescription is determined, our only free parameter is the number density of halos identified at z=2 and this parameter is fixed to match the observed number density at z = 0.3. By tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly compute, from the mock catalog, the distributions of central and satellite LRGs and their internal motions in each host halo at z=0.3. While the SDSS LRGs are galaxies selected by the magnitude and color cuts from the SDSS images and are not necessarily a stellar-mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo occupation distribution for central and satellite LRGs, the projected auto-correlation function of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, in the angle-averaged redshift-space power spectrum.

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.

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.

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 $\Delta 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, $\sigma_z/(1+z)$, are better than 10-20% level.

We present a stacked weak-lensing analysis of an approximately mass-selected sample of 50 galaxy clusters at 0.15<z<0.3, based on observations with Suprime-Cam on the Subaru Telescope. We develop a new method for selecting lensed background galaxies from which we estimate that our sample of red background galaxies suffers just 1% contamination. We detect the stacked tangential shear signal from the full sample of 50 clusters, based on this red sample of background galaxies, at a total signal-to-noise ratio of S/N=32.7. The Navarro-Frenk-White model is an excellent fit to the data, yielding sub-10% statistical precision on mass and concentration: Mvir=7.19^{+0.53}_{-0.50}\times10^{14}h^{-1}Msol, cvir=5.41^{+0.49}_{-0.45} (c_{200}=4.22^{+0.40}_{-0.36}). Tests of a range of possible systematic errors, including shear calibration and stacking-related issues, indicate that they are sub-dominant to the statistical errors. The concentration parameter obtained from stacking our approximately mass-selected cluster sample is broadly in line with theoretical predictions. Moreover, the uncertainty on our measurement is comparable with the differences between the different predictions in the literature. Overall our results highlight the potential for stacked weak-lensing methods to probe the mean mass density profile of cluster-scale dark matter halos with upcoming surveys, including Hyper-Suprime-Cam, Dark Energy Survey, and KIDS.

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{\rho}_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 \Lambda-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.

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.

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.

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.

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.

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.

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.

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.

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.

Using 1000 ray-tracing simulations for a {\Lambda}-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.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

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

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.

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.

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)

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.

We present a joint weak-lensing/X-ray study of galaxy cluster mass-observable scaling relations, motivated by the critical importance of accurate calibration of mass proxies for future X-ray missions, including eROSITA. We use a sample of 12 clusters at z\simeq0.2 that we have observed with Subaru and XMM-Newton to construct relationships between the weak-lensing mass (M), and three X-ray observables: gas temperature (T), gas mass (Mgas), and quasi-integrated gas pressure (Yx) at overdensities of \Delta=2500, 1000, and 500 with respect to the critical density. We find that Mgas at \Delta\le1000 appears to be the most promising mass proxy of the three, because it has the lowest intrinsic scatter in mass at fixed observable: \sigma_lnM\simeq0.1, independent of cluster dynamical state. The scatter in mass at fixed T and Yx is a factor of \sim2-3 larger than at fixed Mgas, which are indicative of the structural segregation that we find in the M-T and M-Yx relationships. Undisturbed clusters are found to be \sim40% and \sim20% more massive than disturbed clusters at fixed T and Yx respectively at \sim2\sigma significance. In particular, A1914 - a well-known merging cluster - significantly increases the scatter and lowers the the normalization of the relation for disturbed clusters. We also investigated the covariance between intrinsic scatter in M-Mgas and M-T relations, finding that they are positively correlated. This contradicts the adaptive mesh refinement simulations that motivated the idea that Yx may be a low scatter mass proxy, and agrees with more recent smoothed particle hydrodynamic simulations based on the Millennium Simulation. We also propose a method to identify a robust mass proxy based on principal component analysis. The statistical precision of our results are limited by the small sample size and the presence of the extreme merging cluster in our sample.

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.

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.

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.

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.

We present an early broad-brush analysis of Herschel/PACS observations of star-forming galaxies in 8 galaxy clusters drawn from our survey of 30 clusters at z~0.2. We define a complete sample of 192 spectroscopically confirmed cluster members down to L_TIR>3x10^10L_sun and L_K>0.25L_sun. The average K-band and bolometric infrared luminosities of these galaxies both fade by a factor of ~2 from clustercentric radii of 2r_200 to ~0.5r_200, indicating that as galaxies enter the clusters ongoing star-formation stops first in the most massive galaxies, and that the specific star-formation rate (SSFR) is conserved. On smaller scales the average SSFR jumps by 25%, suggesting that in cluster cores processes including ram pressure stripping may trigger a final episode of star-formation that presumably exhausts the remaining gas. This picture is consistent with our comparison of the Herschel-detected cluster members with the cluster mass distributions, as measured in our previous weak-lensing study of these clusters. For example, the spatial distribution of the Herschel sources is positively correlated with the structures in the weak-lensing mass maps at 5sigma significance, with the strongest signal seen at intermediate group-like densities. The strong dependence of the total cluster IR luminosity on cluster mass L_TIR propto M_virial^2 is also consistent with accretion of galaxies and groups of galaxies (i.e. the substructure mass function) driving the cluster IR luminosity. The most surprising result is that roughly half of the Herschel-detected cluster members have redder S_100/S_24 flux ratios than expected, based on the Rieke et al. models. On average cluster members are redder than non-members, and the fraction of red galaxies increases towards the cluster centers, both of which indicate that these colors are not attributable to systematic photometric errors. [Abridged]

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\sigma significance level. The mean ellipticity is found to be e = 0.46 \pm 0.04 (1\sigma), 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).

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.

Future or ongoing galaxy redshift surveys can put stringent constraints on neutrinos masses via the high-precision measurements of galaxy power spectrum, when combined with cosmic microwave background (CMB) information. In this paper we develop a method to model galaxy power spectrum in the weakly nonlinear regime for a mixed dark matter (CDM plus finite-mass neutrinos) model, based on perturbation theory (PT) whose validity is well tested by simulations for a CDM model. In doing this we carefully study various aspects of the nonlinear clustering and then arrive at a useful approximation allowing for a quick computation of the nonlinear power spectrum as in the CDM case. The nonlinear galaxy bias is also included in a self-consistent manner within the PT framework. Thus the use of our PT model can give a more robust understanding of the measured galaxy power spectrum as well as allow for higher sensitivity to neutrino masses due to the gain of Fourier modes beyond the linear regime. Based on the Fisher matrix formalism, we find that BOSS or Stage-III type survey, when combined with Planck CMB information, gives a precision of total neutrino mass constraint, sigma(m_nu,tot) 0.1eV, while Stage-IV type survey may achieve sigma(m_nu,tot) 0.05eV, i.e. more than a 1-sigma detection of neutrino masses. We also discuss possible systematic errors on dark energy parameters caused by the neutrino mass uncertainty. The significant correlation between neutrino mass and dark energy parameters is found, if the information on power spectrum amplitude is included. More importantly, for Stage-IV type survey, a best-fit dark energy model may be biased and falsely away from the underlying true model by more than the 1-sigma statistical errors, if neutrino mass is ignored in the model fitting.

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.

We use 5000 cosmological N-body simulations of 1(Gpc/h)^3 box for the concordance LCDM model in order to study the sampling variances of nonlinear matter power spectrum. We show that the non-Gaussian errors can be important even on large length scales relevant for baryon acoustic oscillations (BAO). Our findings are (1) the non-Gaussian errors degrade the cumulative signal-to-noise ratios (S/N) for the power spectrum amplitude by up to a factor of 2 and 4 for redshifts z=1 and 0, respectively. (2) There is little information on the power spectrum amplitudes in the quasi-nonlinear regime, confirming the previous results. (3) The distribution of power spectrum estimators at BAO scales, among the realizations, is well approximated by a Gaussian distribution with variance that is given by the diagonal covariance component. (4) For the redshift-space power spectrum, the degradation in S/N by non-Gaussian errors is mitigated due to nonlinear redshift distortions. (5) For an actual galaxy survey, the additional shot noise contamination compromises the cosmological information inherent in the galaxy power spectrum, but also mitigates the impact of non-Gaussian errors. The S/N is degraded by up to 30% for a WFMOS-type survey. (6) The finite survey volume causes additional non-Gaussian errors via the correlations of long-wavelength fluctuations with the fluctuations we want to measure, further degrading the S/N values by about 30% even at high redshift z=3.

We present the first galaxy evolution results from the Local Cluster Substructure Survey (LoCuSS), a multi-wavelength survey of 100 X-ray selected galaxy clusters at 0.15<z<0.30. LoCuSS combines far-UV through far-IR observations of cluster galaxies with gravitational lensing analysis and X-ray data to investigate the interplay between the hierarchical assembly of clusters and the evolution of cluster galaxies. Here we present new panoramic Spitzer/MIPS 24micron observations of the merging cluster Abell 1758 at z=0.279 spanning 6.5x6.5Mpc and reaching a 90% completeness limit of 400uJy. We estimate a global cluster SFR of 910\pm320 M_sun/yr within 3 Mpc of the cluster centre, originating from 42 galaxies with L_IR > 5x10^10 L_sun. The obscured activity in A1758 is therefore comparable with that in Cl 0024+1654, the most active cluster previously studied at 24um. The obscured galaxies faithfully trace the cluster potential as revealed by the weak-lensing mass map of the cluster, including numerous mass peaks at R~2-3Mpc that are likely associated with infalling galaxy groups and filamentary structures. However the core (R<500kpc) of A1758N is 2x more active in the IR than that of A1758S, likely reflecting differences in the recent dynamical history of the two clusters. The 24micron results from A1758 therefore suggest that dust-obscured cluster galaxies are common in merging clusters and suggests that obscured activity in clusters is triggered by both the details of cluster-cluster mergers and processes that operate at larger radii including those within in-falling groups. Our ongoing far-UV through far-IR observations of a large sample of clusters should allow us to disentangle the different physical processes responsible for triggering obscured star formation in clusters.

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.

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.

We present the first attempt to analytically study the nonlinear matter power spectrum for a mixed dark matter (cold dark matter plus neutrinos of total mass ~0.1eV) model based on cosmological perturbation theory. The suppression in the power spectrum amplitudes due to massive neutrinos is, compared to the linear regime, enhanced in the weakly nonlinear regime where standard linear theory ceases to be accurate. We demonstrate that, thanks to this enhanced effect and the gain in the range of wavenumbers to which the PT prediction is applicable, the use of such a nonlinear model may enable a precision of sigma(m_nu,tot) ~ 0.07eV in constraining the total neutrino mass for the planned galaxy redshift survey, a factor of 2 improvement compared to the linear regime.

(Abridged) Combining cosmic shear power spectra and cluster counts is powerful to improve cosmological parameter constraints and/or test inherent systematics. However they probe the same cosmic mass density field, if the two are drawn from the same survey region, and therefore the combination may be less powerful than first thought. We investigate the cross-covariance between the cosmic shear power spectra and the cluster counts based on the halo model approach, where the cross-covariance arises from the three-point correlations of the underlying mass density field. Fully taking into account the cross-covariance as well as non-Gaussian errors on the lensing power spectrum covariance, we find a significant cross-correlation between the lensing power spectrum signals at multipoles l~10^3 and the cluster counts containing halos with masses M>10^{14}Msun. Including the cross-covariance for the combined measurement degrades and in some cases improves the total signal-to-noise ratios up to plus or minus 20% relative to when the two are independent. For cosmological parameter determination, the cross-covariance has a smaller effect as a result of working in a multi-dimensional parameter space, implying that the two observables can be considered independent to a good approximation. We also discuss that cluster count experiments using lensing-selected mass peaks could be more complementary to cosmic shear tomography than mass-selected cluster counts of the corresponding mass threshold. Using lensing selected clusters with a realistic usable detection threshold (S/N~6 for a ground-based survey), the uncertainty on each dark energy parameter may be roughly halved by the combined experiments, relative to using the power spectra alone.

We present results from a weak lensing analysis of the galaxy cluster A1689 (z=0.183) based on deep wide-field imaging data taken with Suprime-Cam on Subaru telescope. A maximum entropy method has been used to reconstruct directly the projected mass distribution of A1689 from combined lensing distortion and magnification measurements of red background galaxies.The resulting mass distribution is clearly concentrated around the cD galaxy, and mass and light in the cluster are similarly distributed in terms of shape and orientation. The azimuthally-averaged mass profile from the two-dimensional reconstruction is in good agreement with the earlier results from the Subaru one-dimensional analysis of the weak lensing data, supporting the assumption of quasi-circular symmetry in the projected mass distribution of the cluster.

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.

(abridged) A wide-field galaxy redshift survey allows one to probe galaxy clustering at largest spatial scales, which carries an invaluable information on horizon-scale physics complementarily to the cosmic microwave background (CMB). Assuming the planned survey consisting of z~1 and z~3 surveys with areas of 2000 and 300 square degrees, respectively, we study the prospects for probing dark energy clustering from the measured galaxy power spectrum, assuming the dynamical properties of dark energy are specified in terms of the equation of state and the effective sound speed c_e in the context of an adiabatic cold dark matter dominated model. The dark energy clustering adds a power to the galaxy power spectrum amplitude at spatial scales greater than the sound horizon, and the enhancement is sensitive to redshift evolution of the net dark energy density, i.e. the equation of state. We find that the galaxy survey, when combined with Planck, can distinguish dark energy clustering from a smooth dark energy model such as the quintessence model (c_e=1), when c_e<0.04 (0.02) in the case of the constant equation of state w_0=-0.9 (-0.95). An ultimate full-sky survey of z~1 galaxies allows the detection when c_e<0.08 (0.04) for w_0=0.9 (-0.95). We also investigate a degeneracy between the dark energy clustering and the non-relativistic neutrinos implied from the neutrino oscillation experiments, because the two effects both induce a scale-dependent modification in the galaxy power spectrum shape at largest spatial scales accessible from the galaxy survey. It is shown that a wider redshift coverage can efficiently separate the two effects by utilizing the different redshift dependences, where dark energy clustering is apparent only at low redshifts z<1.

(abridged) High-z galaxy redshift surveys open up exciting possibilities for precision determinations of neutrino masses and inflationary models. The high-z surveys are more useful for cosmology than low-z ones owing to much weaker non-linearities in matter clustering, redshift-space distortion and galaxy bias. We can then utilize the two-dimensional information of the linear power spectrum in angular and redshift space to measure the scale-dependent suppression of matter clustering due to neutrino free-streaming as well as the shape of the primordial power spectrum. To illustrate capabilities of high-z surveys for constraining neutrino masses and the primordial power spectrum, we compare three future redshift surveys covering 300 square degrees at 0.5<z<2, 2<z<4, and 3.5<z<6.5. We find that, combined with the cosmic microwave background data expected from the Planck satellite, these surveys allow precision determination of the total neutrino mass with the projected errors of sigma(m_nu)=0.059, 0.043, and 0.025 eV, respectively, thus yielding a positive detection of the neutrino mass rather than an upper limit, as sigma(m_nu) is smaller than the lower limits to the neutrino masses implied from the neutrino oscillation experiments. The accuracies of constraining the tilt and running index of the primordial power spectrum, sigma(n_s)=(3.8, 3.7, 3.0)x10^-3, and sigma(alpha_s)=(5.9, 5.7, 2.4)x10^-3, respectively, are smaller than the current uncertainties by more than an order of magnitude, which will allow us to discriminate between candidate inflationary models. In particular, the error on alpha_s from the highest redshift survey is not very far away from the prediction of a class of simple inflationary models driven by a massive scalar field with self-coupling, alpha_s=-(0.8-1.2)x10^-3.

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.

The steep mass profile of A1689 derived from recent detailed lensing observations is not readily reconciled with the low concentration halos predicted by the standard CDM model. However, halo triaxiality may act to bias the profile constraints derived assuming a spherically symmetric mass distribution, since lensing relates only to the projected mass distribution. A degree of halo triaxiality is inherent to the CDM structure formation, arising from the collision-less nature of the dark matter. Here we compare the CDM-based model predictions of triaxial halo with the precise lensing measurements of A1689 based on the ACS/HST and Subaru data, over a wide range of 10kpc<r<2Mpc. The model lensing profiles cover the intrinsic spread of halo mass and shape, and are projected over all inclinations when comparing with the data. We show that the model parameters are only weakly constrained and strongly degenerate mainly because of the lack of information along the line of sight. In particular, the limits on the concentration parameter become less restrictive with increasing triaxiality. Yet, by comparing the obtained constraints with expected probability distributions for the axis ratio and concentration parameters computed from numerical simulations, we find that ~6% of cluster-size halos in the CDM model can match the A1689 lensing observations at the 2-sigma level, corresponding to cases where the major-axis of the halo is aligned with the line of sight. Thus halo triaxiality could reduce the apparent discrepancy between theory and observation. This hypothesis needs to be further explored by a statistical lensing study for other clusters as well as by complementary three-dimensional information derived using X-ray, kinematics, and SZ effect observations.

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.

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.

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.

We study the effects of halo substructure and a distribution in the concentration parameter of haloes on large-scale structure statistics. The effects on the power spectrum and bispectrum are studied on the smallest scales accessible from future surveys. We compare halo-model predictions with results based on N-body simulations, but also extend our predictions to 10-kpc scales which will be probed by future simulations. We find that weak-lensing surveys proposed for the coming decade can probe the power spectrum on small enough scales to detect substructure in massive haloes. We discuss the prospects of constraining the mass fraction in substructure in view of partial degeneracies with parameters such as the tilt and running of the primordial power spectrum.

Baryonic oscillations in the galaxy power spectrum have been studied as a way of probing dark-energy models. While most studies have focused on spectroscopic surveys at high redshift, large multi-color imaging surveys have already been planned for the near future. In view of this, we study the prospects for measuring baryonic oscillations from angular statistics of galaxies binned using photometric redshifts. We use the galaxy bispectrum in addition to the power spectrum; this allows us to measure and marginalize over possibly complex galaxy bias mechanisms to get robust cosmological constraints. In our parameter estimation we allow for a weakly nonlinear biasing scheme that may evolve with redshift by two bias parameters in each of ten redshift bins. We find that a multi-color imaging survey that probes redshifts beyond one can give interesting constraints on dark-energy parameters. In addition, the shape of the primordial power spectrum can be measured to better accuracy than with the CMB alone. We explore the impact of survey depth, area, and calibration errors in the photometric redshifts on dark-energy constraints.

By obtaining photometric redshift information, tomography allows us to cross-correlate galaxy ellipticities in different source redshift bins. The cross-correlation is non-vanishing because the different bins share much of the foreground mass distribution from which, over Gpc scales, the lensing signal is built. If the redshift bins are thick enough however, the cross-correlations are insensitive to contamination from the intrinsic alignments of galaxies since these fall off rapidly on scales larger than a few tens of Mpc. We forecast how lensing tomography using only the cross-power spectra can constrain cosmological parameters compared to tomography including the auto-spectra. It is shown that the parameter errors are degraded by only O(10%) for 5 or more source redshift bins. Thus, the cross-power spectrum tomography can be a simple, model-independent means of reducing the intrinsic alignment contamination while retaining most of the constraints on cosmology.

We explore the ability of weak lensing surveys to locate massive clusters. We use both analytic models of dark matter halos and mock weak lensing surveys generated from a large cosmological N-body simulation. The analytic models describe average properties of weak lensing halos and predict the number counts, enabling us to compute an effective survey selection function. We test the model prediction for the peak number counts in weak lensing mass maps against the mock numerical data, and find that the noise due to intrinsic galaxy ellipticities causes a systematic effect which increases the peak counts. We develop a correction scheme for the systematic effect in an empirical manner, and show that, after the correction, the model prediction agrees well with the mock data. The mock data is also used to examine the completeness and efficiency of the weak lensing halo search with fully taking into account the noise and the projection effect by large-scale structures. We show that the detection threshold of S/N=4-5 gives an optimal balance between completeness and efficiency. Our results suggest that, for a weak lensing survey with a galaxy number density of ng=30/arcmin^2 with a mean redshift z=1, the mean number of peaks in the 10sq deg area is N_peak=62 for a detection threshold S/N=4. The contamination rate is 42%, and thus, on average, 36 out of 62 peaks (at least) are signals from real halos. Weak lensing surveys thus provide a reasonably efficient way to searching for massive clusters.

The lensing magnification effect due to large-scale structure is statistically measurable by correlation of size fluctuations in distant galaxy images as well as by the QSO-galaxy cross-correlation. We use the halo model formulation of Takada & Jain (2003) to compute these magnification-induced correlations without employing the weak lensing approximation, mu=1+2kappa. Our predictions thus include the full contribution from non-linear magnification that is due to lensing halos. We compare the model prediction with ray-tracing simulations and find excellent agreement over a range of angular scales (0.'5<theta<30'). The non-linear magnification contribution enhances the amplitude of the magnification correlation relative to the weak lensing approximation, and the non-linear correction is more significant on smaller angular scales and for sources at higher redshifts. The enhancement amounts to 10-25% on arcminite scales for the QSO-galaxy cross-correlation, even after inclusion of a realistic model of galaxy clustering within the host halo. Therefore, it is necessary to account for the non-linear contribution in theoretical models in order to make an unbiased, cosmological interpretation of the precise measurements expected from forthcoming massive surveys.

We use the halo model of clustering to compute two- and three-point correlation functions for weak lensing, and apply them in a new statistical technique to measure properties of massive halos. We present analytical results on the eight shear three-point correlation functions constructed using combination of the two shear components at each vertex of a triangle. We compare the amplitude and configuration dependence of the functions with ray-tracing simulations and find excellent agreement for different scales and models. These results are promising, since shear statistics are easier to measure than the convergence. In addition, the symmetry properties of the shear three-point functions provide a new and precise way of disentangling the lensing E-mode from the B-mode due to possible systematic errors. We develop an approach based on correlation functions to measure the properties of galaxy-group and cluster halos from lensing surveys. Shear correlations on small scales arise from the lensing matter within halos of mass M > 10^13 solar masses. Thus the measurement of two- and three-point correlations can be used to extract information on halo density profiles, primarily the inner slope and halo concentration. We demonstrate the feasibility of such an analysis for forthcoming surveys. We include covariances in the correlation functions due to sample variance and intrinsic ellipticity noise to show that 10% accuracy on profile parameters is achievable with surveys like the CFHT Legacy survey, and significantly better with future surveys. Our statistical approach is complementary to the standard approach of identifying individual objects in survey data and measuring their properties.

With the advent of high-quality surveys in cosmology the full three-point correlation function will be a valuable statistic for describing structure formation models. It contains information on cosmological parameters and detailed halo properties that cannot be extracted from the two-point correlation function. We use the halo clustering model to analytically calculate the three-point correlation function (3PCF) for general cosmological fields. We present detailed results for the configuration dependence of the 3-dimensional mass and galaxy distributions and the 2-dimensional cosmic shear field. We work in real space, where higher order correlation functions on small scales are easier to measure and interpret, but halo model calculations get rapidly intractable. Hence we develop techniques for accurate evaluations of the 1, 2 and 3-halo contributions to the 3PCF. The 3PCF violates the hierarchical ansatz in both its scale and configuration dependence. We study the behavior of the coefficient Q in the hierarchical expansion from large, quasilinear scales down to about 20 kpc. We find that the nonlinear 3PCF is sensitive to the halo profile of massive halos, especially its inner slope. We model the distribution of galaxies in halos and show that the 3PCF of red galaxies has a weaker configuration and scale dependence than the mass, while for blue galaxies it is very sensitive to the parameters of the galaxy formation model. The 3PCF from weak lensing on the other hand shows different scalings due to projection effects and a sensitivity to cosmological parameters.

The three-point correlation function (3PCF) of the spin-2 fields, the cosmic shear and cosmic microwave background polarization fields, is a statistical measure of non-Gaussian signals. At each vertex of a triangle, the shear field has two independent components. The resulting eight possible 3PCFs were recently investigated by Schneider & Lombardi (2002) and Zaldarriaga & Scoccimarro (2002). Using rotation and parity transformations they addressed the question: how many non-zero components for the shear 3PCF exist? We address this question using an analytical model and measurements from ray-tracing simulations. We show that all the eight 3PCFs are generally non-zero and have comparable amplitude. These eight 3PCFs can be used to improve the signal-to-noise from survey data and their configuration dependence can be used to separate the contribution from E and B-modes. We estimate the signal-to-noise for measuring the shear 3PCF from weak lensing surveys using simulated maps that include the noise due to intrinsic ellipticities. A deep lensing survey with area exceeding 10 square degrees would allow for the detection of the shear 3PCFs, but a survey area of order 100 square degrees is needed for accurate measurements.

Intracluster magnetic fields with $ \sim \mu 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.

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.

The statistical property of the weak lensing fields is studied quantitatively using the ray-tracing simulations. Motivated by the empirical lognormal model that characterizes the probability distribution function(PDF) of the three-dimensional mass distribution excellently, we critically investigate the validity of lognormal model in the weak lensing statistics. Assuming that the convergence field, $\kappa$, is approximately described by the lognormal distribution, we present analytic formulae of convergence for the one-point PDF and the Minkowski functionals. Comparing those predictions with ray-tracing simulations in various cold dark matter models, we find that the one-point lognormal PDF can describe the non-Gaussian tails of convergence fields accurately up to $\nu\sim10$, where $\nu$ is the level threshold given by $\nu\equiv\kappa/\var^{1/2}$, although the systematic deviation from lognormal prediction becomes manifest at higher source redshift and larger smoothing scales. The lognormal formulae for Minkowski functionals also fit to the simulation results when the source redshift is low. Accuracy of the lognormal-fit remains good even at the small angular scales, where the perturbation formulae by Edgeworth expansion break down. On the other hand, lognormal models does not provide an accurate prediction for the statistics sensitive to the rare events such as the skewness and the kurtosis of convergence. We therefore conclude that the empirical lognormal model of the convergence field is safely applicable as a useful cosmological tool, as long as we are concerned with the non-Gaussianity of $\nu\simlt5$ for low source redshift samples.

The observed magnetic field of microgauss strength in clusters of galaxies should induce the Faraday rotation effect on the linearly polarized cosmic microwave background (CMB) radiation when the CMB radiation propagates through a cluster at low redshift. Employing the Press-Schechter prescription for the cluster abundance under the cold dark matter (CDM) scenario and a plausible isothermal $\beta$-model for the gas distribution, we compute angular power spectra of the CMB polarization fields including the Faraday rotation mixing effect under the simple assumption of uniform magnetic field configuration across a cluster. As a result, we find that a parity-odd B-type polarization pattern is statistically generated on the observed sky, even when the primary polarization only contains the parity-even E-mode component, such as in the case of pure scalar perturbations. The generated B-type polarization has a peak with amplitude of $\sim0.1 \mu{\rm K}(B_0/0.1 \mu{\rm G})(\nu_0/10 {\rm GHz})^{-2}$ at angular scales of $l\approx 1000$ for the currently favored adiabatic $\Lambda$CDM model. This result also implies that, if the magnetic field has a $0.5 \mu{\rm G}$ strength and we observe at lower frequencies such as $\nu_0\simlt 5 {\rm GHz}$, the secondary signal due to the Faraday rotation effect could be comparable to the magnitude of the primary polarization. The frequency dependence of the Faraday rotation can be then used to discriminate the effect from primary and other secondary signals on the CMB polarizations. Our results therefore offer a new empirical opportunity to measure or constrain the intracluster magnetic field in the average sense.

We investigate a cross-correlation between the weak gravitational lensing field of the large-scale structure, $\kappa$, and the secondary temperature fluctuation field, $\Delta$, of cosmic microwave background (CMB) induced by Thomson scattering of CMB photons off the ionized medium in the mildly nonlinear structure. The cross-correlation is expected to observationally unveil the biasing relation between the dark matter and ionized medium distributions in the large-scale structure. We develop a formalism for calculating the cross-correlation function and its angular power spectrum based on the small angle approximation. As a result, we find that leading contribution to the cross-correlation comes from the secondary CMB fluctuation field induced by the quadratic Doppler effects of the bulk velocity field $v$ of ionized medium. The magnitude of the cross-correlation can be estimated as $\kaco{\kappa\Delta}\sim 10^{-10}$ on small angular scales ($l\simgt 1000$) under the currently favored cold CDM model of the structure formation and the simplest scenario of homogeneous reionization after a given redshift of $z_{\rm ion}\simgt 5$. Interestingly, we find that the density-modulated quadratic Doppler effect of $O(\delta_e v^2)$ produces the cross-correlation with the weak lensing field that has {\em linear} dependence on the electron density fluctuation field $\delta_e$ or equivalently on the biasing relation between the dark matter and electron distributions. In other words, the angular power spectrum could be both positive or negative, depending on the positive biasing or antibiasing, respectively. Detection of the cross-correlation thus offers a new opportunity to observationally understand the reionization history of intergalactic medium connected to the dark matter clustering.

We propose a new method for extracting the non-Gaussian signatures on the isotemperature statistics in the cosmic microwave background (CMB) sky, which is induced by the gravitational lensing due to the intervening large-scale structure of the universe. To develop the method, we focus on a specific statistical property of the intrinsic Gaussian CMB field; a field point in the map that has a larger absolute value of the temperature threshold tends to have a larger absolute value of the curvature parameter defined by a trace of second derivative matrix of the temperature field, while the ellipticity parameter similarly defined is uniformly distributed independently of the threshold because of the isotropic nature of the Gaussian field. The weak lensing then causes a stronger distortion effect on the isotemperature contours with higher threshold and especially induces a coherent distribution of the ellipticity parameter correlated with the threshold as a result of the coupling between the CMB curvature parameter and the gravitational tidal shear in the observed map. These characteristic patterns can be statistically picked up by considering three independent characteristic functions, which are obtained from the averages of quadratic combinations of the second derivative fields of CMB over isotemperature contours with each threshold. Consequently, we find that the lensing effect generates non-Gaussian signatures on those functions that have a distinct functional dependence of the threshold. We test the method using numerical simulations of CMB maps and show that the lensing signals can be measured definitely, provided that we use CMB data with sufficiently low noise and high angular resolution.

We apply the morphological descriptions of two-dimensional contour map, the so-called Minkowski functionals (the area fraction, circumference, and Euler characteristics), to the convergence field $\kappa(\bm{\theta})$ of the large-scale structure reconstructed from the shear map produced by the ray-tracing simulations. The perturbation theory of structure formation has suggested that the non-Gaussian features on the Minkowski functionals with respect to the threshold in the weakly nonlinear regime are induced by the three skewness parameters of $\kappa$ that are sensitive to the density parameter of matter, $\Omega_{\rm m}$. We show that, in the absence of noise due to the intrinsic ellipticities of source galaxies with which the perturbation theory results can be recovered, the accuracy of $\Omega_{\rm m}$ determination is improved by $\sim 20%$ using the Minkowski functionals compared to the conventional method of using the direct measure of skewness.

We investigate the non-Gaussian signatures in Cosmic Microwave Background (CMB) maps induced by the intervening large-scale structure through weak lensing. In order to measure the deviation from the Gaussian behavior of the intrinsic temperature anisotropies, we use a family of three morphological descriptors, the so-called Minkowski functionals. We show analytically how these quantities depend on the temperature threshold, and compare the results to numerical experiments including the instrumental effects of Planck. Minkowski functionals can directly measure the statistical properties of the displacement field and hence provide useful constraints on large-scale structure formation in the past.

We investigate the dependence of cosmological parameters on the number count of peaks (local maxima and minima) in the cosmic microwave background (CMB) sky. The peak statistics contains the whole information of acoustic oscillations in the angular power spectrum $C_l$ over $l$-space and thus it can place complementary constraints on the cosmological parameters to those obtained from measurements of $C_l$. Based on the instrumental specifications of Planck, we find that the number count of peaks can provide new constraints on the combination of the matter density $\Omega_{\rm m}$ and the Hubble parameter $h$ approximately scaled as $\Omega_{\rm m} h^{-4.9}$ for a flat $\Lambda$CDM model with $\Omega_{m}=0.3$ and $h=0.7$. Therefore, we suggest that combining it with the constraints from $C_l$ scaled as $\Omega_{m}h^{3.8}$ (or commonly $\Omega_{\rm m}h^2$) can potentially determine $\Omega_{\rm m}$ or equivalently solve the cosmic degeneracy by the CMB data alone.

We present quantitative investigations of the weak lensing effect on the two-point correlation functions of local maxima (hotspots), $\xipk(\theta)$, in the cosmic microwave background (CMB) maps. The lensing effect depends on the projected mass fluctuations between today and the redshift $z_{\rm rec}\approx1100$. If adopting the Gaussian assumption for the primordial temperature fluctuations field, the peak statistics can provide an additional information about the intrinsic distribution of hotspots that those pairs have some characteristic separation angles. The weak lensing then redistributes hotspots in the observed CMB maps from the intrinsic distribution and consequently imprints non-Gaussian signatures onto $\xipk(\theta)$. Especially, since the intrinsic $\xipk(\theta)$ has a pronounced depression feature around the angular scale of $\theta\approx 70'$ for a flat universe, the weak lensing induces a large smoothing at the scale. We show that the lensing signature therefore has an advantage to effectively probe mass fluctuations with large wavelength modes around $\lambda\approx 50 h^{-1}{\rm Mpc}$. To reveal the detectability, we performed numerical experiments with specifications of {\em MAP} and {\em Planck Surveyor} including the instrumental effects of beam smoothing and detector noise. It is then found that our method can successfully provide constraints on amplitude of the mass fluctuations and cosmological parameters in a flat universe with and without cosmological constant, provided that we use maps with 65% sky coverage expected from Planck.

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.

We formulate the Lagrangian perturbation theory to solve the non-linear dynamics of self-gravitating fluid within the framework of the post-Newtonian approximation in general relativity, using the (3+1) formalism. Our formulation coincides with Newtonian Lagrangian perturbation theory developed by Buchert for the scale much smaller than the horizon scale, and with the gauge invariant linearized theory in the longitudinal gauge conditions for the linear regime. These are achieved by using the gauge invariant quantities at the initial time when the linearized theory is valid enough. The post-Newtonian corrections in the solution of the trajectory field of fluid elements are calculated in the explicit forms. Thus our formulation allows us to investigate the evolution of the large-scale fluctuations involving relativistic corrections from the early regime such as the decoupling time of matter and radiation until today. As a result, we are able to show that naive Newtonian cosmology to the structure formation will be a good approximation even for the perturbations with scales not only inside but also beyond the present horizon scale in the longitudinal coordinates. Although the post-Newtonian corrections are small, it is shown that they have a growing transverse mode which is not present in Newtonian theory as well as in the gauge invariant linearized theory. Such post-Newtonian order effects might produce characteristic appearances of the large-scale structure formation, for example, through the observation of anisotropies of the cosmic microwave background radiation (CMB). Furthermore since our approach has a straightforward Newtonian limit, it will be also convenient for numerical implementation based on the presently available Newtonian simulation. ......

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.

We point out the existence of a new type of growing transverse mode in the gravitational instability. This appears as a post-Newtonian effect to Newtonian dynamics. We demonstrate this existence by formulating the Lagrangian perturbation theory in the framework of the cosmological post-Newtonian approximation in general relativity. Such post-Newtonian order effects might produce characteristic appearances of large-scale structure formation, for example, through the observation of anisotropy of the cosmic microwave background radiation (CMB).