The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

Clustering of the large scale structure provides complementary information to the measurements of the cosmic microwave background anisotropies through power spectrum and bispectrum of density perturbations. Extracting the bispectrum information, however, is more challenging than it is from the power spectrum due to the complex models and the computational cost to measure the signal and its covariance. To overcome these problems, we adopt a proxy statistic, skew spectrum which is a cross-spectrum of the density field and its quadratic field. By applying a large smoothing filter to the density field, we show the theory fits the simulations very well. With the spectra and their full covariance estimated from $N$-body simulations as our "mock" Universe, we perform a global fits for the cosmological parameters. The results show that adding skew spectrum to power spectrum the $1\sigma$ marginalized errors for parameters $ b_1^2A_s, n_s$ and $f_{\rm NL}^{\rm loc}$ are reduced by $31\%, 22\%, 44\%$, respectively. This is the answer to the question posed in the title and indicates that the skew spectrum will be a fast and effective method to access complementary information to that enclosed in the power spectrum measurements, especially for the forthcoming generation of wide-field galaxy surveys.

In this paper, we study the evolution of the ionization fraction $x_e(z)$ during the epoch of reionization by using the dispersion measures (DM) of fast radio bursts (FRBs). Different from the previous studies, here we turn to consider the large-scale clustering information of observed DM of FRB catalog, which only needs the rough redshift distribution, instead of the exact redshift information of each FRB. Firstly, we consider the instantaneous "\texttt{tanh}" model for $x_e(z)$ and find that including the auto-correlation information of the mock catalog, about $10^4$ FRBs with the intrinsic DM scatter of 100 $\rm pc/cm^3$ spanning 20\% of all sky, could significantly improve the constraint on the width $\Delta_z$ of the model, when comparing with that from the CMB data alone. The evolution shape of the ionization fraction will be tightly narrowed, namely the duration of the epoch of reionization has been shrunk, $z_{\rm dur}<1.15$ (68\% C.L.). Furthermore, we also use another redshift-asymmetric reionization model and obtain that the FRB mock catalog could measure the ionization fraction at $z=6$ precisely with the $1\sigma$ error $\Delta x_e(z=6)=0.009$, which means that the large-scale clustering information of observed DM of FRB catalog is very sensitive to the ionization fraction of the end of reionization epoch. We conclude that the observation of high-redshift FRBs could be a complementary probe to study the reionization history in the future.

We propose a method based on the process of extracting gravitational wave (GW) parameters from GW signals to simulate the binary neutron-star (BNS) merging events. We simulate 1000 GW standard sirens based on the observation of the Einstein Telescope (ET). Almost all the simulated GW data are in the redshift range of $[0,3]$. The role of the GW standard siren in the inference of the cosmological parameters is investigated. We find that the GW data can help improve the accuracy of cosmological parameters. Moreover, the degeneracy of cosmological parameters is broken by the GW data. The GW standard siren is helpful for the constraint of the cosmological parameters.

In this work, we present a new method to estimate cosmological parameters accurately based on Artificial Neural Network (ANN), and a code called ECoPANN (Estimating Cosmological Parameters with ANN) is developed to achieve parameter inference. We test the ANN method by estimating the basic parameters of the concordance cosmological model using the simulated temperature power spectrum of the cosmic microwave background (CMB). The results show that the ANN performs excellently on best-fit values and errors of parameters, as well as correlations between parameters when compared with that of the Markov chain Monte Carlo (MCMC) method. Besides, for a well trained ANN model, it is capable of estimating parameters for multiple experiments that have different precisions, which can greatly reduce the consumption of time and computing resources for parameter inference. Moreover, the well trained ANN is capable of discovering new potential physics in both the current and future higher precision observations. In addition, we extend the ANN to a multi-branch network to achieve joint constraint on parameters. We test the multi-branch network using the simulated temperature and polarization power spectra of CMB, type Ia supernovae, and baryon acoustic oscillation, and almost obtain the same results as the MCMC method. Therefore, we propose that the ANN can provide an alternative way to accurate and fast estimate cosmological parameters, and ECoPANN can be applied to the research of cosmology and even other broader scientific fields.

In this paper, we estimate the cross-correlation power spectra between the Planck 2018 cosmic microwave background (CMB) temperature anisotropy map and the unresolved $\gamma$-ray background (UGRB) from the 9-years Fermi-Large Area Telescope (LAT) data. In this analysis, we use up to nine energy bins over a wide energy range of [0.631, 1000] GeV from the Fermi-LAT UGRB data. Firstly, we find that the Fermi data with the energy ranges [1.202, 2.290] GeV and [17.38, 36.31] GeV show the positive evidence for the Integrated Sachs-Wolfe (ISW) effect at $1.8\sigma$ confidence level, and the significance would be increased to $2.7\sigma$ when using these two energy bins together. Secondly, we apply the single power-law model to normalize the amplitude and use all the nine Fermi energy bins to measure the significance of the ISW effect, we obtained ${\rm A_{amp}}=0.95 \pm 0.53$ ($68\%$ C.L.). For the robustness test, we implement a null hypothesis by randomizing the Fermi mock maps of nine energy bins and obtain the non-detection of ISW effect, which confirms that the ISW signal comes from the Fermi-LAT diffuse $\gamma$-ray data and is consistent with the standard $\Lambda$CDM model prediction essentially. We use a cross-correlation coefficient to show the relation between different energy bins. Furthermore, we vary the cut ranges $|b|$ of galactic plane on the mask of Fermi map and carefully check the consequent influence on the ISW signal detection.

In this work, we achieve the determination of the cosmic curvature $\Omega_K$ in a cosmological model-independent way, by using the Hubble parameter measurements $H(z)$ and type Ia supernovae (SNe Ia). In our analysis, two nonlinear interpolating tools are used to reconstruct the Hubble parameter, one is the Artificial Neural Network (ANN) method, and the other is the Gaussian process (GP) method. We find that $\Omega_K$ based on the GP method can be greatly influenced by the prior of $H_0$, while the ANN method can overcome this. Therefore, the ANN method may have more advantages than GP in the measurement of the cosmic curvature. Based on the ANN method, we find a spatially open universe is preferred by the current $H(z)$ and SNe Ia data, and the difference between our result and the value inferred from Planck CMB is $1.6\sigma$. In order to test the reliability of the ANN method, and the potentiality of the future gravitational waves (GW) standard sirens in the measurement of the cosmic curvature, we constrain $\Omega_K$ using the simulated Hubble parameter and GW standard sirens in a model-independent way. We find that the ANN method is reliable and unbiased, and the error of $\Omega_K$ is $\sim0.186$ when 100 GW events with electromagnetic counterparts are detected, which is $\sim56\%$ smaller than that constrained from the Pantheon SNe Ia. Therefore, the data-driven method based on ANN has potential in the measurement of the cosmic curvature.

In this work, we propose a new nonparametric approach for reconstructing a function from observational data using an Artificial Neural Network (ANN), which has no assumptions about the data and is a completely data-driven approach. We test the ANN method by reconstructing functions of the Hubble parameter measurements $H(z)$ and the distance-redshift relation $D_L(z)$ of Type Ia supernovae. We find that both $H(z)$ and $D_L(z)$ can be reconstructed with high accuracy. Furthermore, we estimate cosmological parameters using the reconstructed functions of $H(z)$ and $D_L(z)$ and find the results are consistent with those obtained using the observational data directly. Therefore, we propose that the function reconstructed by ANN can represent the actual distribution of observational data and can be used for parameter estimation in further cosmological research. In addition, we present a new strategy for training and evaluating the neural network, and a code for reconstructing functions using ANN has been developed and is available for download.

There is a persistent $H_0$-tension, now at more than $\gtrsim 4\sigma$ level, between the local distance ladder value and the \emph{Planck} cosmic microwave background measurement, in the context of flat $\Lambda$CDM model. We reconstruct $H(z)$ in a cosmological-model-independent way using three low-redshift distance probes including the latest data from baryon acoustic oscillation, Type Ia supernova and four gravitational lensing Time-Delay observations. We adopt general parametric models of $H(z)$ and assume a Gaussian prior on the sound horizon at drag epoch, $r_{\mathrm s}$, from \emph{Planck} measurement. The reconstructed $H_0$ using Pantheon SN Ia and BAO data are consistent with the \emph{Planck} flat $\Lambda$CDM value. When including the GLTD data, $H_0$ increases mildly, yet remaining discrepant with the local measurement at $\sim 2.5\sigma$ level. Our reconstructions being blind to the dark sectors at low redshift, we reaffirm the earlier claims that the Hubble tension is not likely to be solved by modifying the energy budget of the low-redshift universe. We further forecast the constraining ability of future realistic mock BAO data from DESI and GLTD data from LSST, combining which, we anticipate that the uncertainty of the inferred $H_0$ would be improved by $\sim 38\%$, reaching $\sigma_{H_0} \approx 0.56$ uncertainty level.

We report the identification of a positive cross-correlation signal between the unresolved $\gamma$-ray emission, measured by the \emph{Fermi} Large Area Telescope, and four different galaxy cluster catalogues. The selected catalogues peak at low-redshift and span different frequency bands, including infrared, optical and X-rays. The signal-to-noise ratio of the detected cross-correlation amounts to 3.5 in the most significant case. We investigate and comment about its possible origin, in terms of compact $\gamma$-ray emission from AGNs inside clusters or diffuse emission from the intracluster medium. The analysis has been performed by introducing an accurate estimation of the cross-correlation power-spectrum covariance matrix, built with mock realisations of the gamma and galaxy cluster maps. Different methods to produce the mock realizations starting from the data maps have been investigated and compared, identifying suitable techniques which can be generalized to other cross-correlation studies.

The presence of a significant amount of gravitational radiation in the early Universe affects the total energy density and hence the expansion rate in the early epoch. In this work, we develop a physical model to connect the parameter of relativistic degree of freedom $N_\mathrm{eff}$ with the amplitude and shape of primordial tensor power spectrum, and use the CMB temperature and polarization data from {\it Planck} and BICEP2/KECK Array, and the primordial deuterium measurement from damped Lyman-$\alpha$ (DLA) systems to constrain this model. We find that with this extra relation $\Delta N_{\rm eff}(r,n_{\rm t})$, the tensor-to-scalar ratio $r$ is constrained to be $r<0.07$ ($3\sigma$ C.L.) and the tilt of tensor power spectrum is $n_\mathrm{t}=-0.01\pm 0.31$ ($1\sigma$ C.L.) for {\it Planck}+BICEP2+KECK+[D/H] data. This achieves a much tighter constraint on the tensor spectrum and provides a stringent test for cosmic inflation models. In addition, the current constraint on $N_{\rm eff}=3.122 \pm 0.171$ excludes the possibility of fourth neutrino species at more than $5\sigma$ C.L.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

In this letter we present constraints on the scale-dependent "local" type primordial non-Gaussianity, which is described by non-Gaussianity's spectral index $n_{\rm NG}$, from the NRAO VLA Sky Survey and the quasar catalog of the Sloan Digital Sky Survey (SDSS) Data Release 6, together with the SDSS Data Release 12 photo-z sample. Here, we use the auto-correlation analyses of these three probes and their cross-correlation analyses with the cosmic microwave background (CMB) temperature map, and obtain the tight constraint on the spectral index: $n_{\rm NG}=0.2 ^{+0.7}_{-1.0}$ ($1\,\sigma$ C.L.), which shows the first competitive constraint on the running of non-Gaussianity from current large-scale structure clustering data. Furthermore, we also perform the forecast calculations and improve the limit of $n_{\rm NG}$ using the future Euclid mission, and obtain the standard deviation at 68\% confidence level: $\Delta n_{\rm NG}=1.74$ when considering the fiducial value $f_{\rm NL}=3$, which provides the complementary constraining power to those from the CMB bispectrum information.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

Dark energy might interact with dark matter in a direct, non-gravitational way, which can help remedy several theoretical defects. In order to find out the properties of interacting dark energy models, it is necessary to investigate the cosmological perturbations in detail. However, due to the improper use of pressure perturbation of dark energy, a large-scale instability at the early stage occurs occasionally. In recent years, parametrized post-Friedmann approach has been used to calculate the perturbation equations in the interacting dark energy scenario. Under this framework, the dark energy pressure perturbation was replaced by the relationship between the momentum density of dark energy and other components on a large scale. However, this paper shows that if the interaction terms are related to the velocity perturbation of dark energy, the density perturbation of dark matter and the matter power spectrum will diverge when the equation of state parameter of dark energy $w$ is close to $-1$. A simple parameterization steering clear of this problem is proposed in this paper which is a more general form and can be applied to explore the interaction between dark matter and dark energy by using various cosmological data.

