The Nobel Prize winning confirmation in 1998 of the accelerated expansion of our Universe put into sharp focus the need of a consistent theoretical model to explain the origin of this acceleration. As a result over the past two decades there has been a huge theoretical and observational effort into improving our understanding of the Universe. The cosmological equations describing the dynamics of a homogeneous and isotropic Universe are systems of ordinary differential equations, and one of the most elegant ways these can be investigated is by casting them into the form of dynamical systems. This allows the use of powerful analytical and numerical methods to gain a quantitative understanding of the cosmological dynamics derived by the models under study. In this review we apply these techniques to cosmology. We begin with a brief introduction to dynamical systems, fixed points, linear stability theory, Lyapunov stability, centre manifold theory and more advanced topics relating to the global structure of the solutions. Using this machinery we then analyse a large number of cosmological models and show how the stability conditions allow them to be tightly constrained and even ruled out on purely theoretical grounds. We are also able to identify those models which deserve further in depth investigation through comparison with observational data. This review is a comprehensive and detailed study of dynamical systems applications to cosmological models focusing on the late-time behaviour of our Universe, and in particular on its accelerated expansion. In self contained sections we present a large number of models ranging from canonical and non-canonical scalar fields, interacting models and non-scalar field models through to modified gravity scenarios. Selected models are discussed in detail and interpreted in the context of late-time cosmology.

The morphological properties of large scale structure of the Universe can be fully described by four Minkowski functionals (MFs), which provide important complementary information to other statistical observables such as the widely used 2-point statistics in configuration and Fourier spaces. In this work, for the first time, we present the differences in the morphology of large scale structure caused by modifications to general relativity (to address the cosmic acceleration problem), by measuring the MFs from N-body simulations of modified gravity and general relativity. We find strong statistical power when using the MFs to constrain modified theories of gravity: with a galaxy survey that has survey volume $\sim 0.125 (h^{-1}$Gpc$)^3$ and galaxy number density $\sim 1 / (h^{-1}$Mpc$)^{3}$, the two normal-branch DGP models and the F5 $f(R)$ model that we simulated can be discriminated from $\Lambda$CDM at a significance level >~ 5$\sigma$ with an individual MF measurement. Therefore, the MF of large scale structure is potentially a powerful probe of gravity, and its application to real data deserves active explorations.

One of the main goals of modern cosmology is to search for primordial gravitational waves by looking on their imprints in the B-type polarization in the cosmic microwave background radiation. However, this signal is contaminated by various sources, including cosmic weak lensing, foreground radiations, instrumental noises, as well as the E-to-B leakage caused by the partial sky surveys, which should be well understood to avoid the misinterpretation of the observed data. In this paper, we adopt the E/B decomposition method suggested by Smith in 2006, and study the imprints of E-to-B leakage residuals in the constructed B-type polarization maps, $\mathcal{B}(\hat{n})$, by employing various statistical tools. We find that the effects of E-to-B leakage are negligible for the $\mathcal{B}$-mode power spectrum, as well as the skewness and kurtosis analyses of $\mathcal{B}$-maps. However, if employing the morphological statistical tools, including Minkowski functionals and/or Betti numbers, we find the effect of leakage can be detected at very high confidence level, which shows that in the morphological analysis, the leakage can play a significant role as a contaminant for measuring the primordial B-mode signal and must be taken into account for a correct explanation of the data.

We give the three-dimensional dynamical autonomous systems for most of the popular scalar field dark energy models including (phantom) quintessence, (phantom) tachyon, k-essence and general non-canonical scalar field models, change the dynamical variables from variables $(x, y, \lambda)$ to observable related variables $(w_{\phi}, \Omega_{\phi}, \lambda)$, and show the intimate relationships between those scalar fields that the three-dimensional system of k-essence can reduce to (phantom) tachyon, general non-canonical scalar field can reduce to (phantom) quintessence and k-essence can also reduce to (phantom) quintessence for some special cases. For the applications of the three-dimensional dynamical systems, we investigate several special cases and give the exactly dynamical solutions in detail. In the end of this paper, we argue that, it is more convenient and also has more physical meaning to express the differential equations of dynamical systems in $(w_{\phi}, \Omega_{\phi}, \lambda)$ instead of variables $(x, y, \lambda)$ and to investigate the dynamical system in 3-Dimension instead of 2-Dimension. We also raise a question about the possibility of the chaotic behavior in the spatially flat single scalar field FRW cosmological models in the presence of ordinary matter.

