Using redshift space distortion data, we perform model-independent reconstructions of the growth history of matter inhomogeneity in the expanding Universe using two methods: crossing statistics and Gaussian processes. We then reconstruct the corresponding history of the Universe background expansion and fit it to type Ia supernovae data, putting constraints on $(\Omega_\mathrm{m,0},\sigma_{8,0})$. The results obtained are consistent with the concordance flat-$\Lambda$CDM model and General Relativity as the gravity theory given the current quality of the inhomogeneity growth data.

We explore a class of primordial power spectra that can fit the observed anisotropies in the cosmic microwave background well and that predicts a value for the Hubble parameter consistent with the local measurement of $H_0 = 74$ km/s/Mpc. This class of primordial power spectrum consists of a continuous deformation between the best-fit power law primordial power spectrum and the primordial power spectrum derived from the modified Richardson-Lucy deconvolution algorithm applied to the $C_\ell$s of best-fit power law primordial power spectrum. We find that linear interpolation half-way between the power law and modified Richardson-Lucy power spectra fits the Planck data better than the best-fit $\Lambda$CDM by $\Delta\log\mathcal{L} = 2.5$. In effect, this class of deformations of the primordial power spectra offer a new dimension which is correlated with the Hubble parameter. This correlation causes the best-fit value for $H_0$ to shift and the uncertainty to expand to $H_0 = 70.2 \pm 1.2$ km/s/Mpc. When considering the Planck dataset combined with the Cepheid $H_0$ measurement, the best-fit $H_0$ becomes $H_0 = 71.8 \pm 0.9$ km/s/Mpc. We also compute a Bayes factor of $\log K = 5.7$ in favor of the deformation model.

With Planck CMB observations, we established the spectral amplitude and tilt of the primordial power spectrum. Evidence of a red spectral tilt ($n_\mathrm{s}=0.96$) at $8\sigma$ provides strong support for the inflationary mechanism, especially the slow-roll of the effective scalar field in its nearly flat potential as the generator of scalar primordial perturbations. With the next generation of Large Scale Structure surveys, we expect to probe primordial physics beyond the overall shape and amplitude of the main, smooth and slowly-changing part of the inflaton potential. Using the specifications for the upcoming Euclid survey, we investigate to what extent we can constrain the inflation potential beyond its established slow-roll behaviour. We provide robust forecasts with Euclid and Planck mock data from 9 fiducial power spectra that contain suppression and wiggles at different cosmological scales, using the Wiggly Whipped Inflation framework to generate these features in the primordial spectrum. We include both Euclid cosmic shear and galaxy clustering, with a conservative cut-off for non-linear scales. Using MCMC simulations, we obtain an improvement in constraints in the WWI potential, as well an improvement for the background cosmology parameters. We find that apart from improving the constraints on the overall scale of the inflationary potential by 40-50%, we can also identify oscillations in the primordial spectrum that are present within intermediate to small scales ($k\sim0.01-0.2\,\mathrm{Mpc^{-1}}$).

We determine the Hubble constant $H_0$ precisely ($2.3\%$ uncertainty) in a manner independent of cosmological model through Gaussian process regression, using strong lensing and supernova data. Strong gravitational lensing of a variable source can provide a time-delay distance $D_{\Delta t}$ and angular diameter distance to the lens $D_{\rm{d}}$. These absolute distances can anchor Type Ia supernovae, which give an excellent constraint on the shape of the distance-redshift relation. Updating our previous results to use the H0LiCOW program's milestone dataset consisting of six lenses, four of which have both $D_{\Delta t}$ and $D_{\rm{d}}$ measurements, we obtain $H_0=72.8_{-1.7}^{+1.6}\rm{\ km/s/Mpc}$ for a flat universe and $H_0=77.3_{-3.0}^{+2.2}\rm{\ km/s/Mpc}$ for a non-flat universe. We carry out several consistency checks on the data and find no statistically significant tensions, though a noticeable redshift dependence persists in a particular systematic manner that we investigate. Speculating on the possibility that this trend of derived Hubble constant with lens distance is physical, we show how this can arise through modified gravity light propagation, which would also impact the weak lensing $\sigma_8$ tension.

Consistent observations indicate that some of the important cosmological parameters measured through the local observations are in huge tension with their measurements from the global observations (within the minimal $\Lambda$CDM cosmology). The tensions in those cosmological parameters have been found to be either weakened or reconciled with the introduction of new degrees of freedom that effectively increases the underlying parameter space compared to the minimal $\Lambda$CDM cosmology. It might be interesting to investigate the above tensions within the context of an emergent dark energy scenario proposed recently by Li and Shafieloo~\cite{Li:2019yem}. We find that the tension on $H_0$ is clearly alleviated within 68\% confidence level with an improvement of the $\chi^2$ for CMB, for the above emergent dark energy model having only six free parameters similar to the spatially flat $\Lambda$CDM model. The tension on $H_0$ is still alleviated for every combined datasets considered in the work, however, such alleviation occurs by worsening the $\chi^2$ compared to the $\chi^2$ for $\Lambda$CDM model obtained for the same combined dataset.

In this paper, we quantify the non-linear effects from $k$-essence dark energy through an effective parameter $\mu$ that encodes the additional contribution of a dark energy fluid or a modification of gravity to the Poisson equation. This is a first step toward quantifying non-linear effects of dark energy/modified gravity models in a more general approach. We compare our $N$-body simulation results from $k$-evolution with predictions from the linear Boltzmann code $\texttt{CLASS}$, and we show that for the $k$-essence model one can safely neglect the difference between the two potentials, $ \Phi -\Psi$, and short wave corrections appearing as higher order terms in the Poisson equation, which allows us to use single parameter $\mu$ for characterizing this model. We also show that for a large $k$-essence speed of sound the $\texttt{CLASS}$ results are sufficiently accurate, while for a low speed of sound non-linearities in matter and in the $k$-essence field are non-negligible. We propose a $\tanh$-based parameterisation for $\mu$, motivated by the results for two cases with low ($c_s^2=10^{-7}$) and high ($c_s^2=10^{-4}$) speed of sound, to include the non-linear effects based on the simulation results. This parametric form of $\mu$ can be used to improve Fisher forecasts or Newtonian $N$-body simulations for $k$-essence models.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

In this paper, we propose a new approach to determining cosmological distances to active galactic nuclei (AGN) via light travel-time arguments, which can be extended from nearby sources to very high redshift sources. The key assumption is that the variability seen in AGN is constrained by the speed of light and therefore provides an estimate of the linear size of an emitting region. This can then be compared with the angular size measured with very long baseline interferometry (VLBI) in order to derive a distance. We demonstrate this approach on a specific well studied low redshift (z = 0.0178) source 3C84 (NGC 1275), which is the bright radio core of the Perseus Cluster. We derive an angular diameter distance including statistical errors of $D_{A} = 72^{+5}_{-6}$ Mpc for this source, which is consistent with other distance measurements at this redshift. Possible sources of systematic errors and ways to correct for them are discussed.

We reconstruct the expansion history of the Universe using type Ia supernovae (SN Ia) in a manner independent of any cosmological model assumptions. To do so, we implement a non-parametric iterative smoothing method on the Joint Light-curve Analysis (JLA) data while exploring the SN Ia light-curve hyper-parameter space by Monte Carlo Markov Chain sampling. We test to see how the posteriors of these hyper-parameters depend on cosmology, whether using different dark energy models or reconstructions shift these posteriors. Our constraints on the SN Ia light-curve hyper-parameters from our model-independent analysis are very consistent with the constraints from using different parametrizations of the equation of state of dark energy, namely the flat $\Lambda\mathrm{CDM}$ cosmology, the Chevallier-Polarski-Linder (CPL) model, and the Phenomenologically Emergent Dark Energy (PEDE) model. This implies that the distance moduli constructed from the JLA data are mostly independent of the cosmological models. We also studied that the possibility the light-curve parameters evolve with redshift and our results show consistency with no evolution. The reconstructed expansion history of the Universe and dark energy properties also seem to be in good agreement with the expectations of the standard $\Lambda\mathrm{CDM}$ model. However, our results also indicate that the data still allow for considerable flexibility in the expansion history of the Universe.

Accurate estimation of the Hubble constant, and other cosmological parameters, from distances measured by cosmic gravitational wave sirens requires sufficient allowance for the dark energy evolution. We demonstrate how model independent statistical methods, specifically Gaussian process regression, can remove bias in the reconstruction of $H(z)$, and can be combined model independently with supernova distances. This allows stringent tests of both $H_0$ and $\Lambda$CDM, and can detect unrecognized systematics. We also quantify the redshift systematic control necessary for the use of dark sirens, showing that it must approach spectroscopic precision to avoid significant bias.

We introduce a generalised form of an emergent dark energy model with one degree of freedom for the dark energy sector that has the flexibility to include both $\Lambda$CDM as well as the PEDE model (Phenomenologically Emergent Dark Energy) proposed by Li & Shafieloo (2019) as two of its special limits. This allows us to compare statistically these models in a straightforward way and following conventional Bayesian approach. The free parameter for the dark energy sector, namely $\Delta$, has the value of $0$ for the case of the $\Lambda$ and $1$ for the case of PEDE and its posterior fitting the generalised parametric form to different data can directly show the consistency of the models to the data. Fitting the introduced parametric form to Planck CMB data and most recent $H_0$ results from local observations of Cepheids and Supernovae (Riess et al. 2019), we show that the $\Delta=0$ associated with the $\Lambda$CDM model would fall out of 4$\sigma$ confidence limits of the derived posterior of the $\Delta$ parameter. In contrast, PEDE model can satisfy the combination of the observations. This is another support for the case of PEDE model with respect to the standard $\Lambda$CDM model if we trust the reliability of both Planck CMB data and local $H_0$ observations.

