We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05<z_{cl}<1.38$. We discover 641 candidate lens systems, of which 536 are new. From the full sample, 47 are almost certainly bonafide lenses, 181 of them are highly probable lenses and 413 are possible lens systems. Additionally, we present 131 lens candidates at galaxy-scale serendipitously discovered during the inspection. We obtained spectroscopic follow-up of 10 candidates using the X-shooter. With this follow-up, we confirm 8 systems as strong gravitational lenses. Of the remaining two, one of the sources is too faint to detect any emission, and the other has a tentative redshift close to the lens redshift, but additional arcs in this system are yet to be observed spectroscopically. Since the HSC-SSP is an ongoing survey, we expect to find $\sim600$ definite or probable lenses using this procedure and even more if combined with other lens finding methods.

We show that peculiar velocities of Type Ia supernovae can be used to derive constraints on the sum of neutrino masses, $\Sigma m_{\nu}$, and dark energy equation of state, $w = w_0+w_a(1-a)$, from measurements of the magnitude-redshift relation, complementary to galaxy redshift and weak lensing surveys. Light from a supernova propagates through a perturbed Universe so the luminosity distance is modified from its homogeneous prediction. This modification is proportional to the matter density fluctuation and its time derivative due to gravitational lensing and peculiar velocity respectively. At low redshifts, the peculiar velocity signal dominates while at high redshifts lensing does. We show that using lensing and peculiar velocity of supernovae from the upcoming surveys WFIRST and ZTF, without other observations, we can constrain $\Sigma m_{\nu} \lesssim 0.31$ eV, $\sigma(w_0) \lesssim 0.02$, and ${\sigma(w_a)} \lesssim 0.18$ ($1-\sigma$ CL) in the $\Sigma m_{\nu}$-$w_0$-$w_a$ parameter space, where all the other cosmological parameters are fixed. We find that adding peculiar velocity information from low redshifts shrinks the volume of the parameter ellipsoid in this space by $\sim 33$%. We also allow $\Omega_{\text{CDM}}$ to vary as well as $\Sigma m_{\nu}$, $w_0$ and $w_a$, and demonstrate how these constraints degrade as a consequence.

We investigate how the cosmological constraints from SNe Ia are improved by including the effects of weak-lensing convergence. To do so, we introduce the lognormal function as the convergence PDF modeling the lensing scatter of SN Ia magnitude, and apply a sample selection for SNe Ia to avoid strongly lensed samples. Comparing with the contribution of other uncertainties (e.g., the intrinsic magnitude scatter), we find that the lensing effect is dominant at $z > 1$. Then forecasting the parameter constraints for the Wide-Field InfraRed Survey Telescope survey, we show that considering the weak-lensing effect, the constraints on the density parameters $\Omega_m$ or $\Omega_{\Lambda}$, and the dark energy equation of state $w$ are improved, especially for SNe Ia samples at higher redshift $z > 1$. Furthermore, we see that the degeneracy between the total mass of neutrino $\Sigma m_{\nu}$ and the (cold) dark matter density parameter $\Omega_{c}$ can be resolved.

The Hubble constant problem is the discrepancy between different measurements of the Hubble constant in different scales. We show that this problem can be resolved within the general relativistic framework of the perturbation theory in the inhomogeneous universe, with the help of spatial averaging procedure over a finite local domain in the $t=\mbox{const.}$ hypersurface. The idea presented in this paper is unique in the sense that it has all of the following properties. a) It is based on the general relativistic perturbation theory, with ordinary dust matter only. No strange matter nor energy components are required. b) The employment of the spatially invariant averaging procedure on the finite domain is essential. c) The key is the first-order effect of the inhomogeneities in the linear perturbation theory. No non-linear effects are required.

We investigate the possibility to have a constraint on the mass of thermal warm dark matter (WDM) particle from the expected data of the Wide Field Infrared Survey Telescope (WFIRST) survey if all the dark matter is warm. For this purpose we consider the lensing effect of large scale structure based on the warm dark matter scenario on the apparent magnitude of SNe Ia. We use HALOFIT as non-linear matter power spectrum and the variance of PDF. We preform a Fisher matrix analysis and obtain the lower bound of $m_{\rm WDM}>0.167$keV.

Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constraint on the cosmological parameters. For this purpose, the accurate shape measurement of background galaxies is absolutely important in which any systematic error in the measurement should be carefully corrected. One of the main systematic error comes from photon noise which is Poisson noise of flux from the atmosphere(noise bias). We investigate how the photon noise makes a systematic error in shear measurement within the framework of ERA method we developed in earlier papers and gives a practical correction method. The method is tested by simulations with real galaxy images with various conditions and it is confirmed that it can correct $80 \sim 90\%$ of the noise bias except for galaxies with very low signal to noise ratio.

Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constrain on the cosmological parameters. For this purpose any systematic error in the measurement should be corrected. One of the main systematic error comes from Pixel noise which is Poisson noise of flux from atmosphere. We investigate how the pixel noise makes systematic error in shear measurement based on ERA method and develop the correction method. This method is tested by simulations with various conditions and it is confirmed that the correction method can correct $80 \sim 90\%$ of the systematic error except very low signal to noise ratio galaxies.

