Non-parametric reconstruction of the dark energy equation of state (EoS) aims to determine the EoS as a function of redshift without invoking any particular dark energy model, so that the resulting EoS can be free of model-induced biases or artifacts. Without proper regularization, however, such reconstruction is often overwhelmed by the noise of poorly constrained modes. An intuitive regularization scheme is to assume a priori the dark energy EoS to evolve at most slowly with time, which may be enforced by a correlation between the EoS at different epochs. Indeed, studies that impose the correlation prior are able to significantly reduce the uncertainties of the reconstructed EoS and even show hints for dynamical dark energy. In this work, we examine the correlation prior using mock datasets of type Ia supernovae (SNe Ia), baryonic acoustic oscillations (BAOs), age-derived Hubble parameter, Hubble constant, and cosmic microwave background. We find that even though the prior is designed to disfavor evolving equations of state, it can still accommodate spurious oscillating features at high significance. Within the 1000 mock datasets of existing observations that are generated for the concordance cosmological model, i.e., the input dark energy EoS $w=-1$, there are 688 (69) cases recovering an EoS that departs from $-1$ by more than $1\sigma$ ($2\sigma$) in one or more redshift bins. The reconstructed EoS turns up and down markedly in many cases. Moreover, inverting the signs of the randomly assigned errors of the mock data more or less reverses the behavior of the EoS. Spurious results occur even more frequently when idealized SN Ia and BAO data from future surveys are included. Our tests suggest that further studies are needed to ensure accurate reconstruction of the EoS with the correlation prior.

The Chinese Space Station Optical Survey (CSS-OS) is a planned full sky survey operated by the Chinese Space Station Telescope (CSST). It can simultaneously perform the photometric imaging and spectroscopic slitless surveys, and will probe weak and strong gravitational lensing, galaxy clustering, individual galaxies and galaxy clusters, active galactic nucleus (AGNs), and so on. It aims to explore the properties of dark matter and dark energy and other important cosmological problems. In this work, we focus on two main CSS-OS scientific goals, i.e. the weak gravitational lensing (WL) and galaxy clustering surveys. We generate the mock CSS-OS data based on the observational COSMOS and zCOSMOS catalogs. We investigate the constraints on the cosmological parameters from the CSS-OS using the Markov Chain Monte Carlo (MCMC) method. The intrinsic alignments, galaxy bias, velocity dispersion, and systematics from instrumental effects in the CSST WL and galaxy clustering surveys are also included, and their impacts on the constraint results are discussed. We find that the CSS-OS can improve the constraints on the cosmological parameters by a factor of a few (even one order of magnitude in the optimistic case), compared to the current WL and galaxy clustering surveys. The constraints can be further enhanced when performing joint analysis with the WL, galaxy clustering, and galaxy-galaxy lensing data. Therefore, the CSS-OS is expected to be a powerful survey for exploring the Universe. Since some assumptions may be still optimistic and simple, it is possible that the results from the real survey could be worse. We will study these issues in details with the help of simulations in the future.

The Chinese Space Station Optical Survey (CSS-OS) is a major science project of the Space Application System of the China Manned Space Program. This survey is planned to perform both photometric imaging and slitless spectroscopic observations, and it will focus on different cosmological and astronomical goals. Most of these goals are tightly dependent on the accuracy of photometric redshift (photo-z) measurement, especially for the weak gravitational lensing survey as a main science driver. In this work, we assess if the current filter definition can provide accurate photo-z measurement to meet the science requirement. We use the COSMOS galaxy catalog to create a mock catalog for the CSS-OS. We compare different photo-z codes and fitting methods that using the spectral energy distribution (SED) template-fitting technique, and choose to use a modified LePhare code in photo-z fitting process. Then we investigate the CSS-OS photo-z accuracy in certain ranges of filter parameters, such as band position, width, and slope. We find that the current CSS-OS filter definition can achieve reasonably good photo-z results with sigma_z~0.02 and outlier fraction ~3%.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The cosmic expansion history, mapped by the Hubble parameter as a function of redshift, offers the most direct probe of the dark energy equation of state. One way to determine the Hubble parameter at different redshifts is essentially differentiating the cosmic age or distance with respect to redshift, which may incur large numerical errors with observational data. Taking the scenario that the Hubble parameter increases monotonically with redshift as a reasonable prior, we propose to enforce the monotonicity when reconstructing the Hubble parameter at a series of redshifts. Tests with mock type Ia supernova (SN Ia) data show that the monotonicity prior does not introduce significant biases and that errors on the Hubble parameter are greatly reduced compared to those determined with a flat prior at each redshift. Results from real SN Ia data are in good agreement with those based on ages of passively evolving galaxies. Although the Hubble parameter reconstructed from SN Ia distances does not necessarily provide more information than the distances themselves do, it offers a convenient way to compare with constraints from other methods. Moreover, the monotonicity prior is expected to be helpful to other probes that measure the Hubble parameter at multiple redshifts (e.g., baryon acoustic oscillations), and it may be generalized to other cosmological quantities that are reasonably monotonic with redshift.