Motivated by the current status of the cosmological observations and significant tensions in the estimated values of some key parameters assuming the standard $\Lambda$CDM model, we propose a simple but radical phenomenological emergent dark energy model where dark energy has no effective presence in the past and emerges at the later times. Theoretically, in this phenomenological dark energy model with zero degree of freedom (similar to a $\Lambda$CDM model), one can derive that the equation of state of dark energy increases from $-\frac{2}{3 {\rm{ln}}\, 10} -1$ in the past to $-1$ in the future. We show that by setting a hard-cut 2$\sigma$ lower bound prior for the $H_0$ that associates with $97.72\%$ probability from the recent local observations~\citep{riess2019large}, this model can satisfy different combinations of cosmological observations at low and high redshifts (SNe Ia, BAO, Ly{$\alpha$} BAO and CMB) substantially better than the concordance $\Lambda$CDM model with $\Delta \chi^2_{bf} \sim -41.08$ and $\Delta\,{\rm{DIC}} \sim-35.38$. If there are no substantial systematics in SN Ia, BAO or Planck CMB data and assuming reliability of the current local $H_0$ measurements, there is a very high probability that with slightly more precise measurement of the Hubble constant our proposed phenomenological model rules out the cosmological constant with decisive statistical significance and is a strong alternative to explain combination of different cosmological observations. This simple phenomenologically emergent dark energy model can guide theoretically motivated dark energy model building activities.

We investigate constraints on some key cosmological parameters by confronting metastable dark energy models with different combinations of the most recent cosmological observations. Along with the standard $\Lambda$CDM model, two phenomenological metastable dark energy models are considered: (\romannumeral1) DE decays exponentially, (\romannumeral2) DE decays into dark matter. We find that: (1) when considering the most recent supernovae and BAO data, and assuming a fiducial $\Lambda$CDM model, the inconsistency in the estimated value of the $\Omega_{\rm{m,0}}h^2$ parameter obtained by either including or excluding Planck CMB data becomes very much substantial and points to a clear tension~\citep{sahni2014model,zhao2017dynamical}; (2) although the two metastable dark energy models that we study provide greater flexibility in fitting the data, and they indeed fit the SNe Ia+BAO data substantially better than $\Lambda$CDM, they are not able to alleviate this tension significantly when CMB data are included; (3) while local measurements of the Hubble constant are significantly higher relative to the estimated value of $H_0$ in our models (obtained by fitting to SNe Ia and BAO data), the situation seems to be rather complicated with hints of inconsistency among different observational data sets (CMB, SNe Ia+BAO and local $H_0$ measurements). Our results indicate that we might not be able to remove the current tensions among different cosmological observations by considering simple modifications of the standard model or by introducing minimal dark energy models. A complicated form of expansion history, different systematics in different data and/or a non-conventional model of the early Universe might be responsible for these tensions.

Strongly lensed quasar systems with time delay measurements provide "time delay distances", which are a combination of three angular diameter distances and serve as powerful tools to determine the Hubble constant $H_0$. However, current results often rely on the assumption of the $\Lambda$CDM model. Here we use a model-independent method based on Gaussian process to directly constrain the value of $H_0$. By using Gaussian process regression, we can generate posterior samples of unanchored supernova distances independent of any cosmological model and anchor them with strong lens systems. The combination of a supernova sample with large statistics but no sensitivity to $H_0$ with a strong lens sample with small statistics but $H_0$ sensitivity gives a precise $H_0$ measurement without the assumption of any cosmological model. We use four well-analyzed lensing systems from the state-of-art lensing program H0LiCOW and the Pantheon supernova compilation in our analysis. Assuming the Universe is flat, we derive the constraint $H_0=72.2 \pm 2.1\,$km/s/Mpc, a precision of $2.9\%$. Allowing for cosmic curvature with a prior of $\Omega_{k}=[-0.2,0.2]$, the constraint becomes $H_0=73.0_{-3.0}^{+2.8}\,$km/s/Mpc.

Previously we proposed a novel method to inspect the isotropy of the properties of gamma-ray bursts (GRBs) such as their duration, fluences and peak fluxes at various energy bands and different time scales, complementary to existing studies of spatial distribution of GRBs by other authors. The method was then applied on the Fermi GBM Burst Catalog containing 1591 GRBs and except one particular direction where we noticed some hints of violation from statistical isotropy, the rest of the data showed consistency with isotropy. In this work we apply our method with some minor modifications to the updated Fermi/GBM data sample containing 2266 GRBs, thus $\sim 40$ % larger. We also test two other major GRB catalogs, the BATSE Current GRB Catalog of the CGRO satellite containing $\sim 2000$ bursts and the Swift/BAT Gamma-Ray Burst Catalog containing $\sim 1200$ bursts. The new results using the updated data are consistent with our previous findings and no statistically significant anisotropic feature in the observed properties of these samples of all GRBs is found.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

During the next decade, gravitational waves will be observed from hundreds of binary inspiral events. When the redshifts of the host galaxies are known, these events can be used as `standard sirens', sensitive to the expansion rate of the Universe. Measurements of the Hubble constant $H_0$ from standard sirens can be done independently of other cosmological probes, and events occurring at $z<0.1$ will allow us to infer $H_0$ independently of cosmological models. The next generation of spectroscopic galaxy surveys will play a crucial role in reducing systematic uncertainties in $H_0$ from standard sirens, particularly for the numerous `dark sirens' which do not have an electromagnetic counterpart. In combination with large spectroscopic data sets, standard sirens with an EM counterpart are expected to constrain $H_0$ to $\sim 1-2\%$ precision within the next decade. This is competitive with the best estimates of $H_0$ obtained to date and will help illuminate the current tension between existing measurements. Information on the galaxies that host the gravitational wave events will also shed light on the origin and evolution of compact object binaries.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We examine the Pantheon supernovae distance data compilation in a model independent analysis to test the validity of cosmic history reconstructions beyond the concordance $\Lambda$CDM cosmology. Strong deviations are allowed by the data at $z\gtrsim1$ in the reconstructed Hubble parameter, $Om$ diagnostic, and dark energy equation of state. We explore three interpretations: 1) possibility of the true cosmology being far from $\Lambda$CDM, 2) supernovae property evolution, and 3) survey selection effects. The strong (and theoretically problematic) deviations at $z\gtrsim1$ vanish and good consistency with $\Lambda$CDM is found with a simple Malmquist-like linear correction. The adjusted data is robust against the model independent iterative smoothing reconstruction. However, we caution that while by eye the original deviation from $\Lambda$CDM is striking, $\chi^2$ tests do not show the extra linear correction parameter is statistically significant, and a model-independent Gaussian Process regression does not find significant evidence for the need for correction at high-redshifts.

Lack of knowledge about the background expansion history of the Universe from independent observations makes it problematic to obtain a precise and accurate estimation of the Hubble constant $H_0$ from gravitational wave standard sirens, even with electromagnetic counterpart redshifts. Simply fitting simultaneously for the matter density in a flat \lcdm\ model can reduce the precision on $H_0$ from 1\% to 5\%, while not knowing the actual background expansion model of the universe (e.g.\ form of dark energy) can introduce substantial bias in estimation of the Hubble constant. When the statistical precision is at the level of 1\% uncertainty on $H_0$, biases in non-\lcdm\ cosmologies that are consistent with current data could reach the 3$\sigma$ level. To avoid model-dependent biases, statistical techniques that are appropriately agnostic about model assumptions need to be employed.

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We discuss the discordance between the estimated values of the cosmological parameters from Planck assuming the concordance $\Lambda$CDM model and low-redshift measurements. In particular, we consider the Hubble constant mismatch between Planck temperature constraint for the $\Lambda$CDM model and the Riess et. al. local measurements as well as the discordance between the estimated value of $S_8$ from Planck and some weak lensing surveys such as Kilo Degree Survey (KiDS-450) and Dark Energy Survey (DES) observations. The discordance can come from a wide range of non-standard cosmological or astrophysical processes as well as from some particular systematics of the observations. In this paper, without considering any particular astrophysical process or extension to the standard model at the background level, we seek solely to project the effect of these differences in the values of the key cosmological parameters on to the shape of the primordial power spectrum (PPS). In order to realise this goal, we uncover the shape of the PPS by implementing the Modified Richardson-Lucy algorithm (MRL) that fits the Planck temperature data as acceptably as the case of the standard model of cosmology, but with a Hubble constant consistent with local measurements as well as improving the consistency between the derived $S_8$ and $\sigma_8$ parameters with estimations of the weak lensing surveys.

We combine model-independent reconstructions of the expansion history from the latest Pantheon supernovae distance modulus compilation and measurements from baryon acoustic oscillation to test some important aspects of the concordance model of cosmology namely the FLRW metric and flatness of spatial curvature. We then use the reconstructed expansion histories to fit growth measurement from redshift-space distortion and obtain strong constraints on $(\Omega_\mathrm{m},\gamma,\sigma_8)$ in a model independent manner. Our results show consistency with a spatially flat FLRW Universe with general relativity to govern the perturbation in the structure formation and the cosmological constant as dark energy. However, we can also see some hints of tension among different observations within the context of the concordance model related to high redshift observations ($z > 1$) of the expansion history. This supports earlier findings of Sahni et al. (2014) & Zhao et al. (2017) and highlights the importance of precise measurement of expansion history and growth of structure at high redshifts.