We investigate how accurately the total mass of neutrinos is constrained from the magnitude dispersion of SNe Ia due to the effects of gravitational lensing. For this purpose, we use the propagation equation of light bundles in a realistic inhomogeneous universe and propose a sample selection for supernovae to avoid difficulties associated with small-scale effects such as strong lensing or shear effects. With a fitting formula for the nonlinear matter power spectrum taking account of the effects of massive neutrinos, we find that in our model it is possible to obtain the upper limit $\Sigma m_{\nu} \simeq 1.0[{\rm eV}]$ for future optical imaging surveys with the Wide-Field InfraRed Survey Telescope and Large Synoptic Survey Telescope. Furthermore, we discuss how far we need to observe SNe Ia and to what extent we have to reduce the magnitude error except for lensing in order to realize the current tightest limit $\Sigma m_{\nu} < 0.2[{\rm eV}]$.

We improve the ERA(Ellipticity of Re-smeared Artificial image) method of PSF(Point Spread Function) correction in weak lensing shear analysis in order to treat realistic shape of galaxies and PSF. This is done by re-smearing PSF and the observed galaxy image smeared by a RSF(Re-Smearing Function), and allows us to use a new PSF with a simple shape and to correct PSF effect without any approximations and assumptions. We perform numerical test to show that the method applied for galaxies and PSF with some complicated shapes can correct PSF effect with systematic error less than 0.1%. We also apply ERA method for real data of Abell 1689 cluster to confirm that it is able to detect the systematic weak lensing shear pattern. The ERA method requires less than 0.1 or 1 second to correct PSF for each object in numerical test and real data analysis, respectively.

We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.

We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $\beta_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $\beta_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $\beta_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5\sigma$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts.

We generalize ERA method of PSF correction for more realistic situations. The method re-smears the observed galaxy image(galaxy image smeared by PSF) and PSF image by an appropriate function called Re-Smearing Function(RSF) to make new images which have the same ellipticity with the lensed (before smeared by PSF) galaxy image. It has been shown that the method avoids a systematic error arising from an approximation in the usual PSF correction in moment method such as KSB for simple PSF shape. By adopting an idealized PSF we generalize ERA method applicable for arbitrary PSF. This is confirmed with simulated complex PSF shapes. We also consider the effect of pixel noise and found that the effect causes systematic overestimation.

Analysis of strong gravitational lensing data is important in this era of precision cosmology. The objective of the present study is to directly compare the analysis of strong gravitational lens systems using different lens model software and similarly parameterized models to understand the differences and limitations of the resulting models. The software lens model translation tool, HydraLens, was used to generate multiple models for four strong lens systems including COSMOS J095930+023427, SDSS J1320+1644, SDSSJ1430+4105 and J1000+0021. All four lens systems were modeled with PixeLens, Lenstool, glafic, and Lensmodel. The input data and parameterization of each lens model was similar for the four model programs used to highlight differences in the output results. The calculation of the Einstein radius and enclosed mass for each lens model was comparable. The results were more dissimilar if the masses of more than one lens potential were free-parameters. The image tracing algorithms of the software are different, resulting in different output image positions and differences in time delay and magnification calculations, as well as ellipticity and position angle of the resulting lens model. In a comparison of different software versions using identical model input files, results differed significantly when using two versions of the same software. These results further support the need for future lensing studies to include multiple lens models, use of open software, availability of lens model files use in studies, and computer challenges to develop new approaches. Future studies need a standard nomenclature and specification of the software used to allow improved interpretation, reproducibility and transparency of results.

We developed a new method that uses ellipticity defined by 0th order moments (0th-ellipticity) for weak gravitational lensing shear analysis. Although there is a strong correlation between the ellipticity calculated using this approach and the usual ellipticity defined by the 2nd order moment, the ellipticity calculated here has a higher signal-to-noise ratio because it is weighted to the central region of the image. These results were confirmed using data for Abell 1689 from the Subaru telescope. For shear analysis, we adopted the ellipticity of re-smeared artificial image (ERA) method for point spread function (PSF) correction, and we tested the precision of this 0th-ellipticity with simple simulation, then we obtained the same level of precision with the results of ellipticity defined by quadrupole moments. Thus, we can expect that weak lensing analysis using 0 shear will be improved in proportion to the statistical error.

