We present a weak lensing cluster search using Hyper Suprime-Cam Subaru Strategic Program (HSC survey) first-year data. We pay special attention to the dilution effect of cluster member and foreground galaxies on weak lensing signals from clusters of galaxies; we adopt the globally normalized weak lensing estimator which is least affected by cluster member galaxies, and we select source galaxies by using photometric redshift information to mitigate the effect of foreground galaxies. We produce six samples of source galaxies with different low-$z$ galaxy cuts, construct weak lensing mass maps for each of the source sample, and search for high peaks in the mass maps that cover an effective survey area of $\sim$120 deg$^2$. We combine six catalogs of high peaks into a sample of cluster candidates which contains 124 high peaks with signal-to-noise ratios greater than five. We cross-match the peak sample with the public optical cluster catalog constructed from the same HSC survey data to identify cluster counterparts of the peaks. We find that 107 out of 124 peaks have matched clusters within 5 arcmin from peak positions. Among them, we define a sub-sample of 64 secure clusters that we use to examine dilution effects on our weak lensing cluster search. We find that source samples with the low-$z$ galaxy cuts mitigate the dilution effect on weak lensing signals of high-$z$ clusters ($z > 0.3$), and thus combining multiple peak catalogs from different source samples improves the efficiency of weak lensing cluster searches.

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 present the first results of a pilot X-ray study of 37 rich galaxy clusters at $0.1<z<1.1$ in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) field. Diffuse X-ray emissions from these clusters were serendipitously detected in the XMM-Newton fields of view. We systematically analyze X-ray images of 37 clusters and emission spectra of a subsample of 17 clusters with high photon statistics by using the XMM-Newton archive data. The frequency distribution of the offset between the X-ray centroid or peak and the position of the brightest cluster galaxy was derived for the optical cluster sample. The fraction of relaxed clusters estimated from the X-ray peak offsets in 17 clusters is $29\pm11(\pm13)$\%, which is smaller than that of the X-ray cluster samples such as HIFLUGCS. Since the optical cluster search is immune to the physical state of X-ray-emitting gas, it is likely to cover a larger range of the cluster morphology. We also derived the luminosity-temperature relation and found that the slope is marginally shallower than those of X-ray-selected samples and consistent with the self-similar model prediction of 2. Accordingly, our results show that the X-ray properties of the optical clusters are marginally different from those observed in the X-ray samples. The implication of the results and future prospects are briefly discussed.

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.

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.

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 present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong.

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.