We present mosaicked long-slit spectral maps of 18 nearby Active Galactic Nuclei (AGNs), 2 LINERs, and 4 star-forming galaxies. With the resulting data cubes taken using the Kast dual spectrograph on the 3 m Shane telescope of the Lick Observatory, we measure the aperture effects on the spectroscopic classification of AGNs. With more starlight included in a larger aperture, the nuclear spectrum that is Seyfert-like may become contaminated. We generated standard spectroscopic classification diagrams in different observing apertures. These show quantitatively how the ensemble of Seyferts migrates toward the H $\scriptsize{\textrm{II}}$ region classification when being observed with increasing aperture sizes. But the effect ranges widely in individual active galaxies. Some of the less luminous Seyferts shfit by a large amount, while some other barely move or even shift in different directions. We find that those Seyfert galaxies with the fraction of nuclear H$\alpha$ emission lower than 0.2 of the host galaxy, 2-10 keV hard X-ray luminosity lower than $10^{43}$ erg s$^{-1}$, and the observed nuclear [O $\scriptsize{\textrm{III}}$] luminosity lower than $10^{40.5}$ erg s$^{-1}$, are more likely to change activity classification type when the entire host galaxy is included. Overall, 4 of our 24 galaxies (18 Seyferts) change their spectral activity classification type when observed with a very large aperture.

In this paper, we propose to estimate the spatial curvature of the universe and the cosmic opacity in a model-independent way with expansion rate measurements, $H(z)$, and type Ia supernova (SNe Ia). On the one hand, using a nonparametric smoothing method Gaussian process, we reconstruct a function $H(z)$ from opacity-free expansion rate measurements. Then, we integrate the $H(z)$ to obtain distance modulus $\mu_{\rm H}$, which is dependent on the cosmic curvature. On the other hand, distances of SNe Ia can be determined by their photometric observations and thus are opacity-dependent. In our analysis, by confronting distance moduli $\mu_{\rm H}$ with those obtained from SNe Ia, we achieve estimations for both the spatial curvature and the cosmic opacity without any assumptions for the cosmological model. Here, it should be noted that light curve fitting parameters, accounting for the distance estimation of SNe Ia, are determined in a global fit together with the cosmic opacity and spatial curvature to get rid of the dependence of these parameters on cosmology. In addition, we also investigate whether the inclusion of different priors for the present expansion rate ($H_0$: global estimation, $67.74\pm 0.46~\rm km~ s^{-1} ~Mpc^{-1}$, and local measurement, $73.24\pm 1.74~\rm km~ s^{-1} ~Mpc^{-1}$) exert influence on the reconstructed $H(z)$ and the following estimations of the spatial curvature and cosmic opacity. Results show that, in general, a spatially flat and transparent universe is preferred by the observations. Moreover, it is suggested that priors for $H_0$ matter a lot. Finally, we find that there is a strong degeneracy between the curvature and the opacity.

We investigate the nature of the extragalactic unresolved gamma-ray background (UGRB) by cross-correlating several galaxy catalogs with sky-maps of the UGRB built from 78 months of Pass 8 Fermi-Large Area Telescope data. This study updates and improves similar previous analyses in several aspects. Firstly, the use of a larger gamma-ray dataset allows us to investigate the energy dependence of the cross-correlation in more detail, using up to 8 energy bins over a wide energy range of [0.25-500] GeV. Secondly, we consider larger and deeper catalogs (2MASS Photometric-Redshift catalog, 2MPZ; WISE x SuperCOSMOS, WISC; and SDSS-DR12 photometric-redshift dataset) in addition to the ones employed in the previous studies (NVSS and SDSS-QSOs). Thirdly, we exploit the redshift information available for the above catalogs to divide them into redshift bins and perform the cross-correlation separately in each of them. Our results confirm, with higher statistical significance, the detection of cross-correlation signal between the UGRB maps and all the catalogs considered, on angular scales smaller than 1 degree. Significances range from 16.3 sigma for NVSS, 7 sigma for SDSS-DR12 and WISC, 5 sigma for 2MPZ and 4 sigma for SDSS-QSOs. Furthermore, including redshift tomography, the significance of the SDSS-DR12 signal strikingly rises up to 12 sigma and the one of WISC to 10.6 sigma. We offer a simple interpretation of the signal in the framework of the halo model. The precise redshift and energy information allows us to clearly detect a change over redshift in the spectral and clustering behavior of the gamma-ray sources contributing to the UGRB.

In the present paper, we investigate the cosmographic problem using the bias-variance trade-off. We find that both the z-redshift and the $y=z/(1+z)$-redshift can present a small bias estimation. It means that the cosmography can describe the supernova data more accurately. Minimizing risk, it suggests that cosmography up to the second order is the best approximation. Forecasting the constraint from future measurements, we find that future supernova and redshift drift can significantly improve the constraint, thus having the potential to solve the cosmographic problem. We also exploit the values of cosmography on the deceleration parameter and equation of state of dark energy $w(z)$. We find that supernova cosmography cannot give stable estimations on them. However, much useful information was obtained, such as that the cosmography favors a complicated dark energy with varying $w(z)$, and the derivative $dw/dz<0$ for low redshift. The cosmography is helpful to model the dark energy.

Luminosity correlations of long Gamma-ray bursts (GRB) are extensively proposed as an effective complementarity to trace the Hubble diagram of Universe at high redshifts, which is of great importance to explore properties of dark energy. Recently, several empirical luminosity correlations have been statistically proposed from GRB observations. However, to treat GRB as the distance indicator, there are two key issues: the redshift evolution of luminosity correlations and their calibrations. In this paper, we choose the Amati relation, the correlation between the peak spectra energy and the equivalent isotropic energy of GRBs ($E_{\rm p}-E_{\rm iso}$), as an example, and find that the current GRB dataset implies that there could be a evolution of the luminosity correlation with respect to the redshift. Therefore, we propose an extended Amati relation with two extra redshift-dependent terms to correct the redshift evolution of GRB relation. Secondly, we carefully check the reliability of the calibration method using the low-redshift GRB data. Importantly, we find that the low-redshift calibration method does not take whole correlations between $\Omega_{\rm m}$ and coefficients into account. Neglecting these correlation information can break the degeneracies and obtain the biased constraint on $\Omega_{\rm m}$ which is very sensitive to values of parameters for the calibration. A small shift of parameters of "calibrated" relation could significantly change the final constraint on $\Omega_{\rm m}$ in the low-redshift calibration method. Finally, we simulate several GRB samples with different statistical errors and find that, in order to correctly recover the fiducial value of $\Omega_{\rm m}$ using the low-redshift calibration method, we need a large number of GRB samples with high precisions.

The possible slowing down of cosmic acceleration was widely studied. However, the imposition of dark energy parametrization brought some tensions. In our recent paper, we test this possibility using a model-independent method, Gaussian processes. However, the reason of generating these tensions is still closed. In the present paper, we analyse the derivative of deceleration parameter to solve the problems. The reconstruction of the derivative again suggests that no slowing down of acceleration is presented within 95\% C.L. from current observational data. We then deduce its constraint on dark energy. The corresponding constraint clearly reveals the reason of tension between different models in previous work. We also study the essential reason of why current data cannot convincingly measure the slowing down of acceleration. The constraints indicate that most of current data are not in the allowed region.

We investigate the primordial gravitational waves (PGWs) in the general scenario where the inflation is preceded by a pre-inflationary stage with the effective equation of state $w$. Comparing with the results in the usual inflationary models, the power spectrum of PGWs is modified in two aspects: One is the mixture of the perturbation modes caused by he presence of the pre-inflationary period, and the other is the thermal initial state formed at the Planck era of the early Universe. By investigating the observational imprints of these modifications on the B-mode polarization of cosmic microwave background (CMB) radiation, we obtain the constraints on the conformal temperature of the thermal gravitational-wave background $T<5.01\times 10^{-4}$Mpc$^{-1}$ and a tensor-to-scalar ratio $r<0.084$ ($95\%$ confident level), which follows the bounds on total number of e-folds $N>63.5$ for the model with $w=1/3$, and $N>65.7$ for that with $w=1$. By taking into account various noises and the foreground radiations, we forecast the detection possibility of the thermal gravitational-wave background by the future CMBPol mission, and find that if $r>0.01$, the detection is possible as long as $T>1.5\times 10^{-4}$Mpc$^{-1}$. However, the effect of different $w$ is quite small, and it seems impossible to determine its value from the potential observations of CMBPol mission.

In this paper we use a recently compiled data set, which comprises 118 galactic-scale strong gravitational lensing (SGL) systems to constrain the statistic property of SGL system, as well as the curvature of universe without assuming any fiducial cosmological model. Based on the singular isothermal ellipsoid (SIE) model of SGL system, we obtain that the constrained curvature parameter $\Omega_{\rm k}$ is close to zero from the SGL data, which is consistent with the latest result of planck measurement. More interestingly, we find that the parameter $f$ in the SIE model is strongly correlated with the curvature $\Omega_{\rm k}$. Neglecting this correlation in the analysis will significantly overestimate the constraining power of SGL data on the curvature. Furthermore, the obtained constraint on $f$ is different from previous results: $f=1.105\pm0.030$ ($68\%$ C.L.), which means that the standard singular isothermal sphere (SIS) model ($f=1$) is disfavored by the current SGL data at more than $3\sigma$ confidence level. We also divide the whole SGL data into two parts according to the centric stellar velocity dispersion $\sigma_{\rm c}$ and find that the larger value of $\sigma_{\rm c}$ the subsample has, the more favored the standard SIS model is. Finally, we extend the SIE model by assuming the power-law density profiles for the total mass density, $\rho=\rho_0(r/r_0)^{-\alpha}$, and luminosity density, $\nu=\nu_0(r/r_0)^{-\delta}$, and obtain the constraints on the power-law indexes: $\alpha=1.95\pm0.04$ and $\delta=2.40\pm0.13$ at 68\% confidence level. When assuming the power-law index $\alpha=\delta=\gamma$, this scenario is totally disfavored by the current SGL data, $\chi^2_{\rm min,\gamma} - \chi^2_{\rm min,SIE} \simeq 53$.

We report the detection of a cross-correlation signal between {\it Fermi} Large Area Telescope diffuse gamma-ray maps and catalogs of clusters. In our analysis, we considered three different catalogs: WHL12, redMaPPer and PlanckSZ. They all show a positive correlation with different amplitudes, related to the average mass of the objects in each catalog, which also sets the catalog bias. The signal detection is confirmed by the results of a stacking analysis. The cross-correlation signal extends to rather large angular scales, around 1 degree, that correspond, at the typical redshift of the clusters in these catalogs, to a few to tens of Mpc, i.e. the typical scale-length of the large scale structures in the Universe. Most likely this signal is contributed by the cumulative emission from AGNs associated to the filamentary structures that converge toward the high peaks of the matter density field in which galaxy clusters reside. In addition, our analysis reveals the presence of a second component, more compact in size and compatible with a point-like emission from within individual clusters. At present, we cannot distinguish between the two most likely interpretations for such a signal, i.e. whether it is produced by AGNs inside clusters or if it is a diffuse gamma-ray emission from the intra-cluster medium. We argue that this latter, intriguing, hypothesis might be tested by applying this technique to a low redshift large mass cluster sample.

Much focus was on the possible slowing down of cosmic acceleration under the dark energy parametrization. In the present paper, we investigate this subject using the Gaussian processes (GP), without resorting to a particular template of dark energy. The reconstruction is carried out by abundant data including luminosity distance from Union2, Union2.1 compilation and gamma-ray burst, and dynamical Hubble parameter. It suggests that slowing down of cosmic acceleration cannot be presented within 95\% C.L., in considering the influence of spatial curvature and Hubble constant. In order to reveal the reason of tension between our reconstruction and previous parametrization constraint for Union2 data, we compare them and find that slowing down of acceleration in some parametrization is only a "mirage". Although these parameterizations fits well with the observational data, their tension can be revealed by high order derivative of distance $D$. Instead, GP method is able to faithfully model the cosmic expansion history.

