The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

We present spectroscopic confirmation of three new two-image gravitationally lensed quasars, compiled from existing strong lens and X-ray catalogs. Images of HSC J091843.27-022007.5 show a red galaxy with two blue point sources at either side, separated by 2.26 arcsec. This system has a source and a lens redshifts $z_s = 0.804$ and $z_l = 0.459$, respectively, as obtained by our follow-up spectroscopic data. CXCO J100201.50+020330.0 shows two point sources separated by 0.85 arcsec on either side of an early-type galaxy. The follow-up spectroscopic data confirm the fainter quasar has the same redshift with the brighter quasar from the SDSS fiber spectrum at $z_s = 2.016$. The deflecting foreground galaxy is a typical early-type galaxy at a redshift of $z_l = 0.439$. SDSS J135944.21+012809.8 has two point sources with quasar spectra at the same redshift $z_s = 1.096$, separated by 1.05 arcsec, and fits to the HSC images confirm the presence of a galaxy between these. These discoveries demonstrate the power of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)'s deep imaging and wide sky coverage. Combined with existing X-ray source catalogues and follow-up spectroscopy, the HSC-SSP provides us unique opportunities to find multiple-image quasars lensed by a foreground galaxy.

We provide a new fitting formula of the matter bispectrum in the nonlinear regime calibrated by high-resolution cosmological $N$-body simulations of $41$ cold dark matter ($w$CDM, $w=$ constant) models around the Planck 2015 best-fit parameters. As the parameterization in our fitting function is similar to that in Halofit, our fitting is named BiHalofit. The simulation volume is sufficiently large ($> 10 \, {\rm Gpc}^3$) to cover almost all measurable triangle bispectrum configurations in the universe. The function is also calibrated using one-loop perturbation theory at large scales ($k<0.3 \, h \, {\rm Mpc}^{-1}$). Our formula reproduced the matter bispectrum to within $10 \, (15) \, \%$ accuracy in the Planck 2015 model at wavenumber $k< 3 \, (10) \, h \, {\rm Mpc}^{-1}$ and redshifts $z=0-3$. The other $40$ $w$CDM models obtained poorer fits, with accuracy approximating $20 \, \%$ at $k<3 \, h \, {\rm Mpc}^{-1}$ and $z=0-1.5$ (the deviation includes the $10 \, \%$-level sample variance of the simulations). We also provide a fitting formula that corrects the baryonic effects such as radiative cooling and active galactic nucleus feedback, using the latest hydrodynamical simulation IllustrisTNG. We demonstrate that our new formula more accurately predicts the weak-lensing bispectrum than the existing fitting formulas. This formula will assist current and future weak-lensing surveys and cosmic microwave background lensing experiments. Numerical codes of the formula are available, written in Python, C and Fortran.

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 report the serendipitous discovery of HSC J0904$-$0102, a quadruply-lensed Lyman break galaxy (LBG) in the Survey of Gravitationally-lensed Objects in Hyper Suprime-Cam Imaging (SuGOHI). Owing to its point-like appearance, the source was thought to be a lensed active galactic nucleus. We obtained follow-up spectroscopic data with the Gemini Multi-Object Spectrographs on the Gemini South Telescope, which confirmed this to be a lens system. The deflecting foreground galaxy is a typical early-type galaxy at a high redshift of $z_{\ell} = 0.957$ with stellar velocity dispersion $\sigma_v=259\pm56$ km~s$^{-1}$. The lensed source is identified as an LBG at $z_{\rm s} = 3.403$, based on the sharp drop bluewards of Ly$\alpha$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $\theta_{\rm Ein} = 1. 23^{\prime\prime}$ and a total mass within the Einstein radius of $M_{\rm Ein} = (5.55\pm 0.24) \times 10^{11}M_{\odot}$ corresponding to a velocity dispersion of $\sigma_{\rm SIE}= 283\pm 3$ km~s$^{-1}$, which is in good agreement with the value derived spectroscopically. The most isolated lensed LBG image has a magnification of $\sim 6.5$. In comparison with other lensed LBGs and typical $z\sim4$ LBG populations, HSC J0904$-$0102 is unusually compact, an outlier at $>2\sigma$ confidence. Together with a previously discovered SuGOHI lens, HSC J1152$+$0047, that is similarly compact, we believe that the HSC Survey is extending LBG studies down to smaller galaxy sizes.

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

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 results of an imaging observation campaign conducted with the Subaru Telescope adaptive optics system (IRCS+AO188) on 28 gravitationally lensed quasars (23 doubles, 1 quad, and 1 possible triple, and 3 candidates) from the SDSS Quasar Lens Search. We develop a novel modelling technique that fits analytical and hybrid point spread functions (PSFs), while simultaneously measuring the relative astrometry, photometry, as well as the lens galaxy morphology. We account for systematics by simulating the observed systems using separately observed PSF stars. The measured relative astrometry is comparable with that typically achieved with the Hubble Space Telescope, even after marginalizing over the PSF uncertainty. We model for the first time the quasar host galaxies in 5 systems, without a-priory knowledge of the PSF, and show that their luminosities follow the known correlation with the mass of the supermassive black hole. For each system, we obtain mass models far more accurate than those previously published from low-resolution data, and we show that in our sample of lensing galaxies the observed light profile is more elliptical than the mass, for ellipticity > 0.25. We also identify eight doubles for which the sources of external and internal shear are more reliably separated, and should therefore be prioritized in monitoring campaigns aimed at measuring time-delays in order to infer the Hubble constant.