The Large Synoptic Survey Telescope (LSST) is a high etendue imaging facility that is being constructed atop Cerro Pachon in Northern Chile. It is scheduled to begin science operations in 2022. With an 8.4m (6.5m effective) aperture, a three-mirror design achieving a seeing-limited 9.6deg^2 field of view, and a 3.2 Gigapixel camera, the LSST has the deep-wide-fast imaging capability to carry out an 18,000deg^2 survey in six passbands (ugrizy) to a coadded depth of r~27.5 over 10 years using 90% of its observational time. The remaining 10% will be devoted to deeper and faster time-domain observations and smaller surveys. In total, each patch of the sky in the main survey will receive 800 visits allocated across the six passbands with 30s exposure visits. The huge volume of high-quality LSST data will provide a wide range of science opportunities and open a new era of precision cosmology with unprecedented statistical power and tight control of systematic errors. In this review, we give a brief account of the LSST cosmology program with an emphasis on dark energy investigations. The LSST will address dark energy physics and cosmology in general by exploiting diverse precision probes including large-scale structure, weak lensing, type Ia supernovae, galaxy clusters, and strong lensing. Combined with the cosmic microwave background data, these probes form interlocking tests on the cosmological model and the nature of dark energy. The LSST data products will be made available to the U.S. and Chilean scientific communities and to international partners with no proprietary period. Close collaborations with coeval surveys observing at a variety of wavelengths, resolutions, depths, and timescales will be a vital part of the LSST science program, which will not only enhance specific studies but also allow a more complete understanding of the universe through different windows.

Intrapixel nonuniformity is known to exist in CCD and CMOS image sensors, though the effects in backside illuminated (BSI) CCDs are too small to be a concern for most astronomical observations. However, projects like the Large Synoptic Survey Telescope require precise knowledge of the detector characteristics, and intrapixel effects may need more attention. By scanning CCD and CMOS cameras with a small light spot (unresolved by the optics), we find in the images that the spot's flux, centroid displacement, and ellipticity vary periodically on the pixel scale in most cases. The amplitude of variation depends on not only the detector but also how well the spot is sampled by the pixels. With a spot radius of 2 pixels (encircling 80% energy) as measured, the flux and the ellipticity extracted from the BSI CCD camera vary by 0.2-0.3% (rms) and 0.005 (rms), respectively, while the deviation of the centroid position (rms ~ 0.01 pixel) is not correlated with the pixels. The effects are more pronounced for the BSI CMOS camera and even worse for the frontside illuminated CMOS camera. The results suggest that a closer examination of the intrapixel effects is needed for precision astronomy.

The Large Synoptic Survey Telescope (LSST) will survey the southern sky from 2022--2032 with unprecedented detail. Since the observing strategy can lead to artifacts in the data, we investigate the effects of telescope-pointing offsets (called dithers) on the $r$-band coadded 5$\sigma$ depth yielded after the 10-year survey. We analyze this survey depth for several geometric patterns of dithers (e.g., random, hexagonal lattice, spiral) with amplitude as large as the radius of the LSST field-of-view, implemented on different timescales (per season, per night, per visit). Our results illustrate that per night and per visit dither assignments are more effective than per season. Also, we find that some dither geometries (e.g., hexagonal lattice) are particularly sensitive to the timescale on which the dithers are implemented, while others like random dithers perform well on all timescales. We then model the propagation of depth variations to artificial fluctuations in galaxy counts, which are a systematic for large-scale structure studies. We calculate the bias in galaxy counts caused by the observing strategy, accounting for photometric calibration uncertainties, dust extinction, and magnitude cuts; uncertainties in this bias limit our ability to account for structure induced by the observing strategy. We find that after 10 years of the LSST survey, the best dither strategies lead to uncertainties in this bias smaller than the minimum statistical floor for a galaxy catalog as deep as $r$$<$27.5. A few of these strategies bring the uncertainties close to the statistical floor for $r$$<$25.7 after only one year of survey.

