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.

I report the result of a which-way experiment based on Young's double-slit experiment. It reveals which slit photons go through while retaining the (self) interference of all the photons collected. The idea is to image the slits using a lens with a narrow aperture and scan across the area where the interference fringes would be. The aperture is wide enough to separate the slits in the images, i.e., telling which way. The illumination pattern over the pupil is reconstructed from the series of slit intensities. The result matches the double-slit interference pattern well. As such, the photon's wave-like and particle-like behaviors are observed simultaneously in a straightforward and thus unambiguous way. The implication is far reaching. For one, it presses hard, at least philosophically, for a consolidated wave-and-particle description of quantum objects, because we can no longer dismiss such a challenge on the basis that the two behaviors do not manifest at the same time. A bold proposal is to forgo the concept of particles. Then, Heisenberg's uncertainty principle would be purely a consequence of waves without being ordained upon particles.

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.

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