We compare the measured angular cross-correlation between the Fermi-LAT gamma-ray sky and catalogues of extra-galactic objects with the expected signal induced by weakly interacting massive particle (WIMP) dark matter (DM). We include a detailed description of the contribution of astrophysical gamma-ray emitters such as blazars, misaligned AGN and star forming galaxies, and perform a global fit to the measured cross-correlation. Five catalogues are considered: SDSS-DR6 quasars, 2MASS galaxies, NVSS radio galaxies, SDSS-DR8 Luminous Red Galaxies and SDSS-DR8 main galaxy sample. To model the cross-correlation signal we use the halo occupation distribution formalism to estimate the number of galaxies of a given catalogue in DM halos and their spatial correlation properties. We discuss uncertainties in the predicted cross-correlation signal arising from the DM clustering and WIMP microscopic properties, which set the DM gamma-ray emission. The use of different catalogues probing objects at different redshifts reduces significantly, though not completely, the degeneracy among the different gamma-ray components. We find that the presence of a significant WIMP DM signal is allowed by the data but not significantly preferred by the fit, although this is mainly due to a degeneracy with the misaligned AGN component. With modest substructure boost, the sensitivity of this method excludes thermal annihilation cross sections at 95% C.L. for WIMP masses up to few tens of GeV. Constraining the low-redshift properties of astrophysical populations with future data will further improve the sensitivity to DM.

In this paper, we use the model dependent method to revisit the constraint on the well-known cosmic distance duality relation (CDDR). By using the latest SNIa samples, such as Union2.1, JLA and SNLS, we find that the SNIa data alone can not constrain the cosmic opacity parameter $\varepsilon$, which denotes the deviation from the CDDR, $d_{\rm L} = d_{\rm A}(1+z)^{2+\varepsilon}$, very well. The constraining power on $\varepsilon$ from the luminosity distance indicator provided by SNIa and GRB is hardly to be improved at present. When we include other cosmological observations, such as the measurements of Hubble parameter, the baryon acoustic oscillations and the distance information from cosmic microwave background, we obtain the tightest constraint on the cosmic opacity parameter $\varepsilon$, namely the 68\% C.L. limit: $\varepsilon=0.023\pm0.018$. Furthermore, we also consider the evolution of $\varepsilon$ as a function of $z$ using two methods, the parametrization and the principle component analysis, and do not find the evidence for the deviation from zero. Finally, we simulate the future SNIa and Hubble measurements and find the mock data could give very tight constraint on the cosmic opacity $\varepsilon$ and verify the CDDR at high significance.

We develop an efficient method based on the linear regression algorithm to probe the cosmological CPT violation using the CMB polarisation data. We validate this method using simulated CMB data and apply it to recent CMB observations. We find that a combined data sample of BICEP1 and BOOMERanG 2003 favours a nonzero isotropic rotation angle at $2.3\sigma$ confidence level, ie, $\Delta\alpha=-3.3 \pm1.4$ deg (68% CL) with systematics included.

Searching for the signal of primordial gravitational waves in the B-modes (BB) power spectrum is one of the key scientific aims of the cosmic microwave background (CMB) polarization experiments. However, this could be easily contaminated by several foreground issues, such as the thermal dust emission. In this paper we study another mechanism, the cosmic birefringence, which can be introduced by a CPT-violating interaction between CMB photons and an external scalar field. Such kind of interaction could give rise to the rotation of the linear polarization state of CMB photons, and consequently induce the CMB BB power spectrum, which could mimic the signal of primordial gravitational waves at large scales. With the recent polarization data of BICEP2 and the joint analysis data of BICEP2/Keck Array and Planck, we perform a global fitting analysis on constraining the tensor-to-scalar ratio $r$ by considering the polarization rotation angle which can be separated into a background isotropic part and a small anisotropic part. Since the data of BICEP2 and Keck Array experiments have already been corrected by using the "self-calibration" method, here we mainly focus on the effects from the anisotropies of CMB polarization rotation angle. We find that including the anisotropies in the analysis could slightly weaken the constraints on $r$, when using current CMB polarization measurements. We also simulate the mock CMB data with the BICEP3-like sensitivity. Very interestingly, we find that if the effects of the anisotropic polarization rotation angle can not be taken into account properly in the analysis, the constraints on $r$ will be dramatically biased. This implies that we need to break the degeneracy between the anisotropies of the CMB polarization rotation angle and the CMB primordial tensor perturbations, in order to measure the signal of primordial gravitational waves accurately.

If dark matter (DM) is composed by particles which are non-gravitationally coupled to ordinary matter, their annihilations or decays in cosmic structures can result in detectable radiation. We show that the most powerful technique to detect a particle DM signal outside the Local Group is to study the angular cross-correlation of non-gravitational signals with low-redshift gravitational probes. This method allows to enhance signal-to-noise from the regions of the Universe where the DM-induced emission is preferentially generated. We demonstrate the power of this approach by focusing on GeV-TeV DM and on the recent cross-correlation analysis between the 2MASS galaxy catalogue and the Fermi-LAT gamma-ray maps. We show that this technique is more sensitive than other extragalactic gamma-ray probes, such as the energy spectrum and angular autocorrelation of the extragalactic background, and emission from clusters of galaxies. Intriguingly, we find that the measured cross-correlation can be well fitted by a DM component, with thermal annihilation cross section and mass between 10 and 100 GeV, depending on the small-scale DM properties and gamma-ray production mechanism. This solicits further data collection and dedicated analyses.

In this paper, we study the case of single field inflation within the framework of modified gravity theory where the gravity part has an arbitrary form $f(R)$. Via a conformal transformation, this case can be transformed into its Einstein frame where it looks like a two-field inflation model. However, due to the existence of the isocurvature modes in such a multi-degree-of-freedom (m.d.o.f.) system, the (curvature) perturbations are not equivalent in two frames, so in despite of its convenience, it is illegal to treat the perturbations in its Einstein frame as the "real" ones as we always do for pure $f(R)$ theory or single field with nonminimal coupling. Here by pulling the results of curvature perturbations back into its original Jordan frame, we show explicitly the power spectrum and spectral index of the perturbations in the Jordan frame, as well as how it differs from the Einstein frame. We also fit our results with the newest Planck data. Since there are large parameter space in these models, we show that it is easy to fit the data very well.

Building on our previous cross-correlation analysis (Xia et al. 2011) between the isotropic gamma-ray background (IGRB) and different tracers of the large-scale structure of the universe, we update our results using 60-months of data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. We perform a cross-correlation analysis between the IGRB and objects that may trace the astrophysical sources of the IGRB: SDSS-DR6 QSOs, the SDSS-DR8 Main Galaxy Sample, Luminous Red Galaxies (LRGs) in the SDSS catalog, 2MASS galaxies, and radio NVSS galaxies. The benefit of correlating the Fermi-LAT signal with catalogs of objects at various redshifts is to provide tomographic information on the IGRB which is crucial to separate the various contributions and to clarify its origin. We observe a significant (>3.5 sigma) cross-correlation signal on angular scales smaller than 1 deg in the NVSS, 2MASS and QSO cases and, at lower statistical significance (~3.0 sigma), with SDSS galaxies. These results are robust against the choice of the statistical estimator, estimate of errors, map cleaning procedure and instrumental effects. Finally, we test the hypothesis that the IGRB observed by Fermi-LAT originates from the summed contributions of three types of unresolved extragalactic sources: BL Lacs, FSRQs and Star-Forming Galaxies (SFGs). We find that a model in which the IGRB is mainly produced by SFGs ($72^{+23}_{-37}$% with 2 sigma errors), with BL Lacs and FSRQs giving a minor contribution, provides a good fit to the data. We also consider a possible contribution from Misaligned Active Galactic Nuclei, and we find that, depending on the details of the model and its uncertainty, they can also provide a substantial contribution, partly degenerate with the SFG one. (abridged)

In this paper we use the current and future cosmic microwave background (CMB) experiments to test the Charge-Parity-Time Reversal (CPT) symmetry. We consider a CPT-violating interaction in the photon sector $\mathcal{L}_{\rm cs}\sim p_\mu A_\nu \tilde{F}^{\mu\nu}$ which gives rise to a rotation of the polarization vectors of the propagating CMB photons. By combining the nine-year WMAP, BOOMERanG 2003 and BICEP1 observations, we obtain the current constraint on the isotropic rotation angle $\bar{\alpha} = -2.12 \pm 1.14$ ($1\sigma$), indicating an about $2\sigma$ significance of the CPT violation. Here, we particularly take the systematic errors of CMB measurements into account. Then, we study the effects of the anisotropies of the rotation angle [$\Delta{\alpha}({\bf \hat{n}})$] on the CMB polarization power spectra in detail. Due to the small effects, the current CMB polarization data can not constrain the related parameters very well. We obtain the 95\% C.L. upper limit of the variance of the anisotropies of the rotation angle $C^\alpha(0) < 0.035$ from all the CMB datasets. More interestingly, including the anisotropies of rotation angle could lower the best fit value of $r$ and relax the tension on the constraints of $r$ between BICEP2 and Planck. Finally, we investigate the capabilities of future Planck polarization measurements on $\bar{\alpha}$ and $\Delta{\alpha}({\bf \hat{n}})$. Benefited from the high precision of Planck data, the constraints of the rotation angle can be significantly improved.

We investigate the cosmological implications of the recently constructed 5-dimensional braneworld cosmology with gravitating Nambu-Goto matching conditions. Inserting both matter and radiation sectors, we first extract the analytical cosmological solutions. Additionally, we use observational data from Type Ia Supernovae (SNIa) and Baryon Acoustic Oscillations (BAO), along with requirements of Big Bang Nucleosynthesis (BBN), in order to impose constraints on the parameters of the model. We find that the scenario at hand is in very good agreement with observations, and thus a small departure from the standard Randall-Sundrum scenario is allowed.

In this paper, we set the new limits on the equation of state parameter (EoS) of dark energy with the observations of cosmic microwave background radiation (CMB) from Planck satellite, the type Ia supernovae from Pan-STARRS and the baryon acoustic oscillation (BAO). We consider two parametrization forms of EoS: a constant $w$ and time evolving $w(a)=w_0+w_a(1-a)$. The results show that with a constant EoS, $w=-1.141\pm{0.075}$ ($68\%~C.L.$), which is consistent with $\Lambda$CDM at about $2\sigma$ confidence level. For a time evolving $w(a)$ model, we get $w_0=-1.09^{+0.16}_{-0.18}$ ($1\sigma~C.L.$), $w_a=-0.34^{+0.87}_{-0.51}$ ($1\sigma~C.L.$), and in this case $\Lambda$CDM can be comparable with our observational data at $1\sigma$ confidence level. In order to do the parametrization independent analysis, additionally we adopt the so called principal component analysis (PCA) method, in which we divide redshift range into several bins and assume $w$ as a constant in each redshift bin (bin-w). In such bin-w scenario, we find that for most of the bins cosmological constant can be comparable with the data, however, there exists few bins which give $w$ deviating from $\Lambda$CDM at more than $2\sigma$ confidence level, which shows a weak hint for the time evolving behavior of dark energy. To further confirm this hint, we need more data with higher precision.

The BICEP2 collaboration reports a detection of primordial cosmic microwave background (CMB) B-mode with a tensor-scalar ratio $r=0.20^{+0.07}_{-0.05}$ (68% C.L.). However, this result is in tension with the recent Planck limit, $r<0.11$ (95% C.L.), on constraining inflation models. In this Letter we consider an inflationary cosmology with a preceding nonsingular bounce which gives rise to observable signatures on primordial perturbations. One interesting phenomenon is that both the primordial scalar and tensor modes can have a step feature on their power spectra, which nicely cancels the tensor excess power on the CMB temperature power spectrum. By performing a global analysis, we obtain the 68% C.L. constraints on the parameters of the model from the Planck+WP and BICEP2 data together: the jump scale $\log_{10}(k_{\rm b}/{\rm Mpc}^{-1})=-2.4\pm0.2$ and the spectrum amplitude ratio of bounce-to-inflation $r_B\equiv P_{\rm m} / A_{\rm s} = 0.71\pm0.09$. Our result reveals that the bounce inflation scenario can simultaneously explain the Planck and BICEP2 observations better than the standard $\Lambda$CDM model, and can be verified by the future CMB polarization measurements.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