The form of the primordial power spectrum (PPS) of cosmological scalar (matter density) perturbations is not yet constrained satisfactorily in spite of the tremendous amount of information from the Cosmic Microwave Background (CMB) data. While a smooth power-law-like form of the PPS is consistent with the CMB data, some PPS with small non-smooth features at large scales can also fit the CMB temperature and polarization data with similar statistical evidence. Future CMB surveys cannot help distinguish all such models due to the cosmic variance at large angular scales. In this paper, we study how well we can differentiate be- tween such featured forms of the PPS not otherwise distinguishable using CMB data. We ran 15 N-body DESI-like simulations of these models to explore this approach. Showing that statistics such as the halo mass function and the two-point correlation function are not able to distinguish these models in a DESI-like survey, we advocate to avoid reducing the dimensionality of the problem by demonstrating that the use of a simple three-dimensional count-in-cell density field can be much more effective for the purpose of model distinction.

Reconstructing the expansion history of the Universe from type Ia supernovae data, we fit the growth rate measurements and put model-independent constraints on some key cosmological parameters, namely, $\Omega_\mathrm{m},\gamma$, and $\sigma_8$. The constraints are consistent with those from the concordance model within the framework of general relativity, but the current quality of the data is not sufficient to rule out modified gravity models. Adding the condition that dark energy density should be positive at all redshifts, independently of its equation of state, further constrains the parameters and interestingly supports the concordance model.

We consider the prospects of probing features in the primordial power spectrum with future Cosmic Microwave Background (CMB) polarization measurements. In the scope of the inflationary scenario, such features in the spectrum can be produced by local non-smooth pieces in an inflaton potential (smooth and quasi-flat in general) which in turn may originate from fast phase transitions during inflation in other quantum fields interacting with the inflaton. They can fit some outliers in the CMB temperature power spectrum which are unaddressed within the standard inflationary ${\mathrm{\Lambda}}$CDM model. We consider Wiggly Whipped Inflation (WWI) as a theoretical framework leading to improvements in the fit to the Planck 2015 temperature and polarization data in comparison with the standard inflationary models, although not at a statistically significant level. We show that some type of features in the potential within the WWI models, leading to oscillations in the primordial power spectrum that extend to intermediate and small scales can be constrained with high confidence (at 3$\sigma$ or higher confidence level) by an instrument as the Cosmic ORigins Explorer (CORE). In order to investigate the possible confusion between inflationary features and footprints from the reionization era, we consider an extended reionization history with monotonic increase of free electrons with decrease in redshift. We discuss the present constraints on this model of extended reionization and future predictions with CORE. We also project, to what extent, this extended reionization can create confusion in identifying inflationary features in the data.

We perform an anisotropic clustering analysis of 1,133,326 galaxies from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR) 12 covering the redshift range $0.15<z<0.69$. The geometrical distortions of the galaxy positions, caused by incorrect cosmological model assumptions, are captured in the anisotropic two-point correlation function on scales 6 -- 40 $h^{-1}\rm Mpc$. The redshift evolution of this anisotropic clustering is used to place constraints on the cosmological parameters. We improve the methodology of Li et al. 2016, to enable efficient exploration of high dimensional cosmological parameter spaces, and apply it to the Chevallier-Polarski-Linder parametrization of dark energy, $w=w_0+w_a{z}/({1+z})$. In combination with the CMB, BAO, SNIa and $H_0$ from Cepheid data, we obtain $\Omega_m = 0.301 \pm 0.008,\ w_0 = -1.042 \pm 0.067,\ $ and $w_a = -0.07 \pm 0.29$ (68.3\% CL). Adding our new AP measurements to the aforementioned results reduces the error bars by $\sim$30 -- 40\% and improves the dark energy figure of merit by a factor of $\sim$2. We check the robustness of the results using realistic mock galaxy catalogues.

Cosmic spatial curvature is a fundamental geometric quantity of the Universe. We investigate a model independent, geometric approach to measure spatial curvature directly from observations, without any derivatives of data. This employs strong lensing time delays and supernova distance measurements to measure the curvature itself, rather than just testing consistency with flatness. We define two curvature estimators, with differing error propagation characteristics, that can crosscheck each other, and also show how they can be used to map the curvature in redshift slices, to test constancy of curvature as required by the Robertson-Walker metric. Simulating realizations of redshift distributions and distance measurements of lenses and sources, we estimate uncertainties on the curvature enabled by next generation measurements. The results indicate that the model independent methods, using only geometry without assuming forms for the energy density constituents, can determine the curvature at the $\sim6\times10^{-3}$ level.

We propose a new class of metastable dark energy (DE) phenomenological models in which the DE decay rate does not depend on external parameters such as the scale factor or the curvature of the Universe. Instead, the DE decay rate is assumed to be a constant depending only on intrinsic properties of DE and the type of a decay channel, similar to case of the radioactive decay of unstable particles and nuclei. As a consequence, the DE energy density becomes a function of the proper time elapsed since its formation, presumably in the very early Universe. Such a natural type of DE decay can profoundly affect the expansion history of the Universe and its age. Metastable DE can decay in three distinct ways: (i) exponentially, (ii) into dark matter, (iii) into dark radiation. Testing metastable DE models with observational data we find that the decay half-life must be many times larger than the age of the Universe. Models in which dark energy decays into dark matter lead to lower values of the Hubble parameter at large redshifts relative to $\Lambda$CDM. Consequently these models provide a better fit to cosmological BAO data (especially data from from high redshift quasars) than concordance ($\Lambda$CDM) cosmology.

We present a measurement of the anisotropic and isotropic Baryon Acoustic Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample with optimal redshift weights. Applying the redshift weights improves the constraint on the BAO dilation parameter $\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

The sky distribution of Gamma-Ray Bursts (GRBs) has been intensively studied by various groups for more than two decades. Most of these studies test the isotropy of GRBs based on their sky number density distribution. In this work we propose an approach to test the isotropy of the Universe through inspecting the isotropy of the properties of GRBs such as their duration, fluences and peak fluxes at various energy bands and different time scales. We apply this method on the {\em Fermi} / Gamma-ray Burst Monitor (GBM) data sample containing 1591 GRBs. The most noticeable feature we found is near the Galactic coordinates $l\approx 30^\circ$, $b\approx 15^\circ$ and radius $r\approx 20^\circ-40^\circ$. The inferred probability for the occurrence of such an anisotropic signal (in a random isotropic sample) is derived to be less than a percent in some of the tests while the other tests give results consistent with isotropy. These are based on the comparison of the results from the real data with the randomly shuffled data samples. Considering large number of statistics we used in this work (which some of them are correlated to each other) we can anticipate that the detected feature could be result of statistical fluctuations. Moreover, we noticed a considerably low number of GRBs in this particular patch which might be due to some instrumentation or observational effects that can consequently affect our statistics through some systematics. Further investigation is highly desirable in order clarify about this result, e.g. utilizing a larger future {\em Fermi} / GBM data sample as well as data samples of other GRB missions and also looking for possible systematics.

A persistent signal of power asymmetry on opposite hemispheres of CMB sky was seen in full-sky temperature measurements made so far. This asymmetry was seen in microwave sky from WMAP as well as PLANCK satellites, and calls for attention the larger question of \emph{statistical isotropy}, one of the foundational principles of modern cosmology. In this work we present an analysis of polarized CMB maps from PLANCK 2015 full mission data. We apply the local variance estimator on low resolution $E-$mode maps from PLANCK 2015 polarization \texttt{Commander} solution. We find a significant hemispherical power asymmetry in polarization data on large angular scales, at the level of $\sim 2.6-3.9\%$ depending on the galactic mask, and the circular disc radius used for computing local variance maps. However the direction is found to be pointing broadly towards CMB kinetic dipole direction. Precise measurements of CMB polarization in future will shed light on this apparent discrepancy in the anisotropy axis seen in temperature and polarized CMB sky, and likely influence of systematics on our findings.

Non-parametric reconstruction methods, such as Gaussian process (GP) regression, provide a model-independent way of estimating an underlying function and its uncertainty from noisy data. We demonstrate how GP-reconstruction can be used as a consistency test between a given data set and a specific model by looking for structures in the residuals of the data with respect to the model's best-fit. Applying this formalism to the Planck temperature and polarisation power spectrum measurements, we test their global consistency with the predictions of the base $\Lambda$CDM model. Our results do not show any serious inconsistencies, lending further support to the interpretation of the base $\Lambda$CDM model as cosmology's gold standard.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

We test the consistency between Planck temperature and polarization power spectra and the concordance model of $\Lambda$ Cold Dark Matter cosmology ($\Lambda$CDM) within the framework of Crossing statistics. We find that Planck TT best fit $\Lambda$CDM power spectrum is completely consistent with EE power spectrum data while EE best fit $\Lambda$CDM power spectrum is not consistent with TT data. However, this does not point to any systematic or model-data discrepancy since in the Planck EE data, uncertainties are much larger compared to the TT data. We also investigate the possibility of any deviation from $\Lambda$CDM model analyzing the Planck 2015 data. Results from TT, TE and EE data analysis indicate that no deviation is required beyond the flexibility of the concordance $\Lambda$CDM model. Our analysis thus rules out any strong evidence for beyond the concordance model in the Planck spectra data. We also report a mild amplitude difference comparing temperature and polarization data, where temperature data seems to have slightly lower amplitude than expected (consistently at all multiples), as we assume both temperature and polarization data are realizations of the same underlying cosmology.

We modify the algorithm we proposed before in Aghamousa & Shafieloo (2015) on time delay estimation of the strong lens systems incorporating weighted cross-correlation and weighted summation of correlation coefficients. We show the high performance of this algorithm by applying it on Time Delay Challenge (TDC1) simulated data. We apply then our proposed method on the light curves of the lensed quasar SDSS J1001+5027 since this system has been well studied by other groups to compare our results with their findings. In this work we propose a new estimator namely "mirror" estimator along with a list of criteria for reliability test of estimation. Our mirror estimator results to $-117.1^{+7.1}_{-3.7}$ and $-117.1^{+7.2}_{-8.8}$ using simple Monte Carlo simulations and simulated light curves provided by Rathna Kumar et al. (2013) respectively. Although TDC1 simulations do not reflect the properties of SDSS J1001+5027 light curves, using these simulations results to smaller uncertainty which shows the higher quality observations can lead to substantially more precise time delay estimation. Our time delay estimation is in agreement with findings of the other groups for this strong lens system and the difference in the size of the error bars reflects the importance of appropriate light curve simulations.