The control of photometric redshift (photo-$z$) errors is a crucial and challenging task for precision weak lensing cosmology. The spacial cross-correlations (equivalently, the angular cross power spectra) of galaxies between tomographic photo-$z$ bins are sensitive to the true redshift distribution $n_i(z)$ of each bin and hence can help calibrate the photo-$z$ error distribution for weak lensing surveys. Using Fisher matrix analysis, we investigate the contributions of various components of the angular power spectra to the constraints of $n_i(z)$ parameters and demonstrate the importance of the cross power spectra therein, especially when catastrophic photo-$z$ errors are present. We further study the feasibility of reconstructing $n_i(z)$ from galaxy angular power spectra using Markov Chain Monte Carlo estimation. Considering an LSST-like survey with $10$ photo-$z$ bins, we find that the underlying redshift distribution can be determined with a fractional precision ($\sigma(\theta)/\theta$ for parameter $\theta$) of roughly $1\%$ and $10\%$ for the mean redshift and width of $n_i(z)$, respectively.

The likelihood function is a crucial element of parameter estimation. In analyses of galaxy overdensities and weak lensing shear, one often approximates the likelihood of the power spectrum with a Gaussian distribution. The posterior probability derived from such a likelihood deviates considerably from the exact posterior on the largest scales probed by any survey, where the central limit theorem does not apply. We show that various forms of Gaussian likelihoods can have a significant impact on the estimation of the primordial non-Gaussianity parameter f_NL from the galaxy angular power spectrum. The Gaussian plus log-normal likelihood, which has been applied successfully in analyses of the cosmic microwave background, outperforms the Gaussian likelihoods. Nevertheless, even if the exact likelihood of the power spectrum is used, the estimated parameters may be still biased. As such, the likelihoods and estimators need to be thoroughly examined for potential systematic errors.

We characterize the baryon acoustic oscillations (BAO) feature in halo two-point statistics using N-body simulations. We find that nonlinear damping of the BAO signal is less severe for halos in the mass range we investigate than for dark matter. The amount of damping depends weakly on the halo mass. The correlation functions show a mass-dependent drop of the halo clustering bias below roughly 90 Mpc/h, which coincides with the scale of the BAO trough. The drop of bias is 4% for halos with mass M>10^{14} Msun/h and reduces to roughly 2% for halos with mass M>10^{13} Msun/h. In contrast, halo biases in simulations without BAO change more smoothly around 90 Mpc/h. In Fourier space, the bias of M>10^{14} Msun/h halos decreases smoothly by 11% from wavenumber k = 0.012 h/Mpc to 0.2 h/Mpc, whereas that of M>10^{13} Msun/h halos decreases by less than 4% over the same range. By comparing the halo biases in pairs of otherwise identical simulations, one with and the other without BAO, we also observe a modulation of the halo bias. These results suggest that precise calibrations of the mass-dependent BAO signal and scale-dependent bias on large scales would be needed for interpreting precise measurements of the two-point statistics of clusters or massive galaxies in the future.

We calculate the correlation function of 13,904 galaxy clusters of z \leq 0.4 selected from the cluster catalog of Wen, Han & Liu. The correlation function can be fitted with a power-law model \xi(r)=(r/R_0)^{\gamma} on the scales of 10h^(-1)Mpc \leq r \leq 50h^(-1)Mpc, with a larger correlation length of R_0=18.84\pm0.27 h^(-1)Mpc for clusters with a richness of R \geq 15 and a smaller length of R_0=16.15 \pm 0.13 h^(-1)Mpc for clusters with a richness of R \geq 5. The power law index of \gamma=2.1 is found to be almost the same for all cluster subsamples. A pronounced baryon acoustic oscillations (BAO) peak is detected at r ~ 110 h^(-1)Mpc with a significance of ~ 1.9\sigma. By analyzing the correlation function in the range of 20 h^(-1)Mpc \leq r \leq 200 h^(-1)Mpc, we find the constraints on distance parameters are D_v(0.276)=1077\pm55(1\sigma) Mpc and h=0.73 \pm 0.039 (1\sigma), which are consistent with the WMAP 7-year cosmology. However, the BAO signal from the cluster sample is stronger than expected and leads to a rather low matter density \Omega_m h^2=0.093\pm0.0077 (1\sigma), which deviates from the WMAP 7-year result by more than 3 \sigma. The correlation function of the GMBCG cluster sample is also calculated and our detection of the BAO feature is confirmed.