We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological $N$-body simulations. The data sets include full-sky convergence and shear maps from redshifts $z=0.05$ to $5.3$ at intervals of $150 \, h^{-1}{\rm Mpc}$ comoving radial distance (corresponding to a redshift interval of $\Delta z \simeq 0.05$ at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to $z=5.3$. The dark matter halos are identified from the same $N$-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and Luminous Red Galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy-galaxy and cluster-galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy-galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within $5\%$ up to $\ell = 3000$ (or at an angular scale $\theta > 0.5$ arcsin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

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 the results of clustering analyses of Lyman break galaxies (LBGs) at $z\sim3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. Mean halo masses, calculated as $\langle M_{h} \rangle=10^{11.7}-10^{12.8}h^{-1}M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy shows higher values at $z=3-5$ than $z=0.5-1.5$ over the entire stellar mass range. Satellite fractions of dropout galaxies, even at less massive haloes, are found to drop sharply from $z=2$ down to less than $0.04$ at $z=3-5$. These results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star}<10^{10}M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We derive observationally, for the first time, $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass at a peak in the star-formation efficiency, at $3<z<5$, and it shows a little increasing trend with cosmic time at $z>3$. In addition, $M_{{\rm h}}^{{\rm pivot}}$ and its normalization are found to be almost unchanged during $0<z<5$. Our study shows an observational evidence that galaxy formation is ubiquitously most efficient near a halo mass of $M_{{\rm h}}\sim10^{12}M_{\odot}$ over cosmic time.

We present the results of an halo occupation distribution (HOD) analysis of star-forming galaxies at $z \sim 2$. We obtained high-quality angular correlation functions based on a large sgzK sample, which enabled us to carry out the HOD analysis. The mean halo mass and the HOD mass parameters are found to increase monotonically with increasing $K$-band magnitude, suggesting that more luminous galaxies reside in more massive dark haloes. The luminosity dependence of the HOD mass parameters was found to be the same as in the local Universe; however, the masses were larger than in the local Universe over all ranges of magnitude. This implies that galaxies at $z \sim 2$ tend to form in more massive dark haloes than in the local Universe, a process known as downsizing. By analysing the dark halo mass evolution using the extended Press--Schechter formalism and the number evolution of satellite galaxies in a dark halo, we find that faint Lyman break galaxies at $z \sim 4$ could evolve into the faintest sgzKs $(22.0 < K \leq 23.0)$ at $z \sim 2$ and into the Milky-Way-like galaxies or elliptical galaxies in the local Universe, whereas the most luminous sgzKs $(18.0 \leq K \leq 21.0)$ could evolve into the most massive systems in the local Universe. The stellar-to-halo mass ratio (SHMR) of the sgzKs was found to be consistent with the prediction of the model, except that the SHMR of the faintest sgzKs was smaller than the prediction at $z \sim 2$. This discrepancy may be attributed that our samples are confined to star-forming galaxies.

Ongoing and future wide-field galaxy surveys can be used to locate a number of clusters of galaxies with cosmic shear measurement alone. We study constraints on cosmological models using statistics of weak lensing selected galaxy clusters. We extend our previous theoretical framework to model the statistical properties of clusters in variants of cosmological models as well as in the standard $\Lambda$CDM model. Weak lensing selection of clusters does not rely on the conventional assumption such as the relation between luminosity and mass and/or hydrostatic equilibrium, but a number of observational effects compromise robust identification. We use a large set of realistic mock weak-lensing catalogs as well as analytic models to perform a Fisher analysis and make forecast for constraining two competing cosmological models, $w$CDM model and $f(R)$ model proposed by Hu & Sawicki, with our lensing statistics. We show that weak lensing selected clusters are excellent probe of cosmology when combined with cosmic shear power spectrum even in presence of galaxy shape noise and masked regions. With the information of weak lensing selected clusters, the precision of cosmological parameter estimate can be improved by a factor of $\sim1.6$ and $\sim8$ for $w$CDM model and $f(R)$ model, respectively. Hyper Suprime-Cam survey with sky coverage of $1250$ squared degrees can constrain the equation of state of dark energy $w_{0}$ with a level of $\Delta w_0 \sim0.1$. It can also constrain the additional scalar degree of freedom in $f(R)$ model with a level of $|f_{R0}| \sim5\times10^{-6}$, when constraints from cosmic microwave background measurements are incorporated. Future weak lensing surveys with sky coverage of $20,000$ squared degrees will place tighter constraints on $w_{0}$ and $|f_{R0}|$ even without cosmic microwave background measurements.

We explore a variety of statistics of clusters selected with cosmic shear measurement by utilizing both analytic models and large numerical simulations. We first develop a halo model to predict the abundance and the clustering of weak lensing selected clusters. Observational effects such as galaxy shape noise are included in our model. We then generate realistic mock weak lensing catalogs to test the accuracy of our analytic model. To this end, we perform full-sky ray-tracing simulations that allow us to have multiple realizations of a large continuous area. We model the masked regions on the sky using the actual positions of bright stars, and generate 200 mock weak lensing catalogs with sky coverage of ~1000 squared degrees. We show that our theoretical model agrees well with the ensemble average of statistics and their covariances calculated directly from the mock catalogues. With a typical selection threshold, ignoring shape noise correction causes overestimation of the clustering of weak lensing selected clusters with a level of about 10%, and shape noise correction boosts the cluster abundance by a factor of a few. We calculate the cross-covariances using the halo model with accounting for the effective reduction of the survey area due to masks. The covariance of the cosmic shear auto power spectrum is affected by the mode-coupling effect that originates from sky masking. Our model and the results can be readily used for cosmological analysis with ongoing and future weak lensing surveys.

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 present results of weak lensing cluster counts obtained from 11 sq.deg SuprimeCam data. Although the area is much smaller than previous work dealing with weak lensing peak statistics, the number density of galaxies usable for weak lensing analysis is about twice as large as those. The higher galaxy number density reduces the noise in the weak lensing mass maps, and thus increases the signal-to-noise ratio of peaks of the lensing signal due to massive clusters. This enables us to construct a weak lensing selected cluster sample by adopting a high threshold S/N, such that the contamination rate due to false signals is small. We find 6 peaks with S/N>5. For all the peaks, previously identified clusters of galaxies are matched within a separation of 1 arcmin, demonstrating good correspondence between the peaks and clusters of galaxies. We evaluate the statistical error using mock weak lensing data, and find Npeak=6+/-3.1 in an effective area of 9.0 sq.deg. We compare the measured weak lensing cluster counts with the theoretical model prediction based on halo models and place the constraint on Omega_m-sigma_8 plane which is found to be consistent with currently standard LCDM models. It is demonstrated that the weak lensing cluster counts can place a unique constraint on sigma_8-c_0 plane, where c_0 is the normalization of the dark matter halo mass-concentration relationship. Finally we discuss prospects for ongoing/future wide field optical galaxy surveys.

We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density.

We present the results of an X-ray follow-up campaign targeting 10 Weak Lensing (WL) selected galaxy clusters from a Subaru WL survey. Eight clusters were studied with dedicated Chandra pointings, whereas archival X-ray data were used for the remaining two clusters. The WL clusters appear to fit the same scaling relation between X-ray luminosity and temperature as X-ray selected clusters. However, when we consider the luminosity-mass relation, the WL selected clusters appear underluminous by a factor 3.9$\pm$0.9 (or, alternatively, more massive by 2.9$\pm$0.2), compared to X-ray selected clusters. Only by considering various observational effects that could potentially bias WL masses, can this difference be reconciled. We used X-ray imaging data to quantify the dynamical state of the clusters and found that one of the clusters appears dynamically relaxed, and two of the clusters host a cool core, consistent with Sunyaev-Zel'dovich effect selected clusters. Our results suggest that regular, cool core clusters may be over-represented in X-ray selected samples.

We examined the anisotropic point spread function (PSF) of Suprime-Cam data utilizing dense star field data. We decomposed the PSF ellipticities into three components, the optical aberration, atmospheric turbulence, and chip-misalignment in an empirical manner, and evaluated the amplitude of each component. We found that, for long-exposure data, the optical aberration has the largest contribution to the PSF ellipticities, which could be modeled well by a simple analytic function based on the lowest-order aberration theory. The statistical properties of PSF ellipticities resulting from the atmospheric turbulence were investigated by using the numerical simulations. The simulation results are in a reasonable agreement with the observed data. It is also found that the optical PSF can be well corrected by the standard correction method with a polynomial fitting function. However, for the atmospheric PSF, its correction is affected by the common limitation caused by sparse sampling of PSFs due to a limited number of stars. We also examined the effects of the residual PSF anisotropies on Suprime-Cam cosmic shear data. We found that the shape and amplitude of the B-mode shear variance are broadly consistent with those of the residual PSF ellipticities measured from the dense star field data. This indicates that most of the sources of residual systematic are understood, which is an important step for cosmic shear statistics to be a practical tool of the precision cosmology.

Sky masking is unavoidable in wide-field weak lensing observations. We study how masks affect the measurement of statistics of matter distribution probed by weak gravitational lensing. We first use 1000 cosmological ray-tracing simulations to examine in detail the impact of masked regions on the weak lensing Minkowski Functionals (MFs). We consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks increase the variance of the convergence field and the expected values of the MFs are biased. The bias then affects the non-Gaussian signals induced by the gravitational growth of structure. We then explore how masks affect cosmological parameter estimation. We calculate the cumulative signal-to-noise ratio S/N for masked maps to study the information content of lensing MFs. We show that the degradation of S/N for masked maps is mainly determined by the effective survey area. We also perform simple chi^2 analysis to show the impact of lensing MF bias due to masked regions. Finally, we compare ray-tracing simulations with data from a Subaru 2 deg^2 survey in order to address if the observed lensing MFs are consistent with those of the standard cosmology. The resulting chi^2/n_dof = 29.6/30 for three combined MFs, obtained with the mask effects taken into account, suggests that the observational data are indeed consistent with the standard LambdaCDM model. We conclude that the lensing MFs are powerful probe of cosmology only if mask effects are correctly taken into account.

Weak lensing provides an important route toward collecting samples of clusters of galaxies selected by mass. Subtle systematic errors in image reduction can compromise the power of this technique. We use the B-mode signal to quantify this systematic error and to test methods for reducing this error. We show that two procedures are efficient in suppressing systematic error in the B-mode: (1) refinement of the mosaic CCD warping procedure to conform to absolute celestial coordinates and (2) truncation of the smoothing procedure on a scale of 10$^{\prime}$. Application of these procedures reduces the systematic error to 20% of its original amplitude. We provide an analytic expression for the distribution of the highest peaks in noise maps that can be used to estimate the fraction of false peaks in the weak lensing $\kappa$-S/N maps as a function of the detection threshold. Based on this analysis we select a threshold S/N = 4.56 for identifying an uncontaminated set of weak lensing peaks in two test fields covering a total area of $\sim 3$deg$^2$. Taken together these fields contain seven peaks above the threshold. Among these, six are probable systems of galaxies and one is a superposition. We confirm the reliability of these peaks with dense redshift surveys, x-ray and imaging observations. The systematic error reduction procedures we apply are general and can be applied to future large-area weak lensing surveys. Our high peak analysis suggests that with a S/N threshold of 4.5, there should be only 2.7 spurious weak lensing peaks even in an area of 1000 deg$^2$ where we expect $\sim$ 2000 peaks based on our Subaru fields.

We study the prospects for measuring the dark matter distribution of voids with stacked weak lensing. We select voids from a large set of $N$-body simulations, and explore their lensing signals with the full ray-tracing simulations including the effect of the large-scale structure along the line-of-sight. The lensing signals are compared with simple void model predictions to infer the three-dimensional mass distribution of voids. We show that the stacked weak lensing signals are detected at significant level (S/N$\geq5$) for a 5000 degree$^2$ survey area, for a wide range of void radii up to $\sim50$ Mpc. The error from the galaxy shape noise little affects lensing signals at large scale. It is also found that dense ridges around voids have a great impact on the weak lensing signals, suggesting that proper modeling of the void density profile including surrounding ridges is essential for extracting the average total underdens mass of voids.

We study the cosmological information contained in the Minkowski Functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local type primordial non-Gaussianity parameter f_NL. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing, to generate 100 independent convergence maps of 25 deg^2 field-of-view for f_NL = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude sigma_8. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam survey with the proposed survey area of 1500 deg^2, the primordial non-Gaussianity can be constrained with a level of f_NL ~ 80 and w ~ 0.036 by weak lensing MFs. If simply scaled by the effective survey area, a 20000 deg^2 lensing survey using Large Synoptic Survey Telescope will give constraints of f_NL ~ 25 and w ~ 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.

We examine scatter and bias in weak lensing selected clusters, employing both an analytic model of dark matter haloes and numerical mock data of weak lensing cluster surveys. We pay special attention to effects of the diversity of dark matter distributions within clusters. We find that peak heights of the lensing convergence map correlates rather poorly with the virial mass of haloes. The correlation is tighter for the spherical overdensity mass with a higher mean interior density. We examine the dependence of the halo shape on the peak heights, and find that the rms scatter caused by the halo diversity scales linearly with the peak heights with the proportionality factor of 0.1-0.2. The noise originated from the halo shape is found to be comparable to the source galaxy shape noise and the cosmic shear noise. We find the significant halo orientation bias, i.e., weak lensing selected clusters on average have their major axes aligned with the line-of-sight direction. We compute the orientation bias using an analytic triaxial halo model and obtain results quite consistent with the ray-tracing results. We develop a prescription to analytically compute the number count of weak lensing peaks taking into account all the main sources of scatters in peak heights. We find that the improved analytic predictions agree well with the simulation results for high S/N peaks. We also compare the expected number count with our weak lensing analysis results for 4 sq deg of Subaru/Suprime-Cam observations and find a good agreement.

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

Using a large set of ray-tracing in N-body simulations, we examine lensing profiles around massive dark haloes in detail, with a particular emphasis on the profile at around the virial radii. We compare radial convergence profiles, which are measured accurately in the ray-tracing simulations by stacking many dark haloes, with our simple analytic model predictions. Our analytic models consist of a main halo, which is modelled by the Navarro-Frenk-White (NFW) density profile with three different forms of the truncation, plus the correlated matter (2-halo term) around the main halo. We find that the smoothly truncated NFW profile best reproduces the simulated lensing profiles, out to more than 10 times the virial radius. We then use this analytic model to investigate potential biases in cluster weak lensing studies in which a single, untruncated NFW component is usually assumed in interpreting observed signals. We find that cluster masses, inferred by fitting reduced tangential shear profiles with the NFW profile, tend to be underestimated by ~5-10% if fitting is performed out to ~10'-30'. In contrast, the concentration parameter is overestimated typically by ~20% for the same fitting range. We also investigate biases in computing the signal-to-noise ratio of weak lensing mass peaks, finding them to be <4% for significant mass peaks. In the Appendices, we provide useful formulae for the smoothly truncated NFW profile.

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.

Galaxy clusters are usually detected in blind optical surveys via suitable filtering methods. We present an optimal matched filter which maximizes their signal-to-noise ratio by taking advantage of the knowledge we have of their intrinsic physical properties and of the data noise properties. In this paper we restrict our application to galaxy magnitudes, positions and photometric redshifts if available, and we also apply the filter separately to weak lensing data. The method is suitable to be naturally extended to a multi-band approach which could include not only additional optical bands but also observables with different nature such as X-rays. For each detection, the filter provides its significance, an estimate for the richness and for the redshift even if photo-z are not given. The provided analytical error estimate is tested against numerical simulations. We finally apply our method to the COSMOS field and compare the results with previous cluster detections obtained with different methods. Our catalogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including both optical and weak-lensing information.

We present deep J-, H-, and Ks-band imaging data of the MOIRCS Deep Survey (MODS), which was carried out with Multi-Object Infrared Camera and Spectrograph (MOIRCS) mounted on the Subaru telescope in the GOODS-North region. The data reach 5sigma total limiting magnitudes for point sources of J=23.9, H=22.8, and Ks=22.8 (Vega magnitude) over 103 arcmin^2 (wide field). In 28 arcmin^2 of the survey area, which is ultra deep field of the MODS (deep field), the data reach the 5sigma depths of J=24.8, H=23.4, and Ks=23.8. The spatial resolutions of the combined images are FWHM ~ 0.6 arcsec and ~ 0.5 arcsec for the wide and deep fields in all bands, respectively. Combining the MODS data with the multi-wavelength public data taken with the HST, Spitzer, and other ground-based telescopes in the GOODS field, we construct a multi-wavelength photometric catalog of Ks-selected sources. Using the catalog, we present Ks-band number counts and near-infrared color distribution of the detected objects, and demonstrate some selection techniques with the NIR colors for high redshift galaxies. These data and catalog are publicly available via internet.

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

(Abridged) We present results of Suzaku observations of the intracluster medium (ICM) in Abell 1689, combined with complementary analysis of the SDSS data and weak and strong lensing analysis of Subaru/Suprime-Cam and HST/ACS observations. Faint X-ray emission from the ICM around the virial radius is detected at 4.0 sigma significance. We find anisotropic gas temperature and entropy distributions in cluster outskirts correlated with large-scale structure of galaxies. The high temperature and entropy region in the northeastern (NE) outskirts is connected to an overdense filamentary structure. The outskirt regions in contact with low density void environments have low gas temperatures and entropies, deviating from hydrostatic equilibrium. These results suggest that thermalization of the ICM occurs faster along the filamentary structures than the void regions. A joint X-ray and lensing analysis shows that the hydrostatic mass is $\sim60-90%$ of spherical lensing one but comparable to a triaxial halo mass within errors in $0.6r_{2500} \simlt r \simlt 0.8r_{500}$, and that it is significantly biased as low as $\simlt60%$ within $0.4r_{2500}$, irrespective of mass models. The thermal gas pressure within $r_{500}$ is, at most, $\sim50$--60% of the total pressure to balance fully the gravity of the spherical lensing mass, and $\sim30$--40% around the virial radius. Although these constitute lower limits when one considers the possible halo triaxiality, these small relative contributions of thermal pressure would require additional sources of pressure, such as bulk and/or turbulent motions.

We use very deep near-infrared (NIR) imaging data obtained in MOIRCS Deep Survey (MODS) to investigate the evolution of the galaxy stellar mass function back to z~3. The MODS data reach J=24.2, H=23.1, K=23.1 (5sigma, Vega magnitude) over 103 arcmin^2 (wide) and J=25.1, H=23.7, K=24.1 over 28 arcmin^2 (deep) in the GOODS-North region. The wide and very deep NIR data allow us to measure the number density of galaxies down to low stellar mass (10^9-10^10 Msun) even at high redshift with high statistical accuracy. The normalization of the mass function decreases with redshift and the integrated stellar mass density becomes ~ 8-18% of the local value at z~2 and ~ 4-9% at z~3, which are consistent with results of previous studies in general fields. Furthermore, we found that the low-mass slope becomes steeper with redshift from alpha ~- 1.3 at z~1 to alpha ~- 1.6 at z~3, and that the evolution of the number density of low-mass (10^9-10^10 Msun) galaxies is weaker than that of M* (~10^11 Msun) galaxies. This indicates that the contribution of low-mass galaxies to the total stellar mass density has been significant at high redshift. The steepening of the low-mass slope with redshift is opposite trend expected from the stellar mass dependence of the specific star formation rate reported in previous studies. The present result suggests that the hierarchical merging process overwhelmed the effect of the stellar mass growth by star formation and was very important for the stellar mass assembly of these galaxies at 1<~z<~3.

We combine measurements of weak gravitational lensing from the CFHTLS-Wide survey, supernovae Ia from CFHT SNLS and CMB anisotropies from WMAP5 to obtain joint constraints on cosmological parameters, in particular, the dark energy equation of state parameter w. We assess the influence of systematics in the data on the results and look for possible correlations with cosmological parameters. We implement an MCMC algorithm to sample the parameter space of a flat CDM model with a dark-energy component of constant w. Systematics in the data are parametrised and included in the analysis. We determine the influence of photometric calibration of SNIa data on cosmological results by calculating the response of the distance modulus to photometric zero-point variations. The weak lensing data set is tested for anomalous field-to-field variations and a systematic shape measurement bias for high-z galaxies. Ignoring photometric uncertainties for SNLS biases cosmological parameters by at most 20% of the statistical errors, using supernovae only; the parameter uncertainties are underestimated by 10%. The weak lensing field-to-field variance pointings is 5%-15% higher than that predicted from N-body simulations. We find no bias of the lensing signal at high redshift, within the framework of a simple model. Assuming a systematic underestimation of the lensing signal at high redshift, the normalisation sigma_8 increases by up to 8%. Combining all three probes we obtain -0.10<1+w<0.06 at 68% confidence (-0.18<1+w<0.12 at 95%), including systematic errors. Systematics in the data increase the error bars by up to 35%; the best-fit values change by less than 0.15sigma. [Abridged]

We present the first results of a MOS campaign to follow up cluster candidates located via weak lensing. Our main goals are to search for spatial concentrations of galaxies that are plausible optical counterparts of the weak lensing signals, and to determine the cluster redshifts from those of member galaxies. Around each of 36 targeted cluster candidates, we obtain 15-32 galaxy redshifts. For 28 of these targets, we confirm a secure cluster identification. This includes three cases where two clusters at different redshifts are projected along the same line-of-sight. In 6 of the 8 unconfirmed targets, we find multiple small galaxy concentrations at different redshifts. In both the remaining two targets, a single small galaxy concentration is found. We evaluate the weak lensing mass of confirmed clusters. For a subsample of our most cleanly measured clusters, we investigate the statistical relation between their weak lensing mass and the velocity dispersion of their member galaxies, comparing our sample with optically and X-ray selected samples from the literature. Our lensing-selected clusters are consistent with sigma_v=sigma_sis, with a similar scatter to the optically and X-ray selected clusters. We thus find no evidence of selection bias compared to these other techniques. We also derive an empirical relation between the cluster mass and the galaxy velocity dispersion, which is in reasonable agreement with the prediction of N-body simulations in the LCDM cosmology.

We qualitatively examine properties of artificial deformation in shapes of objects (galaxies and stars) induced by the pixelization effects (also called as the aliasing effects) using toy mock simulation images. We pay a special attention to the the second pixelization because it might be a potential source of a systematic noise in a weak lensing analysis. In particular, it is found that resampling with rotation induces artificial ellipticities in object shapes having a periodic concentric-circle-shaped pattern. Our major findings are as follows. (1) Root-mean-square (RMS) of artificial ellipticities in object shapes induced by the first pixelization effect can be as large as RMS > 10^-2 if a characteristic size of objects (e.g., the FWHM) is smaller than twice of the pixel size. While for larger objects, it quickly becomes very small (RMS < 10^-5). (2) The amplitude of the shape deformation induced by the second pixelization effect depends on the object size. It also depends strongly on an interpolation scheme adopted to carry out resampling and on the grid size of the output pixels. The RMS of ellipticities in object shapes induced by the second pixelization effect can be suppressed to well below 10^{-2} if one adopts a proper interpolation scheme (implemented in popular image processing softwares). We also discuss an impact of the pixelization effects on a weak lensing analysis.

We consider the extent to which future imaging surveys of galaxies can distinguish between dark energy and modified gravity models for the origin of the cosmic acceleration. Dynamical dark energy models may have similar expansion rates as models of modified gravity, yet predict different growth of structure histories. We parameterize the cosmic expansion by the two parameters, $w_0$ and $w_a$, and the linear growth rate of density fluctuations by Linder's $\gamma$, independently. Dark energy models generically predict $\gamma \approx 0.55$, while the DGP model $\gamma \approx 0.68$. To determine if future imaging surveys can constrain $\gamma$ within 20 percent (or $\Delta\gamma<0.1$), we perform the Fisher matrix analysis for a weak lensing survey such as the on-going Hyper Suprime-Cam (HSC) project. Under the condition that the total observation time is fixed, we compute the Figure of Merit (FoM) as a function of the exposure time $\texp$. We find that the tomography technique effectively improves the FoM, which has a broad peak around $\texp\simeq {\rm several}\sim 10$ minutes; a shallow and wide survey is preferred to constrain the $\gamma$ parameter. While $\Delta\gamma < 0.1$ cannot be achieved by the HSC weak-lensing survey alone, one can improve the constraints by combining with a follow-up spectroscopic survey like WFMOS and/or future CMB observations.

We present the results of an ongoing weak lensing survey conducted with the Subaru telescope whose initial goal is to locate and study the distribution of shear-selected structures or halos. Using a Suprime-cam imaging survey spanning 21.82 square degree, we present a catalog of 100 candidate halos located from lensing convergence maps. Our sample is reliably drawn from that subset of our survey area, (totaling 16.72 square degree) uncontaminated by bright stars and edge effects and limited at a convergence signal to noise ratio of 3.69. To validate the sample detailed spectroscopic measures have been made for 26 candidates using the Subaru multi-object spectrograph, FOCAS. All are confirmed as clusters of galaxies but two arise as the superposition of multiple clusters viewed along the line of sight. Including data available in the literature and an ongoing Keck spectroscopic campaign, a total of 41 halos now have reliable redshifts. For one of our survey fields, the XMM LSS (Pierre et al. 2004) field, we compare our lensing-selected halo catalog with its X-ray equivalent. Of 15 halos detected in the XMM-LSS field, 10 match with published X-ray selected clusters and a further 2 are newly-detected and spectroscopically confirmed in this work. Although three halos have not yet been confirmed, the high success rate within the XMM-LSS field (12/15) confirms that weak lensing provides a reliable method for constructing cluster catalogs, irrespective of the nature of the constituent galaxies or the intracluster medium.

We present the first measurement of clustering properties of low mass galaxies with a stellar mass down to M_*~10^9 Msun at 1<z<4 in 24.4 arcmin^2 of the GOODS-North region with a depth of K_{AB}~25, based on the near infrared observations performed with MOIRCS at the Subaru Telescope. The correlation amplitude strongly depends on the K-band flux, color, and stellar mass of the galaxies. We find that K-band luminous galaxies have a larger correlation length than K-band faint galaxies. For color selected samples at 2<z<4, distant red galaxies with J-K>1.3 show a large bias of b~7.2+-1.3 on scales up to \theta~100" or 3.1 comoving Mpc, while blue galaxies with 0.5<J-K<1.3, in which most Lyman break galaxies are populated, have a weak clustering signal on large scales, but a possible strong small scale excess at \theta<10". For massive galaxies with M_*>~10^{10} Msun, we estimate a correlation length and bias to be r_0~4.5 h^{-1} Mpc and b=1.9-3.5, which are much larger than those of low mass (M_*~10^9-10^{10} Msun) galaxies. The comparison of our measurements with analytic CDM models constrains the properties of hosting dark halos, and indicates that the low mass galaxies would be progenitors of galaxies with a typical luminosity of L<~L_* in the local Universe. The blue galaxies in low mass samples are more strongly clustered in more massive halos with higher occupation numbers than low mass red galaxies. This fact suggests an environment effect due to the halo mass on star formation activity at high-z.

The results from weak gravitational lensing analyses are subject to a cosmic variance error term that has previously been estimated assuming Gaussian statistics. In this letter we address the issue of estimating cosmic variance errors for weak lensing surveys in the non-Gaussian regime. Using standard cold dark matter model ray-tracing simulations characterized by Omega_m=0.3, Omega_Lambda=0.7, h=0.7, sigma_8=1.0 for different survey redshifts z_s, we determine the variance of the two-point shear correlation function measured across 64 independent lines of sight. We compare the measured variance to the variance expected from a random Gaussian field and derive a redshift-dependent non-Gaussian calibration relation. We find that the ratio can be as high as ~30 for a survey with source redshift z_s ~ 0.5 and ~10 for z_s ~ 1. The transition scale theta_c above which the ratio is consistent with unity, is found to be theta_c ~ 20 arcmin for z_s ~ 0.5 and theta_c ~ 10 arcmin for z_s ~ 1. We provide fitting formula to our results permitting the estimation of non-Gaussian cosmic variance errors for any weak lensing analysis, and discuss the impact on current and future surveys. A more extensive set of simulations will however be required to investigate the dependence of our results on cosmology, specifically on the amplitude of clustering.

We report angular correlation function (ACF) of Lyman Break Galaxies (LBGs) with unprecedented statistical quality on the basis of 16,920 LBGs at z=4 detected in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field. The ACF significantly departs from a power law, and shows an excess on small scale. Particularly, the ACF of LBGs with i'<27.5 have a clear break between the small and large-scale regimes at the angular separation of ~7'' whose projected length corresponds to the virial radius of dark halos with a mass of 10^11-12 Mo, indicating multiple LBGs residing in a single dark halo. Both on small (2''<theta<3'') and large (40''<theta<400'') scales, clustering amplitudes monotonically increase with luminosity for the magnitude range of i'=24.5-27.5, and the small-scale clustering shows a stronger luminosity dependence than the large-scale clustering. The small-scale bias reaches b~10-50, and the outskirts of small-scale excess extend to a larger angular separation for brighter LBGs. The ACF and number density of LBGs can be explained by the cold dark matter model.

We examine masses of hosting haloes of two photometrically-selected high-z galaxy samples: the old passively-evolving galaxies (OPEGs) and Lyman Break Galaxies (LBGs) both taken from the Subaru/XMM-Newton Deep Survey (SXDS). The large survey area of the SXDS (1sq deg) allows us to measure the angular two-point correlation functions to a wide separation of >10 arcmin with a good statistical quality. We utilize the halo model prescription for estimating characteristic masses of hosting haloes from the measured large-scale clustering amplitudes. It is found that the hosting halo mass positively correlates with the luminosity of galaxies. Then, adopting the extended Press-Schechter model (EPS), we compute the predictions for the mass evolution of the hosting haloes in the framework of the cold dark matter (CDM) cosmology in order to make an evolutionary link between the two galaxy samples at different redshifts and to identify their present-day descendants by letting their haloes evolve forward in time. It is found that, in the view of the mass evolution of hosting haloes in the CDM model, bright LBGs are consistent with being the progenitor of the OPEGs, whereas it is less likely that the LBG population, as a whole, have evolved into the OPEG population. It is also found that the present-day descendants of both the bright LBGs and OPEGs are likely to be located in massive systems such as groups of galaxies or clusters of galaxies. Finally, we estimate the hosting halo mass of local early-type galaxy samples from the 2dF and SDSS based on the halo model and it turns out that their expected characteristic mass of hosting haloes is in good agreement with the EPS predictions for the descendant's mass of both the bright LBGs and OPEGs.

We present Subaru spectroscopy of the early-type lensing galaxies in the lens systems HST14113+5211 (z=0.464) and B2045+265 (z=0.868), being aimed at measuring the velocity dispersion of the lensing galaxies as an important component to their mass distributions and internal dynamics. For HST14113+5211 we have obtained 174+/-20km/s (1-sigma) inside an aperture of 0.9", and for B2045+265 we have obtained 213+/-23km/s inside an aperture of 1.2". To extract the significance of these information on the mass distributions and stellar dynamics of the lensing galaxies, we construct two alternative models for both the mass and velocity distributions. It is found that the mass-to-light ratios as derived from a joint lensing and stellar dynamics analysis are virtually consistent with those obtained from the use of the Fundamental Plane. Also, for HST14113+5211, we find that the total mass distribution is well reproduced by a power-law density profile with an index of 2, thereby suggesting that a singular isothermal model is a good fit to the lensing galaxy. In contrast, the corresponding slope for B2045+265 is flatter than isothermal, suggesting that additional contributions to lensing mass from surrounding structures (possibly the group of galaxies) provide the observed angular separations of multiple images.

We present a simple model that accurately describes various statistical properties of peculiar velocities of dark matter halos. We pay particular attention to the following two effects; first, the evolution of the halo peculiar velocity depends on the local matter density, instead of the global density. Second, dark matter halos are biased tracers of the underlying mass distribution, thus halos tend to be located preferentially at high density regions. For the former, we develop an empirical model calibrated with N-body simulations, while for the latter, we use a conventional halo bias models based on the extended Press-Schechter model combined with an empirical log-normal probability distribution function of the mass density distribution. We find that compared with linear theory, the present model significantly improves the accuracy of predictions of statistical properties of the halo peculiar velocity field including the velocity dispersion, the probability distribution function, and the pairwise velocity dispersion at large separations. Thus our model predictions may be useful in analyzing future observations of the peculiar velocities of galaxy clusters.

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, QSOs and also the CMB temperature map), is known to consist of two components; a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a Lambda-CDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating those two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive halos with M>10^14h^-1Msun at z<1 and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arc minute, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at sub-arcminute scales.

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 explore empirical relations between three different populations of high-redshift galaxies and their hosting dark halos employing the halo model approach. Specifically we consider Lyman-break galaxies at z\sim4 and at z\sim5, and Lyman-Alpha emitters at z\simeq 4.86. We adopt a halo occupation function prescription to parameterize the properties of their hosting halos and the efficiency of halo-dependent star formation. We find that the two LBG samples are well described by the halo model with an appropriate HOF. We obtain constraints on properties of their hosting halos. A typical mass of hosting halos for LBGs is 5\times10^{11}h^{-1}M_\odot and the expected number of LBGs per halo is \sim0.5, therefore there is an approximate one-to-one correspondence between halos and LBGs. We also find a sign of the minimum mass of LBG hosting halos decreasing with time. We discuss implications of these findings on the star formation history of LBGs. On the other hand, for LAEs, our simple HOF prescription fails to reproduce simultaneously the observed angular correlation function and the number density. This might imply either that the distribution of LAEs within hosting halos differs from that of dark matter, or that the strong large-scale correlation is due to the existence of an unusual, large overdense region, and so the survey region is not a representative of the z\sim5 universe, the definite answer should wait for a much wider survey of LAEs at high redshifts.