We re-analyze the nonlocal gravity model of Deser and Woodard which was proposed to account for the current phase of cosmic acceleration. We show that the growth of perturbations predicted by this nonlocal gravity model when its background evolution is fixed by some particular non-$\Lambda$CDM models can be substantially lower than when its background is fixed by $\Lambda$CDM. This can be seen when we consider the background expansion by a dark energy model with a slightly less negative equation of state with respect to cosmological constant. Our results hints towards a fact that the choice of the background expansion can play a crucial role how this nonlocal gravity model can fit the growth history data. While the growth data might show better consistency to GR models (among the background models we studied so far), it seems the nonlocal gravity model studied in this work is able to show comparable consistency to the growth data as well. Showing this consistency can be considered as a significant result since this model can stand as a distinguishable alternative to the standard model of cosmology.

Using measurements of $H(z)$ and $d_\mathrm{A}(z)$ from the Baryon Oscillation Spectroscopic Survey (BOSS) DR12 and luminosity distances from the Joint Lightcurve Analysis (JLA) compilation of supernovae (SN), we measure $H_0 r_\mathrm{d}$ without any model assumption. Our measurement of $H_0 r_\mathrm{d} = $($10033.20^{+333.10}_{-371.81}$ (SN) $\pm$ 128.19 (BAO)) km s$^{-1}$ is consistent with Planck constrains for the flat {\Lambda}CDM model. We also report that higher expansion history rates $h(z)$ (among the possibilities) as well as lower-bound values of $H_0 r\mathrm{d}$ result in better internal consistency among the independent data (H(z)rd and $d_\mathrm{A}(z)/r_\mathrm{d}$ from BAO at $z=0.32$ and z=0.57 and $d_\mathrm{L}$ from JLA) we used in this work. This can be interpreted as an interesting and independent support of Planck cosmology without using any cosmic microwave background data. We then combine these observables to test the Friedmann-Lema\^itre-Robertson-Walker (FLRW) metric and the flatness of the Universe in a model-independent way at two redshifts, namely 0.32 and 0.57, by introducing a new diagnostic for flat-FLRW, $\Theta(z)$, which only depends on observables of BAO and SN data. Our results are consistent with a flat-FLRW Universe within $2{\sigma}$.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

With the Planck 2015 Cosmic Microwave Background (CMB) temperature and polarization data, we search for possible features in the primordial power spectrum (PPS). We revisit the Wiggly Whipped Inflation (WWI) framework and demonstrate how generation of some particular primordial features can improve the fit to Planck data. WWI potential allows the scalar field to transit from a steeper potential to a nearly flat potential through a discontinuity either in potential or in its derivatives. WWI offers the inflaton potential parametrizations that generate a wide variety of features in the primordial power spectra incorporating most of the localized and non-local inflationary features that are obtained upon reconstruction from temperature and polarization angular power spectrum. At the same time, in a single framework it allows us to have a background parameter estimation with a nearly free-form primordial spectrum. Using Planck 2015 data, we constrain the primordial features in the context of Wiggly Whipped Inflation and present the features that are supported both by temperature and polarization. WWI model provides more than $13$ improvement in $\chi^2$ fit to the data with respect to the best fit power law model considering combined temperature and polarization data from Planck and B-mode polarization data from BICEP and Planck dust map. We use 2-4 extra parameters in the WWI model compared to the featureless strict slow roll inflaton potential. We find that the differences between the temperature and polarization data in constraining background cosmological parameters such as baryon density, cold dark matter density are reduced to a good extent if we use primordial power spectra from WWI. We also discuss the extent of bispectra obtained from the best potentials in arbitrary triangular configurations using the BI-spectra and Non-Gaussianity Operator (BINGO).

We study features in the bispectrum of the primordial curvature perturbation correlated with the reconstructed primordial power spectrum from the observed cosmic microwave background temperature data. We first show how the bispectrum can be completely specified in terms of the power spectrum and its first two derivatives, valid for any configuration of interest. Then using a model-independent reconstruction of the primordial power spectrum from the Planck angular power spectrum of temperature anisotropies, we compute the bispectrum in different triangular configurations. We find that in the squeezed limit at k ~ 0.06/Mpc and k ~ 0.014/Mpc there are marginal 2sigma deviations from the standard featureless bispectrum, which meanwhile is consistent with the reconstructed bispectrum in the equilateral configuration.

We report the first test of isotropy of the Universe in the matter dominated epoch using the Lyman-$\alpha$ forest data from the high redshift quasars ($z>2$) from SDSS-III BOSS-DR9 datasets. Using some specified data cuts, we obtain the probability distribution function (PDF) of the Lyman-$\alpha$ forest transmitted flux and use the statistical moments of the PDF to address the isotropy of the Universe. In an isotropic Universe one would expect the transmitted flux to have consistent statistical characteristics in different parts of the sky. We trisect the total survey area of 3275 ${\rm deg}^2$ along the galactic latitude and using quadrant convention. We also make three subdivisions in the data for three different signal-to-noise-ratios (SNR). Finally we obtain and compare the statistical moments in the mean redshifts of 2.3, 2.6 and 2.9. We find, that the moments from all patches agree at all redshifts and at all SNRs, within 3$\sigma$ uncertainties. Since Lyman-$\alpha$ transmitted flux directly maps the neutral hydrogen distribution in the inter galactic medium (IGM), our results indicate, within the limited survey area and sensitivity of the data, the distribution of the neutral hydrogen in the Universe is consistent with isotropic distribution. We should mention that we report few deviations from isotropy in the data with low statistical significance. Increase in survey area and larger amount of data are needed to make any strong conclusion about these deviations.

We test the isotropy of the local Universe using low redshift supernova data from various catalogs and the non-parametric method of smoothed residuals. Using a recently developed catalog which combines supernova data from various surveys, we show that the isotropic hypothesis of a Universe with zero velocity perturbation can be rejected with moderate significance, with $p$-value $\sim 0.07$ out to redshift $z < 0.045$. We estimate the direction of maximal anisotropy on the sky for various pre-existing catalogs and show that it remains relatively unaffected by the light curve fitting procedure. However the recovered direction is biased by the underlying distribution of data points on the sky. We estimate both the uncertainty and bias in the direction by creating mock data containing a randomly oriented bulk flow and using our method to reconstruct its direction. We conclude that the inhomogeneous nature of the data introduces a directional bias in galactic latitude of approximately $|\Delta b_{\rm max}| \sim 18^{\circ}$ for the supernova catalog considered in this work, and after correcting for this effect we infer the direction of maximum anisotropy as $(b,\ell) = (20^{\circ}, 276^{\circ}) \pm (12^{\circ},29^{\circ})$ in galactic coordinates. Finally we compare the anisotropic signal in the data to mock realisations in which large scale velocity perturbations are consistently accounted for at the level of linear perturbation theory. We show that including the effect of the velocity perturbation in our mock catalogs degrades the significance of the anisotropy considerably, with $p$-value increasing to $\sim 0.29$. One can conclude from our analysis that there is a moderate deviation from isotropy in the supernova data, but the signal is consistent with a large scale bulk velocity expected within $\Lambda$CDM.

The observable time delays between the multiple images of strong lensing systems with time variable sources can provide us with some valuable information to probe the expansion history of the Universe. Estimation of these time delays can be very challenging due to complexities of the observed data where there are seasonal gaps, various noises and systematics such as unknown microlensing effects. In this paper we introduce a novel approach to estimate the time delays for strong lensing systems implementing various statistical methods of data analysis including the method of smoothing and cross-correlation. The method we introduce in this paper has been recently used in TDC0 and TDC1 Strong Lens Time Delay Challenges and has shown its power in reliable and precise estimation of time delays dealing with data with different complexities.

We present a novel approach to test the consistency of the cosmological models with multiband CMB data using a nonparametric approach. In our analysis we calibrate the REACT (Risk Estimation and Adaptation after Coordinate Transformation) confidence levels associated with distances in function space (\textit{confidence distances}) based on the Monte Carlo simulations in order to test the consistency of an assumed cosmological model with observation. To show the applicability of our algorithm, we confront Planck 2013 temperature data with concordance model of cosmology considering two different Planck spectra combination. In order to have an accurate quantitative statistical measure to compare between the data and the theoretical expectations, we calibrate REACT confidence distances and perform a bias control using many realizations of the data. Our results in this work using Planck 2013 temperature data put the best fit $\Lambda$CDM model at $95\% (\sim 2\sigma)$ confidence distance from the center of the nonparametric confidence set while repeating the analysis excluding the Planck $217 \times 217$ GHz spectrum data, the best fit $\Lambda$CDM model shifts to $70\% (\sim 1\sigma)$ confidence distance. The most prominent features in the data deviating from the best fit $\Lambda$CDM model seems to be at low multipoles $ 18 < \ell < 26$ at greater than $2\sigma$, $\ell \sim 750$ at $\sim1$ to $2 \sigma$ and $\ell \sim 1800$ at greater than $2\sigma$ level. Excluding the $217\times217$ GHz spectrum the feature at $\ell \sim 1800$ becomes substantially less significance at $\sim1$ to $2 \sigma$ confidence level. Results of our analysis based on the new approach we propose in this work are in agreement with other analysis done using alternative methods.