Aiming at comparing different morphological models of galaxy clusters, we use two new methods to make a cosmological model-independent test of the distance-duality (DD) relation. The luminosity distances come from Union2 compilation of Supernovae Type Ia. The angular diameter distances are given by two cluster models (De Filippis et al. and Bonamente et al.). The advantage of our methods is that it can reduce statistical errors. Concerning the morphological hypotheses for cluster models, it is mainly focused on the comparison between elliptical $\beta$-model and spherical $\beta$-model. The spherical $\beta$-model is divided into two groups in terms of different reduction methods of angular diameter distances, i.e. conservative spherical $\beta$-model and corrected spherical $\beta$-model. Our results show that the DD relation is consistent with the elliptical $\beta$-model at $1\sigma$ confidence level (CL) for both methods, whereas for almost all spherical $\beta$-model parameterizations, the DD relation can only be accommodated at $3\sigma$ CL, particularly for the conservative spherical $\beta$-model. In order to minimize systematic uncertainties, we also apply the test to the overlap sample, i.e. the same set of clusters modeled by both De Filippis et al. and Bonamente et al.. It is found that the DD relation is compatible with the elliptically modeled overlap sample at $1\sigma$ CL, however for most of the parameterizations, the DD relation can not be accommodated even at $3\sigma$ CL for any of the two spherical $\beta$-models. Therefore it is reasonable that the marked triaxial ellipsoidal model is a better geometrical hypothesis describing the structure of the galaxy cluster compared with the spherical $\beta$-model if the DD relation is valid in cosmological observations.

We estimate the Rees--Sciama (RS) effect of super structures on the cosmic microwave background (CMB) temperature fluctuations and identify a related effect on galaxy redshifts. By numerically solving the geodesic equation, we find that both superclusters and supervoids can decrease the temperature of the CMB by several micro Kelvin in the central region and increase the temperature slightly in the surrounding area due to the RS effect. The two components of the RS effect, redshift and gravitational time delay, largely cancel each other, leaving an equivalent but much smaller effect on the CMB photons that started out at the same time from the distorted last scattering surface. For galaxies, the time delay effect is separable from the redshift effect, and the slight change to the redshift induced by super structures can be at the percent level of large-scale rms bulk velocities, which might only be detected statistically. On much smaller scales, a tiny redshift difference between two images of a strongly lensed source should exist in general, which is related to the Hubble expansion rate at the source redshift. However, as Loeb (1998) pointed out, observational issues and the proper motion of the structure would make such a measurement impossible.

The uncertainty in the redshift distributions of galaxies has a significant potential impact on the cosmological parameter values inferred from multi-band imaging surveys. The accuracy of the photometric redshifts measured in these surveys depends not only on the quality of the flux data, but also on a number of modeling assumptions that enter into both the training set and SED fitting methods of photometric redshift estimation. In this work we focus on the latter, considering two types of modeling uncertainties: uncertainties in the SED template set and uncertainties in the magnitude and type priors used in a Bayesian photometric redshift estimation method. We find that SED template selection effects dominate over magnitude prior errors. We introduce a method for parameterizing the resulting ignorance of the redshift distributions, and for propagating these uncertainties to uncertainties in cosmological parameters.

Dark energy is an important science driver of many upcoming large-scale surveys. With small, stable seeing and low thermal infrared background, Dome A, Antarctica, offers a unique opportunity for shedding light on fundamental questions about the universe. We show that a deep, high-resolution imaging survey of 10,000 square degrees in \emph{ugrizyJH} bands can provide competitive constraints on dark energy equation of state parameters using type Ia supernovae, baryon acoustic oscillations, and weak lensing techniques. Such a survey may be partially achieved with a coordinated effort of the Kunlun Dark Universe Survey Telescope (KDUST) in \emph{yJH} bands over 5000--10,000 deg$^2$ and the Large Synoptic Survey Telescope in \emph{ugrizy} bands over the same area. Moreover, the joint survey can take advantage of the high-resolution imaging at Dome A to further tighten the constraints on dark energy and to measure dark matter properties with strong lensing as well as galaxy--galaxy weak lensing.

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

Much of the science case for the next generation of deep, wide-field optical/infrared surveys has been driven by the further study of dark energy. This is a laudable goal (and the subject of a companion white paper by Zhan et al.). However, one of the most important lessons of the current generation of surveys is that the interesting science questions at the end of the survey are quite different than they were when the surveys were being planned. The current surveys succeeded in this evolving terrain by being very general tools that could be applied to a number of very fundamental measurements. Likewise, the accessibility of the data enabled the broader cosmological and astronomical community to generate more science than the survey collaborations could alone. With that in mind, we should consider some of the basic physical and cosmological questions that surveys like LSST and JDEM-Wide will be able to address.