We report the discovery of 13 confirmed two-image quasar lenses from a systematic search for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search (SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular separations ($\lesssim$2") confirmed from our spectroscopy, which are either projected pairs, physical binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. The newly discovered quasar lens system, SDSS J1452+4224 at zs$\approx$4.8 is one of the highest redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the lens sample presented here does not have a well-defined selection function.

The most widely used Galactic extinction map (Schlegel, Finkbeiner, & Davis 1998, SFD) is constructed assuming that the observed FIR fluxes entirely come from the Galactic dust. According to the earlier suggestion by Yahata et al. (2007), we consider how far-infrared (FIR) emission of galaxies affects the SFD map. We first compute the surface number density of SDSS DR7 galaxies as a function of the $r$-band extinction, $A_{r,\rm SFD}$. We confirm that the surface densities of those galaxies positively correlate with $A_{r,\rm SFD}$ for $A_{r,\rm SFD}<0.1$, as first discovered by Yahata et al. (2007) for SDSS DR4 galaxies. Next we construct an analytic model to compute the surface density of galaxies taking account of the contamination of their FIR emission. We adopt a log-normal probability distribution for the ratio of $100\mu {\rm m}$ and $r$-band luminosities of each galaxy, $y \equiv (\nu L)_{100\mu {\rm m}}/(\nu L)_r$. Then we search for the mean and r.m.s values of $y$ that fit the observed anomaly using the analytic model. The required values to reproduce the anomaly are roughly consistent with those measured from the stacking analysis of SDSS galaxies (Kashiwagi, Yahata, & Suto 2013). Due to the limitation of our statistical modeling, we are not yet able to remove the FIR contamination of galaxies from the extinction map. Nevertheless the agreement with the model prediction suggests that the FIR emission of galaxies is mainly responsible for the observed anomaly. While the corresponding systematic error in the Galactic extinction map is 0.1 to 1mmag, it is directly correlated with galaxy clustering, and thus needs to be carefully examined in precision cosmology.

Upcoming wide-area weak lensing surveys are expensive both in time and cost and require an optimal survey design in order to attain maximum scientific returns from a fixed amount of available telescope time. The super-sample covariance (SSC), which arises from unobservable modes that are larger than the survey size, significantly degrades the statistical precision of weak lensing power spectrum measurement even for a wide-area survey. Using the 1000 mock realizations of the log-normal model, which approximates the weak lensing field for a $\Lambda$-dominated cold dark matter model, we study an optimal survey geometry to minimize the impact of SSC contamination. For a continuous survey geometry with a fixed survey area, a more elongated geometry such as a rectangular shape of 1:400 side-length ratio reduces the SSC effect and allows for a factor 2 improvement in the cumulative signal-to-noise ratio ($S/N$) of power spectrum measurement up to $\ell_{\rm max}\simeq $ a few $10^3$, compared to compact geometries such as squares or circles. When we allow the survey geometry to be disconnected but with a fixed total area, assuming $1\times 1$ sq. degrees patches as the fundamental building blocks of survey footprints, the best strategy is to locate the patches with $\sim 15$ degrees separation. This separation angle corresponds to the scale at which the two-point correlation function has a negative minimum. The best configuration allows for a factor 100 gain in the effective area coverage as well as a factor 2.5 improvement in the $S/N$ at high multipoles, yielding a much wider coverage of multipoles than in the compact geometry.

We report the discovery of four doubly imaged quasar lenses. All the four systems are selected as lensed quasar candidates from the Sloan Digital Sky Survey data. We confirm their lensing hypothesis with additional imaging and spectroscopic follow-up observations. The discovered lenses are SDSS J0743+2457 with the source redshift z_s=2.165, the lens redshift z_l=0.381, and the image separation theta=1.034", SDSS J1128+2402 with z_s=1.608 and theta=0.844", SDSS J1405+0959 with z_s=1.810, z_l~0.66, and theta=1.978", and SDSS J1515+1511 with z_s=2.054, z_l=0.742, and theta=1.989". It is difficult to estimate the lens redshift of SDSS J1128+2402 from the current data. Two of the four systems (SDSS J1405+0959 and SDSS J1515+1511) are included in our final statistical lens sample to derive constraints on dark energy and the evolution of massive galaxies.

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

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

We report the discovery of a pair of quasars at $z=1.487$, with a separation of $8\farcs585\pm0\farcs002$. Subaru Telescope infrared imaging reveals the presence of an elliptical and a disk-like galaxy located almost symmetrically between the quasars, creating a cross-like configuration. Based on absorption lines in the quasar spectra and the colors of the galaxies, we estimate that both galaxies are located at redshift $z=0.899$. This, as well as the similarity of the quasar spectra, suggests that the system is a single quasar multiply imaged by a galaxy group or cluster acting as a gravitational lens, although the possibility of a binary quasar cannot be fully excluded. We show that the gravitational lensing hypothesis implies these galaxies are not isolated, but must be embedded in a dark matter halo of virial mass $\sim 4 \times 10^{14}\ h_{70}^{-1}\ {M}_\odot$ assuming an NFW model with a concentration parameter of $c_{vir}=6$, or a singular isothermal sphere profile with a velocity dispersion of $\sim 670$ km s$^{-1}$. We place constraints on the location of the dark matter halo, as well as the velocity dispersions of the galaxies. In addition, we discuss the influence of differential reddening, microlensing and intrinsic variability on the quasar spectra and broadband photometry.

The cosmological power spectrum of the coherent matter flow is measured exploiting an improved prescription for the apparent anisotropic clustering pattern in redshift space. New statistical analysis is presented to provide an optimal observational platform to link the improved redshift distortion theoretical model to future real datasets. The statistical power as well as robustness of our method are tested against 60 realizations of 8 Gpc/h^3 dark matter simulation maps mocking the precision level of upcoming wide--deep surveys. We showed that we can accurately extract the velocity power spectrum up to quasi linear scales of k~0.1 h/Mpc at z = 0.35 and up to k~0.15 h/Mpc at higher redshifts within a couple of percentage precision level. Our understanding of redshift space distortion is proved to be appropriate for precision cosmology, and our statistical method will guide us to righteous path to meet the real world.

We measure the small-scale (comoving separation 10 kpc/h < r_p < 200 kpc/h) two-point correlation function of quasars using a sample of 26 spectroscopically confirmed binary quasars at 0.6<z<2.2 from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). Thanks to careful candidate selections and extensive follow-up observations of the SQLS, which is aimed at constructing a complete quasar lens sample, our sample of binary quasars is also expected to be nearly complete within a specified range of angular separations and redshifts. The measured small-scale correlation function rises steeply toward smaller scales, which is consistent with earlier studies based on incomplete or smaller binary quasar samples. We find that the quasar correlation function can be fitted by a power-law reasonably well over 4 order of magnitudes, with the best-fit slope of xi(r)\propto r^{-1.92}. We interpret the measured correlation function within the framework of the Halo Occupation Distribution (HOD). We propose a simple model which assumes a constant fraction of quasars that appear as satellites in dark matter haloes, and find that measured small-scale clustering signals constrain the satellite fraction to f_sat=0.054_{-0.016}^{+0.017} for a singular isothermal sphere number density profile of satellites. We note that the HOD modelling appears to underpredict clustering signals at the smallest separations of r_p ~ 10 kpc/h unless we assume very steep number density profiles (such as an NFW profile with the concentration parameter c_vir > 30), which may be suggestive of enhanced quasar activities by direct interactions.

We present a statistical analysis of the final lens sample from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of a complete subsample of 19 lensed quasars selected from 50,836 source quasars is compared with theoretical expectations, with particular attention to the selection function. Assuming that the velocity function of galaxies does not evolve with redshift, the SQLS sample constrains the cosmological constant to \Omega_\Lambda=0.79^{+0.06}_{-0.07}(stat.)^{+0.06}_{-0.06}(syst.) for a flat universe. The dark energy equation of state is found to be consistent with w=-1 when the SQLS is combined with constraints from baryon acoustic oscillation (BAO) measurements or results from the Wilkinson Microwave Anisotropy Probe (WMAP). We also obtain simultaneous constraints on cosmological parameters and redshift evolution of the galaxy velocity function, finding no evidence for redshift evolution at z<1 in any combinations of constraints. For instance, number density evolution quantified as \nu_n=d\ln\phi_*/d\ln(1+z) and the velocity dispersion evolution \nu_\sigma=d\ln\sigma_*/d\ln(1+z) are constrained to \nu_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and \nu_\sigma=-0.05^{+0.19}_{-0.16}(stat.)^{+0.03}_{-0.03}(syst.) respectively when the SQLS result is combined with BAO and WMAP for flat models with a cosmological constant. We find that a significant amount of dark energy is preferred even after fully marginalizing over the galaxy evolution parameters. Thus the statistics of lensed quasars robustly confirm the accelerated cosmic expansion.

We present the final statistical sample of lensed quasars from the Sloan Digital Sky Survey (SDSS) Quasar Lens Search (SQLS). The well-defined statistical lens sample consists of 26 lensed quasars brighter than i=19.1 and in the redshift range of 0.6<z<2.2 selected from 50,836 spectroscopically confirmed quasars in the SDSS Data Release 7 (DR7), where we restrict the image separation range to 1"<\theta<20" and the i-band magnitude differences in two image lenses to be smaller than 1.25 mag. The SDSS DR7 quasar catalog also contains 36 additional lenses identified with various techniques. In addition to these lensed quasars, we have identified 81 pairs of quasars from follow-up spectroscopy, 26 of which are physically associated binary quasars. The statistical lens sample covers a wide range of image separations, redshifts, and magnitudes, and therefore is suitable for systematic studies of cosmological parameters and surveys of the structure and evolution of galaxies and quasars.

The quasar SDSS J133401.39+331534.3 at z = 2.426 is found to be a two-image gravitationally lensed quasar with the image separation of 0.833. The object is first identified as a lensed quasar candidate in the Sloan Digital Sky Survey Quasar Lens Search, and then confirmed as a lensed system from follow-up observations at the Subaru and University of Hawaii 2.2-meter telescopes. We estimate the redshift of the lensing galaxy to be 0.557 based on absorption lines in the quasar spectra as well as the color of the galaxy. In particular, we observe the system with the Subaru Telescope AO188 adaptive optics with laser guide star, in order to derive accurate astrometry, which well demonstrates the usefulness of the laser guide star adaptive optics imaging for studying strong lens systems. Our mass modeling with improved astrometry implies that a nearby bright galaxy $\sim 4"$ apart from the lensing galaxy is likely to affect the lens potential.

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

We present new measurements of the coherent motion of galaxies based on observations of the large-scale redshift-space distortions seen in the two-dimensional two-point correlation function of Luminous Red Galaxies in Data Release Seven of the Sloan Digital Sky Survey. We have developed a new methodology for estimating these coherent motions, which is less dependent on the details of galaxy bias and of the cosmological model to explain the late-time acceleration of the expansion of the Universe. We measure a one-dimensional velocity dispersion of galaxies on large-scales of \sigma_v=3.01^{+0.45}_{-0.46} Mpc/h and \sigma_v=3.69^{+0.47}_{-0.47} \mpcoh$ at a mean redshift of z=0.25 and 0.38 respectively. These values are fully consistent with predictions for a WMAP7-normalised LCDM Universe and inconsistent at confidence of 3.8\sigma with a Dvali-Gabadadze-Porrati (DGP) model for the Universe. We can convert the units of these $\sigma_v$ measurements to 270^{+40}_{-41} km/s and 320^{+41}_{-41} km/s respectively (assuming a $\Lambda$CDM universe), which are lower that expected based on recent low redshift (z<0.2) measurements of the peculiar velocity field (or "bulk flows"). It is difficult to directly compare these measurements as they cover different redshift ranges and different areas of the sky. However, one possible cosmological explanation for this discrepancy is that our Galaxy is located in unusually over, or under, dense region of the Universe.

We study the sample variance of the matter power spectrum for the standard Lambda Cold Dark Matter universe. We use a total of 5000 cosmological N-body simulations to study in detail the distribution of best-fit cosmological parameters and the baryon acoustic peak positions. The obtained distribution is compared with the results from the Fisher matrix analysis with and without including non-Gaussian errors. For the Fisher matrix analysis, we compute the derivatives of the matter power spectrum with respect to cosmological parameters using directly full nonlinear simulations. We show that the non-Gaussian errors increase the unmarginalized errors by up to a factor 5 for k_{max}=0.4h/Mpc if there is only one free parameter provided other parameters are well determined by external information. On the other hand, for multi-parameter fitting, the impact of the non-Gaussian errors is significantly mitigated due to severe parameter degeneracies in the power spectrum. The distribution of the acoustic peak positions is well described by a Gaussian distribution, with its width being consistent with the statistical interval predicted from the Fisher matrix. We also examine systematic bias in the best-fit parameter due to the non-Gaussian errors. The bias is found to be smaller than the 1 sigma statistical error for both the cosmological parameters and the acoustic scale positions.

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<\theta<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be \Omega_\Lambda=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

We report the discovery and confirmation of eight new two-image lensed quasars by the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. The lenses are SDSSJ0904+1512 (image separation \theta=1"13, source redshift z_s=1.826), SDSSJ1054+2733 (\theta=1"27, z_s=1.452), SDSSJ1055+4628 (\theta=1"15, z_s=1.249), SDSSJ1131+1915 (\theta=1"46, z_s=2.915), SDSSJ1304+2001 (\theta=1"87, z_s=2.175), SDSSJ1349+1227 (\theta=3"00, z_s=1.722), SDSSJ1455+1447 (\theta=1"73, z_s=1.424), and SDSSJ1620+1203 (\theta=2"77, z_s=1.158). Three of them, SDSSJ1055+4628, SDSSJ1455+1447, and SDSSJ1620+1203, satisfy the criteria for constructing our statistical sample for studying the cosmological model. Based on galactic absorption lines of the lens galaxies, we also derive lens redshifts of z_l=0.398 and z_l=0.513 for SDSSJ1620+1203 and the previously discovered lens SDSSJ0746+4403, respectively.

We develop an optimized technique to extract density--density and velocity--velocity spectra out of observed spectra in redshift space. The measured spectra of the distribution of halos from redshift distorted mock map are binned into 2--dimensional coordinates in Fourier space so as to be decomposed into both spectra using angular projection dependence. With the threshold limit introduced to minimize nonlinear suppression, the decomposed velocity--velocity spectra are reasonably well measured up to scale k=0.07 h/Mpc, and the measured variances using our method are consistent with errors predicted from a Fisher matrix analysis. The detectability is extendable to k\sim 0.1 h/Mpc with more conservative bounds at the cost of weakened constraint.

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

We report the discovery of five gravitationally lensed quasars from the Sloan Digital Sky Survey (SDSS). All five systems are selected as two-image lensed quasar candidates from a sample of high-redshift (z>2.2) SDSS quasars. We confirmed their lensing nature with additional imaging and spectroscopic observations. The new systems are SDSS J0819+5356 (source redshift z_s=2.237, lens redshift z_l=0.294, and image separation \theta=4.04"), SDSS J1254+2235 (z_s=3.626, \theta=1.56"), SDSS J1258+1657 (z_s=2.702, \theta=1.28"), SDSS J1339+1310 (z_s=2.243, \theta=1.69"), and SDSS J1400+3134 (z_s=3.317, \theta=1.74"). We estimate the lens redshifts of the latter four systems to be z_l=0.2-0.8 from the colors and magnitudes of the lensing galaxies. We find that the image configurations of all systems are well reproduced by standard mass models. Although these lenses will not be included in our statistical sample of z_s<2.2 lenses, they expand the number of lensed quasars which can be used for high-redshift galaxy and quasar studies.

We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a \Lambda Cold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

We present the discovery of a new quadruply lensed quasar. The lens system, SDSS J1330+1810 at z_s=1.393, was identified as a lens candidate from the spectroscopic sample of the Sloan Digital Sky Survey. Optical and near-infrared images clearly show four quasar images with a maximum image separation of 1.76", as well as a bright lensing galaxy. We measure a redshift of the lensing galaxy of z_l=0.373 from absorption features in the spectrum. We find a foreground group of galaxies at z=0.31 centred ~120" southwest of the lens system. Simple mass models fit the data quite well, including the flux ratios between images, although the lens galaxy appears to be ~1 mag brighter than expected by the Faber-Jackson relation. Our mass modelling suggests that shear from nearby structure is affecting the lens potential.

We critically examine how well the evolution of large-scale density perturbations is followed in cosmological $N$-body simulations. We first run a large volume simulation and perform a mode-by-mode analysis in three-dimensional Fourier space. We show that the growth of large-scale fluctuations significantly deviates from linear theory predictions. The deviations are caused by {\it nonlinear} coupling with a small number of modes at largest scales owing to finiteness of the simulation volume. We then develop an analytic model based on second-order perturbation theory to quantify the effect. Our model accurately reproduces the simulation results. For a single realization, the second-order effect appears typically as ``zig-zag'' patterns around the linear-theory prediction, which imprints artificial ``oscillations'' that lie on the real baryon-acoustic oscillations. Although an ensemble average of a number of realizations approaches the linear theory prediction, the dispersions of the realizations remain large even for a large simulation volume of several hundred megaparsecs on a side. For the standard $\Lambda$CDM model, the deviations from linear growth rate are as large as 10 percent for a simulation volume with $L = 500h^{-1}$Mpc and for a bin width in wavenumber of $\Delta k = 0.005h$Mpc$^{-1}$, which are comparable to the intrinsic variance of Gaussian random realizations. We find that the dispersions scales as $\propto L^{-3/2} \Delta k^{-1/2}$ and that the mean dispersion amplitude can be made smaller than a percent only if we use a very large volume of $L > 2h^{-1}$Gpc. The finite box size effect needs to be appropriately taken into account when interpreting results from large-scale structure simulations for future dark energy surveys using baryon acoustic oscillations.

We study the large-scale anisotropic two-point correlation function using 46,760 luminous red galaxies at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We measure the correlation function as a function of separations parallel and perpendicular to the line-of-sight in order to take account of anisotropy of the large-scale structure in redshift space. We find a slight signal of baryonic features in the anisotropic correlation function, i.e., a ``baryon ridge'' which corresponds to a baryon acoustic peak in the spherically averaged correlation function which has already been reported using the same sample. The baryon ridge has primarily a spherical structure with a known radius in comoving coordinates. It enables us to divide the redshift distortion effects into dynamical and geometrical components and provides further constraints on cosmological parameters, including the dark energy equation-of-state. With an assumption of a flat $\Lambda$ cosmology, we find the best-fit values of $\Omega_{\rm m} = 0.218^{+0.047}_{-0.037}$ and $\Omega_{\rm b} = 0.047^{+0.016}_{-0.016}$ (68% C.L.) when we use the overall shape of the anisotropic correlation function of $40<s<200\himpc$ including a scale of baryon acoustic oscillations. When an additional assumption $\Omega_{\rm b}h^2=0.024$ is adopted, we obtain $\Omega_{\rm DE}=0.770^{+0.051}_{-0.040}$ and $w=-0.93^{+0.45}_{-0.35}$. These constraints are estimated only from our data of the anisotropic correlation function, and they agree quite well with values both from the cosmic microwave background (CMB) anisotropies and from other complementary statistics using the LRG sample. With the CMB prior from the 3 year WMAP results, we give stronger constraints on those parameters.

We identify a third image in the unique quasar lens SDSS J1029+2623, the second known quasar lens produced by a massive cluster of galaxies. The spectrum of the third image shows similar emission and absorption features, but has a redder continuum than the other two images which can be explained by differential extinction or microlensing. We also identify several lensed arcs. Our observations suggest a complicated structure of the lens cluster at z~0.6. We argue that the three lensed images are produced by a naked cusp on the basis of successful mass models, the distribution of cluster member galaxies, and the shapes and locations of the lensed arcs. Lensing by a naked cusp is quite rare among galaxy-scale lenses but is predicted to be common among large-separation lensed quasars. Thus the discovery can be viewed as support for an important theoretical prediction of the standard cold dark matter model.

We present the discovery of four gravitationally lensed quasars selected from the spectroscopic quasar catalog of the Sloan Digital Sky Survey. We describe imaging and spectroscopic follow-up observations that support the lensing interpretation of the following four quasars: SDSS J0832+0404 (image separation \theta=1.98", source redshift z_s=1.115, lens redshift z_l=0.659); SDSS J1216+3529 (\theta=1.49", z_s=2.012); SDSS J1322+1052 (\theta=2.00", z_s=1.716); and SDSS J1524+4409 (\theta=1.67", z_s=1.210, z_l=0.320). Each system has two lensed images. We find that the fainter image component of SDSS J0832+0404 is significantly redder than the brighter component, perhaps because of differential reddening by the lensing galaxy. The lens potential of SDSS J1216+3529 might be complicated by the presence of a secondary galaxy near the main lensing galaxy.

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We report the discovery of a new gravitationally lensed quasar, SDSS J131339.98+515128.3, at a redshift of 1.875 with an image separation of 1.24". The lensing galaxy is clearly detected in visible-light follow-up observations. We also identify three absorption-line doublets in the spectra of the lensed quasar images, from which we measure the lens redshift to be 0.194. Like several other known lenses, the lensed quasar images have different continuum slopes. This difference is probably the result of reddening and microlensing in the lensing galaxy. The lensed quasar was selected by correlating Sloan Digital Sky Survey (SDSS) spectroscopic quasars with Two Micron All Sky Survey (2MASS) sources and choosing quasars that show near-infrared (IR) excess. The near-IR excess can originate, for example, from the contribution of the lensing galaxy at near-IR wavelengths. We show that the near-IR excess technique is indeed an efficient method to identify lensed systems from a large sample of quasars.

We report the discovery of a quadruply imaged quasar, SDSSJ125107.57+293540.5, selected from the Sloan Digital Sky Survey. Follow-up imaging reveals that the system consists of four blue point-like components in a typical cusp lens geometry surrounding a central red galaxy. The source redshift is 0.802 and the lens redshift is 0.410. The maximum image separation between the lensed components is 1"79. While the image configuration is well reproduced by standard mass models with reasonable parameter values, the flux ratios predicted by these models differ from the observed ratios in all bands. This is suggestive of small-scale structures in this lens, although the definitive identification of the anomaly requires more accurate photometry and astrometry.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We report the discovery of a cluster-scale lensed quasar, SDSS J1029+2623, selected from the Sloan Digital Sky Survey. The lens system exhibits two lensed images of a quasar at z_s=2.197. The image separation of 22.5" makes it the largest separation lensed quasar discovered to date. The similarity of the optical spectra and the radio loudnesses of the two components support the lensing hypothesis. Images of the field show a cluster of galaxies at z_l~0.55 that is responsible for the large image separation. The lensed images and the cluster light center are not collinear, which implies that the lensing cluster has a complex structure.

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm.

We report the discoveries of the two-image gravitationally lensed quasars, SDSS J0746+4403 and SDSS J1406+6126, selected from the Sloan Digital Sky Survey (SDSS). SDSS J0746+4403, which will be included in our lens sample for statistics and cosmology, has a source redshift of z_s=2.00, an estimated lens redshift of z_l~0.3, and an image separation of 1.08". SDSS J1406+6126 has a source redshift of z_s=2.13, a spectroscopically measured lens redshift of z_l=0.27, and an image separation of 1.98". We find that the two quasar images of SDSS J1406+6126 have different intervening MgII absorption strengths, which are suggestive of large variations of absorbers on kpc scales. The positions and fluxes of both the lensed quasar systems are easily reproduced by simple mass models with reasonable parameter values. These objects bring to 18 the number of lensed quasars that have been discovered from the SDSS data.

We report the discovery of the two-image gravitationally lensed quasar SDSS J133222.62+034739.9 (SDSS J1332+0347) with an image separation of Delta_theta=1.14". This system consists of a source quasar at z_s=1.445 and a lens galaxy at z_l=0.191. The agreement of the luminosity, ellipticity and position angle of the lens galaxy with those expected from lens model confirms the lensing hypothesis.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

We present a first morphometric investigation of a preliminary sample from the SDSS of 154287 galaxies with apparent magnitude 14.5<m_r<17.5 and redshift 0.001<z<0.4. We measure the Minkowski Functionals, which are a complete set of morphological descriptors. To account for the complicated wedge--like geometry of the present survey data, we construct isodensity contour surfaces from the galaxy positions in redshift space and employ two complementary methods of computing the Minkowski Functionals. We find that the observed Minkowski Functionals for SDSS galaxies are consistent with the prediction of a Lambda--dominated spatially--flat Cold Dark Matter model with random--Gaussian initial conditions, within the cosmic variance estimated from the corresponding mock catalogue. We expect that future releases of the SDSS survey will allow us to distinguish morphological differences in the galaxy distribution with regard to different morphological type and luminosity ranges.

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

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

We calculate the one-point probability distribution function (PDF) for cosmic density in non-linear regime of the gravitational evolution. Under the local approximation that the evolution of cosmic fluid fields can be characterized by the Lagrangian local dynamics with finite degrees of freedom, the analytic expressions of PDF are derived taking account of the smoothing effect. The validity and the usefulness of the local approximation are then discussed comparing those results with N-body simulations in a Gaussian initial condition. Adopting the ellipsoidal collapse model (ECM) and the spherical collapse model (SCM) as Lagrangian local dynamics, we found that the PDFs from the local approximation excellently match the simulation results in the case of the cold dark matter initial spectrum. As for the scale-free initial spectra given by $P(k)\propto k^n$, N-body result suffers from spurious numerical effects, which prevent us to give a detailed comparison. Nevertheless, at the quality of N-body data, the model predictions based on the ECM and the SCM quantitatively agree with N-body results in cases with spectral index n<0. For the index n>=0, choice of the Lagrangian local dynamics becomes crucial for an accurate prediction and a more delicate modeling is required, however, we find that the model prediction based on the ECM provides a better approximation to the N-body results of cumulants and PDFs.

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 a general framework to treat the evolution of one-point probability distribution function (PDF) for cosmic density $\delta$ and velocity-divergence fields $\theta$. In particular, we derive an evolution equation for the one-point PDFs and consider the stochastic nature associated with these quantities. Under the local approximation that the evolution of cosmic fluid fields can be characterized by the Lagrangian local dynamics with finite degrees of freedom, evolution equation for PDFs becomes a closed form and consistent formal solutions are constructed. Adopting this local approximation, we explicitly evaluate the one-point PDFs $P(\delta)$ and $P(\theta)$ from the spherical and the ellipsoidal collapse models as the representative Lagrangian local dynamics. In a Gaussian initial condition, while the local density PDF from the ellipsoidal model almost coincides with the that of the spherical model, differences between spherical and ellipsoidal collapse model are found in the velocity-divergence PDF. Importantly, the joint PDF of local density, $P(\delta,t;\delta',t')$, evaluated at the same Lagrangian position but at the different times $t$ and $t'$ from the ellipsoidal collapse model exhibits a large amount of scatter. The mean relation between $\delta$ and $\delta'$ does fail to match the one-to-one mapping obtained from spherical collapse model. Moreover, the joint PDF $P(\delta;\theta)$ from the ellipsoidal collapse model shows a similar stochastic feature, both of which are indeed consistent with the recent result from N-body simulations.

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 present the first analysis of three-dimensional genus statistics for the SDSS EDR galaxy sample. Due to the complicated survey volume and the selection function, analytic predictions of the genus statistics for this sample are not feasible, therefore we construct extensive mock catalogs from N-body simulations in order to compare the observed data with model predictions. This comparison allows us to evaluate the effects of a variety of observational systematics on the estimated genus for the SDSS sample, including the shape of the survey volume, the redshift distortion effect, and the radial selection function due to the magnitude limit. The observed genus for the SDSS EDR galaxy sample is consistent with that predicted by simulations of a $\Lambda$-dominated spatially-flat cold dark matter model. Standard ($\Omega_0=1$) cold dark matter model predictions do not match the observations. We discuss how future SDSS galaxy samples will yield improved estimates of the genus.

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 quantitatively study the probability distribution function (PDF) of cosmological nonlinear density fluctuations from N-body simulations with Gaussian initial condition. In particular, we examine the validity and limitations of one-point and two-point log-normal PDF models against those directly estimated from the simulations. We find that the one-point log-normal PDF describes very accurately the cosmological density distribution even in the nonlinear regime (the rms variance \sigma_{nl} \simlt 4 and the over-density \delta \simlt 100). Furthermore the two-point log-normal PDFs are also in good agreement with the simulation data from linear to fairly nonlinear regime, while slightly deviate from them for \delta \simlt -0.5. Thus the log-normal PDF can be used as a useful empirical model for the cosmological density fluctuations. While this conclusion is fairly insensitive to the shape of the underlying power spectrum of density fluctuations P(k), models with substantial power on large scales, i.e., n\equiv d\ln P(k)/d \ln k \simlt -1, are better described by the log-normal PDF. On the other hand, we note that the one-to-one mapping of the initial and the evolved density fields consistent with the log-normal model does not approximate the broad distribution of their mutual correlation even on average. Thus the origin of the phenomenological log-normal PDF approximation still remains to be understood.