Using modified Richardson-Lucy algorithm we reconstruct the primordial power spectrum (PPS) from Planck Cosmic Microwave Background (CMB) temperature anisotropy data. In our analysis we use different combinations of angular power spectra from Planck to reconstruct the shape of the primordial power spectrum and locate possible features. Performing an extensive error analysis we found the dip near $\ell\sim750-850$ represents the most prominent feature in the data. Feature near $\ell\sim1800-2000$ is detectable with high confidence only in 217 GHz spectrum and is apparently consequence of a small systematic as described in the revised Planck 2013 papers. Fixing the background cosmological parameters and the foreground nuisance parameters to their best fit baseline values, we report that the best fit power law primordial power spectrum is consistent with the reconstructed form of the PPS at 2$\sigma$ C.L. of the estimated errors (apart from the local features mentioned above). As a consistency test, we found the reconstructed primordial power spectrum from Planck temperature data can also substantially improve the fit to WMAP-9 angular power spectrum data (with respect to power-law form of the PPS) allowing an overall amplitude shift of $\sim2.5\%$. In this context low-$\ell$ and 100 GHz spectrum from Planck which have proper overlap in the multipole range with WMAP data found to be completely consistent with WMAP-9 (allowing amplitude shift). As another important result of our analysis we do report the evidence of gravitational lensing through the reconstruction analysis. Finally we present two smooth form of the PPS containing only the important features. These smooth forms of PPS can provide significant improvements in fitting the data (with respect to the power law PPS) and can be helpful to give hints for inflationary model building.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

Estimation of the angular power spectrum is one of the important steps in Cosmic Microwave Background (CMB) data analysis. Here, we present a nonparametric estimate of the temperature angular power spectrum for the Planck 2013 CMB data. The method implemented in this work is model-independent, and allows the data, rather than the model, to dictate the fit. Since one of the main targets of our analysis is to test the consistency of the $\Lambda$CDM model with Planck 2013 data, we use the nuisance parameters associated with the best-fit $\Lambda$CDM angular power spectrum to remove foreground contributions from the data at multipoles $\ell \geq50$. We thus obtain a combined angular power spectrum data set together with the full covariance matrix, appropriately weighted over frequency channels. Our subsequent nonparametric analysis resolves six peaks (and five dips) up to $\ell \sim1850$ in the temperature angular power spectrum. We present uncertainties in the peak/dip locations and heights at the $95\%$ confidence level. We further show how these reflect the harmonicity of acoustic peaks, and can be used for acoustic scale estimation. Based on this nonparametric formalism, we found the best-fit $\Lambda$CDM model to be at $36\%$ confidence distance from the center of the nonparametric confidence set -- this is considerably larger than the confidence distance ($9\%$) derived earlier from a similar analysis of the WMAP 7-year data. Another interesting result of our analysis is that at low multipoles, the Planck data do not suggest any upturn, contrary to the expectation based on the integrated Sachs-Wolfe contribution in the best-fit $\Lambda$CDM cosmology.

We test the isotropy of the local distribution of galaxies using the 2MASS extended source catalogue. By decomposing the full sky survey into distinct patches and using a combination of photometric and spectroscopic redshift data, we use both parametric and non-parametric methods to obtain the shape of the luminosity function in each patch. We use the shape of the luminosity function to test the statistical isotropy of the underlying galaxy distribution. The parametric estimator shows some evidence of a hemispherical asymmetry in the north/south Galactic plane. However the non-parametric estimator exhibits no significant anisotropy, with the galaxy distribution being consistent with the assumption of isotropy in all regions considered. The parametric asymmetry is attributed to the relatively poor fit of the functional form to the underlying data. When using the non-parametric estimator, we do find a dipole in the shape of the luminosity function, with maximal deviation from isotropy at galactic coordinate $(b,l)=(30^{\circ},315^{\circ})$. However we can ascribe no strong statistical significance to this observation.

Baryon Acoustic Oscillations (BAO) allow us to determine the expansion history of the Universe, thereby shedding light on the nature of dark energy. Recent observations of BAO's in the SDSS DR9 and DR11 have provided us with statistically independent measurements of $H(z)$ at redshifts of 0.57 and 2.34, respectively. We show that these measurements can be used to test the cosmological constant hypothesis in a model independent manner by means of an improved version of the $Om$ diagnostic. Our results indicate that the SDSS DR11 measurement of $H(z) = 222 \pm 7$ km/sec/Mpc at $z = 2.34$, when taken in tandem with measurements of $H(z)$ at lower redshifts, imply considerable tension with the standard $\Lambda$CDM model. Our estimation of the new diagnostic $Omh^2$ from SDSS DR9 and DR11 data, namely $Omh^2 \approx 0.122 \pm 0.01$, which is equivalent to $\Omega_{0m}h^2$ for the spatially flat $\Lambda$CDM model, is in tension with the value $\Omega_{0m}h^2 = 0.1426 \pm 0.0025$ determined for $\Lambda$CDM from Planck+WP. This tension is alleviated in models in which the cosmological constant was dynamically screened (compensated) in the past. Such evolving dark energy models display a pole in the effective equation of state of dark energy at high redshifts, which emerges as a smoking gun test for these theories.

Motivated by the idea that inflation occurs at the GUT symmetry breaking scale, in this paper we construct a new class of large field inflaton potentials where the inflaton starts with a power law potential; after initial period of relative fast roll that lasts until after a few e-folds inside the horizon, it transits to the attractor of the slow roll part of the potential with a lower power. Due to the initial fast roll stages of inflation, we find a suppression in scalar primordial power at large scales and at the same time the choice of the potential can provide us a tensor primordial spectrum with high amplitude. This suppression in scalar power with a large tensor-to-scalar ratio helps us to reconcile the Planck and BICEP2 data in a single framework. We find that a transition from a cubic to quadratic form of inflaton potential generates an appropriate suppression in power of scalar primordial spectrum that provides significant improvement in fit compared to power law model when compared with Planck and BICEP2 data together. We calculate the extent of non-Gaussianity, specifically, the bispectrum for the best fit potential and show that it is consistent with Planck bispectrum constraints.

Motivated by BICEP2 results on the CMB polarization B-mode which imply primordial gravitational waves are produced when the Universe has the expansion rate of about $H \approx 10^{14}$ GeV, and by deviations from a smooth power-law behaviour for multipoles $\ell <50$ in the CMB temperature anisotropy power spectrum found in the WMAP and Planck experiments, we have expanded our class of large field inflationary models that fit both the BICEP2 and Planck CMB observations consistently. These best-fitted large field models are found to have a transition from a faster roll to the slow roll $V(\phi)=m^2 \phi^2/2$ inflation at a field value around 14.6~${\rm M_{Pl}}$ and thus a potential energy of $V(\phi) \sim (10^{16}\,{\rm GeV})^4$. In general this transition with sharp features in the inflaton potential produces not only suppression of scalars relative to tensor modes at small $k$ but also introduces wiggles in the primordial perturbation spectrum. These wiggles are shown to be useful to explain some localized features in the CMB angular power spectrum and can also have other observational consequences. Thus, primordial GW can be used now to make a tomography of inflation determining its fine structure. The resulting Wiggly Whipped Inflation scenario is described in details and the anticipated perturbation power spectra, CMB power spectra, non-Gaussianity and other observational consequences are calculated and compared to existing and forthcoming observations.

Combining Planck CMB temperature [1] and BICEP2 B-mode polarization data [2,3] we show qualitatively that, assuming inflationary consistency relation, the power-law form of the scalar primordial spectrum is ruled out at more than $3\sigma$ CL. This is an important finding, since the power-law form of the scalar primordial spectrum is one of the main assumptions of concordance model of cosmology and also a direct prediction of many inflationary scenarios. We show that a break or step in the form of the primordial scalar perturbation spectrum, similar to what we studied recently analyzing Planck data [4] can address both Planck and BICEP2 results simultaneously. Our findings also indicate that the data may require more flexibilities than what running of scalar spectral index can provide. Finally we show that an inflaton potential, originally appeared in [5] can generate both the step and the break model of scalar primordial spectrum in two different limits. The discussed potential is found to be favored by Planck data but marginally disfavored by BICEP2 results as it produces slightly lower amplitude of tensor primordial spectrum. Hence, if the tensor-to-scalar ratio ($r$) quoted by BICEP2 persists, it is of importance that we generate inflationary models with large $r$ and at the same time provide suppression in scalar primordial spectrum at large scales.

We discuss some details regarding the method of smoothed residuals, which has recently been used to search for anisotropic signals in low-redshift distance measurements (Supernovae). In this short note we focus on some details regarding the implementation of the method, particularly the issue of effectively detecting signals in data that are inhomogeneously distributed on the sky. Using simulated data, we argue that the original method proposed in Colin et al. [1] will not detect spurious signals due to incomplete sky coverage, and that introducing additional Gaussian weighting to the statistic as in [2] can hinder its ability to detect a signal. Issues related to the width of the Gaussian smoothing are also discussed.

We revisit the question of hemispherical power asymmetry in the WMAP and Planck temperature sky maps by measuring the local variance over the sky and on disks of various sizes. For the 2013 Planck sky map we find that none of the 1000 available isotropic Planck "Full Focal Plane" simulations have a larger variance asymmetry than that estimated from the data, suggesting the presence of an anisotropic signature formally significant at least at the $3.3\sigma$ level. For the WMAP 9-year data we find that 5 out of 1000 simulations have a larger asymmetry. The preferred direction for the asymmetry from the Planck data is $(l,b)=(212^\circ,-13^\circ)$, in good agreement with previous reports of the same hemispherical power asymmetry.