Recently, BICEP2 collaboration has released their results on the measurements of the CMB polarizations. In the framework of the $\Lambda$CDM with a power law form of the scalar primordial power spectrum, this new measurement on the B-mode puts a tight constraint on tensor/scalar ratio $r=0.2^{+0.07}_{-0.05}~(1\sigma)$, which however is in tension with the Planck limit, $r<0.11$ (95% C.L.). In this paper, we consider various extensions of the $\Lambda$CDM model by introducing extra cosmological parameters such as the equation state of dark energy $w$, the curvature of the universe $\Omega_k$, the running of the scalar power spectrum index $\alpha_s$, the sum of the neutrino mass $\Sigma m_{\nu}$, the effective number of neutrinos $N_{eff}$, the tensor power spectrum index $n_{t}$, then perform the global fit to the BICEP2, Planck and also the SN and BAO data. Our results show the cosmological parameters $\alpha_s$ and $n_t$ are highly and directly correlated with the tensor/scalar ratio $r$. However indirectly the parameters $\Omega_k$ and $N_{eff}$ are also correlated with $r$. We will in this paper give the numerical values on the parameters introduced and show explicitly how the tension between the BICEP2 and Planck is effectively alleviated by the inclusion of the parameters $\Omega_k$, $\alpha_s$, $N_{eff}$ and $n_{t}$ separately.

In this paper, we study the cosmological constraints from the measurements of Hubble parameters---$H(z)$ data. Here, we consider two kinds of $H(z)$ data: the direct $H_0$ probe from the Hubble Space Telescope (HST) observations of Cepheid variables with $H_0=73.8\pm2.4$ ${\rm km\,s^{-1}\,Mpc^{-1}}$ and several measurements on the Hubble parameter at high redshifts $H(z)$. Employing Markov Chain Monte Carlo method, we also combine the WMAP nine-year data (WMAP9), the baryon acoustic oscillations (BAO) and type Ia supernovae (SNIa) "Union2.1" compilation to determine the cosmological parameters, such as the equation of state (EoS) of dark energy $w$, the curvature of the universe $\Omega_k$, the total neutrino mass $\sum{m_\nu}$, the effective number of neutrinos $N_{\rm eff}$, and the parameters associated with the power spectrum of primordial fluctuations. These $H(z)$ data provide extra information on the accelerate rate of our Universe at high redshifts. Therefore, adding these $H(z)$ data significantly improves the constraints on cosmological parameters, such as the number of relativistic species. Moreover, we find that direct prior on $H_0$ from HST can also give good constraints on some parameters, due to the degeneracies between these parameters and $H_0$.

In the context of the study of the Integrated Sachs Wolfe effect (ISW), we construct a template of the projected density distribution up to $z\simeq 0.7$ by using the Luminous Galaxies (LGs) from the Sloan Digital Sky Survey DR8. We use a photo-z catalogue trained with more than a hundred thousand galaxies from BOSS in the SDSS DR8 imaging area. We consider two different LG samples whose selection matches that of SDSS-III/BOSS: the LOWZ sample ($z\in [0.15,0.5]$) and the CMASS sample ($z\in[0.4,0.7]$). When building the LG density maps we use the information from star density, survey footprint, seeing conditions, sky emission, dust extinction and airmass to explore the impact of these artifacts on the two LG samples. In agreement with previous studies, we find that the CMASS sample is particularly sensitive to Galactic stars, which dominate the contribution to the auto-angular power spectrum below $\ell=7$. Other potential systematics affect mostly the low multipole range ($\ell\in[2,7]$), but leave fluctuations on smaller scales practically unchanged. The resulting power spectra in the multipole range $\ell\in[2,100]$ for the LOWZ, CMASS and LOWZ+CMASS samples are compatible with linear $\Lambda$CDM expectations and constant bias values of $b=1.98 \pm 0.11$, $2.08\pm0.14$ and $1.88\pm 0.11$, respectively, with no traces of non-Gaussianity: $f_{\rm NL}^{\rm local}=59\pm 75$ at 95% confidence level for the full LOWZ+CMASS sample in the range $\ell\in[4,100]$. After cross-correlating WMAP-9yr data with the LOWZ+CMASS LG density field, the ISW signal is detected at the level of 1.62--1.69$\,\sigma$. While this result is in close agreement with predictions from Monte Carlo simulations in the concordance $\Lambda$CDM model, it cannot rule out by itself an Einstein-de Sitter scenario, and has a moderately low signal compared to previous studies conducted on subsets of this LG sample.

We use a suite of N-body simulations that incorporate massive neutrinos as an extra-set of particles to investigate their effect on the halo mass function. We show that for cosmologies with massive neutrinos the mass function of dark matter haloes selected using the spherical overdensity (SO) criterion is well reproduced by the fitting formula of Tinker et al. (2008) once the cold dark matter power spectrum is considered instead of the total matter power, as it is usually done. The differences between the two implementations, i.e. using $P_{\rm cdm}(k)$ instead of $P_{\rm m}(k)$, are more pronounced for large values of the neutrino masses and in the high end of the halo mass function: in particular, the number of massive haloes is higher when $P_{\rm cdm}(k)$ is considered rather than $P_{\rm m}(k)$. As a quantitative application of our findings we consider a Planck-like SZ-clusters survey and show that the differences in predicted number counts can be as large as $30\%$ for $\sum m_\nu = 0.4$ eV. Finally, we use the Planck-SZ clusters sample, with an approximate likelihood calculation, to derive Planck-like constraints on cosmological parameters. We find that, in a massive neutrino cosmology, our correction to the halo mass function produces a shift in the $\sigma_8(\Omega_{\rm m}/0.27)^\gamma$ relation which can be quantified as $\Delta \gamma \sim 0.05$ and $\Delta \gamma \sim 0.14$ assuming one ($N_\nu=1$) or three ($N_\nu=3$) degenerate massive neutrino, respectively. The shift results in a lower mean value of $\sigma_8$ with $\Delta \sigma_8 = 0.01$ for $N_\nu=1$ and $\Delta \sigma_8 = 0.02$ for $N_\nu=3$, respectively. Such difference, in a cosmology with massive neutrinos, would increase the tension between cluster abundance and Planck CMB measurements.

We perform a combined perturbation and observational investigation of the scenario of non-minimal derivative coupling between a scalar field and curvature. First we extract the necessary condition that ensures the absence of instabilities, which is fulfilled more sufficiently for smaller coupling values. Then using Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) observations, we show that, contrary to its significant effects on inflation, the non-minimal derivative coupling term has a negligible effect on the universe acceleration, since it is driven solely by the usual scalar-field potential. Therefore, the scenario can provide a unified picture of early and late time cosmology, with the non-minimal derivative coupling term responsible for inflation, and the usual potential responsible for late-time acceleration. Additionally, the fact that the necessary coupling term does not need to be large, improves the model behavior against instabilities.

Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.

Recently, the Planck collaboration has released the first cosmological papers providing the high resolution, full sky, maps of the cosmic microwave background (CMB) temperature anisotropies. It is crucial to understand that whether the accelerating expansion of our universe at present is driven by an unknown energy component (Dark Energy) or a modification to general relativity (Modified Gravity). In this paper we study the coupled dark energy models, in which the quintessence scalar field nontrivially couples to the cold dark matter, with the strength parameter of interaction $beta$. Using the Planck data alone, we obtain that the strength of interaction between dark sectors is constrained as $beta < 0.102$ at $95%$ confidence level, which is tighter than that from the WMAP9 data alone. Combining the Planck data with other probes, like the Baryon Acoustic Oscillation (BAO), Type-Ia supernovae ``Union2.1 compilation'' and the CMB lensing data from Planck measurement, we find the tight constraint on the strength of interaction $beta < 0.052$ ($95%$ C.L.). Interestingly, we also find a non-zero coupling $beta = 0.078 pm 0.022$ ($68%$ C.L.) when we use the Planck, the ``SNLS'' supernovae samples, and the prior on the Hubble constant from the Hubble Space Telescope (HST) together. This evidence for the coupled dark energy models mainly comes from a tension between constraints on the Hubble constant from the Planck measurement and the local direct $H_0$ probes from HST.

We perform a forecast analysis on how well a Euclid-like photometric galaxy cluster survey will constrain the total neutrino mass and effective number of neutrino species. We base our analysis on the Monte Carlo Markov Chains technique by combining information from cluster number counts and cluster power spectrum. We find that combining cluster data with CMB measurements from Planck improves by more than an order of magnitude the constraint on neutrino masses compared to each probe used independently. For the LCDM+m_nu model the 2 sigma upper limit on total neutrino mass shifts from M_nu < 0.35 eV using cluster data alone to M_nu < 0.031 eV when combined with CMB data. When a non-standard model with N_eff number of neutrino species is considered, we estimate N_eff<3.14 (95% CL), while the bounds on neutrino mass are relaxed to M_nu < 0.040 eV. This accuracy would be sufficient for a 2 sigma detection of neutrino mass even in the minimal normal hierarchy scenario. We also consider scenarios with a constant dark energy equation of state and a non-vanishing curvature. When these models are considered the error on M_nu is only slightly affected, while there is a larger impact of the order of ~ 15 % and ~ 20% respectively on the 2 sigma error bar of N_eff with respect to the standard case. We also treat the LCDM+m_nu+N_eff case with free nuisance parameters, which parameterize the uncertainties on the cluster mass determination. In this case, the upper bounds on M_nu are relaxed by a factor larger than two, M_nu < 0.083 eV (95% CL), hence compromising the possibility of detecting the total neutrino mass with good significance. We thus confirm the potential that a large optical/near-IR cluster survey, like that to be carried out by Euclid, could have in constraining neutrino properties [abridged].

Recently, the Planck collaboration has released the first cosmological papers providing the highest resolution, full sky, maps of the cosmic microwave background (CMB) temperature anisotropies. In this paper we study a phenomenological model which interpolates between the pure $\Lambda$CDM model and the Dvali-Gabadadze-Porrati (DGP) braneworld model with an additional parameter $\alpha$. Firstly, we calculate the "distance information" of Planck data which includes the "shift parameter" $R$, the "acoustic scale" $l_A$, and the photon decoupling epoch $z_\ast$ in different cosmological models and find that this information is almost independent on the input models we use. Then, we compare the constraints on the free parameter $\alpha$ of the DGP model from the "distance information" of Planck and WMAP data and find that the Planck data with high precision do not improve the constraint on $\alpha$, but give the higher median value and the better limit on the current matter density fraction $\Omega_m$. Then, combining the "distance information" of Planck measurement, baryon acoustic oscillations (BAO), type Ia supernovae (SNIa) and the prior on the current Hubble constant (HST), we obtain the tight constraint on the parameter $\alpha < 0.20$ at $95\%$ confidence level, which implies that the flat DGP model has been ruled out by the current cosmological data. Finally, we allow the additional parameter $\alpha < 0$ in our calculations and interestingly obtain $\alpha=-0.29\pm0.20$ ($68\%$ C.L.), which means the current data slightly favor the effective equation of state $w_{\rm eff}<-1$. More importantly, the tension between constraints on $H_0$ from different observational data has been eased.

The Planck collaboration has recently published maps of the Cosmic Microwave Background radiation with the highest precision. In the standard flat $\Lambda$CDM framework, Planck data show that the Hubble constant $H_0$ is in tension with that measured by the several direct probes on $H_0$. In this paper, we perform a global analysis from the current observational data in the general dark energy models and find that resolving this tension on $H_0$ requires the dark energy model with its equation of state (EoS) $w\neq-1$. Firstly, assuming the $w$ to be a constant, the Planck data favor $w < -1$ at about $2\,\sigma$ confidence level when combining with the supernovae "SNLS" compilation. And consequently the value derived on $H_0$, $H_0=71.3\pm2.0$ ${\rm km\,s^{-1}\,Mpc^{-1}}$ (68% C.L.), is consistent with that from direct $H_0$ probes. We then investigate the dark energy model with a time-evolving $w$, and obtain the 68% C.L. constraints $w_0=-0.81\pm0.19$ and $w_a=-1.9\pm1.1$ from the Planck data and the "SNLS" compilation. Current data still slightly favor the Quintom dark energy scenario with EoS across the cosmological constant boundary $w\equiv-1$.

Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds. This paper describes the component separation framework adopted by Planck. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to l = 2000. The parameter constraints on LambdaCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations. For the low-frequency component, the spectral index varies widely over the sky, ranging from about beta = -4 to -2. Considering both morphology and prior knowledge of the low frequency components, the index map allows us to associate a steep spectral index (beta < -3.2) with strong anomalous microwave emission, corresponding to a spinning dust spectrum peaking below 20 GHz, a flat index of beta > -2.3 with strong free-free emission, and intermediate values with synchrotron emission.