We present measurements of the cosmic shear correlation in the shapes of galaxies in the Suprime-Cam 2.1 deg^2 R_c-band imaging data. As an estimator of the shear correlation originated from the gravitational lensing, we adopt the aperture mass variance. We detect a non-zero E mode variance on scales between 2 and 40arcmin. We also detect a small but non-zero B mode variance on scales larger than 5arcmin. We compare the measured E mode variance to the model predictions in CDM cosmologies using maximum likelihood analysis. A four-dimensional space is explored, which examines sigma_8, Omega_m, Gamma and zs (a mean redshift of galaxies). We include three possible sources of error: statistical noise, the cosmic variance estimated using numerical experiments, and a residual systematic effect estimated from the B mode variance. We derive joint constraints on two parameters by marginalizing over the two remaining parameters. We obtain an upper limit of Gamma<0.5 for zs>0.9 (68% confidence). For a prior Gamma\in[0.1,0.4] and zs\in[0.6,1.4], we find sigma_8=(0.50_{-0.16}^{+0.35})Omega_m^{-0.37} for flat cosmologies and sigma_8=(0.51_{-0.16}^{+0.29})Omega_m^{-0.34}$ for open cosmologies (95% confidence). If we take the currently popular LCDM model, we obtain a one-dimensional confidence interval on sigma_8 for the 95.4% level, 0.62<\sigma_8<1.32 for zs\in[0.6,1.4]. Information on the redshift distribution of galaxies is key to obtaining a correct cosmological constraint. An independent constraint on Gamma from other observations is useful to tighten the constraint.