Within the context of the concordance model of cosmology we test the consistency of the angular power spectrum data from WMAP and Planck looking for possible systematics. The best fit concordance model to each observation is used as a mean function along with a Crossing function with an orthogonal basis to fit the data from the other observation searching for any possible deviation. We report that allowing an overall amplitude shift in the observed angular power spectra of the two observations, the best fit mean function from Planck data is consistent with WMAP 9 year data but the best fit mean function generated from WMAP-9 data is not consistent with Planck data at the $3\sigma$ level. This is an expected result when there is no clear systematic/tension between two observations and one of them has a considerably higher precision. We conclude that there is no clear tension between Planck and WMAP 9 year angular power spectrum data from a statistical point of view (allowing the overall amplitude shift). Our result highlights the fact that while the angular power spectrum from cosmic microwave background observations is a function of various cosmological parameters, comparing individual parameters might be misleading in the presence of cosmographic degeneracies. Another main result of our analysis is the importance of the overall amplitudes of the observed spectra from Planck and WMAP observations. Fixing the amplitudes at their reported values results in an unresolvable tension between the two observations at more than $3\sigma$ level which can be a hint towards a serious systematic.

We confront the concordance (standard) model of cosmology, the spatially flat $\Lambda$CDM Universe with power-law form of the primordial spectrum with Planck CMB angular power spectrum data searching for possible smooth deviations beyond the flexibility of the standard model. The departure from the concordance cosmology is modeled in the context of Crossing statistic and statistical significance of this deviation is used as a measure to test the consistency of the standard model to the Planck data. Derived Crossing functions suggest the presence of some broad features in angular spectrum beyond the expectations of the concordance model. Our results indicate that the concordance model of cosmology is consistent to the Planck data only at 2 to 3$\sigma$ confidence level if we allow smooth deviations from the angular power spectrum given by the concordance model. This might be due to random fluctuations or may hint towards smooth features in the primordial spectrum or departure from another aspect of the standard model. Best fit Crossing functions indicate that there are lack of power in the data at both low-$\ell$ and high-$\ell$ with respect to the concordance model. This hints that we may need some modifications in the foreground modeling to resolve the significant inconsistency at high-$\ell$. However, presence of some systematics at high-$\ell$ might be another reason for the deviation we found in our analysis.

One of the main goals of physical cosmology is to reconstruct the expansion history of the universe and finding the actual model of dark energy. In this article I review the difficulties of understanding dark energy and discuss about two strategic approaches, 'reconstructing dark energy' and 'falsifying dark energy models'. While one can use the data to reconstruct the expansion history of the universe and the properties of dark energy using novel approaches, considering the data limitations and its uncertainties we have to deal with cosmographic degeneracies that makes it difficult to distinguish between different dark energy models. On the other hand one can use the power of the data to falsify an assumed model using advanced statistical techniques. Within all these issues, focusing on falsification of cosmological constant has a particular importance since finding any significant deviation from $\Lambda$ would result to a break through in theoretical physics, ruling out the standard concordance model of cosmology.

With the recently published Cosmic Microwave Background data from Planck we address the optimized binning of the primordial power spectrum. As an important modification to the usual binning of the primordial spectrum, along with the spectral amplitude of the bins, we allow the position of the bins also to vary. This technique enables us to address the location of the possible broad physical features in the primordial spectrum with relatively smaller number of bins compared to the analysis performed earlier. This approach is in fact a reconstruction method looking for broad features in the primordial spectrum and avoiding fitting noise in the data. Performing Markov Chain Monte Carlo analysis we present samples of the allowed primordial spectra with broad features consistent with Planck data. To test how realistic it is to have step-like features in primordial spectrum we revisit an inflationary model, proposed by A. A. Starobinsky which can address the similar features obtained from the binning of the spectrum. Using the publicly available code BINGO, we numerically calculate the local $f_{\rm NL}$ for this model in equilateral and arbitrary triangular configurations of wavevectors and show that the obtained non-Gaussianity for this model is consistent with Planck results. In this paper we have also considered different spectral tilts at different bins to identify the cosmological scale that the spectral index needs to have a red tilt and it is interesting to report that spectral index cannot be well constrained up to $k \sim 0.01~ {\rm Mpc}^{-1}$.

In this paper, using and implementing a new line of sight CMB code, called CMBAns [1], that allows us to modify H(z) for any given feature at any redshift we study the effect of changes in the expansion history of the Universe on the CMB power spectrum. Motivated by the detailed analytical calculations of the effects of the changes in H(z) on ISW plateau and CMB low multipoles, we study two phenomenological parametric form of the expansion history using WMAP data and through MCMC analysis. Our MCMC analysis shows that the standard LCDM cosmological model is consistent with the CMB data allowing the expansion history of the Universe vary around this model at different redshifts. However, our analysis also shows that a decaying dark energy model proposed in [2] has in fact a marginally better fit than the standard cosmological constant model to CMB data. Concordance of our studies here with the previous analysis showing that Baryon Acoustic Oscillation (BAO) and supernovae data (SN Ia) also prefer mildly this decaying dark energy model to LCDM, makes this finding interesting and worth further investigation.

Constraints on the main cosmological parameters using CMB or large scale structure data are usually based on power-law assumption of the primordial power spectrum (PPS). However, in the absence of a preferred model for the early universe, this raises a concern that current cosmological parameter estimates are strongly prejudiced by the assumed power-law form of PPS. In this paper, for the first time, we perform cosmological parameter estimation allowing the free form of the primordial spectrum. This is in fact the most general approach to estimate cosmological parameters without assuming any particular form for the primordial spectrum. We use direct reconstruction of the PPS for any point in the cosmological parameter space using recently modified Richardson-Lucy algorithm however other alternative reconstruction methods could be used for this purpose as well. We use WMAP 9 year data in our analysis considering CMB lensing effect and we report, for the first time, that the flat spatial universe with no cosmological constant is ruled out by more than 4\sigma confidence limit without assuming any particular form of the primordial spectrum. This would be probably the most robust indication for dark energy using CMB data alone. Our results on the estimated cosmological parameters show that higher values of baryonic and matter density and lower value of Hubble parameter (in comparison to the estimated values by assuming power-law PPS) is preferred by the data. However, the estimated cosmological parameters by assuming free form of the PPS have overlap at 1\sigma confidence level with the estimated values assuming power-law form of PPS.

We have improved further the error sensitive Richardson-Lucy deconvolution algorithm making it applicable directly on the un-binned measured angular power spectrum of Cosmic Microwave Background observations to reconstruct the form of the primordial power spectrum. This improvement makes the application of the method significantly more straight forward by removing some intermediate stages of analysis allowing a reconstruction of the primordial spectrum with higher efficiency and precision and with lower computational expenses. Applying the modified algorithm we fit the WMAP 9 year data using the optimized reconstructed form of the primordial spectrum with more than 300 improvement in \chi^2 with respect to the best fit power-law. This is clearly beyond the reach of other alternative approaches and reflects the efficiency of the proposed method in the reconstruction process and allow us to look for any possible feature in the primordial spectrum projected in the CMB data. Though the proposed method allow us to look at various possibilities for the form of the primordial spectrum, all having good fit to the data, proper error-analysis is needed to test for consistency of theoretical models since, along with possible physical artefacts, most of the features in the reconstructed spectrum might be arising from fitting noises in the CMB data. Reconstructed error-band for the form of the primordial spectrum using many realizations of the data, all bootstrapped and based on WMAP 9 year data, shows proper consistency of power-law form of the primordial spectrum with the WMAP 9 data at all wave numbers. Including WMAP polarization data in to the analysis have not improved much our results due to its low quality but we expect Planck data will allow us to make a full analysis on CMB observations on both temperature and polarization separately and in combination.

We compare two different probes of the expansion history of the universe, namely, luminosity distances from type Ia supernovae and angular diameter distances from galaxy clusters, using the Bayesian interpretation of Crossing statistic [1, 2] in conjunction with the assumption of cosmic duality relation. Our analysis is conducted independently of any a-priori assumptions about the nature of dark energy. The model independent method which we invoke searches for inconsistencies between SNIa and galaxy cluster data sets. If detected such an inconsistency would imply the presence of systematics in either of the two data sets. Simulating observations based on expected WFIRST supernovae data and X-ray eROSITA + SZ Planck cluster data, we show that our method allows one to detect systematics with high precision and without advancing any hypothesis about the nature of dark energy.

We use Gaussian Processes to map the expansion history of the universe in a model independent manner from the Union2.1 supernovae data and then apply these reconstructed results to solve for the growth history. By comparing this to BOSS and WiggleZ large scale structure data we examine whether growth is determined wholly by expansion, with the measured gravitational growth index testing gravity without assuming a model for dark energy. A further model independent analysis looks for redshift dependent deviations of growth from the general relativity solution without assuming the growth index form. Both approaches give results consistent with general relativity.

Baryon Acoustic Oscillations (BAO) provide an important standard ruler which can be used to probe the recent expansion history of our universe. We show how a simple extension of the Om diagnostic, which we call Om3, can combine standard ruler information from BAO with standard candle information from type Ia supernovae (SNIa) to yield a powerful novel null diagnostic of the cosmological constant hypothesis. A unique feature of Om3 is that it requires minimal cosmological assumptions since its determination does not rely upon prior knowledge of either the current value of the matter density and the Hubble constant, or the distance to the last scattering surface. Observational uncertainties in these quantities therefore do not affect the reconstruction of Om3. We reconstruct Om3 using the Union 2.1 SNIa data set and BAO data from SDSS, WiggleZ and 6dFGS. Our results are consistent with dark energy being the cosmological constant. We show how Om and Om3 can be used to obtain accurate model independent constraints on the properties of dark energy from future data sets such as BigBOSS.

We present that by combining Crossing Statistic and Smoothing method one can reconstruct the expansion history of the universe with a very high precision without considering any prior on the cosmological quantities such as the equation of state of dark energy. We show that the presented method performs very well in reconstruction of the expansion history of the universe independent of the underlying models and it works well even for non-trivial dark energy models with fast or slow changes in the equation of state of dark energy. Accuracy of the reconstructed quantities along with independence of the method to any prior or assumption gives the proposed method advantages to the other non-parametric methods proposed before in the literature. Applying on the Union 2.1 supernovae combined with WiggleZ BAO data we present the reconstructed results and test the consistency of the two data sets in a model independent manner. Results show that latest available supernovae and BAO data are in good agreement with each other and spatially flat LCDM model is in concordance with the current data.