The light rays from a source are subject to a local inhomogeneous geometry generated by inhomogeneous matter distribution as well as the existence of collapsed objects. In this paper we investigate the effect of inhomogeneities and the existence of collapsed objects on the propagation of light rays and evaluate changes in the magnitude-redshift relation from the standard relationship found in a homogeneous FRW universe. We give the expression of the correlation function and the variance for the perturbation of apparent magnitude, and calculate it numerically by using the non-linear matter power spectrum. We use the lognormal probability distribution function for the density contrast and spherical collapse model to truncate the power spectrum in order to estimate the blocking effect by collapsed objects. We find that the uncertainties in $\Omega_m$ is $\sim 0.02$, and that of $w$ is $\sim 0.04$. We also discuss a possible method to extract these effects from real data which contains intrinsic ambiguities associated with the absolute magnitude.

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 propose a new method for Point Spread Function (PSF) correction in weak gravitational lensing shear analysis using an artificial image with the same ellipticity as the lensed image. This avoids the systematic error associated with the approximation in PSF correction used in previous approaches. We test the new method with simulated objects which have Gaussian or Cersic profiles smeared by a Gaussian PSF, and confirm that there is no systematic error.

We present a 4 deg^2 weak gravitational lensing survey of subhalos in the very nearby Coma cluster using the Subaru/Suprime-Cam. The large apparent size of cluster subhalos allows us to measure the mass of 32 subhalos detected in a model-independent manner, down to the order of 10^-3 of the virial mass of the cluster. Weak-lensing mass measurements of these shear-selected subhalos enable us to investigate subhalo properties and the correlation between subhalo masses and galaxy luminosities for the first time. The mean distortion profiles stacked over subhalos show a sharply truncated feature which is well-fitted by a Navarro-Frenk-White (NFW) mass model with the truncation radius, as expected due to tidal destruction by the main cluster. We also found that subhalo masses, truncation radii, and mass-to-light ratios decrease toward the cluster center. The subhalo mass function, dn/dln M_sub, in the range of 2 orders of magnitude in mass, is well described by a single power law or a Schechter function. Best-fit power indices of 1.09_-0.32^+0.42 for the former model and 0.99_-0.23^+0.34 for the latter, are in remarkable agreement with slopes of ~0.9-1.0 predicted by the cold dark matter paradigm. The tangential distortion signals in the radial range of 0.02-2Mpc/h from the cluster center show a complex structure which is well described by a composition of three mass components of subhalos, the NFW mass distribution as a smooth component of the main cluster, and a lensing model from a large scale structure behind the cluster. Although the lensing signals are 1 order of magnitude lower than those for clusters at z~0.2, the total signal-to-noise ratio, S/N=13.3, is comparable to, or higher, because the enormous number of background source galaxies compensates for the low lensing efficiency of the low lensing efficiency of the nearby cluster.

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.

Analysis of strong gravitational lensing depends on software analysis of observational data. The purpose of this study was to evaluate the behavior of strong gravitational lens modeling software with changes in redshift. Four different strong gravitational lens software modeling codes were directly compared (Lenstool / glafic, two light traces mass codes, and GRALE / PixeLens, two non-light traces mass codes) in the analysis of model data as well as analysis of the giant gravitational quasar SDSSJ1004+4112. A generalized model for time delay calculation shows that calculated time delay is proportional to $D_{d}D_{s}/D_{ds}$.The percent change in time delays calculated for each system at each redshift tested were compared with percent change in the value of $D_{d}D_{s}/D_{ds}$. A simple point mass model was tested with each code. Five models were used with a constant $z_{lens}$ and a varying $z_{source}$, and five models with a constant $z_{source}$ and a varying $z_{lens}$. The effects of changing geometry were similarly investigated for SDSSJ1004+4112. In general, the changes in time delay were of a similar magnitude and direction although, some calculated time delays varied by as much as 30 percent from changes in $D_{d}D_{s}/D_{ds}$. Changes in calculated mass for the point mass model with a constant $z_{source}$ were almost identical to changes in $D_{d}D_{s}/D_{ds}$ for three of the four codes tested. These data demonstrate the effect of changes in redshift on parameters calculated by each of the codes as compared to changes in $D_{d}D_{s}/D_{ds}$. The paucity of existing direct comparison studies of strong gravitational lensing supports the need for more studies of this kind. These results show that even small changes in redshift affects the calculation of time delay and mass, and that the effect on the calculations is dependent on the particular software used.

Despite expanding research activity in gravitational lens modeling, there is no particular software which is considered a standard. Much of the gravitational lens modeling software is written by individual investigators for their own use. Some gravitational lens modeling software is freely available for download but is widely variable with regard to ease of use and quality of documentation. This review of 13 software packages was undertaken to provide a single source of information. Gravitational lens models are classified as parametric models or non-parametric models, and can be further divided into research and educational software. Software used in research includes the GRAVLENS package (with both gravlens and lensmodel), Lenstool, LensPerfect, glafic, PixeLens, SimpLens, Lensview, and GRALE. In this review, GravLensHD, G-Lens, Gravitational Lensing, lens and MOWGLI are categorized as educational programs that are useful for demonstrating various aspects of lensing. Each of the 13 software packages is reviewed with regard to software features (installation, documentation, files provided, etc.) and lensing features (type of model, input data, output data, etc.) as well as a brief review of studies where they have been used. Recent studies have demonstrated the utility of strong gravitational lensing data for mass mapping, and suggest increased use of these techniques in the future. Coupled with the advent of greatly improved imaging, new approaches to modeling of strong gravitational lens systems are needed. This is the first systematic review of strong gravitational lens modeling software, providing investigators with a starting point for future software development to further advance gravitational lens modeling research.