Recently, the nine-year data release of the Wilkinson Microwave Anisotropy Probe (WMAP9) found that the inflationary models with the scalar spectral index n_s \geq 1 are excluded at about 5\sigma confidence level. In this paper, we set the new limits on the scalar spectral index in different cosmological models combining the WMAP9 data with the small-scale cosmic microwave background measurement from the South Pole Telescope, baryon acoustic oscillation data, Hubble Telescope measurements of the Hubble constant, and supernovae luminosity distance data. In most of extended cosmological models, e.g. with a dark energy equation of state, the constraints on n_s do not change significantly. The Harrison-Zel'dovich-Peebles (HZ) scale invariant spectrum is still disfavored at more than 4\sigma confidence level. However, when considering the model with a number of relativistic species N_{eff}, we obtain the limit on the spectral index of n_s=0.980\pm0.011 (1\sigma), due to the strong degeneracy between n_s and N_{eff}. The HZ spectrum now is still consistent with the current data at 95% confidence level.

[Abridged] We present a comprehensive investigation of the cosmological evolution of the luminosity function (LF) of galaxies and active galactic nuclei (AGN) in the infrared (IR). Based on the observed dichotomy in the ages of stellar populations of early-type galaxies on one side and late-type galaxies on the other, the model interprets the epoch-dependent LFs at z \geq 1.5 using a physical model for the evolution of proto-spheroidal galaxies and of the associated AGNs, while IR galaxies at z<1.5 are interpreted as being mostly late-type 'cold' (normal) and 'warm' (starburst) galaxies. As for proto-spheroids, in addition to the epoch-dependent LFs of stellar and AGN components separately, we have worked out the evolving LFs of these objects as a whole (stellar plus AGN component). The model provides a physical explanation for the observed positive evolution of both galaxies and AGNs up to z \simeq 2.5 and for the negative evolution at higher redshifts, for the sharp transition from Euclidean to extremely steep counts at (sub-)mm wavelengths, as well as the (sub-)mm counts of strongly lensed galaxies, that are hard to account for by alternative, physical or phenomenological, approaches. The evolution of late-type galaxies and of z<1.5 AGNs is described using a parametric phenomenological approach. The modeled AGN contributions to the counts and to the cosmic infrared background (CIB) are always subdominant with maximal at mid-IR wavelengths. The model provides a good fit to the multi-wavelength (from the mid-IR to millimeter waves) data on LFs at different redshifts and on number counts (both global and per redshift slices). A prediction of the present model is a systematic variation with wavelength of the populations dominating the counts and the contributions to the CIB intensity. The implied specific trend for cross-wavelength CIB power spectra is found to be in good agreement with the data.

Cosmography is the tool that makes possible to untie the interpretation of cosmological observations from the definition of any dynamical prior. We review the constraints on the cosmographic parameter obtained using the most thorough data set ensemble available. We focus on some specific topics about the statistically based selection of the most stringent fitting expansion.

In this paper, we study the degeneracies among several cosmological parameters in detail and discuss their impacts on the determinations of these parameters from the current and future observations. By combining the latest data sets, including type-Ia supernovae "Union2.1" compilation, WMAP seven-year data and the baryon acoustic oscillations from the SDSS Data Release Seven, we perform a global analysis to determine the cosmological parameters, such as the equation of state of dark energy w, the curvature of the universe \Omega_k, the total neutrino mass \sum{m_\nu}, and the parameters associated with the power spectrum of primordial fluctuations (n_s, \alpha_s and r). We pay particular attention on the degeneracies among these parameters and the influences on their constraints, by with or without including these degeneracies, respectively. We find that $w$ and \Omega_k or \sum{m_\nu} are strongly correlated. Including the degeneracies will significantly weaken the constraints. Furthermore, we study the capabilities of future observations and find these degeneracies can be broken very well. Consequently, the constraints of cosmological parameters can be improved dramatically.

We have investigated the Cosmic Infrared Background (CIB) anisotropies in the framework of the physical evolutionary model for proto-spheroidal galaxies by Granato et al. (2004). After having re-calibrated the cumulative flux function $dS/dz$ at $\lambda \ge 850\,\mu$m using the available determinations of the shot noise amplitude (the original model already correctly reproduces it at shorter wavelengths) the CIB power spectra at wavelengths from $250\,\mu$m to $2\,$mm measured by {\it Planck}, {\it Herschel}, SPT and ACT experiments have been fitted using the halo model with only 2 free parameters, the minimum halo mass and the power-law index of the mean occupation function of satellite galaxies. The best-fit {\it minimum} halo mass is $\log(M_{\rm min}/M_\odot) = 12.24 \pm 0.06$, higher than, but consistent within the errors, with the estimate by Amblard et al. (2011) and close to the estimate by Planck Collaboration (2011). The redshift evolution of the volume emissivity of galaxies yielded by the model is found to be consistent with that inferred from the data. The derived {\it effective} halo mass, $M_{\rm eff} \simeq 5\times 10^{12}\,M_\odot$, of $z\simeq 2$ sub-millimeter galaxies is close to that estimated for the most efficient star-formers at the same redshift. The effective bias factor and the comoving clustering radius at $z\simeq 2$ yielded by the model are substantially lower than those found for a model whereby the star formation is fueled by steady gas accretion, but substantially higher than those found for a merging-driven galaxy evolution with a top-heavy initial mass function.

We use the galaxy angular power spectrum at $z\sim0.5-1.2$ from the Canada-France-Hawaii-Telescope Legacy Survey Wide fields (CFHTLS-Wide) to constrain separately the total neutrino mass $\sum{m_\nu}$ and the effective number of neutrino species $N_{\rm{eff}}$. This survey has recently benefited from an accurate calibration of the redshift distribution, allowing new measurements of the (non-linear) matter power spectrum in a unique range of scales and redshifts sensitive to neutrino free streaming. Our analysis makes use of a recent model for the effect of neutrinos on the weakly non-linear matter power spectrum derived from accurate N-body simulations. We show that CFHTLS, combined with WMAP7 and a prior on the Hubble constant provides an upper limit of $\sum{m_\nu}<0.29\,$eV and $N_{\rm{eff}} =4.17^{+1.62}_{-1.26}$ (2$\,\sigma$ confidence levels). If we omit smaller scales which may be affected by non-linearities, these constraints become $\sum{m_\nu}<0.41\,$eV and $N_{\rm{eff}} =3.98^{+2.02}_{-1.20}$ (2$\,\sigma$ confidence levels). Finally we show that the addition of other large scale structures probes can further improve these constraints, demonstrating that high redshift large volumes surveys such as CFHTLS are complementary to other cosmological probes of the neutrino mass.

We perform a cosmographic analysis using several cosmological observables such as the luminosity distance moduli, the volume distance, the angular diameter distance and the Hubble parameter. These quantities are determined using different data sets: Supernovae type Ia and Gamma Ray Bursts, the Baryonic Acoustic Oscillations, the cosmic microwave background power spectrum and the Hubble parameter as measured from surveys of galaxies. This data set allows to put constraints on the cosmographic expansion with unprecedented precision. We also present forecasts for the coefficients of the kinematic expansion using future but realistic data sets: constraints on the coefficients of the expansions are likely to improve by a factor ten with the upcoming large scale structure probes. Finally, we derive the set of the cosmographic parameters for several cosmological models (including $\Lambda$CDM) and compare them with our best fit set. While distance measurements are unable to discriminate among these models, we show that the inclusion of the Hubble data set leads to strong constraints on the lowest order coefficients and in particular it is incompatible with $\Lambda$CDM at 3-$\sigma$ confidence level. We discuss the reliability of this determination and suggest further observations which might be of crucial importance for the viability of cosmographic tests in the next future.

In this paper we study the possibility of testing Charge-Parity-Time Reversal (CPT) symmetry with cosmic microwave background (CMB) experiments. We consider two kinds of Chern-Simons (CS) term, electromagnetic CS term and gravitational CS term, and study their effects on the CMB polarization power spectra in detail. By combining current CMB polarization measurements, the seven-year WMAP, BOOMERanG 2003 and BICEP observations, we obtain a tight constraint on the rotation angle $\Delta\alpha=-2.28\pm1.02$ deg ($1\,\sigma$), indicating a $2.2\,\sigma$ detection of the CPT violation. Here, we particularly take the systematic errors of CMB measurements into account. After adding the QUaD polarization data, the constraint becomes $-1.34<\Delta\alpha<0.82$ deg at 95% confidence level. When comparing with the effect of electromagnetic CS term, the gravitational CS term could only generate TB and EB power spectra with much smaller amplitude. Therefore, the induced parameter $\epsilon$ can not be constrained from the current polarization data. Furthermore, we study the capabilities of future CMB measurements, Planck and CMBPol, on the constraints of $\Delta\alpha$ and $\epsilon$. We find that the constraint of $\Delta\alpha$ can be significantly improved by a factor of 15. Therefore, if this rotation angle effect can not be taken into account properly, the constraints of cosmological parameters will be biased obviously. For the gravitational CS term, the future Planck data still can not constrain $\epsilon$ very well, if the primordial tensor perturbations are small, $r <0.1$. We need the more accurate CMBPol experiment to give better constraint on $\epsilon$.

In this work, starting from 21 months of data from the Fermi-Large Area Telescope, we derive maps of the residual isotropic gamma-ray emission, a relevant fraction of which is expected to be contributed by the extragalactic diffuse gamma-ray background. We compute the angular two-point auto-correlation function of the residual Fermi-LAT maps at energies E>1GeV, E>3GeV and E>30GeV well above the Galactic plane and find no significant correlation signal. This is, indeed, what is expected if the EGB were contributed by BL Lacertae, Flat Spectrum Radio Quasars or star-forming galaxies, since, in this case, the predicted signal is very weak. Then, we search for the Integrated Sachs-Wolfe signature by cross-correlating the Fermi-LAT maps with the WMAP7-Cosmic Microwave Background map. We find a cross-correlation consistent with zero, even though the expected signal is larger than that of the EGB auto-correlation. Finally, in an attempt to constrain the nature of the gamma-ray background we cross-correlate the Fermi-LAT maps with the angular distributions of objects that may contribute to the EGB: QSOs in the SDSS-DR6 catalog, NVSS galaxies, 2MASS galaxies and LRG in the SDSS catalog. The cross-correlation is always consistent with zero, in agreement with theoretical expectations, but we find (with low statistical significance) some interesting features that may indicate that some specific classes of objects contribute to the EGB. A chi2 analysis confirms that the correlation properties of the 21-month data do not provide strong constraints of the EGB origin. However, the results suggest that the situation will significantly improve with the 5- and 10-year Fermi-LAT data. The future EGB analysis will then allow placing significant constraints on the nature of the EGB and might provide in addition a detection of the ISW signal. (Abridged)

In this paper we measure the angular power spectra $C_\ell$ of three high-redshift large-scale structure probes: the radio sources from the NRAO VLA Sky Survey (NVSS), the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs) and the MegaZ-LRG (DR7), the final SDSS II Luminous Red Galaxy (LRG) photometric redshift survey. We perform a global analysis of the constraints on the amplitude of primordial non-Gaussianity from these angular power spectra, as well as from their cross-correlation power spectra with the cosmic microwave background (CMB) temperature map. In particular, we include non-Gaussianity of the type arising from single-field slow roll, multifields, curvaton (local type), and those which effects on the halo clustering can be described by the equilateral template (related to higher-order derivative type non-Gaussianity) and by the enfolded template (related to modified initial state or higher-derivative interactions). When combining all data sets, we obtain limits of $f_{\rm NL}=48\pm20$, $f_{\rm NL}=50\pm265$ and $f_{\rm NL}=183\pm95$ at 68% confidence level for local, equilateral and enfolded templates, respectively. Furthermore, we explore the constraint on the cubic correction $g_{\rm NL}\phi^3$ on the bias of dark matter haloes and obtain a limit of $-1.2\times10^5<g_{\rm NL}<11.3\times10^5$ at 95% confidence level.