In this Letter, we report the observational constraints on the Hu-Sawicki $f(R)$ theory derived from weak lensing peak abundances, which are closely related to the mass function of massive halos. In comparison with studies using optical or x-ray clusters of galaxies, weak lensing peak analyses have the advantages of not relying on mass-baryonic observable calibrations. With observations from the Canada-France-Hawaii-Telescope Lensing Survey, our peak analyses give rise to a tight constraint on the model parameter $|f_{R0}|$ for $n=1$. The $95\%$ CL limit is $\log_{10}|f_{R0}| < -4.82$ given WMAP9 priors on $(\Omega_{\rm m}, A_{\rm s})$. With Planck15 priors, the corresponding result is $\log_{10}|f_{R0}| < -5.16$.

Two of the most commonly used tools to constrain the primordial non-Gaussianity are the bispectrum and the Minkowski functionals of CMB temperature anisotropies. These two measures of non-Gaussianity in principle provide distinct (though correlated) information, but in the past constraints from them have only been loosely compared, and not statistically combined. In this work we evaluate, for the first time, the covariance matrix between the local non-Gaussianity coefficient fnl estimated through the bispectrum and Minkowski functionals. We find that the estimators are positively correlated, with corerlation coefficient r ~ 0.3. Using the WMAP7 data to combine the two measures and accounting for the point-source systematics, we find the combined constraint fnl=37+/-28, which has a ~20% smaller error than either of the individual constraints.

We give a full investigation on the dynamics of power-law kinetic quintessence $L(X, \phi)=V(\phi)(-X+X^2)$ by considering the potential related parameter $\Gamma$($=\frac{V V''}{V'^2}$) as a function of another potential parameter $\lambda$($=\frac{V'}{\kappa V^{3/2}}$), which correspondingly extends the analysis of the dynamical system of our universe from two-dimension to three-dimension. Beside the critical points found in previous papers, we find a new de-Sitter-like dominant attractor(\textit{cp$6$}) and give its stable condition using the center manifold theorem. For the dark energy dominant solution(\textit{cp$6$} and \textit{cp$7$}), it could be distinguished from canonical quintessence and tachyon models since the sound speed $c_s^2=0$ or $c_s^2\ll 1$. For the scaling solution (\textit{cp$8$}), it is very interesting that the sound speed $c_s^2=1/5$ while it behaves as ordinary matter. We therefore point out that the power-law kinetic quintessence should have different signatures on cold dark matter power spectrum and cosmic microwave background both at early time when this scalar field is an early dark energy with $\Omega_\phi$ being non-negligible at high redshift and at late time when it drives the accelerating expansion. We even do not know whether there are any degeneracies of the impacts between these two epoches. They are expected to be investigated in future.

Imperfect photometric calibration of galaxy surveys due to either astrophysical or instrumental effects leads to biases in measuring galaxy clustering and in the resulting cosmological parameter measurements. More interestingly (and disturbingly), the spatially varying calibration also generically leads to violations of statistical isotropy of the galaxy clustering signal. Here we develop, for the first time, a formalism to propagate the effects of photometric calibration variations with arbitrary spatial dependence across the sky to the observed power spectra and to the cosmological parameter constraints. We develop an end-to-end pipeline to study the effects of calibration, and illustrate our results using specific examples including Galactic dust extinction and survey-dependent magnitude limits as a function of zenith angle of the telescope. We establish requirements on the control of calibration so that it doesn't significantly bias constraints on dark energy and primordial non-Gaussianity. Two principal findings are: 1) largest-angle photometric calibration variations (dipole, quadrupole and a few more modes, though not the monopole) are the most damaging, and 2) calibration will need to be understood at the 0.1%-1% level (i.e. rms variations mapped out to accuracy between 0.001 and 0.01 mag), though the precise requirement strongly depends on the faint-end slope of the luminosity function and the redshift distribution of galaxies in the survey.