The coming decade will be an exciting period for dark energy research, during which astronomers will address the question of what drives the accelerated cosmic expansion as first revealed by type Ia supernova (SN) distances, and confirmed by later observations. The mystery of dark energy poses a challenge of such magnitude that, as stated by the Dark Energy Task Force (DETF), "nothing short of a revolution in our understanding of fundamental physics will be required to achieve a full understanding of the cosmic acceleration." The lack of multiple complementary precision observations is a major obstacle in developing lines of attack for dark energy theory. This lack is precisely what next-generation surveys will address via the powerful techniques of weak lensing (WL) and baryon acoustic oscillations (BAO) -- galaxy correlations more generally -- in addition to SNe, cluster counts, and other probes of geometry and growth of structure. Because of their unprecedented statistical power, these surveys demand an accurate understanding of the observables and tight control of systematics. This white paper highlights the opportunities, approaches, prospects, and challenges relevant to dark energy studies with wide-deep multiwavelength photometric redshift surveys. Quantitative predictions are presented for a 20000 sq. deg. ground-based 6-band (ugrizy) survey with 5-sigma depth of r~27.5, i.e., a Stage 4 survey as defined by the DETF.

Baryon acoustic oscillations (BAOs) and weak lensing (WL) are complementary probes of cosmology. We explore the distance and growth factor measurements from photometric BAO and WL techniques and investigate the roles of the distance and growth factor in constraining dark energy. We find for WL that the growth factor has a great impact on dark energy constraints but is much less powerful than the distance. Dark energy constraints from WL are concentrated in considerably fewer distance eigenmodes than those from BAO, with the largest contributions from modes that are sensitive to the absolute distance. Both techniques have some well determined distance eigenmodes that are not very sensitive to the dark energy equation of state parameters w_0 and w_a, suggesting that they can accommodate additional parameters for dark energy and for the control of systematic uncertainties. A joint analysis of BAO and WL is far more powerful than either technique alone, and the resulting constraints on the distance and growth factor will be useful for distinguishing dark energy and modified gravity models. The Large Synoptic Survey Telescope (LSST) will yield both WL and angular BAO over a sample of several billion galaxies. Joint LSST BAO and WL can yield 0.5% level precision on ten comoving distances evenly spaced in log(1+z) between redshift 0.3 and 3 with cosmic microwave background priors from Planck. In addition, since the angular diameter distance, which directly affects the observables, is linked to the comoving distance solely by the curvature radius in the Friedmann-Robertson-Walker metric solution, LSST can achieve a pure metric constraint of 0.017 on the mean curvature parameter Omega_k of the universe simultaneously with the constraints on the comoving distances.

Since type Ia Supernovae (SNe) explode in galaxies, they can, in principle, be used as the same tracer of the large-scale structure as their hosts to measure baryon acoustic oscillations (BAOs). To realize this, one must obtain a dense integrated sampling of SNe over a large fraction of the sky, which may only be achievable photometrically with future projects such as the Large Synoptic Survey Telescope. The advantage of SN BAOs is that SNe have more uniform luminosities and more accurate photometric redshifts than galaxies, but the disadvantage is that they are transitory and hard to obtain in large number at high redshift. We find that a half-sky photometric SN survey to redshift z = 0.8 is able to measure the baryon signature in the SN spatial power spectrum. Although dark energy constraints from SN BAOs are weak, they can significantly improve the results from SN luminosity distances of the same data, and the combination of the two is no longer sensitive to cosmic microwave background priors.

Using a number of numerical tests and analytic arguments we investigate how measurements of cosmic shear lead to constraints on dark energy. We find that, in contrast to the case with galaxy number density correlation functions, standard rulers in the matter power spectrum play no significant role. Sensitivity to distance ratios is provided by the ratios in the lensing kernel. An absolute distance scale can only be established by breaking a potential degeneracy between growth and distance which can be done if the growth-redshift relation and distance-redshift relations are parameterized with sufficiently few parameters. For the quality of dark energy determination, growth determination is primarily important because it improves the distance reconstructions. The information about dark energy in the growth-redshift relation is always of secondary importance though the amount it contributes is highly dependent on what priors are taken in the cosmological parameter space. We also explain the dependence of dark energy constraints from cosmic shear, relative distance measures (supernovae) and absolute distance measures (baryon acoustic oscillations) on assumptions about the mean curvature.