In this paper we present the constraints on cold dark matter (CDM) isocurvature contributions to the cosmological perturbations. By employing Markov Chain Monte Carlo method (MCMC), we perform a global analysis for cosmological parameters using the latest astronomical data, such as 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) observations, matter power spectrum from the Sloan Digital Sky Survey (SDSS) luminous red galaxies (LRG), and "Union2" type Ia Supernovae (SNIa) sample. We find that the correlated mixture of adiabatic and isocurvature modes are mildly better fitting to the current data than the pure adiabatic ones, with the minimal $\chi^2$ given by the likelihood analysis being reduced by 3.5. We also obtain a tight limit on the fraction of the CDM isocurvature contributions, which should be less than 14.6% at 95% confidence level. With the presence of the isocurvature modes, the adiabatic spectral index becomes slightly bigger, n_s^{\rm adi}=0.972\pm0.014~(1\,\sigma), and the tilt for isocurvature spectrum could be large, namely, the best fit value is n_s^{\rm iso}=3.020. Finally, we discuss the effect on WMAP normalization priors, shift parameter R, acoustic scale l_A and z_{*}, from the CDM isocurvaure perturbation. By fitting the mixed initial condition to the combined data, we find the mean values of R, l_A and z_{*} can be changed about 2.9\sigma, 2.8\sigma and 1.5\sigma respectively, comparing with those obtained in the pure adiabatic condition.

We present an analysis of the constraints on the amplitude of primordial non-Gaussianity of local type described by the dimensionless parameter $f_{\rm NL}$. These constraints are set by the auto-correlation functions (ACFs) of two large scale structure probes, the radio sources from NRAO VLA Sky Survey (NVSS) and the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs), as well as by their cross-correlation functions (CCFs) with the cosmic microwave background (CMB) temperature map (Integrated Sachs-Wolfe effect). Several systematic effects that may affect the observational estimates of the ACFs and of the CCFs are investigated and conservatively accounted for. Our approach exploits the large-scale scale-dependence of the non-Gaussian halo bias. The derived constraints on {$f_{\rm NL}$} coming from the NVSS CCF and from the QSO ACF and CCF are weaker than those previously obtained from the NVSS ACF, but still consistent with them. Finally, we obtain the constraints on $f_{\rm NL}=53\pm25$ ($1\,\sigma$) and $f_{\rm NL}=58\pm24$ ($1\,\sigma$) from NVSS data and SDSS DR6 QSO data, respectively.

The NRAO VLA Sky Survey (NVSS) is the only dataset that allows an accurate determination of the auto-correlation function (ACF) on angular scales of several degrees for Active Galactic Nuclei (AGNs) at typical redshifts $z \simeq 1$. Surprisingly, the ACF is found to be positive on such large scales while, in the framework of the standard hierarchical clustering scenario with Gaussian primordial perturbations it should be negative for a redshift-independent effective halo mass of order of that found for optically-selected quasars. We show that a small primordial non-Gaussianity can add sufficient power on very large scales to account for the observed NVSS ACF. The best-fit value of the parameter $f_{\rm NL}$, quantifying the amplitude of primordial non-Gaussianity of local type is $f_{\rm NL}=62 \pm 27$ ($1\,\sigma$ error bar) and $25<f_{\rm NL}<117$ ($2\,\sigma$ confidence level), corresponding to a detection of non-Gaussianity significant at the $\sim 3\,\sigma$ confidence level. The minimal halo mass of NVSS sources is found to be $M_{\rm min}=10^{12.47\pm0.26}h^{-1}M_{\odot}$ ($1\,\sigma$) strikingly close to that found for optically selected quasars. We discuss caveats and possible physical and systematic effects that can impact on the results.

The electrodynamics modified by the Chern-Simons term $\mathcal{L}_{cs}\sim p_{\mu}A_{\nu}\tilde F^{\mu\nu}$ with a non-vanishing $p_{\mu}$ violates the {\it Charge-Parity-Time Reversal} symmetry (CPT) and rotates the linear polarizations of the propagating {\it Cosmic Microwave Background} (CMB) photons. In this paper we measure the rotation angle $\Delta\alpha$ by performing a global analysis on the current CMB polarization measurements from the {\it seven-year Wilkinson Microwave Anisotropy Probe} (WMAP7), {\it BOOMERanG 2003} (B03), BICEP and QUaD using a Markov Chain Monte Carlo method. Neglecting the systematic errors of these experiments, we find that the results from WMAP7, B03 and BICEP all are consistent and their combination gives $\Delta\alpha=-2.33\pm0.72~{\rm deg}~(68\%~C.L.)$, indicating a $3\,\sigma$ detection of the CPT violation. The QUaD data alone gives $\Delta\alpha=0.59\pm0.42~{\rm deg}~(68\%~C.L.)$ which has an opposite sign for the central value and smaller error bar compared to that obtained from WMAP7, B03 and BICEP. When combining all the polarization data together, we find $\Delta\alpha=-0.04\pm0.35~{\rm deg}~(68\%~C.L.)$ which significantly improves the previous constraint on $\Delta\alpha$ and test the validity of the fundamental CPT symmetry at a higher level.

We review the paradigm of quintom cosmology. This scenario is motivated by the observational indications that the equation of state of dark energy across the cosmological constant boundary is mildly favored, although the data are still far from being conclusive. As a theoretical setup we introduce a no-go theorem existing in quintom cosmology, and based on it we discuss the conditions for the equation of state of dark energy realizing the quintom scenario. The simplest quintom model can be achieved by introducing two scalar fields with one being quintessence and the other phantom. Based on the double-field quintom model we perform a detailed analysis of dark energy perturbations and we discuss their effects on current observations. This type of scenarios usually suffer from a manifest problem due to the existence of a ghost degree of freedom, and thus we review various alternative realizations of the quintom paradigm. The developments in particle physics and string theory provide potential clues indicating that a quintom scenario may be obtained from scalar systems with higher derivative terms, as well as from non-scalar systems. Additionally, we construct a quintom realization in the framework of braneworld cosmology, where the cosmic acceleration and the phantom divide crossing result from the combined effects of the field evolution on the brane and the competition between four and five dimensional gravity. Finally, we study the outsets and fates of a universe in quintom cosmology. In a scenario with null energy condition violation one may obtain a bouncing solution at early times and therefore avoid the Big Bang singularity. Furthermore, if this occurs periodically, we obtain a realization of an oscillating universe. Lastly, we comment on several open issues in quintom cosmology and their connection to future investigations.

In this paper, we combine the the latest observational data, including the WMAP five-year data (WMAP5), the baryon acoustic oscillations (BAO) and type Ia supernovae (SN) "union" compilation, and use the Markov Chain Monte Carlo method to determine the dark energy parameters. We pay particular attention to the Integrated Sache-Wolfe (ISW) data from the cross-correlations of cosmic microwave background (CMB) and large scale structure (LSS). In the \Lambda CDM model, we find that the ISW data, as a complement to the WMAP data, could significantly improve the constraint of curvature \Omega_k. We also check the improvement of constraints from the new prior on the Hubble constant and find this new prior could improve the constraint of \Omega_k by a factor of 2. Finally, we study the dynamical evolving EoS of dark energy from the current observational data. Based on the dynamical dark energy model, parameterizing as w(a)=w_0+w_a(1-a), we find that the \Lambda CDM model remains a good fit to the current data. When taking into account the ISW data, the error bars of w_0 and w_a could be shrunk slightly. Current constraints on the dynamical dark energy model are not conclusive. The future precision measurements are needed.

We constrain the parameters describing the kinematical state of the universe using a cosmographic approach, which is fundamental in that it requires a very minimal set of assumptions (namely to specify a metric) and does not rely on the dynamical equations for gravity. On the data side, we consider the most recent compilations of Supernovae and Gamma Ray Bursts catalogues. This allows to further extend the cosmographic fit up to $z = 6.6$, i.e. up to redshift for which one could start to resolve the low z degeneracy among competing cosmological models. In order to reliably control the cosmographic approach at high redshifts, we adopt the expansion in the improved parameter $y = z/(1+z)$. This series has the great advantage to hold also for $z > 1$ and hence it is the appropriate tool for handling data including non-nearby distance indicators. We find that Gamma Ray Bursts, probing higher redshifts than Supernovae, have constraining power and do require (and statistically allow) a cosmographic expansion at higher order than Supernovae alone. Exploiting the set of data from Union and GRBs catalogues, we show (for the first time in a purely cosmographic approach parametrized by deceleration $q_0$, jerk $j_0$, snap $s_0$) a definitively negative deceleration parameter $q_0$ up to the 3$\sigma$ confidence level. We present also forecasts for realistic data sets that are likely to be obtained in the next few years.

The coupled dark energy models, in which the quintessence scalar field nontrivially couples to the cold dark matter, have been proposed to explain the coincidence problem. In this paper we study the perturbations of coupled dark energy models and the effects of this interaction on the current observations. Here, we pay particular attention to its imprint on the late-time Integrated Sachs-Wolfe (ISW) effect. We perform a global analysis of the constraints on this interaction from the current observational data. Considering the typical exponential form as the interaction form, we obtain that the strength of interaction between dark sectors is constrained as $\beta<0.085$ at 95% confidence level. Furthermore, we find that future measurements with smaller error bars could improve the constraint on the strength of coupling by a factor two, when compared to the present constraints.

In this paper we rely on the quasar (QSO) catalog of the Sloan Digital Sky Survey Data Release Six (SDSS DR6) of about one million photometrically selected QSOs to compute the Integrated Sachs-Wolfe (ISW) effect at high redshift, aiming at constraining the behavior of the expansion rate and thus the behaviour of dark energy at those epochs. This unique sample significantly extends previous catalogs to higher redshifts while retaining high efficiency in the selection algorithm. We compute the auto-correlation function (ACF) of QSO number density from which we extract the bias and the stellar contamination. We then calculate the cross-correlation function (CCF) between QSO number density and Cosmic Microwave Background (CMB) temperature fluctuations in different subsamples: at high z>1.5 and low z<1.5 redshifts and for two different choices of QSO in a conservative and in a more speculative analysis. We find an overall evidence for a cross-correlation different from zero at the 2.7\sigma level, while this evidence drops to 1.5\sigma at z>1.5. We focus on the capabilities of the ISW to constrain the behaviour of the dark energy component at high redshift both in the \LambdaCDM and Early Dark Energy cosmologies, when the dark energy is substantially unconstrained by observations. At present, the inclusion of the ISW data results in a poor improvement compared to the obtained constraints from other cosmological datasets. We study the capabilities of future high-redshift QSO survey and find that the ISW signal can improve the constraints on the most important cosmological parameters derived from Planck CMB data, including the high redshift dark energy abundance, by a factor \sim 1.5.

The accelerating expansion of our universe at present could be driven by an unknown energy component (Dark Energy) or a modification of general relativity (Modified Gravity). In this note we revisit the constraints on a phenomenological model which interpolates between the pure \LambdaCDM model and the Dvali-Gabadadze-Porrati (DGP) braneworld model with an additional parameter \alpha. Combining the cosmic microwave background (CMB), baryon acoustic oscillations (BAO) and type Ia supernovae (SNIa), as well as some high-redshift observations, such as the gamma-ray bursts (GRB) and the measurements of linear growth factors (LGF), we obtain the tight constraint on the parameter \alpha=0.254\pm0.153 (68% C.L.), which implies that the flat DGP model is incompatible with the current observations, while the pure \LambdaCDM model still fits the data very well. Finally, we simulate the future measurements with higher precisions and find that the constraint on \alpha can be improved by a factor two, when compared to the present constraints.

Early dark energy models, for which the contribution to the dark energy density at high redshifts is not negligible, influence the growth of cosmic structures and could leave observable signatures that are different from the standard cosmological constant cold dark matter ($\Lambda$CDM) model. In this paper, we present updated constraints on early dark energy using geometrical and dynamical probes. From WMAP five-year data, baryon acoustic oscillations and type Ia supernovae luminosity distances, we obtain an upper limit of the dark energy density at the last scattering surface (lss), $\Omega_{\rm EDE}(z_{\rm lss})<2.3\times10^{-2}$ (95% C.L.). When we include higher redshift observational probes, such as measurements of the linear growth factors, Gamma-Ray Bursts (GRBs) and Lyman-$\alpha$ forest (\lya), this limit improves significantly and becomes $\Omega_{\rm EDE}(z_{\rm lss})<1.4\times10^{-3}$ (95% C.L.). Furthermore, we find that future measurements, based on the Alcock-Paczy\'nski test using the 21cm neutral hydrogen line, on GRBs and on the \lya forest, could constrain the behavior of the dark energy component and distinguish at a high confidence level between early dark energy models and pure $\Lambda$CDM. In this case, the constraints on the amount of early dark energy at the last scattering surface improve by a factor ten, when compared to present constraints. We also discuss the impact on the parameter $\gamma$, the growth rate index, which describes the growth of structures in standard and in modified gravity models.