We present explicit expressions for the lens equation for a cored isothermal ellipsoidal gravitational lens as a single real sixth-order algebraic equation in two approaches; 2-dimensional Cartesian coordinates and 3-dimensional polar ones. We find a condition for physical solutions which correspond to at most five images. For a singular isothermal ellipsoid, the sixth-order equation is reduced to fourth-order one for which analytic solutions are well-known. Furthermore, we derive analytic criteria for determining the number of images for the singular lens, which give us simple expressions for the caustics and critical curves. The present formulation offers a useful way for studying galaxy lenses frequently modeled as isothermal ellipsoids.

We present a simple model based on the dark halo approach which provides a useful way to understand key points determining the shape of the non-Gaussian tails of the dark matter one-point probability distribution function(PDF). In particular, using the scale-free models with power-law profile of dark halos, we derive a simple analytic expression for the one-point PDF. It is found that the shape of the PDF changes at the characteristic value of $\delta_*$ which is defined by the smoothed density of a halo with the characteristic mass $M_*$ at the epoch. In cold dark matter models with top-hat smoothing filters, the characteristic smoothed density at present time typically takes the value $\delta_*\gg 1$ for a small smoothing scale $\rth\sim 1$Mpc$/h$ and conversely $\delta_*\ll 1$ for a large smoothing scale $\rth > 10$Mpc$/h$. On the range $\delta/\delta_*<1$, the shape of the PDF is almost solely determined by the outer slope of halos and scales as a power-law. The resultant non-Gaussian tails of PDF then resemble the log-normal PDFs in that range and show a good agreement with N-body simulations, which can be ascribed to the universality of the outer slope of the halo profile. In contrast, tails of one-point PDF in the range $\delta/\delta_*>1$ basically follow the steep exponential tails of the halo mass function, which exhibit a strong sensitivity to both the outer slope of the halo profile and the initial power spectrum. Based on these results, the discussion on the PDF of galaxy distribution and the application to the weak lensing statistics are also presented.