We explore the complementarity between two tomographic probes of the universe: baryon acoustic oscillations (BAOs) and weak lensing (WL). We find that photometric measurements of BAOs alone do not provide very tight constraints on the dark energy equation of state parameters, partially due to our uncertain knowledge of the galaxy bias. WL, on the other hand, is adversely impacted by the uncertainties of the probability distribution of photometric redshift (photo-z) errors. A joint analysis of the two, however, is more robust to these uncertainties and leads to a remarkable improvement over the results of either probe alone. Forecasts of cosmological constraints are provided for various BAO and WL surveys in combination with CMB measurements from Planck. In particular, we find that the joint analysis of BAO and WL power spectra with the Large Synoptic Survey Telescope can tighten the errors of the dark energy equation of state (at a = 0.63) and its rate of change, respectively, to 0.016 and 0.16 (marginalized over 131 other parameters). With supernovae and cluster counting as well as higher-order statistics of the same galaxy and shear data, one can further improve the constraints.

We investigate how well the redshift distributions of galaxies sorted into photometric redshift bins can be determined from the galaxy angular two-point correlation functions. We find that the uncertainty in the reconstructed redshift distributions depends critically on the number of parameters used in each redshift bin and the range of angular scales used, but not on the number of photometric redshift bins. Using six parameters for each photometric redshift bin, and restricting ourselves to angular scales over which the galaxy number counts are normally distributed, we find that errors in the reconstructed redshift distributions are large; i.e., they would be the dominant source of uncertainty in cosmological parameters estimated from otherwise ideal weak lensing or baryon acoustic oscillation data. However, either by reducing the number of free parameters in each redshift bin, or by (unjustifiably) applying our Gaussian analysis into the non-Gaussian regime, we find that the correlation functions can be used to reconstruct the redshift distributions with moderate precision; e.g., with mean redshifts determined to around 0.01. We also find that dividing the galaxies into two spectral types, and thereby doubling the number of redshift distribution parameters, can result in a reduction in the errors in the combined redshift distributions.

With a wariness of Occam's razor awakened by the discovery of cosmic acceleration, we abandon the usual assumption of zero mean curvature and ask how well it can be determined by planned surveys. We also explore the impact of uncertain mean curvature on forecasts for the performance of planned dark energy probes. We find that weak lensing and photometric baryon acoustic oscillation data, in combination with CMB data, can determine the mean curvature well enough that the residual uncertainty does not degrade constraints on dark energy. We also find that determinations of curvature are highly tolerant of photometric redshift errors.

We consider cosmological applications of galaxy number density correlations to be inferred from future deep and wide multi-band optical surveys. We mostly focus on very large scales as a probe of possible features in the primordial power spectrum. We find the proposed survey of the Large Synoptic Survey Telescope may be competitive with future all-sky CMB experiments over a broad range of scales. On very large scales the inferred power spectrum is robust to photometric redshift errors, and, given a sufficient number density of galaxies, to angular variations in dust extinction and photometric calibration errors. We also consider other applications, such as constraining dark energy with the two CMB-calibrated standard rulers in the matter power spectrum, and controlling the effect of photometric redshift errors to facilitate the interpretation of cosmic shear data. We find that deep photometric surveys over wide area can provide constraints that are competitive with spectroscopic surveys in small volumes.

We propose to determine the optimal softening length in N-body halo simulations by minimizing the ensemble-average acceleration error at a small radius r0. This strategy ensures that the error never exceeds the optimal value beyond r0. Furthermore, we derive semi-analytic formulae for calculating the acceleration error due to the discreteness of particles and softened gravity, which are validated by direct N-body force calculations. We estimate that current state-of-the-art halo simulations suffer >~6% acceleration error at 1% of the halo virial radius. The error grows rapidly toward the center and could contribute significantly to the uncertainties of inner halo properties.

Weak lensing surveys that can potentially place strong constraints on dark energy parameters can only do so if the source redshift means and error distributions are very well known. We investigate prospects for controlling errors in these quantities by exploiting their influence on the power spectra of the galaxies. Although, from the galaxy power spectra alone, sufficiently precise and simultaneous determination of redshift biases and variances is not possible, a strong consistency test is. Given the redshift error rms, galaxy power spectra can be used to determine the mean redshift of a group of galaxies to subpercent accuracy. Although galaxy power spectra cannot be used to determine the redshift error rms, they can be used to determine this rms divided by the Hubble parameter, a quantity that may be even more valuable for interpretation of cosmic shear data than the rms itself. We also show that galaxy power spectra, due to the baryonic acoustic oscillations, can potentially lead to constraints on dark energy that are competitive with those due to the cosmic shear power spectra from the same survey.