We investigate the relation between the regularized multi-propagator method, called "Reg PT", and the standard perturbation theory. Reg PT is one of the most successful models to describe nonlinear evolution of dark matter fluctuations. However, Reg PT is a mathematically unproven interpolation formula between the large-scale solution calculated by the standard perturbation theory and the limiting solution in the small scale calculated by the multi-point propagator method. In this paper, we give an alternative explanation for Reg PT in the context of the standard perturbation theory, showing that Reg PT does not ever have more effective information on nonlinear matter evolution than the standard perturbation theory. In other words, the solutions of the standard perturbation theory reproduce the results of $N$-body simulations better than those of Reg PT, especially at the high-$k$ region. This fact means that the standard perturbation theory at the two-loop level is still one of the best predictions of the nonlinear power spectrum to date. Nevertheless, the standard perturbation theory has not been preferred because of the divergent behavior of the solution at small scales. To solve this problem, we also propose a modified standard perturbation theory which avoids the divergence.

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

Precise understanding of nonlinear evolution of cosmological perturbations during inflation is necessary for the correct interpretation of measurements of non-Gaussian correlations in the cosmic microwave background and the large-scale structure of the universe. The "{\delta}N formalism" is a popular and powerful technique for computing non-linear evolution of cosmological perturbations on large scales. In particular, it enables us to compute the curvature perturbation, {\zeta}, on large scales without actually solving perturbed field equations. However, people often wonder why this is the case. In order for this approach to be valid, the perturbed Hamiltonian constraint and matter-field equations on large scales must, with a suitable choice of coordinates, take on the same forms as the corresponding unperturbed equations. We find that this is possible when (1) the unperturbed metric is given by a homogeneous and isotropic Friedmann-Lema\^itre-Robertson-Walker metric; and (2) on large scales and with a suitable choice of coordinates, one can ignore the shift vector (g0i) as well as time-dependence of tensor perturbations to gij/a2(t) of the perturbed metric. While the first condition has to be assumed a priori, the second condition can be met when (3) the anisotropic stress becomes negligible on large scales. However, in order to explicitly show that the second condition follows from the third condition, one has to use gravitational field equations, and thus this statement may depend on the details of theory of gravitation. Finally, as the {\delta}N formalism uses only the Hamiltonian constraint and matter-field equations, it does not a priori respect the momentum constraint. We show that the violation of the momentum constraint only yields a decaying mode solution for {\zeta}, and the violation vanishes when the slow-roll conditions are satisfied.

We propose a new approximated expression for non-linear Dark Matter power spectrum much beyond BAO scales. The proposed expression agrees with the result of N-body simulation with the accuracy better than 2 % up to k=1.0 [h/Mpc] and k=0.7 [h/Mpc] at z=3.0 and z=1.0, respectively. Even at z=0.35, the accuracy remains within 10 % up to k=0.8 [h/Mpc]. In doing so, we used an approximation for the kernel functions used in the Standard Perturbation Theory (SPT) which is also used to prove the Reg PT proposed by Bernardeau et al. (2011).

We apply the Wiener Hermite (WH) expansion to the non-linear evolution of Large-Scale Structure, and obtain an approximate expression for the matter power spectrum in full order of the expansion. This method allows us to expand any random function in terms of an orthonormal bases in space of random functions in such a way that the first order of the expansion expresses the Gaussian distribution, and others are the deviation from the Gaussianity. It is proved that the Wiener Hermite expansion is mathematically equivalent with the $\Gamma$-expansion approach in the Renormalized Perturbation Theory (RPT). While an exponential behavior in the high-$k$ limit has been proved for the mass density and velocity fluctuations of Dark Matter in the RPT, we reprove the behavior in the context of the Wiener Hermite expansion using the result of Standard Perturbation Theory (SPT). We propose a new approximate expression for the matter power spectrum which interpolates the low-$k$ expression corresponding to 1-loop level in SPT and the high-$k$ expression obtained by taking a high-$k$ limit of the WH expansion. The validity of our prescription is specifically verified by comparing with the 2-loop solutions of the SPT. The proposed power spectrum agrees with the result of $N$-body simulation with the accuracy better than 1% or 2% in a range of the BAO scales, where the wave number is about $k= 0.2 \sim 0.4$ $h{\rm Mpc^{-1}}$ at $z=0.5-3.0$. This accuracy is comparable to or slightly less than the ones in the Closure Theory whose the fractional difference from N body result is within 1%. A merit of our method is that the computational time is very short because only single and double integrals are involved in our solution.