We report the first result of weak gravitational lensing survey on a 2.1 sq deg Rc-band image taken with a wide field camera (Suprime-Cam) on the prime focus of 8.2 m Subaru Telescope. The weak lensing mass reconstruction is applied to the data to search for dark matter halos of cluster scale; M >= 10^14 solar mass. The reconstructed convergence field is divided by 1-sigma noise to obtain the signal-to-noise ratio map (S/N-map) of the detection. Local maxima and minima are searched on the S/N-map and the probability distribution function (PDF) of the peaks are created to compare with model predictions. We found excess over noise PDF created from the randomized realization on both positive and negative sides. Negative peaks imply the presence of voids in the dark matter distribution and this is the first report of the detection. Positive peaks, on the other hand, represent the dark matter halos and the number count of the halos on the 2.1 sq deg image is 4.9 +- 2.3 for S/N > 5 where the Gaussian smoothing radius of the convergence map is 1'. The result is consistent with the prediction assuming the Press-Schechter mass function and the NFW halo profile under the cluster normalized CDM cosmology. Although the present statistics is not enough due to the limited field of view, this work demonstrates that dark matter halo count via weak lensing is a promising way to test the paradigm of structure formation and cosmological model.

The systematic magnification of background sources by the weak gravitational-lensing effects of foreground matter, also called cosmic magnification, is becoming an efficient tool both for measuring cosmological parameters and for exploring the distribution of galaxies relative to the dark matter. We extend here the formalism of magnification statistics by estimating the contribution of second-order terms in the Taylor expansion of the magnification and show that the effect of these terms was previously underestimated. We test our analytical predictions against numerical simulations and demonstrate that including second-order terms allows the accuracy of magnification-related statistics to be substantially improved. We also show, however, that both numerical and analytical estimates can provide only lower bounds to real correlation functions, even in the weak lensing regime. We propose to use count-in-cells estimators rather than correlation functions for measuring cosmic magnification since they can more easily be related to correlations measured in numerical simulations.