One of the characteristics of the "Matter Bounce" scenario, an alternative to cosmological inflation for producing a scale-invariant spectrum of primordial adiabatic fluctuations on large scales, is a break in the power spectrum at a characteristic scale, below which the spectral index changes from n = 1 to n = 3. We study the constraints which current cosmological data place on the location of such a break, and more generally on the position of the break and the slope at length scales smaller than the break. The observational data we use include the WMAP five-year data set (WMAP5), other CMB data from BOOMERanG, CBI, VSA, and ACBAR, large-scale structure data from the Sloan Digital Sky Survey (SDSS, their luminous red galaxies sample), Type Ia Supernovae data (the "Union" compilation), and the Sloan Digital Sky Survey Lyman-\alpha forest power spectrum (Ly\alpha) data. We employ the Markov Chain Monte Carlo method to constrain the features in the primordial power spectrum which are motivated by the matter bounce model. We give an upper limit on the length scale where the break in the spectrum occurs.

In this paper we revisit the issue of determining the oscillating primordial scalar power spectrum and oscillating equation of state of dark energy from the astronomical observations. By performing a global analysis with the Markov Chain Monte Carlo method, we find that the current observations from five-year WMAP and SDSS-LRG matter power spectrum, as well as the "union" supernovae sample, constrain the oscillating index of primordial spectrum and oscillating equation of state of dark energy with the amplitude less than $|n_{\rm amp}|<0.116$ and $|w_{\rm amp}|<0.232$ at 95% confidence level, respectively. This result shows that the oscillatory structures on the primordial scalar spectrum and the equation of state of dark energy are still allowed by the current data. Furthermore, we point out that these kinds of modulation effects will be detectable (or gotten a stronger constraint) in the near future astronomical observations, such as the PLANCK satellite, LAMOST telescope and the currently ongoing supernovae projects SNLS.

In this paper, we study the cosmological implications of the 100 square degree Weak Lensing survey (the CFHTLS-Wide, RCS, VIRMOS-DESCART and GaBoDS surveys). We combine these weak lensing data with the cosmic microwave background (CMB) measurements from the WMAP5, BOOMERanG, CBI, VSA, ACBAR, the SDSS LRG matter power spectrum and the Type Ia Supernoave (SNIa) data with the "Union" compilation (307 sample), using the Markov Chain Monte Carlo method to determine the cosmological parameters. Our results show that the \Lambda CDM model remains a good fit to all of these data. For the dynamical dark energy model with time evolving EoS parameterized as w_{\DE}(a) = w_0 + w_a (1-a), we find that the best-fit model implying the mildly preference of Quintom model whose EoS gets across the cosmological constant boundary during evolution. Regarding the total neutrino mass limit, we obtain the upper limit, \sum m_{\nu}< 0.471 eV (95% C.L.) within the framework of the flat \Lambda CDM model. Due to the obvious degeneracies between the neutrino mass and the EoS of dark energy model, this upper limit will be relaxed by a factor of 2 in the framework of dynamical dark energy models. For the constraints on the inflation parameters, we find that the upper limit on the ratio of the tensor to scalar is r<0.35 (95% C.L.) and the inflationary models with the slope n_s\geq1 are excluded at more than 2 \sigma confidence level. In this paper we pay particular attention to the contribution from the weak lensing data and find that the current weak lensing data do improve the constraints on matter density \Omega_m, \sigma_8, \sum{m_{\nu}}, and the EoS of dark energy.

In this paper we study the sensitivity of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) project to the determination of cosmological parameters, employing the Monte Carlo Markov Chains (MCMC) method. For comparison, we first analyze the constraints on cosmological parameters from current observational data, including WMAP, SDSS and SN Ia. We then simulate the 3D matter power spectrum data expected from LAMOST, together with the simulated CMB data for PLANCK and the SN Ia from 5-year Supernovae Legacy Survey (SNLS). With the simulated data, we investigate the future improvement on cosmological parameter constraints, emphasizing the role of LAMOST. Our results show the potential of LAMOST in probing for the cosmological parameters, especially in constraining the equation-of-state (EoS) of the dark energy and the neutrino mass.

In this paper, we combine the latest observational data, including the WMAP five-year data (WMAP5), BOOMERanG, CBI, VSA, ACBAR, as well as the Baryon Acoustic Oscillations (BAO) and Type Ia Supernoave (SN) "Union" compilation (307 sample) to determine the cosmological parameters. Our results show that the $\Lambda$CDM model remains a good fit to the current data. In a flat universe, we obtain the tight limit on the constant EoS of dark energy as, $w=-0.977\pm0.056$ ($1 \sigma$). For the dynamical dark energy models with time evolving EoS, we find that the best-fit values are $w_0=-1.08$ and $w_1=0.368$, implying the preference of Quintom model whose EoS gets across the cosmological constant boundary. For the curvature of universe, our results give $-0.012<\Omega_k<0.009$ (95% C.L.) when fixing $w_{\DE}=-1$. When considering the dynamics of dark energy, the flat universe is still a good fit to the current data. Regarding the neutrino mass limit, we obtain the upper limits, $\sum m_{\nu}<0.533$ eV (95% C.L.) within the framework of the flat $\Lambda$CDM model. When adding the SDSS Lyman-$\alpha$ forest power spectrum data, the constraint on $\sum m_{\nu}$ can be significantly improved, $\sum m_{\nu}<0.161$ eV (95% C.L.). Assuming that the primordial fluctuations are adiabatic with a power law spectrum, within the $\Lambda$CDM model, we find that the upper limit on the ratio of the tensor to scalar is $r<0.200$ (95% C.L.) and the inflationary models with the slope $n_s\geq1$ are excluded at more than $2 \sigma$ confidence level. However, in the framework of dynamical dark energy models, the allowed region in the parameter space of ($n_s$,$r$) is enlarged significantly. Finally, we find no evidence for the large running of the spectral index. (Abridged)

Recently, the WMAP group has published their five-year data and considered the constraints on the time evolving equation of state of dark energy for the first time from the WMAP distance information. In this paper, we study the effectiveness of the usage of these distance information and find that these compressed CMB information can give similar constraints on dark energy parameters compared with the full CMB power spectrum if dark energy perturbations are included, however, once incorrectly neglecting the dark energy perturbations, the difference of the results are sizable.

In this letter, we propose a model of inflationary cosmology with a bounce preceded and study its primordial curvature perturbations. Our model gives rise to a primordial power spectrum with a feature of oscillation on large scales compared with the nearly scale-invariant spectrum generated by the traditional slow rolling inflation model. We will show this effect changes the Cosmic Microwave Background (CMB) temperature power spectrum and the Large Scale Structure (LSS) matter power spectrum. And further with a detailed simulation we will point out this signal is detectable to the forthcoming observations, such as PLANCK and LAMOST.

In this Letter we update our previous result on the test of CPT symmetry with Cosmic Microwave Background (CMB) measurements. A CPT violating interaction in the photon sector ${\cal L}_{cs} \sim p_{\mu}A_{\nu}\tilde F^{\mu\nu}$ gives rise to a rotation of the polarization vectors of the propagating CMB photons. Recently the WMAP group used the newly released polarization data of WMAP5 to measure this rotation angle $\Delta\alpha$ and obtained $\Delta\alpha=-1.7\pm2.1$ deg ($1\sigma$). However, in their analysis the BOOMERanG 2003 data is not included. Here we revisit this issue by combining the full data of WMAP5 and BOOMERanG 2003 angular power spectra for the measurement of this rotation angle $\Delta\alpha$ and find that $\Delta\alpha=-2.6\pm1.9$ deg at a 68% confidence level.

In this paper we study the possibility of testing \emph{CPT} symmetry with Cosmic Microwave Background (CMB) measurements. Working with an effective lagrangian of the photon with \emph{CPT} violation ${\cal L} \sim p_{\mu}A_{\nu}\tilde F^{\mu\nu}$ which causes the polarization vectors of the propagating CMB photons rotated, we determine the rotation angle $\Delta\alpha$ using the BOOMERanG 2003 and the WMAP3 angular power spectra. In this analysis we have included the newly released $TC$ and $GC$ ($l<450$) information of WMAP3 and found $\Delta\alpha=-6.2\pm3.8$ deg at $1\sigma$ confidence level. This result increases slightly the significance for the \emph{CPT} violation obtained in our previous paper (Feng \emph{et al.}, 2006) $\Delta\alpha=-6.0 \pm 4.0$ deg (1$\sigma$). Furthermore we examine the constraint on the rotation angle with the simulated Planck data. Our results show that the Planck will be sensitive to $\Delta \alpha$ at the level of 0.057 deg and able to test the \emph{CPT} symmetry with a higher precision.

In this paper we perform a global analysis of the constraints on the inflationary parameters in the presence of dynamical dark energy models from the current observations, including the three-year Wilkinson Microwave Anisotropy Probe (WMAP3) data, Boomerang-2K2, CBI, VSA, ACBAR, SDSS LRG, 2dFGRS and ESSENCE (192 sample). We use the analytic description of the inflationary power spectra in terms of the Horizon-flow parameters $\{\epsilon_i\}$. With the first order approximation in the slow-roll expansion, we find that the constraints on the Horizon-flow parameters are $\epsilon_1<0.014 (95% C.L.)$ and $\epsilon_2=0.034\pm0.024 (1\sigma)$ in the $\Lambda$CDM model. In the framework of dynamical dark energy models, the constraints become obviously weak, $\epsilon_1<0.022 (95% C.L.)$ and $\epsilon_2=-0.006\pm0.039 (1\sigma)$, and the inflation models with a "blue" tilt, which are excluded about $2\sigma$ in the $\Lambda$CDM model, are allowed now. With the second order approximation, the constraints on the Horizon-flow parameters are significantly relaxed further. If considering the non-zero $\epsilon_3$, the large running of the scalar spectral index is found for the $\Lambda$CDM model, as well as the dynamical dark energy models.

We investigate the constraints on cosmological parameters especially for EoS of dark energy, inflationary parameters, neutrino mass and curvature of universe using simulated Planck data. Firstly we determine cosmological parameters with current observations including ESSENCE, WMAP3, Boomerang-2K2, CBI, VSA, ACBAR, SDSS LRG and 2dFGRS, and take best-fit model as the fiducial model in simulations. In simulations we pay attention to the effects of dynamical dark energy in determination of cosmological parameters. We add simulated SNAP data to do all the simulations. Using present data, we find Quintom dark energy model is mildly favored while \LambdaCDM remains a good fit. In the framework of dynamical dark energy, the constraints on inflationary parameters, m_{\nu} and \Omega_{K} become weak compared with the constraints in \LambdaCDM. Intriguingly, we find that the inflationary models with a "blue" tilt, which are excluded about 2\sigma in \LambdaCDM model, are well within 2\sigma region with the presence of the dynamics of dark energy. The upper limits of neutrino mass are weakened by a factor of 2 (95% C.L.), say, m_{\nu}<1.59 eV and m_{\nu}<1.53 eV for two forms of parametrization of the equation of state of dark energy. The flat universe is a good fit to the current data, namely, |\Omega_{K}|<0.03 (95% C.L.). With the simulated Planck and SNAP data, dynamical dark energy and \LambdaCDM might be distinguished at 4\sigma. And uncertainties of inflationary parameters, m_{\nu} and \Omega_{K} can be reduced obviously. We also constrain the rotation angle \Delta\alpha, denoting possible cosmological CPT violation, with simulated Planck and CMBpol data and find that our results are much more stringent than current constraint and will verify cosmological CPT symmetry with a higher precision. (Abridged)