We analyze the flux power spectrum and its covariance using simulated Lyman alpha forests. We find that pseudo-hydro techniques are good approximations of hydrodynamical simulations at high redshift. However, the pseudo-hydro techniques fail at low redshift because they are insufficient for characterizing some components of the low-redshift intergalactic medium, notably the warm-hot intergalactic medium. Hence, to use the low-redshift Lyman alpha flux power spectrum to constrain cosmology, one would need realistic hydrodynamical simulations. By comparing one-dimensional mass statistics with flux statistics, we show that the nonlinear transform between density and flux quenches the fluctuations so that the flux power spectrum is much less sensitive to cosmological parameters than the one-dimensional mass power spectrum. The covariance of the flux power spectrum is nearly Gaussian. As such, the uncertainties of the underlying mass power spectrum could still be large, even though the flux power spectrum can be precisely determined from a small number of lines of sight.

We analyse the covariance of the one-dimensional mass power spectrum along lines of sight. The covariance reveals the correlation between different modes of fluctuations in the cosmic density field and gives the sample variance error for measurements of the mass power spectrum. For Gaussian random fields, the covariance matrix is diagonal. As expected, the variance of the measured one-dimensional mass power spectrum is inversely proportional to the number of lines of sight that are sampled from each random field. The correlation between lines of sight in a single field may alter the covariance. However, lines of sight that are sampled far apart are only weakly correlated, so that they can be treated as independent samples. Using N-body simulations, we find that the covariance matrix of the one-dimensional mass power spectrum is not diagonal for the cosmic density field due to the non-Gaussianity and that the variance is much higher than that of Gaussian random fields. From the covariance, one will be able to determine the cosmic variance in the measured one-dimensional mass power spectrum as well as to estimate how many lines of sight are needed to achieve a target precision.

We investigate the impact of the intracluster medium on the weak-lensing shear power spectrum (PS). Using a halo model we find that, compared to the dark matter only case, baryonic pressure leads to a suppression of the shear PS on the order of a few percent or more for $l \gtrsim 1000$. Cooling/cooled baryons and the intergalactic medium can further alter the shear PS. Therefore, the interpretation of future precision weak lensing data at high multipoles must take into account the effects of baryons.

(Abridged) I investigate statistical properties of one-dimensional fields in the universe such as the Lyman alpha forest and inverted line-of-sight densities. Because of gravitational clustering, the cosmic density field is already quite non-Gaussian on scales below 10 Mpc/h at redshift z=3. The covariance of the one-dimensional mass power spectrum indicates that Fourier modes in the cosmic density field are highly correlated and that the variance of the measured one-dimensional mass power spectrum is much higher than the expectation for Gaussian random fields. It is found that rare high-density structures contribute significantly to the covariance. The window function due to the length of lines of sight introduces additional correlations between different Fourier modes. In practice, one observes quasar spectra instead of one-dimensional density fields. I show that the nonlinear transform between density and flux quenches the fluctuations so that the flux power spectrum is less sensitive to cosmological parameters than the one-dimensional mass power spectrum. Finally, I provide a method for inverting Lyman alpha forests and obtaining line-of-sight densities, so that statistics can be measured from one-dimensional density fields directly.

The case is made that the intracluster medium (ICM) in spiral-rich galaxy groups today probably has undergone much slower evolution than that in elliptical-rich groups and clusters. The environments of proto-clusters and proto-groups at z > 2 are likely similar to spiral-rich group environments at lower redshift. Therefore, like the ICM in spiral-rich groups today, the ICM in proto-groups and proto-clusters at z > 2 is predicted to be significantly multiphased. The QSO Lyman alpha forest in the vicinity of galaxies is an effective probe of the ICM at a wide range of redshift. Two recent observations of Lyman alpha absorption around galaxies by Adelberger et al. and by Pascarelle et al are reconciled, and it is shown that observations support the multiphase ICM scenario. Galaxy redshifts must be very accurate for such studies to succeed. This scenario can also explain the lower metallicity and lower hot gas fraction in groups.

We report a scaling of the mean matter density with the width of the saturated Lyman alpha absorptions. This property is established using the ``pseudo-hydro'' technique (Croft et al. 1998). It provides a constraint for the inversion of the Lyman alpha forest, which encounters difficulty in the saturated region. With a Gaussian density profile and the scaling relation, a simple inversion of the simulated Lyman alpha forests shows that the one-dimensional mass power spectrum is well recovered on scales above 2 Mpc/h, or roughly k < 0.03 s/km, at z=3. The recovery underestimates the power on small scales, but improvement is possible with a more sophisticated algorithm.