Gaussian processes provide a method for extracting cosmological information from observations without assuming a cosmological model. We carry out cosmography -- mapping the time evolution of the cosmic expansion -- in a model-independent manner using kinematic variables and a geometric probe of cosmology. Using the state of the art supernova distance data from the Union2.1 compilation, we constrain, without any assumptions about dark energy parametrization or matter density, the Hubble parameter and deceleration parameter as a function of redshift. Extraction of these relations is tested successfully against models with features on various coherence scales, subject to certain statistical cautions.

By introducing Crossing functions and hyper-parameters I show that the Bayesian interpretation of the Crossing Statistics [1] can be used trivially for the purpose of model selection among cosmological models. In this approach to falsify a cosmological model there is no need to compare it with other models or assume any particular form of parametrization for the cosmological quantities like luminosity distance, Hubble parameter or equation of state of dark energy. Instead, hyper-parameters of Crossing functions perform as discriminators between correct and wrong models. Using this approach one can falsify any assumed cosmological model without putting priors on the underlying actual model of the universe and its parameters, hence the issue of dark energy parametrization is resolved. It will be also shown that the sensitivity of the method to the intrinsic dispersion of the data is small that is another important characteristic of the method in testing cosmological models dealing with data with high uncertainties.

The observables of the perturbed universe, CMB anisotropy and large structures, depend on a set of cosmological parameters, as well as, the assumed nature of primordial perturbations. In particular, the shape of the primordial power spectrum (PPS) is, at best, a well motivated assumption. It is known that the assumed functional form of the PPS in cosmological parameter estimation can affect the best fit parameters and their relative confidence limits. In this paper, we demonstrate that a specific assumed form actually drives the best fit parameters into distinct basins of likelihood in the space of cosmological parameters where the likelihood resists improvement via modifications to the PPS. The regions where considerably better likelihoods are obtained allowing free form PPS lie outside these basins. In the absence of a preferred model of inflation, this raises a concern that current cosmological parameters estimates are strongly prejudiced by the assumed form of PPS. Our results strongly motivate approaches toward simultaneous estimation of the cosmological parameters and the shape of the primordial spectrum from upcoming cosmological data. It is equally important for theorists to keep an open mind towards early universe scenarios that produce features in the PPS.

Recent detection of high-redshift, massive clusters through Sunyaev-Zel'dovich observations has opened up a new way to test cosmological models. It is known that detection of a single supermassive cluster at a very high redshift can rule out many cosmological models all together. However, since dealing with different observational biases makes it difficult to test the likeliness of the data assuming a cosmological model, most of the cluster data (except those with high mass-redshift) stays untouched in confronting cosmological models with cluster observations. We propose here that one can use the relative abundance of the clusters with different masses at different redshifts to test the likeliness of the data in the context of cosmological models. For this purpose we propose a simple parametric form for the efficiency of observing clusters at different mass-redshift and we test if the standard LCDM model can explain the observed abundance of the clusters using this efficiency parameterization. We argue that one cannot expect an unusual and highly parametric form of the efficiency function to fit the observed data assuming a theoretical model. Using many realizations of Monte Carlo simulations we show that the standard spatially flat LCDM model is barely consistent with the SPT cluster data using a simple and plausible two-dimensional efficiency function for detection of the clusters. More cluster data are needed to make any strong conclusion.

We examine the dark energy and matter densities allowed by precision measurements of distances out to various redshifts, in the presence of spatial curvature and (near) arbitrary behavior of the dark energy equation of state. Degeneracies among the parameters permit a remarkably large variation in their values when using only distance measurements of the late time universe and making no assumptions about the dark energy or curvature. Going beyond distance measurements to a lower limit on the growth of structure bounds the allowed region significantly but still leaves considerable freedom to match a flat LCDM model with dark energy very different from a cosmological constant. The combination of distances with Hubble parameter, gravitational lensing or other large scale structure data is essential to determining robustly the cosmological model.

We propose a new statistic that has been designed to be used in situations where the intrinsic dispersion of a data set is not well known: The Crossing Statistic. This statistic is in general less sensitive than `chi^2' to the intrinsic dispersion of the data, and hence allows us to make progress in distinguishing between different models using goodness of fit to the data even when the errors involved are poorly understood. The proposed statistic makes use of the shape and trends of a model's predictions in a quantifiable manner. It is applicable to a variety of circumstances, although we consider it to be especially well suited to the task of distinguishing between different cosmological models using type Ia supernovae. We show that this statistic can easily distinguish between different models in cases where the `chi^2' statistic fails. We also show that the last mode of the Crossing Statistic is identical to `chi^2', so that it can be considered as a generalization of `chi^2'.

The question of the transition to global isotropy from our anisotropic local Universe is studied using the Union 2 catalogue of Type Ia supernovae (SNe Ia). We construct a "residual" statistic sensitive to systematic shifts in their brightness in different directions and use this to search in different redshift bins for a preferred direction on the sky in which the SNe Ia are brighter or fainter relative to the 'standard' LCDM cosmology. At low redshift (z<0.05) we find that an isotropic model such as LCDM is barely consistent with the SNe Ia data at 2-3 sigma. A complementary maximum likelihood analysis of peculiar velocities confirms this finding -- there is a bulk flow of around 260 km/sec at z \sim 0.06, which disagrees with LCDM at 1-2 sigma. Since the Shapley concentration is believed to be largely responsible for this bulk flow, we make a detailed study of the infall region: the SNe Ia falling away from the Local Group towards Shapley are indeed significantly dimmer than those falling towards us and on to Shapley. Convergence to the CMB rest frame must occur well beyond Shapley (z>0.06) so the low redshift bulk flow can systematically bias any reconstruction of the expansion history of the Universe. At high redshifts z>0.15 the agreement between the SNe Ia data and the isotropic LCDM model does improve, however, the sparseness and low quality of the data means that LCDM cannot be singled out as the preferred cosmological model.

We use the Om statistic and the Genetic Algorithms (GA) in order to derive a null test on the spatially flat cosmological constant model $\Lambda$CDM. This is done in two steps: first, we apply the GA to the Constitution SNIa data in order to acquire a model independent reconstruction of the expansion history of the Universe $H(z)$ and second, we use the reconstructed $H(z)$ in conjunction with the Om statistic, which is constant only for the $\Lambda$CDM model, to derive our constraints. We find that while $\Lambda$CDM is consistent with the data at the $2\sigma$ level, some deviations from $\Lambda$CDM model at low redshifts can be accommodated.

In this paper we consider a unique model of inflation where the universe undergoes rapid asymmetric oscillations, each cycle lasting millions of Planck time. Over many-many cycles the space-time expands to mimic the standard inflationary scenario. Moreover, these rapid oscillations leave a distinctive periodic signature in ln(k) in the primordial power spectrum, where k denotes the comoving scale. The best fit parameters of the cyclic-inflation model provides a very good fit to the 7-year WMAP data.

The dark energy problem has led to speculation that not only may LCDM be wrong, but that the FLRW models themselves may not even provide the correct family of background models. We discuss how direct measurements of H(z) can be used to formulate tests of the standard paradigm in cosmology. On their own, such measurements can be used to test for deviations from flat LCDM. When combined with supernovae distances, Hubble rate measurements provide a test of the Copernican principle and the homogeneity assumption of the standard model, which is independent of dark energy or metric based theory of gravity. A modification of this test also provides a model independent observable for flatness which decorrelates curvature determination from dark energy. We investigate these tests using Hubble rate measurements from age data, as well as from a Hubble rate inferred from recent measurements of the baryon acoustic oscillations. While the current data is too weak to say anything significant, these tests are exciting prospects for the future.

Features in the primordial power spectrum have been suggested as an explanation for glitches in the angular power spectrum of temperature anisotropies measured by the WMAP satellite. However, these glitches might just as well be artifacts of noise or cosmic variance. Using the effective Delta chi^2 between the best-fit power-law spectrum and a deconvolved primordial spectrum as a measure of "featureness" of the data, we perform a full Monte-Carlo analysis to address the question of how significant the recovered features are. We find that in 26% of the simulated data sets the reconstructed spectrum yields a greater improvement in the likelihood than for the actually observed data. While features cannot be categorically ruled out by this analysis, and the possibility remains that simple theoretical models which predict some of the observed features might stand up to rigorous statistical testing, our results suggest that WMAP data are consistent with the assumption of a featureless power-law primordial spectrum.

We investigate the course of cosmic expansion in its `recent past' using the Constitution SN Ia sample (which includes CfA data at low redshifts), jointly with signatures of baryon acoustic oscillations (BAO) in the galaxy distribution and fluctuations in the cosmic microwave background (CMB). Earlier SN Ia data sets could not address this issue because of a paucity of data at low redshifts. Allowing the equation of state of dark energy (DE) to vary, we find that a coasting model of the universe (q_0=0) fits the data about as well as LCDM. This effect, which is most clearly seen using the recently introduced `Om' diagnostic, corresponds to an increase of Om(z) and q(z) at redshifts z \lleq 0.3. In geometrical terms, this suggests that cosmic acceleration may have already peaked and that we are currently witnessing its slowing down. The case for evolving DE strengthens if a subsample of the Constitution set consisting of SNLS+ESSENCE+CfA SN Ia data is analysed in combination with BAO+CMB using the same statistical methods. The effect we observe could correspond to DE decaying into dark matter (or something else). A toy model which mimics this process agrees well with the combined SN Ia+BAO+CMB data.