In this paper we study the sound speed - $c_s^2$, which is directly related to the classical perturbations - of the dynamical dark energy (DE), especially with an equation of state crossing the cosmological constant boundary in details and show its implications on Cosmic Microwave Background (CMB) Anisotropy. With the present observational data of CMB, Type Ia Supernova (SNIa) and galaxy clustering, we perform a global analysis to constrain the sound speed of dark energy, using the Markov Chain Monte Carlo method. We find that the sound speed of dark energy is weakly constrained by current observations thus the futuristic precision measurements of CMB on a very large angular scale (low multipoles) are necessary.

Due to the lack of low redshift long Gamma-Ray Bursts (GRBs), the circular problem has been a severe obstacle for using GRBs as cosmological candles. In this paper, we present a new method to deal with such a problem in MCMC global fitting analysis. Assuming that a certain type of correlations between different observables exists in a subsample of GRBs, for the parameters involved in the correlation relation, we treat them as free parameters and determine them simultaneously with cosmological parameters through MCMC analysis on GRB data together with other observational data. Then the circular problem is naturally eliminated in this procedure. We take the Ghirlanda relation as an example while keeping in mind the debate about its physical validity. Together with SNe Ia, WMAP and SDSS data, we include 27 GRBs with the reported Ghirlanda relation in our study, and perform MCMC global fitting. We consider the $\Lambda$CDM model and dynamical dark energy models. In each case, in addition to the constraints on the relevant cosmological parameters, we obtain the best fit values as well as the distributions of the correlation parameters $A$ and $C$. We find that the observational data sets other than GRBs can affect $A$ and $C$ considerably through their degeneracies with the cosmological parameters. The results on $A$ and $C$ for different cosmological models are in well agreement within $1\sigma$ range. The best fit value of $A$ in all models being analyzed is $A\sim 1.53$ with $\sigma \sim 0.08$. For $C$, we have the best value in the range of $0.94-0.98$ with $\sigma\sim 0.1$. It is also noted that the distributions of $A$ and $C$ are generally broader than the priors used in many studies in literature. (Abriged)

In this paper we study the cosmological CPT-violation and its implications in baryo/leptogenesis and CMB polarization. We propose specifically a variant of the models of gravitational leptogenesis. By performing a global analysis with the Markov Chain Monte Carlo (MCMC) method, we find the current CMB polarization observations from the three-year WMAP (WMAP3) and the 2003 flight of BOOMERANG (B03) data provide a weak evidence for our model. However to verify and especially exclude this type of mechanism for baryo/leptogenesis with cosmological CPT-violation, the future measurements on CMB polarization from PLANCK and CMBpol are necessary.

With the latest astronomical data including Cosmic Microwave Background (WMAP three year, CBI, ACBAR, VSA), Type Ia Supernova ("gold sample"), Galaxy Clustering (SDSS 3-D matter power, Lyman-$\alpha$ forest and Baryon Acoustic Oscillating (BAO)), we make a global fitting to constrain the mass varying neutrinos. We find that the parameter $\delta$, denoting time evolving of neutrino mass, is weakly constrained and the neutrino mass limit today can be relaxed at least by a factor of two. Adding data of $0\nu2\beta$ decay of Heidelberg-Moscow experiment to our analysis, we find that $\delta$ can be constrained tightly and mass varying neutrinos are favored at about 99.7% confidence level.

We investigate the correlation between the neutrino mass limit and dark energy with time evolving equation of state. Parameterizing dark energy as w=w_0+w_1*z/(1+z), we make a global fit using Markov Chain Monte Carlo technique to determine w_0, w_1, neutrino mass as well as other cosmological parameters simultaneously. We pay particular attention to the correlation between neutrino mass \Sigma m_{\nu} and w_1 using current cosmological observations as well as the future simulated datasets such as PLANCK, SNAP and LAMOST.

We use the newly released 182 Type Ia supernova data combined with the third-year Wilkinson Microwave Anisotropic Probe data (WMAP3) and large scale structure (LSS) information including SDSS and 2dFGRS to constrain the dark energy equation of state (EoS) as well as the curvature of universe $\Omega_K$. Using the full dataset of Cosmic Microwave Background (CMB) and LSS rather than the shift parameter and linear growth factor, we make a Markov Chain Monte Carlo (MCMC) global fit, while paying particular attention to the dark energy perturbation. Parameterizing the EoS as $w_{DE}(a) = w_{0} + w_{1}(1-a)$, we find the best fit of ($w_0,w_1$) is ($-1.053,0.944$) and for $w_{DE}(a) = w_{0} + w_{1}\sin({3/2}\pi \ln(a))$, the best fit for ($w_0,w_1$) is ($-1.614,-1.046$). We find that a flat universe is a good approximation, namely, $|\Omega_K|>0.06$ has been excluded by 2$\sigma$ yet the inclusion of $\Omega_K$ can affect the measurement of DE parameters owing to their correlation and the present systematic effects of SNIa measurements.

The three-year WMAP(WMAP3), combined with other cosmological observations from galaxy clustering and Type Ia Supernova (SNIa), prefers a non-vanishing running of the primordial spectral index independent of the low CMB multipoles. Motivated by this feature we study cosmological constraint on the neutrino mass, which severely depends on what prior we adopt for the spectral shape of primordial fluctuations, taking possible running into account. As a result we find a more stringent constraint on the sum of the three neutrino masses, m_\nu < 0.76eV (2 \sigma), compared with m_\nu < 0.90eV (2 \sigma) for the case where power-law prior is adopted to the primordial spectral shape.

We probe the time dependence of the dark energy equation of state (EOS) in light of three-year WMAP (WMAP3) and the combination with other tentative cosmological observations from galaxy clustering (SDSS) and Type Ia Supernova (SNIa). We mainly focus on cases where the EOS is oscillating or with local bumps. By performing a global analysis with the Markov Chain Monte Carlo (MCMC) method, we find the current observations, in particular the WMAP3 + SDSS data combination, allow large oscillations of the EOS which can leave oscillating features on the (residual) Hubble diagram, and such oscillations are potentially detectable by future observations like SNAP, or even by the CURRENTLY ONGOING SNIa observations. Local bumps of dark energy EOS can also leave imprints on CMB, LSS and SNIa. In cases where the bumps take place at low redshifts and the effective EOS are close to -1, CMB and LSS observations cannot give constraints on such possibilities. However, geometrical observations like (future) SNIa can possibly detect such features. On the other hand when the local bumps take place at higher redshifts beyond the detectability of SNIa, future precise observations like Gamma-ray bursts, CMB and LSS may possibly detect such features. In particular, we find that bump-like dark energy EOS on high redshifts might be responsible for the features of WMAP on ranges l \sim 30-50, which is interesting and deserves addressing further.

It is commonly believed that our Universe has experienced two different stages of accelerated expansion. The early stage is known as inflation and the current acceleration is driven by dark energy. Observing inflation and dark energy dynamics are among the most important aspects of the current cosmological research. In this paper we make a first detailed probe of possible degeneracies between dynamical inflation and dark energy in light of current cosmological observations. We find the inclusion of inflation and dark energy parameters together make the parameter spaces broader, but interestingly deviation from scale-invariant primordial spectrum is sightly more favored by current cosmological observations when one marginalizes over dynamical dark energy models.

We probe the scale dependence of the primordial spectrum in the light of the three-year WMAP (WMAP3) alone and WMAP3 in combination with the other cosmological observations such as galaxy clustering and Type Ia Supernova (SNIa). We pay particular attention to the combination with the Lyman $\alpha$ (Ly$\alpha$) forest. Different from the first-year WMAP (WMAP1), WMAP3's preference on the running of the scalar spectral index on the large scales is now fairly independent of the low CMB multipoles $\ell$. A combination with the galaxy power spectrum from the Sloan Digital Sky Survey (SDSS) prefers a negative running to larger than 2$\sigma$, regardless the presence of low $\ell$ CMB ($2\le \ell \le 23$) or not. On the other hand if we focus on the Power Law $\Lambda$CDM cosmology with only six parameters (matter density $\Omega_m h^2$, baryon density $\Omega_b h^2$, Hubble Constant $H_0$, optical depth $\tau$, the spectral index, $n_s$, and the amplitude, $A_s$, of the scalar perturbation spectrum) when we drop the low $\ell$ CMB contributions WMAP3 is consistent with the Harrison-Zel'dovich-Peebles scale-invariant spectrum ($n_s=1$ and no tensor contributions) at $\sim 1\sigma$. When assuming a simple power law primordial spectral index or a constant running, in case one drops the low $\ell$ contributions ($2\le \ell \le 23$) WMAP3 is consistent with the other observations better, such as the inferred value of $\sigma_8$. We also find, using a spectral shape with a minimal extension of the running spectral index model, LUQAS$+$ CROFT Ly$\alpha$ and SDSS Ly$\alpha$ exhibit somewhat different preference on the spectral shape.

We search for signatures of Lorentz and $CPT$ violations in the cosmic microwave background (CMB) temperature and polarization anisotropies by using the WMAP and the 2003 flight of BOOMERANG (B03) data. We note that if the Lorentz and $CPT$ symmetries are broken by a Chern-Simons term in the effective lagrangian, which couples the dual electromagnetic field strength tensor to an external four-vector, the polarization vectors of propagating CMB photons will get rotated. Using the WMAP data alone, one could put an interesting constraint on the size of such a term. Combined with the B03 data, we found that a nonzero rotation angle of the photons is mildly favored: $\Delta \alpha= -6.0^{+4.0}_{-4.0}$ $^{+3.9}_{-3.7}$ deg (1, 2 $\sigma$).

In this paper we use the recently released Type Ia Supernova (SNIa) data to constrain the interactions between the neutrinos and the dark energy scalar fields. In the analysis we take the dark energy scalars to be either Quintessence-like or Phantom-like. Our results show the data mildly favor a model where the neutrinos couple to a phantom-like dark energy scalar, which implies the equation of state of the coupled system behaves like Quintom scenario in the sense of parameter degeneracy. We find future observations like SNAP are potentially promising to measure the couplings between neutrino and dark energy.

Using the Markov chain Monte Carlo (MCMC) method we perform a global analysis constraining the dynamics of dark energy in light of the supernova (Riess "Gold" samples), galaxy clustering (SDSS 3D power spectra and SDSS lyman-\alpha forest information) and the latest three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations. We have allowed the dark energy equation of state to get across -1 and pay particular attention to the effects when incorrectly neglecting dark energy perturbations. We find the parameter space of dynamical dark energy is now well constrained and neglecting dark energy perturbations will make the parameter space significantly smaller. Dynamical dark energy model where the equation of state crosses -1 is mildly favored and the standard \LambdaCDM model is still a good fit to the current data.

Observing dark energy dynamics is the most important aspect of the dark energy research. In this paper we perform an analysis of the constraints on the property of dark energy from the current astronomical observations. We pay particular attention to the effects of dark energy perturbations. Using the data from SNIa (157 "gold" sample), WMAP and SDSS we find that the best fitting dark energy model is given by the dynamical model with the equation of state across -1. Nevertheless the standard \LambdaCDM models are still a very good fit to the current data and evidence for dynamics is currently not very strong. We also consider the constraints from SNLS.

We study in this paper the perturbations of the quintom dark energy model and the effects of quintom perturbations on the current observations. Quintom describes a scenario of dark energy where the equation of state gets across the cosmological constant boundary $w = -1$ during evolution. We present a new method to show that the conventional dark energy models based on single k-essence field and perfect fluid cannot act as quintom due to the singularities and classical instabilities of perturbations around $w = -1$. One needs to add extra degrees of freedom for successful quintom model buildings. There are no singularities or classical instabilities in perturbations of realistic quintom models and they are potentially distinguishable from the cosmological constant. Basing on the realistic quintom models in this paper we provide one way to include the perturbations for dark energy models with parametrized equation of state across -1. Compare with those assuming no dark energy perturbations, we find that the parameter space which allows the equation of state to get across -1 will be enlarged in general when including the perturbations.

We consider in this paper a simple oscillating Quintom model of dark energy which has two free parameters and an equation of state oscillating and crossing -1. For low redshifts the equation of state of this model resembles itself similar to the linearly parameterized dark energy, however differ substantially at large redshifts. We fit our model to the observational data separately from the new high redshift supernova observations from the HST/GOODS program and previous supernova, CMB and galaxy clustering. Our results show that because of the oscillating feature of our model the constraints from observations at large redshifts such as CMB become less stringent.