We present the analytical formulas, derived based on the halo model, to compute the cross-correlation between the thermal Sunyaev-Zel'dovich (SZ) effect and the distribution of galaxy clusters. By binning the clusters according to their redshifts and masses, this cross-correlation, the so-called stacked SZ signal, reveals the average SZ profile around the clusters. The stacked SZ signal is obtainable from a joint analysis of an arcminute-resolution cosmic microwave background (CMB) experiment and an overlapping optical survey, which allows for detection of the SZ signals for clusters whose masses are below the individual cluster detection threshold. We derive the error covariance matrix for measuring the stacked SZ signal, and then forecast for its detection from ongoing and forthcoming combined CMB-optical surveys. We find that, over a wide range of mass and redshift, the stacked SZ signal can be detected with a significant signal to noise ratio (total S/N \gsim 10), whose value peaks for the clusters with intermediate masses and redshifts. Our calculation also shows that the stacking method allows for probing the clusters' SZ profiles over a wide range of scales, even out to projected radii as large as the virial radius, thereby providing a promising way to study gas physics at the outskirts of galaxy clusters.

Similar to the magnification of the galaxies' fluxes by gravitational lensing, the extinction of the fluxes by comic dust, whose existence is recently detected by Menard et al (2009), also modify the distribution of a flux-selected galaxy sample. We study the anisotropic distortion by dust extinction to the 3D galaxy correlation function, including magnification bias and redshift distortion at the same time. We find the extinction distortion is most significant along the line of sight and at large separations, similar to that by magnification bias. The correction from dust extinction is negative except at sufficiently large transverse separations, which is almost always opposite to that from magnification bias (we consider a number count slope s > 0.4). Hence, the distortions from these two effects tend to reduce each other. At low z (~<1), the distortion by extinction is stronger than that by magnification bias, but at high z, the reverse holds. We also study how dust extinction affects probes in real space of the baryon acoustic oscillations (BAO) and the linear redshift distortion parameter beta. We find its effect on BAO is negligible. However, it introduces a positive scale-dependent correction to beta that can be as large as a few percent. At the same time, we also find a negative scale-dependent correction from magnification bias, which is up to percent level at low z, but to ~40% at high z. These corrections are non-negligible for precision cosmology, and should be considered when testing General Relativity through the scale-dependence of beta.

Suppression in primordial power on the Universe's largest observable scales has been invoked as a possible explanation for large-angle observations in the cosmic microwave background, and is allowed or predicted by some inflationary models. Here we investigate the extent to which such a suppression could be confirmed by the upcoming large-volume redshift surveys. For definiteness, we study a simple parametric model of suppression that improves the fit of the vanilla LCDM model to the angular correlation function measured by WMAP in cut-sky maps, and at the same time improves the fit to the angular power spectrum inferred from the maximum-likelihood analysis presented by the WMAP team. We find that the missing power at large scales, favored by WMAP observations within the context of this model, will be difficult but not impossible to rule out with a galaxy redshift survey with large volume (~100 Gpc^3). A key requirement for success in ruling out power suppression will be having redshifts of most galaxies detected in the imaging survey.