It has recently been pointed out by Kowalski et. al. (arxiv:0804.4142) that there is `an unexpected brightness of the SnIa data at z>1'. We quantify this statement by constructing a new statistic which is applicable directly on the Type Ia Supernova (SnIa) distance moduli. This statistic is designed to pick up systematic brightness trends of SnIa datapoints with respect to a best fit cosmological model at high redshifts. It is based on binning the normalized differences between the SnIa distance moduli and the corresponding best fit values in the context of a specific cosmological model (eg LCDM). We then focus on the highest redshift bin and extend its size towards lower redshifts until the Binned Normalized Difference (BND) changes sign (crosses 0) at a redshift z_c (bin size N_c). The bin size N_c of this crossing (the statistical variable) is then compared with the corresponding crossing bin size N_{mc} for Monte Carlo data realizations based on the best fit model. We find that the crossing bin size N_c obtained from the Union08 and Gold06 data with respect to the best fit LCDM model is anomalously large compared to N_{mc} of the corresponding Monte Carlo datasets obtained from the best fit LCDM in each case. In particular, only 2.2% of the Monte Carlo LCDM datasets are consistent with the Gold06 value of N_c while the corresponding probability for the Union08 value of N_c is 5.3%. Thus, according to this statistic, the probability that the high redshift brightness bias of the Union08 and Gold06 datasets is realized in the context of a (w_0,w_1)=(-1,0) model (LCDM cosmology) is less than 6%. The corresponding realization probability in the context of a (w_0,w_1)=(-1.4,2) model is more than 30% for both the Union08 and the Gold06 datasets.

A sudden small change in the second derivative of the inflaton potential can result in a universal local feature in the spectrum of primordial perturbations generated during inflation. The exact solution describing this feature \cite{minu} is characterized by a step in the spectral index modulated by characteristic oscillations and results in a large running of the spectral index localized over a few e-folds of scale. In this paper we confront this step-like feature with the 5 year WMAP results and demonstrate that it provides a better fit to this data than a featureless initial spectrum. If such a feature exists at all, then it should lie at sufficiently large scales $k_0 \lesssim 0.003 {\rm Mpc}^{-1}$ corresponding to $l \lesssim 40$. The sign of the effect is shown to correspond to the negative running of $n_s$ localized near this scale. This feature could arise as a result of a `mini-waterfall'-type fast second order phase transition experienced by an auxiliary heavy field during inflation, in a model similar to hybrid inflation (though for a different choice of parameters). If this is the case, then the auxiliary field should be positively coupled to the inflaton.

We investigate the cosmological properties of an "induced gravity" brane scenario in the absence of mirror symmetry with respect to the brane. We find that brane evolution can proceed along one of four distinct branches. By contrast, when mirror symmetry is imposed, only two branches exist, one of which represents the self-accelerating brane, while the other is the so-called normal branch. This model incorporates many of the well-known possibilities of brane cosmology including phantom acceleration (w < -1), self-acceleration, transient acceleration, quiescent singularities, and cosmic mimicry. Significantly, the absence of mirror symmetry also provides an interesting way of inducing a sufficiently small cosmological constant on the brane. A small (positive) Lambda-term in this case is induced by a small asymmetry in the values of bulk fundamental constants on the two sides of the brane.

We introduce two new diagnostics of dark energy (DE). The first, Om, is a combination of the Hubble parameter and the cosmological redshift and provides a "null test" of dark energy being a cosmological constant. Namely, if the value of Om(z) is the same at different redshifts, then DE is exactly cosmological constant. The slope of Om(z) can differentiate between different models of dark energy even if the value of the matter density is not accurately known. For DE with an unevolving equation of state, a positive slope of Om(z) is suggestive of Phantom (w < -1) while a negative slope indicates Quintessence (w > -1). The second diagnostic, "acceleration probe"(q-probe), is the mean value of the deceleration parameter over a small redshift range. It can be used to determine the cosmological redshift at which the universe began to accelerate, again without reference to the current value of the matter density. We apply the "Om" and "q-probe" diagnostics to the Union data set of type Ia supernovae combined with recent data from the cosmic microwave background (WMAP5) and baryon acoustic oscillations.

In this paper we implement an improved (error sensitive) Richardson-Lucy deconvolution algorithm on the measured angular power spectrum from the WMAP 3 year data to determine the primordial power spectrum assuming different points in the cosmological parameter space for a flat LCDM cosmological model. We also present the preliminary results of the cosmological parameter estimation by assuming a free form of the primordial spectrum, for a reasonably large volume of the parameter space. The recovered spectrum for a considerably large number of the points in the cosmological parameter space has a likelihood far better than a `best fit' power law spectrum up to \Delta \chi^2_{eff} \approx -30. We use Discrete Wavelet Transform (DWT) for smoothing the raw recovered spectrum from the binned data. The results obtained here reconfirm and sharpen the conclusion drawn from our previous analysis of the WMAP 1st year data. A sharp cut off around the horizon scale and a bump after the horizon scale seem to be a common feature for all of these reconstructed primordial spectra. We have shown that although the WMAP 3 year data prefers a lower value of matter density for a power law form of the primordial spectrum, for a free form of the spectrum, we can get a very good likelihood to the data for higher values of matter density. We have also shown that even a flat CDM model, allowing a free form of the primordial spectrum, can give a very high likelihood fit to the data. Theoretical interpretation of the results is open to the cosmology community. However, this work provides strong evidence that the data retains discriminatory power in the cosmological parameter space even when there is full freedom in choosing the primordial spectrum.

We have improved upon the method of smoothing supernovae data to reconstruct the expansion history of the universe, h(z), using two latest datasets, Gold and SNLS. The reconstruction process does not employ any parameterization and is independent of any dark energy model. The reconstructed h(z) is used to derive the distance factor "A" up to redshift 0.35 and the results are compared with the given value of "A" from detection of baryon acoustic oscillation peak (BAO). We find very good agreement between supernovae observations and the results from BAO for \Omega_{0m} \approx 0.276 \pm 0.023. The estimated values of \Omega_{0m} are completely model-independent and are only based on observational data. The derived values of \Omega_{0m} are then used to reconstruct the equation of state of dark energy, w(z). Using our smoothing method we can demonstrate that while SNLS data are in very good agreement with LCDM, the Gold sample slightly prefers evolving dark energy. We also show that proper estimation of the equation of state of dark energy at the high redshifts would be impossible at the current status of observations.

Precise measurements of the anisotropies in the cosmic microwave background enable us to do an accurate study on the form of the primordial power spectrum for a given set of cosmological parameters. In a previous paper (Shafieloo and Souradeep 2004), we implemented an improved (error sensitive) Richardson-Lucy deconvolution algorithm on the measured angular power spectrum from the first year of WMAP data to determine the primordial power spectrum assuming a concordance cosmological model. This recovered spectrum has a likelihood far better than a scale invariant, or, `best fit' scale free spectra (\Delta ln L = 25 w.r.t. Harrison Zeldovich, and, \Delta ln L = 11 w.r.t. power law with n_s=0.95). In this paper we use Discrete Wavelet Transform (DWT) to decompose the local features of the recovered spectrum individually to study their effect and significance on the recovered angular power spectrum and hence the likelihood. We show that besides the infra-red cut off at the horizon scale, the associated features of the primordial power spectrum around the horizon have a significant effect on improving the likelihood. The strong features are localised at the horizon scale.

We propose a non-parametric method of smoothing supernova data over redshift using a Gaussian kernel in order to reconstruct important cosmological quantities including H(z) and w(z) in a model independent manner. This method is shown to be successful in discriminating between different models of dark energy when the quality of data is commensurate with that expected from the future SuperNova Acceleration Probe (SNAP). We find that the Hubble parameter is especially well-determined and useful for this purpose. The look back time of the universe may also be determined to a very high degree of accuracy (\lleq 0.2 %) in this method. By refining the method, it is also possible to obtain reasonable bounds on the equation of state of dark energy. We explore a new diagnostic of dark energy-- the `w-probe'-- which can be calculated from the first derivative of the data. We find that this diagnostic is reconstructed extremely accurately for different reconstruction methods even if \Omega_m is marginalized over. The w-probe can be used to successfully distinguish between $\Lambda$CDM and other models of dark energy to a high degree of accuracy.

The observed angular power spectrum of the cosmic microwave background temperature anisotropy, $C_l$, is a convolution of a cosmological radiative transport kernel with an assumed primordial power spectrum of inhomogeneities. Exquisite measurements of $C_l$ over a wide range of multipoles from the Wilkinson Microwave Anisotropy Probe (WMAP) has opened up the possibility to deconvolve the primordial power spectrum for a given set of cosmological parameters (base model). We implement an improved (error sensitive) Richardson-Lucy deconvolution algorithm on the measured angular power spectrum from WMAP assuming a concordance cosmological model. The most prominent feature of the recovered $P(k)$ is a sharp, infra-red cut off on the horizon scale. The resultant $C_l$ spectrum using the recovered spectrum has a likelihood far better than a scale invariant, or, `best fit' scale free spectra ($\Delta\ln{\cal L}=25$ {\it w.r.t.} Harrison Zeldovich, and, $\Delta\ln{\cal L}=11$ {\it w.r.t.} power law with $n_s=0.95$). The recovered $P(k)$ has a localized excess just below the cut-off which leads to great improvement of likelihood over the simple monotonic forms of model infra-red cut-off spectra considered in the post WMAP literature. The recovered $P(k)$, in particular, the form of infra-red cut-off is robust to small changes in the cosmological parameters. We show that remarkably similar form of infra-red cutoff is known to arise in very reasonable extensions and refinements of the predictions from simple inflationary scenarios. Our method can be extended to other cosmological observations such as the measured matter power spectrum and, in particular, the much awaited polarization spectrum from WMAP.