This is the third paper on the improvements of systematic errors in our weak lensing analysis using an elliptical weight function, called E-HOLICs. In the previous papers we have succeeded in avoiding error which depends on ellipticity of background image. In this paper, we investigate the systematic error which depends on signal to noise ratio of background image. We find that the origin of the error is the random count noise which comes from Poisson noise of sky counts. Random count noise makes additional moments and centroid shift error, and those 1st orders are canceled in averaging, but 2nd orders are not canceled. We derived the equations which corrects these effects in measuring moments and ellipticity of the image and test their validity using simulation image. We find that the systematic error becomes less than 1% in the measured ellipticity for objects with $S/N>3$.

We have developed a new method(E-HOLICs) of estimating gravitational shear by adopting an elliptical weight function to measure background galaxy images in our previous paper. Following to the previous paper where isotropic Point Spread Function(PSF) correction is calculated, we consider an anisotropic PSF correction in this paper in order to apply E-HOLICs for real data. A an example, E-HOLICs is applied to Subaru data of massive and compact galaxy clusters A370, and is able to detect double peaks in the central region of the cluster consistent with the analysis of strong lensing. We also study the systematic error in E-HOLICs using STEP II simulation. In particular we consider the dependences of signal to noise ratio "S/N" of background galaxies in the shear estimation. Although E-HOLICs does improve systematic error due to the ellipticity dependence as shown in part 1, a systematic error due to the S/N dependence remains, namely E-HOLICs underestimates shear when background galaxies with low S/N objects are used. We discuss a possible improvement of the S/N dependence.

While detection of the "local form" bispectrum of primordial perturbations would rule out all single-field inflation models, multi-field models would still be allowed. We show that multi-field models described by the $\delta N$ formalism obey an inequality between $f_{\rm NL}$ and one of the local-form {\it trispectrum} amplitudes, $\tau_{\rm NL}$, such that $\tau_{\rm NL}>\frac12(\frac65f_{\rm NL})^2$ with a possible logarithmic scale dependence, provided that 2-loop terms are small. Detection of a violation of this inequality would rule out most of multi-field models, challenging inflation as a mechanism for generating the primoridal perturbations.

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

We develop a new method to estimate gravitational shear by adopting an elliptical weight function to measure background galaxy images. In doing so, we introduce a new concept of "zero plane" which is an imaginal source plane where shapes of all sources are perfect circles, and regard the intrinsic shear as the result of an imaginal lensing distortion. This makes the relation between the observed shear, the intrinsic shear and lensing distortion more simple and thus higher-order calculation more easy. The elliptical weight function allows us to measure the mutiplemoment of shape of background galaxies more precisely by weighting highly to brighter parts of image and moreover to reduce systematic error due to insufficient expansion of the weight function in the original approach of KSB. Point Spread Function(PSF) correction in E-HOLICs methods becomes more complicated than those in KSB methods. In this paper we studied isotropic PSF correction in detail. By adopting the lensing distortion as the ellipticity of the weight function, we are able to show that the shear estimation in E-HOLICs method reduces to solve a polynomial in the absolute magnitude of the distortion. We compare the systematic errors between our approach and KSB using STEP2 simulation. It is confirmed that KSB method overestimate the input shear for images with large ellipticities, and E-HOLICs correctly estimate the input shear even for such images. Anisotropic PSF correction and analysis of real data will be presented in forthcoming paper.

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