We investigate the cosmological evolution of the tachyon and phantom-tachyon scalar field by considering the potential parameter $\Gamma$($=\frac{V V"}{V'^2}$) as a function of another potential parameter $\lambda$($=\frac{V'}{\kappa V^{3/2}}$), which correspondingly extends the analysis of the evolution of our universe from two-dimensional autonomous dynamical system to the three-dimension. It allows us to investigate the more general situation where the potential is not restricted to inverse square potential and .One result is that, apart from the inverse square potential, there are a large number of potentials which can give the scaling and dominant solution when the function $\Gamma(\lambda)$ equals $3/2$ for one or some values of $\lambda_{*}$ as well as the parameter $\lambda_{*}$ satisfies condition Eq.(18) or Eq.(19). We also find that for a class of different potentials the dynamics evolution of the universe are actually the same and therefore undistinguishable.

We perform a Markov Chain Monte Carlo analysis of the self-accelerating and normal branch of Dvali-Gabadadze-Porrati braneworld gravity. By adopting a parameterized post-Friedmann description of gravity, we utilize all of the cosmic microwave background data, including the largest scales, and its correlation with galaxies in addition to the geometrical constraints from supernovae distances and the Hubble constant. We find that on both branches brane tension or a cosmological constant is required at high significance with no evidence for the unique Dvali-Gabadadze-Porrati modifications. The cross-over scale must therefore be substantially greater than the Hubble scale H_0 r_c > 3 and 3.5 at the 95% CL with and without uncertainties from spatial curvature. With spatial curvature, the limit from the normal branch is substantially assisted by the galaxy cross-correlation which highlights its importance in constraining infrared modifications to gravity.

By considering the potential parameter $\Gamma$ as a function of another potential parameter $\lambda$[47], We successfully extend the analysis of two-dimensional autonomous dynamical system of quintessence scalar field model to the analysis of three-dimension, which makes us be able to research the critical points of a large number of potentials beyond the exponential potential exactly. We find that there are ten critical points in all, three points $P_{3, 5, 6}$} are general points which are possessed by all quintessence models regardless of the form of potentials and the rest points are closely connected to the concrete potentials. It is quite surprising that, apart from the exponential potential, there are a large number of potentials which can give the scaling solution when the function $f(\lambda)(=\Gamma(\lambda)-1)$ equals zero for one or some values of $\lambda_{*}$ and if the parameter $\lambda_{*}$ also satisfies the condition Eq.(16) or Eq.(17) at the same time. We give the differential equations to derive these potentials $V(\phi)$ from $f(\lambda)$. We also find that, if some conditions are satisfied, the de-Sitter-like dominant point $P_4$ and the scaling solution point $P_9$(or $P_{10}$) can be stable simultaneously but $P_9$ and $P_{10}$ can not be stable simultaneity. Although we survey scaling solutions beyond the exponential potential for ordinary quintessence models in standard general relativity, this method can be applied to other extensively scaling solution models studied in literature[46] including coupled quintessence, (coupled-)phantom scalar field, k-essence and even beyond the general relativity case $H^2 \propto\rho_T^n$. we also discuss the disadvantage of our approach.

Dark energy models with a single scalar field cannot cross the equation of state divide set by a cosmological constant. More general models that allow crossing require additional degrees of freedom to ensure gravitational stability. We show that a parameterized post-Friedmann description of cosmic acceleration provides a simple but accurate description of multiple scalar field crossing models. Moreover the prescription provides a well controlled approximation for a wide range of "smooth" dark energy models. It conserves energy and momentum and is exact in the metric evolution on scales well above and below the transition scale to relative smoothness. Standard linear perturbation tools have been altered to include this description and made publicly available for studies of the dark energy involving cosmological structure out to the horizon scale.

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

Non-gravitational processes, such as feedback from galaxies and their active nuclei, are believed to have injected excess entropy into the intracluster gas, and therefore to have modified the density profiles in galaxy clusters during their formation. Here we study a simple model for this so-called preheating scenario, and ask (i) whether it can simultaneously explain both global X-ray scaling relations and number counts of galaxy clusters, and (ii) whether the amount of entropy required evolves with redshift. We adopt a baseline entropy profile that fits recent hydrodynamic simulations, modify the hydrostatic equilibrium condition for the gas by including approx. 20% non-thermal pressure support, and add an entropy floor K_0 that is allowed to vary with redshift. We find that the observed luminosity-temperature (L-T) relations of low-redshift (z=0.05) HIFLUGCS clusters and high-redshift (z=0.8) WARPS clusters are best simultaneously reproduced with an evolving entropy floor of K_0(z)=341(1+z)^{-0.83}h^{-1/3} keV cm^2. If we restrict our analysis to the subset of bright (kT > 3 keV) clusters, we find that the evolving entropy floor can mimic a self-similar evolution in the L-T scaling relation. This degeneracy with self-similar evolution is, however, broken when (0.5 < kT < 3 keV) clusters are also included. The approx. 60% entropy increase we find from z=0.8 to z=0.05 is roughly consistent with that expected if the heating is provided by the evolving global quasar population. Using the cosmological parameters from the WMAP 3-year data with sigma_8=0.76, our best-fit model underpredicts the number counts of the X-ray galaxy clusters compared to those derived from the 158 deg^2 ROSAT PSPC survey. Treating sigma_8 as a free parameter, we find a best-fit value of sigma_8=0.80+/- 0.02.

The cosmological evolution in Nonlinear Born-Infeld(hereafter NLBI) scalar field theory with negative potentials was investigated. The cosmological solutions in some important evolutive epoches were obtained. The different evolutional behaviors between NLBI and linear(canonical) scalar field theory have been presented. A notable characteristic is that NLBI scalar field behaves as ordinary matter nearly the singularity while the linear scalar field behaves as "stiff" matter. We find that in order to accommodate current observational accelerating expanding universe the value of potential parameters $|m|$ and $|V_0|$ must have an {\it upper bound}. We compare different cosmological evolutions for different potential parameters $m, V_0$.

We present the scalar-tensor gravitational theory with an exponential potential in which pauli metric is regarded as the physical space-time metric. We show that it is essentially equivalent to coupled quintessence(CQ) model. However for baryotropic fluid being radiation there are in fact no coupling between dilatonic scalar field and radiation. We present the critical points for baryotropic fluid and investigate the properties of critical points when the baryotropic matter is specified to ordinary matter. It is possible for all the critical points to be attractors as long as the parameters $\lambda$ and $\beta$ satisfy certain conditions. To demonstrate the attractor behaviors of these critical points, We numerically plot the phase plane for each critical point. Finally with the bound on $\beta$ from the observation and the fact that our universe is undergoing an accelerating expansion, we conclude that present accelerating expansion is not the eventual stage of universe. Moreover, we numerically describe the evolution of the density parameters $\Omega$ and the decelerating factor $q$, and computer the present values of some cosmological parameters, which are consistent with current observational data.

We generally investigate the scalar field model with the lagrangian $L=F(X)-V(\phi)$, which we call it {\it General Non-Canonical Scalar Field Model}. We find that it is a special square potential(with a negative minimum) that drives the linear field solution($\phi=\phi_0t$) while in K-essence model(with the lagrangian $L=-V(\phi)F(X)$) the potential should be taken as an inverse square form. Hence their cosmological evolution are totally different. We further find that this linear field solutions are highly degenerate, and their cosmological evolutions are actually equivalent to the divergent model where its sound speed diverges. We also study the stability of the linear field solution. With a simple form of $F(X)=1-\sqrt{1-2X}$ we indicate that our model may be considered as a unified model of dark matter and dark energy. Finally we study the case when the baryotropic index $\gamma$ is constant. It shows that, unlike the K-essence, the detailed form of F(X) depends on the potential $V(\phi)$. We analyze the stability of this constant $\gamma_0$ solution and find that they are stable for $\gamma_0\leq1$. Finally we simply consider the constant c_s^2 case and get an exact solution for F(X)

We study the potential of a large future weak-lensing survey to constrain dark energy properties by using both the number counts of detected galaxy clusters (sensitive primarily to density fluctuations on small scales) and tomographic shear-shear correlations (restricted to large scales). We use the Fisher matrix formalism, assume a flat universe and parameterize the equation of state of dark energy by w(a)=w_0+w_a(1-a), to forecast the expected statistical errors from either observable, and from their combination. We show that the covariance between these two observables is small, and argue that therefore they can be regarded as independent constraints. We find that when the number counts and the shear-shear correlations (on angular scales l < 1000) are combined, an LSST (Large Synoptic Survey Telescope)-like survey can yield statistical errors on (Omega_DE, w_0, w_a) as tight as (0.003, 0.03, 0.1). These values are a factor of 2-25 better than using either observable alone. The results are also about a factor of two better than those from combining number counts of galaxy clusters and their power spectrum.

In this paper, we regard dilaton in Weyl-scaled induced gravitational theory as a coupled quintessence. Based on this consideration, we investigate the dilaton coupled quintessence(DCQ) model in $\omega-\omega'$ plane, which is defined by the equation of state parameter for the dark energy and its derivative with respect to $N$(the logarithm of the scale factor $a$). We find the scalar field equation of motion in $\omega-\omega'$ plane, and show mathematically the property of attractor solutions which correspond to $\omega_\sigma\sim-1$, $\Omega_\sigma=1$. Finally, we find that our model is a tracking one which belongs to "freezing" type model classified in $\omega-\omega'$ plane.

Based on dilatonic dark energy model, we consider two cases: dilaton field with positive kinetic energy(coupled quintessence) and with negative kinetic energy(phantom). In the two cases, we investigate the existence of attractor solutions which correspond to an equation of state parameter $\omega=-1$ and a cosmic density parameter $\Omega_\sigma=1$. We find that the coupled term between matter and dilaton can't affect the existence of attractor solutions. In the Mexican hat potential, the attractor behaviors, the evolution of state parameter $\omega$ and cosmic density parameter $\Omega$, are shown mathematically. Finally, we show the effect of coupling term on the evolution of $X(\frac{\sigma}{\sigma_0})$ and $Y(\frac{\dot{\sigma}}{\sigma^2_0})$ with respect to $N(lna)$ numerically.

We study the Non-Linear Born-Infeld(NLBI) scalar field model and quintessence model with two different potentials($V(\phi)=-s\phi$ and ${1/2}m^2\phi^2$). We investigate the differences between those two models. We explore the equation of state parameter w and the evolution of scale factor $a(t)$ in both NLBI scalar field and quintessence model. The present age of universe and the transition redshift are also obtained. We use the Gold dataset of 157 SN-Ia to constrain the parameters of the two models. All the results show that NLBI model is slightly superior to quintessence model.

We considered the phantom cosmology with a lagrangian $\displaystyle L=\frac{1}{\eta}[1-\sqrt{1+\eta g^{\mu\nu}\phi_{, \mu}\phi_{, \nu}}]-u(\phi)$, which is original from the nonlinear Born-Infeld type scalar field with the lagrangian $\displaystyle L=\frac{1}{\eta}[1-\sqrt{1-\eta g^{\mu\nu}\phi_{, \mu}\phi_{, \nu}}]-u(\phi)$. This cosmological model can explain the accelerated expansion of the universe with the equation of state parameter $w\leq-1$. We get a sufficient condition for a arbitrary potential to admit a late time attractor solution: the value of potential $u(X_c)$ at the critical point $(X_c,0)$ should be maximum and large than zero. We study a specific potential with the form of $u(\phi)=V_0(1+\frac{\phi}{\phi_0})e^{(-\frac{\phi}{\phi_0})}$ via phase plane analysis and compute the cosmological evolution by numerical analysis in detail. The result shows that the phantom field survive till today (to account for the observed late time accelerated expansion) without interfering with the nucleosynthesis of the standard model(the density parameter $\Omega_{\phi}\simeq10^{-12}$ at the equipartition epoch), and also avoid the future collapse of the universe.