We critically examine the reliability of the dark matter clustering in high-resolution cosmological $N$-body simulations on scales below the mean separation length of particles. The particle discreteness effect imposes the two fundamental limitations on those scales; the lack of the initial fluctuation power and the finite mass resolution. We address this problem applying the dark halo approach and are able to discuss separately how those two limitations affect the dark matter clustering in $N$-body simulations at early epochs. We find that limitations of the dark matter clustering are primarily determined by the mass of particles. By a detailed comparison with three major cosmological simulations, we also find that in order to reproduce a proper amplitude of the dark matter clustering on small scales, halos with a characteristic nonlinear mass, $M_{\rm NL}(z)$ defined by $\sigma_R(M_{\rm NL}; z)=1$, must be resolved in the simulation. This leads to a critical redshift $z_{\rm crit}$ determined by $M_{\rm NL}(z_{\rm crit})=n_{\rm halo} m_{\rm part}$ where $n_{\rm halo}$ is the number of particles necessary to resolve the typical nonlinear mass halo ($\sim 10$). We conclude that, at least as far as the two-point correlation functions are concerned, the dark matter clustering in high-resolution $N$-body simulations on scales below the mean particle separation is reliable down to the gravitational force resolution length only for $z<z_{\rm crit}$, while it is strongly affected by the finite mass resolution for $z>z_{\rm crit}$.