In this paper, we develop a method of performing the one-point statistics of a perturbed density field with a multiresolutional decomposition based on the discrete wavelet transform (DWT). We establish the algorithm of the one-point variable and its moments in considering the effects of Poisson sampling and selection function. We also establish the mapping between the DWT one-point statistics in redshift space and real space, i.e. the algorithm for recovering the DWT one-point statistics from the redshift distortion of bulk velocity, velocity dispersion, and selection function. Numerical tests on N-body simulation samples show that this algorithm works well on scales from a few hundreds to a few Mpc/h for four popular cold dark matter models. Taking the advantage that the DWT one-point variable is dependent on both the scale and the shape (configuration) of decomposition modes, one can design estimators of the redshift distortion parameter (beta) from combinations of DWT modes. When the non-linear redshift distortion is not negligible, the beta estimator from quadrupole-to-monopole ratio is a function of scale. This estimator would not work without adding information about the scale-dependence, such as the power-spectrum index or the real-space correlation function of the random field. The DWT beta estimators, however, do not need such extra information. Numerical tests show that the proposed DWT estimators are able to determine beta robustly with less than 15% uncertainty in the redshift range 0 < z < 3.

We propose a method for estimating the large-scale rms bulk velocity of the cosmic mass field from the transmitted fluxes of Ly$\alpha$ forests. It is based on two linear relationships on large scales: 1) the relation between the fluctuations of the transmission and the underlying density field, and 2) the relation between the density fluctuations and the peculiar velocity field. We show that, with a multiscale decomposition, the two relations can be effectively employed for predicting the rms bulk velocity. Since QSO's Ly$\alpha$ forest is due to the absorptions of diffusely distributed and photoionized IGM, this method provides an independent estimate of the rms bulk velocity at high redshifts, on large scales, and free from the bias of galaxies. Using the transmitted flux of 60 moderate-resolution QSO spectra, the rms bulk velocity is found to be 230$\pm$50 km s$^{-1}$ around redshift $z=2.25$ on scale 23 h$^{-1}$Mpc, and down to 110$\pm$45 km s$^{-1}$ around $z=3.25$ on scale 92 h$^{-1}$Mpc for an LCDM universe ($\Omega=0.3$ and $\Lambda=0.7$). The results are basically consistent with the linear evolution theory.

We analyze the power spectrum of a QSO's Ly_alpha transmitted flux in the discrete wavelet transform (DWT) representation. Although the mean DWT power spectrum is consistent with its counterpart in Fourier representation, the spatial distribution of the local power varies greatly, i.e. the local DWT power spectra show remarkably spiky structures on small scales. To measure these spiky features, we introduce the quantities, roughness, of the local power spectrum, and the correlation between spikes on different scales. We then test the predictions made by the correlation hierarchy model on the roughness and the scale-scale correlations of the local power spectrum. Using the Ly_alpha transmitted flux of the QSO HS1700, we find that the underlying cosmic mass field of the transmitted flux at redshift around z ~ 2.2 can be described by the hierarchical clustering model on physical scales from 2.5 h^{-1} Mpc to few tens h^{-1} kpc in an Einstein-de Sitter universe. However, the non-linear features of the clustering show differences on different scale ranges; 1. On physical scales larger than ~ 1.3 h^{-1} Mpc, the field is almost Gaussian. 2. On scales 1.3 h^{-1} Mpc - 0.3 h^{-1} Mpc, the field is consistent with the correlation hierarchy with a constant value for the coefficient Q_4. 3. On scales less than 300 h^{-1} kpc, the field is no longer Gaussian, but essentially intermittent. In this case, the field can still be fitted by the correlation hierarchy, but the coefficient, Q_4, should be scale-dependent. These three points are strongly supported by the following result: the scale dependencies of Q_4 given by two statistically independent measures, i.e. Q_4^R by the roughness and Q_4^C by scale-scale correlation, are the same in the entire scale range considered.

The intermittent behavior of the space-scale distribution of Ly$\alpha$ transmitted flux of QSO HS1700+64 has been analyzed via a discrete wavelet transform. We found that there are strong indications of intermittency on scales down to about 10 $h^{-1}$ kpc. These are: 1.) the probability distribution function of the local fluctuations of the flux is significantly long-tailed on small scales, and 2.) the local power spectrum of the flux shows prominent spiky structures on small scales. Moreover, the local power spectrum averaged on regions with different sizes shows similar spiky structures. Therefore, the random mass density field traced by the Ly$\alpha$ forests is rougher on smaller scales, consistent with singular clustering.