We obtain the projected mass distributions for two Subaru/Suprime-Cam fields in the southwest region (r\simlt 60') of the Coma cluster (z=0.0236) by weak lensing analysis and detect eight subclump candidates. We quantify the contribution of background large-scale structure (LSS) on the projected mass distributions using SDSS multi-bands and photometric data, under the assumption of mass-to-light ratio for field galaxies. We find that one of eight subclump candidates, which is not associated with any member galaxies, is significantly affected by LSS lensing. The mean projected mass for seven subclumps extracted from the main cluster potential is <M_2D^(corr)> = (5.06\pm1.30)10^12h^-1 M_sun after a LSS correction. A tangential distortion profile over an ensemble of subclumps is well described by a truncated singular-isothermal sphere model and a truncated NFW model. A typical truncated radius of subclumps, r_t\simeq 35 h^-1 kpc, is derived without assuming any relations between mass and light for member galaxies. The radius coincides well with the tidal radius, \sim42 h^-1 kpc, of the gravitational force of the main cluster. Taking into account the incompleteness of data area, a projection effect and spurious lensing peaks, it is expected that mass of cluster substructures account for 19 percent of the virial mass, with 13 percent statistical error. The mass fraction of cluster substructures is in rough agreement with numerical simulations.

We investigate the distance-redshift relation in a realistic inhomogeneous universe where the mass distribution is described by the mass function of Sheth and Tormen. It is found that the derived distance deviates systematically from the standard distance up to 10% depending on the choice of the lowest halo mass in which baryonic matter condensed to form luminous object such as galaxies. Remarkably the derived distance is well approximated by the Dyer-Roeder distance if we choose the clumpiness parameter \alpha calculated by our model.We also discuss the effect of inhomogeneities in the determination of dark energy parameter in the supernovae observation, and find that this effect must be taken into account for the future high redshift supernovae observation.

We investigate the possibility to use higher order moments of gravitational lensed images in the weak lensing analysis. For this purpose we employ spin-2 component of HOLICs(Higher Order Lensing Characteristics) developed by us. We test the weak lensing analysis with spin-2 HOLICs using actual, ground based Subaru observations of the massive galaxy cluster A1689(z=0.183). It turns out that spin-2 HOLICs of order up to 8 are sufficiently applicable for weak lensing analysis after correcting PSF anisotropy as well as isotropic PSF smearing.

We investigate the weak gravitational lensing effect on the two-point correlation function of local maxima (hot spots) in the cosmic 21 cm fluctuation map. The intrinsic two-point function has a pronounced depression feature around the angular scale of $\theta \sim 40$ arcmin, which depends on the observed frequency and corresponds to the scale of the acoustic oscillation of cosmic plasma before the recombination. It is found that the weak lensing induces a large w-dependent smoothing at that scale where w is the equation of state parameter of dark energy and thus provides a useful constraints on the dark energy property combined with the depression angular scales on the two-point correlation function.

We have developed a method for measuring higher-order weak lensing distortions of faint background galaxies, namely the weak gravitational flexion, by fully extending the Kaiser, Squires & Broadhurst method to include higher-order lensing image characteristics (HOLICs) introduced by Okura, Umetsu, & Futamase. We take into account explicitly the weight function in calculations of noisy shape moments and the effect of higher-order PSF anisotropy, as well as isotropic PSF smearing. Our HOLICs formalism allows accurate measurements of flexion from practical observational data in the presence of non-circular, anisotropic PSF. We test our method using mock observations of simulated galaxy images and actual, ground-based Subaru observations of the massive galaxy cluster A1689 ($z=0.183$). From the high-precision measurements of spin-1 first flexion, we obtain a high-resolution mass map in the central region of A1689. The reconstructed mass map shows a bimodal feature in the central $4'\times 4'$ region of the cluster. The major, pronounced peak is associated with the brightest cluster galaxy and central cluster members, while the secondary mass peak is associated with a local concentration of bright galaxies. The refined, high-resolution mass map of A1689 demonstrates the power of the generalized weak lensing analysis techniques for quantitative and accurate measurements of the weak gravitational lensing signal.

We examine the usefulness of the strong lensing statistics to constrain the evolution of the number density of lensing galaxies by adopting the values of the cosmological parameters determined by recent WMAP observation. For this purpose, we employ the lens-redshift test proposed by Kochanek (1992) and constrain the parameters in two evolution models, simple power-law model characterized by the power law indexes $nu_{n}$ and $\nu_{v}$ and the evolution model by Mitchell et al. (2005) based on CDM structure formation scenario. We use the well-defined lens sample from the Sloan Digital Sky Survey (SDSS) and this is similarly sized samples used in the previous studies. Furthermore, we adopt the velocity dispersion function of early-type galaxies based on SDSS DR1 and DR5. It turns out that the indexes of power-law model are consistent with the previous studies, thus our results indicate the mild evolution in the number and velocity dispersion of early-type galaxies out to z = 1. However we found that the values for p and q used by Mitchell et al. are inconsistent with the presently available observational data. More complete sample is necessary to withdraw more realistic determination on these parameters.

We study a possibility to use the octopole moment of gravitationally lensed images as a direct measure of the third-order weak gravitational lensing effect, or the gravitational flexion. It turns out that there is a natural relation between flexion and certain combinations of octopole/higher-multipole moments which we call the Higher Order Lensing Image's Characteristics (HOLICs). This will allow one to measure directly flexion from observable octopole and higher-multipole moments of background images. We show based on simulated observations how the use of HOLICs can improve the accuracy and resolution of a reconstructed mass map, in which we assume Gaussian uncertainties in the shape measurements estimated using deep i'-band data of blank fields observed with Suprime-Cam on the Subaru telescope.

The backreaction of nonlinear inhomogeneities to the cosmic expansion is analyzed in the framework of general relativity with a cosmological constant. By defining the spatially averaged matter energy density, we find that the cosmological constant induces a new type of backreaction whose equation of state is $p=-4/3 \rho$, where $\rho$ and $p$ are the effective energy density and effective pressure of the backreaction in the averaged Friedmann universe. However, the effective density is negative, and thus it decreases the acceleration caused by the cosmological constant.

We give a straightforward and divergence free derivation of the equation of motion for a small but finite object in an arbitrary background using strong field point particle limit. The resulting equation becomes a generalized geodesic for a non-rotating spherical object which is consisitent with previous studies.

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

The backreaction of nonlinear inhomogeneities to the cosmic expansion is re-analyzed in the framework of general relativity. Apparent discrepancies regarding the effect of the nonlinear backreaction, which exist among the results of previous works in different gauges, are resolved. By defining the spatially averaged matter energy density as a conserved quantity in the large comoving volume, it is shown that the nonlinear backreaction neither accelerates nor decelerates the cosmic expansion in a matter-dominated universe. The present result in the Newtonian gauge is consistent with the previous results obtained in the comoving synchronous gauge. Although our work does not give a complete proof, it strongly suggests the following no-go theorem: No cosmic acceleration occurs as a result of the nonlinear backreaction via averaging.

We present a weak-lensing analysis of a region around the galaxy cluster Cl 1604+4304 (z=0.897) on the basis of the deep observations with the HST/WFPC2. We apply a variant of Schneider's aperture mass technique to the observed WFPC2 field and obtain the distribution of weak-lensing signal-to-noise (S/N) ratio within the field. The resulting S/N map reveals a clear pronounced peak located about 1.7 arcmin (850h_{50}^{-1} kpc at z=0.897) southwest of the second peak associated with the optical cluster center determined from the dynamical analysis of Postman et al. A non-linear finite-field inversion method has been used to reconstruct the projected mass distribution from the observed shear field. The reconstructed mass map shows a super-critical feature at the location of the S/N peak as well as in the cluster central region. Assuming the redshift distribution of field galaxies, we obtain the total mass in the observed field to be 1.0 h_{50}^{-1} 10^{15} M_sun for <z>=1.0. The estimated mass within a circular aperture of radius 280h_{50}^{-1} kpc centered on the dark clump is 2.4h_{50}^{-1} 10^{14} M_sun. We have confirmed the existence of the ` dark ' mass concentration from another deep HST observation with a slightly different ~20 arcsec pointing.

We report the discovery and spectroscopic confirmation of three gravitationally-lensed images of a galaxy at z=3.9 in the background of a distant, rich cluster of galaxies at z=0.83, on the basis of observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estimates from deep, high-resolution images by Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The lens model can account simultaneously for the observed image configuration and the flux ratio of the closer pair located close to the critical curve. The parities of the three images are also consistent with the lensing hypothesis. Since this galaxy is apparently bright (i'_{AB} \sim 23.7) for its redshift due to the magnification, it serves as a good high redshift target on which we can make extensive and detailed studies based on multi-wavelength imaging and spectroscopy.

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

We report the first measurement of genus curves for the two-dimensional mass map in the neighborhood of rich, X-ray luminous galaxy cluster MS1054-03 at z=0.83, reconstructed from weak lensing data obtained by Suprime-Cam on the prime focus of 8.2m Subaru telescope. We find that the genus curve measured in the whole survey field deviates from that expected from a random Gaussian field. We show that this non-Gaussianity is induced by the rich cluster in this region, and that the genus curve for the region without the cluster is consistent with the prediction for a random Gaussian field. The results demonstrate the feasibility of the genus statistics to examine the non-Gaussianity due to the large scale structures and to probe the statistical properties of the large scale clustering.

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.

We assess the feasibility of probing dark energy with strong gravitational lensing systems. The capability of the method, which depends on the accuracy with which the lensing systems are modeled, is quantitatively investigated using the Fisher-matrix formalism. We show that this method might place useful constraints on the density parameter and the redshift evolution of the dark energy by combining it with a constraint from supernova measurements. For this purpose, the lens potential needs to be precisely reconstructed. We determine the required quality of data. We also briefly discuss the optimal strategy to constrain the cosmological parameters using gravitational lensing systems.

We investigate properties of $r$-mode instability in slowly rotating relativistic polytropes. Inside the star slow rotation and low frequency formalism that was mainly developed by Kojima is employed to study axial oscillations restored by Coriolis force. At the stellar surface, in order to take account of gravitational radiation reaction effect, we use a near-zone boundary condition instead of the usually imposed boundary condition for asymptotically flat spacetime. Due to the boundary condition, complex frequencies whose imaginary part represents secular instability are obtained for discrete $r$-mode oscillations in some polytropic models. It is found that such discrete $r$-mode solutions can be obtained only for some restricted polytropic models. Basic properties of the solutions are similar to those obtained by imposing the boundary condition for asymptotically flat spacetime. Our results suggest that existence of a continuous part of spectrum cannot be avoided even when its frequency becomes complex due to the emission of gravitational radiation.

The redshift dependence of the abundance of collapsed objects places strong constraints on cosmological models of structure formation. We apply a recently proposed model describing the anisotropic collapse of inhomogeneous spatial domains. Compared with the spherical top-hat model, this generic model leads to significantly more collapsed objects at high redshifts: at redshift one and on the scale of rich clusters a factor of 65. Furthermore, for a fixed normalization of the initial fluctuation spectrum (\sigma_8=1), we predict four times as much presently collapsed objects on the mass-scale of rich clusters within the standard CDM cosmogony, compared to the spherical collapse.

A possible method to reconstruct the cosmic equation of state using strong gravitational lensing systems is proposed. The feasibility of the method is investigated by carrying out the reconstruction on the basis of a simple Monte-Carlo simulation. We show that the method can work and that the cosmic equation of state $w(z)$ can be determined within errors of $\Delta w\sim \pm0.1$ -- $\pm0.2$ when a sufficiently large number of lensing systems ($N\sim 20$) for $z\simlt 1$ are precisely measured. Statistics of lensed sources in a wide and deep survey like the SDSS are also briefly discussed.

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

We propose a possible measurement of the time variability of the vacuum energy using strong gravitational lensing events. As an example we take an Einstein cross lens HST 14176+5226 and demonstrate that the measurement of the velocity dispersion with the accuracy of $\pm$ 5 km/sec will have a chance to determine the time dependence of the vacuum energy as well as the density parameter with the accuracy of order 0.1 if one fixes the lens model.

Using a multiple-lens plane algorithm, we study light propagation in inhomogeneous universes, for 43 different COBE-normalized Cold Dark Matter models, with various values of the density parameter Omega_0, cosmological constant lambda_0, Hubble constant H_0, and rms density fluctuation sigma_8. This is the largest cosmological parameter survey ever done in this field. We performed a total of 3,798 experiments, each experiment consisting of propagating a square beam of angular size 21.9'' x 21.9'' composed of 116,281 light rays from the observer up to redshift z=3. These experiments provide statistics of the magnification, shear, and multiple imaging of distant sources. The results of these experiments can be compared with observations, and eventually help constraining the possible values of the cosmological parameters. Additionally, they provide insight into the gravitational lensing process and its complex relationship with the various cosmological parameters.

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

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

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

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

The tidal gravitational field of galaxy clusters causes a coherent distortion of the images of background sources. Since the distribution of image distortions, namely the shear field, traces the local gravitational potential of a deflector, it can be used to reconstruct the two-dimensional mass distribution of clusters of galaxies. Moreover, the shear field can provide unique information on the redshift distribution of high-redshift galaxies. In this review we summarize recently-developed parameter-free methods of cluster-mass reconstruction based on the shear field, and we apply a mass-reconstruction method to the cluster Abell 370 at redshift 0.375.

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.

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.

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

We present complete constraints imposed from observations of the cosmic microwave background radiation (CMBR) on the chaotic inflationary scenario with a nonminimally coupled inflaton field proposed by Fakir and Unruh (FU). Our constraints are complete in the sense that we investigate both the scalar density perturbation and the tensor gravitational wave in the Jordan frame, as well as in the Einstein frame. This makes the constraints extremely strong without any ambiguities due to the choice of frames. We find that the FU scenario generates tiny tensor contributions to the CMBR relative to chaotic models in minimal coupling theory, in spite of its spectral index of scalar perturbation being slightly tilted. This means that the FU scenario will be excluded if any tensor contributions to CMBR are detected by the forthcoming satellite missions. Conversely, if no tensor nature is detected despite the tilted spectrum, a minimal chaotic scenario will be hard to explain and the FU scenario will be supported.

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

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

We investigate a possibility to restrict the chaotic inflationary scenario with "large" nonminimally coupled inflaton("phi") field considered by Fakir and Unruh by means of the observed CMBR anisotropy. We calculate the contributions of the long wavelength gravitational waves generated in the inflationary period to CMBR anisotropy quadrupole moment. We obtain the constraint lambda/(xi^2) = 1.8 X 10^(-8)/psi_i, where "lambda" is the self coupling, "xi" is the nonminimal coupling constant and "psi_i" means the initial value of psi = 8piG(phi^2). Combining this with previously obtained constraint that 2.1 times squre root of lambda/(xi^2) = dT/T = 1.1 X 10^(-5), we conclude that the initial value has to be psi_i = 1.6 X 10^3. If the self-coupling has a reasonably values of order 10^(-2), then xi = 10^4 and phi_i = 10^(-1) times Planck mass.

We have found a successful model of chaotic inflation with an inflaton coupled nonminimally with gravity. The nonminimal coupling constant $\xi$ runs with the evolution of the inflaton. The running nature of the coupling leads naturally to the situations where the coupling becomes small enough to have sufficient period of the inflation to resolve the cosmological puzzles.

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.

Using the (3+1) formalism in general relativity, we perform the post-Newtonian(PN) approximation to clarify what sort of gauge condition is suitable for numerical analysis of coalescing compact binary neutron stars and gravitational waves from them. We adopt a kind of transverse gauge condition to determine the shift vector. On the other hand, for determination of the time slice, we adopt three slice conditions(conformal slice, maximal slice and harmonic slice) and discuss their properties. Using these conditions, the PN hydrodynamic equations are obtained up through the 2.5PN order including the quadrupole gravitational radiation reaction. In particular, we describe methods to solve the 2PN tensor potential which arises from the spatial 3-metric. It is found that the conformal slice seems appropriate for analysis of gravitational waves in the wave zone and the maximal slice will be useful for describing the equilibrium configurations. The PN approximation in the (3+1) formalism will be also useful to perform numerical simulations using various slice conditions and, as a result, to provide an initial data for the final merging phase of coalescing binary neutron stars which can be treated only by fully general relativistic simulations.