We present new optical spectroscopy of the lens elliptical galax in the ``Einstein Cross'' lens system HST 14176+5226, using the Faint Object Camera and Spectrograph (FOCAS) of the Subaru t Our spectroscopic observations are aimed at measuring the stella dispersion of the lens galaxy, located at high redshift of z_L= as an important component to lens models. We have measured this 230 +- 14 km s^{-1} (1 sigma) inside 0.35 effective radi based on the comparison between the observed galaxy spectrum and templates of three G-K giants by means of the Fourier cross-corr To extract the significance of this information on the geometry universe which also affects the lensing of the background image, to fit three different lens models to the available data of the Provided that the lens galaxy has the structural and dynamical p (i.e., its radial density profile, core radius, and velocity ani similar to those of local elliptical galaxies, we calculate the function for the simultaneous reproduction of both the observed and newly measured velocity dispersion of the lens. Although the interval depends rather sensitively on the adopted lens models o parameters, our experiments suggest the larger likelihood for a cosmological constant, Omega_Lambda: formal 1 sigma lower Omega_Lambda in the flat universe ranges 0.73 to 0.97, where lower limit is basically unavailable. This method for determinin model is thus dependent on lens models but is insensitive to oth ambiguities, such as the dust absorption or the evolutionary eff galaxies. Exploring spectroscopic observations of more lens gala redshift may minimize the model uncertainties and thus place a m constraint on Omega_Lambda.

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 correlation between source galaxies and lensing potentials causes a systematic effect on measurements of cosmic shear statistics, known as the source-lens clustering (SLC) effect. The SLC effect on the skewness of lensing convergence, $S_3$, is examined using a nonlinear semi-analytic approach and is checked against numerical simulations. The semi-analytic calculations have been performed in a wide variety of generic models for the redshift distribution of source galaxies and power-law models for the bias parameter between the galaxy and dark matter distributions. The semi-analytic predictions are tested successfully against numerical simulations. We find the relative amplitude of the SLC effect on $S_3$ to be of the order of five to forty per cent. It depends significantly on the redshift distribution of sources and on the way the bias parameter evolves. We discuss possible measurement strategies to minimize the SLC effects.

We develop a phenomenological model to predict the clustering of dark matter halos on the light-cone by combining several existing theoretical models. Assuming that the velocity field of halos on large scales is approximated by linear theory, we propose an empirical prescription of a scale-, mass-, and time-dependence of halo biasing. We test our model against the Hubble Volume $N$-body simulation and examine its validity and limitations. We find a good agreement in two-point correlation functions of dark matter halos between the phenomenological model predictions and measurements from the simulation for $R>5h^{-1}$Mpc both in the real and redshift spaces. Although calibrated on the mass scale of groups and clusters and for redshifts up to $z\sim2$, the model is quite general and can be applied to a wider range of astrophysical objects, such as galaxies and quasars, if the relation between dark halos and visible objects is specified.

We present ray tracing simulations combined with sets of large N-body simulations. Experiments were performed to explore, for the first time, statistical properties of fluctuations in angular separations of nearby light ray pairs (the so-called lensing excursion angle) induced by weak lensing by large-scale structures. We found that the probability distribution function (PDF) of the lensing excursion angles is not simply Gaussian but has an exponential tail. It is, however, found that the tail, or more generally non-Gaussian nature in the PDF has no significant impact on the weak lensing of the CMB. Moreover, we found that the variance in the lensing excursion angles predicted by the power spectrum approach is in good agreement with our numerical results. These results demonstrate a validity of using the power spectrum approach to compute lensing effects on the CMB.

Gravitational lensing causes a correlation between a population of foreground large-scale structures and the observed number density of the background distant galaxies as a consequence of the flux magnification and the lensing area distortion. This correlation has not been taken into account in calculations of the theoretical predictions of the cosmic shear statistics but may cause a systematic error in a cosmic shear measurement. We examine its impact on the cosmic shear statistics using the semi-analytic approach. We find that the lensing magnification has no practical influence on the cosmic shear variance. Exploring possible shapes of redshift distribution of source galaxies, we find that the relative amplitude of the effect on the convergence skewness is 3% at most

We examine the light-cone effect on the two-point correlation functions using numerical simulations for the first time. Specifically, we generate several sets of dark matter particle distributions on the light-cone up to z=0.4 and z=2 over the field-of-view of \pi degree^2 from cosmological N-body simulations. Then we apply the selection function to the dark matter distribution according to the galaxy and QSO luminosity functions. Finally we compute the two-point correlation functions on the light-cone both in real and in redshift spaces using the pair-count estimator and compare with the theoretical predictions. We find that the previous theoretical modeling for nonlinear gravitational evolution, linear and nonlinear redshift-distortion, and the light-cone effect including the selection function is in good agreement with our numerical results, and thus is an accurate and reliable description of the clustering in the universe on the light-cone.

We study cosmic shear statistics using the ray-tracing simulation combined with a set of large $N$-body simulations. In this contribution, we first describe our method. Then we show some selected results especially focusing on effects of the deflection of light rays and the lens-lens coupling which are neglected in making the theoretical predictions of the cosmic shear statistics.

As pointed out in previous studies, the measurement of the skewness of the convergence field $\kappa$ will be useful in breaking the degeneracy among the cosmological parameters constrained from weak lensing observations. The combination of shot noise and finite survey volume implies that such a measurement is likely to be done in a range of intermediate scales ($0.5'$ to 20') where neither perturbation theory nor the hierarchical ansatz apply. Here we explore the behavior of the skewness of $\kappa$ at these intermediate scales, based on results for the non-linear evolution of the mass bispectrum. We combined different ray-tracing simulations to test our predictions, and we find that our calculations describe accurately the transition from the weakly non-linear to the strongly non-linear regime. We show that the single lens-plane approximation remains accurate even in the non-linear regime, and we explicitly calculate the corrections to this approximation. We also discuss the prospects of measuring the skewness in upcoming weak lensing surveys.

The gravitational lensing magnification or demagnification due to large-scale structures induces a scatter in peak magnitudes of high redshift type Ia supernovae (SNe Ia). The amplitude of the lensing dispersion strongly depends on that of density fluctuations characterized by the $\sigma_8$ parameter. Therefore the value of $\sigma_8$ is constrained by measuring the dispersion in the peak magnitudes. We examine how well SN Ia data will provide a constraint on the value of $\sigma_8$ using a likelihood analysis method. It is found that the number and quality of SN Ia data needed for placing a useful constraint on $\sigma_8$ is attainable with Next Generation Space Telescope.

Distances play important roles in cosmological observations, especially in gravitational lens systems, but there is a problem in determining distances because they are defined in terms of light propagation, which is influenced gravitationally by the inhomogeneities in the universe. In this paper we first give the basic optical relations and the definitions of different distances in inhomogeneous universes. Next we show how the observational relations depend quantitatively on the distances. Finally, we give results for the frequency distribution of different distances and the shear effect on distances obtained using various methods of numerical simulation.

We compute statistical properties of weak gravitational lensing by large-scale structure in three Cold Dark Matter models. We use a P$^3$M $N$-body code to simulate the formation and evolution of large-scale structure in the universe. We perform $1.1\times10^7$ ray-tracing experiments for each model using the multiple lens-plane algorithm. From the results of these experiments, we calculate the probability distribution functions (PDF) of the convergences, shears, and magnifications, and their root-mean-square (rms) values. We find that the rms values of the convergence and shear agree with the predictions of a nonlinear analytical model. We also find that the PDFs of the magnifications $\mu$ have a peak at values slightly smaller than $\mu=1$, and are strongly skewed toward large magnifications. In particular, for the high-density model, a power-law tail appears in the magnification distribution at large magnifications for sources at redshifts $z_s>2$. The rms values of the magnifications essentially agree with the nonlinear analytical predictions for sources at low redshift, but exceed these predictions for high redshift sources, once the power-law tail appears. We study the effect of magnification bias on the luminosity functions of high-redshift quasars, using the calculated PDFs of the magnifications. We show that the magnification bias is moderate in the absence of the power-law tail in the magnification distribution, but depends strongly on the value of the density parameter. In presence of the power-law tail, the bias becomes considerable, especially at the bright end of the luminosity functions.

Weyl focusing effects on the image magnifications are investigated by using the multiple gravitational lens theory. We focus on the gravitational lensing effects due to the small scale virialized objects, such as galaxies and clusters of galaxies. We consider a simple model of an inhomogeneous universe. The matter distribution in the universe is modeled by randomly distributed isothermal objects. We found that, for the majority of the random lines of sight, the Weyl focusing has no significant effect and the image magnification of a point like source within redshift of 5 is dominated by the Ricci focusing.

We investigate the performance of the optimized Post-Zel'dovich approximation in three cold dark matter cosmologies. We consider two flat models with $\Omega_0=1$ (SCDM) and with $\Omega_0=0.3$ ($\Lambda$CDM) and an open model with $\Omega_0=0.3$ (OCDM). We find that the optimization scheme proposed by Wei{\ss}, Gottl\"ober & Buchert (1996), in which the performance of the Lagrangian perturbation theory was optimized only for the Einstein-de Sitter cosmology, shows the excellent performances not only for SCDM model but also for both OCDM and $\Lambda$CDM models. This universality of the excellent performance of the optimized Post-Zel'dovich approximation is explained by the fact that a relation between the Post-Zel'dovich order's growth factor $E(a)$ and Zel'dovich order's one $D(a)$, $E(a)/D^2(a)$, is insensitive to the background cosmologies.

We rederive a relation between gravitational lensing magnification relative to the standard Friedmann distance and one relative to the Dyer-Roeder distance by investigating the null geodesic deviation equation. We show that the relation comes from a natural consequence of the definition of the lensing magnification matrices and is not based on the averaging of the magnifications, which has conventionally been used to derive it. We therefore conclude that the relation is true for each individual ray bundle.

The amplification of the protogalaxy candidate cB58 due to gravitational lensing by the foreground cluster of galaxies MS1512.4+3647 is quantified based on recent ROSAT and ASCA X-ray observations. It is found that the amplification is at most 25 for any reasonable cosmological model with or without cosmological constant. It is also argued that the system may be used to place new constraints on the value of the cosmological constant. The gas mass fraction for this cluster is found to be about 0.2.