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 provide an update on the constraints on extended reionization histories with the Planck 2018 cosmic microwave background anisotropy data. The Planck 2018 data on large angular scales improve the measurement of the $E$-mode polarization reionization bump at low multipoles providing the possibility to improve our previous results. Using a minor modification to the original Poly-reion model for the reionization history, we find that the Planck 2018 data significantly improve all our previous results: we find as optical depth of $\tau=0.0572_{-0.0075}^{+0.0064}$ at 68% CL, that early onsets of reionization are strongly disfavoured, i.e. redshift when the reionization begins, $z_{xe=0}=18.18_{-10.89}^{+1.61}$ at 68% CL,and that reionization duration (defined between 10% and 99% reionization) is significantly reduced, i.e. $\Delta_z^{Reion}=4.59_{-2.45}^{+1.67}$ at 68% CL. We explore possible correlations between reionization histories and cosmological parameters, including important extensions beyond $\Lambda$CDM. We find that the degeneracy between reionization and scalar spectral index,neutrino mass sum, spatial curvature, dark matter annihilation and other non-standard models are significantly reduced.The reduction of the error bars and the degeneracies, together with the shift towards lower values of the optical depth that we observe in the Poly-reion model are mainly driven by the new low-$\ell$ polarization likelihood of Planck 2018 baseline based on the HFI data. This is confirmed also by the results derived without this likelihood and the ones with different alternatives to the baseline that are presented for a subset of models.

Primordial black holes (PBHs) generated by gravitational collapse of large primordial overdensities can be a fraction of the observed dark matter. In this paper, we introduce a mechanism to produce a large peak in the primordial power spectrum (PPS) in two-field inflationary models characterized by two stages of inflation based on a large non-canonical kinetic coupling. This mechanism is generic to several two-field inflationary models, due to a temporary tachyonic instability of the isocurvature perturbations at the transition between the two stages of inflation. We numerically compute the primordial perturbations from largest scales to the small scales corresponding to that of PBHs using an extension of BINGO (BI-spectra and Non-Gaussianity Operator). Moreover we numerically compute the stochastic background of gravitational waves (SBGW) produced by second order scalar perturbations within frequencies ranging from nano-Hz to KHz that covers the observational scales corresponding to Pulsar Timing Arrays, Square Kilometer Array to that of Einstein telescope. We discuss the prospect of its detection by these proposed and upcoming gravitational waves experiments.

We investigate the generation of features at large scales in the primordial power spectrum (PPS) when inflation is driven by two scalar fields. In canonical single field models of inflation, these features are often generated due to deviations from the slow-roll regime. While deviations from slow-roll can be naturally achieved in two field models due to a sharp turn in the trajectory in the field space, features at the largest scales of the types suggested by CMB temperature anisotropies are more difficult to achieve in models involving two canonical scalar fields due to the presence of isocurvature fluctuations. We show instead that a coupling between the kinetic terms of the scalar fields can easily produce such features. We discuss models whose theoretical predictions are consistent with current observations and highlight the implications of our results.

Cosmic Microwave Background temperature and polarization anisotropies from Planck have estimated a lower value of the optical depth to reionization ($\tau$) compared to WMAP. A significant period in the reionization history would then fall within $6<z< 10$, where detection of galaxies with Hubble Frontier Fields program and independent estimation of neutral hydrogen in the inter galactic medium by Lyman-$\alpha$ observations are also available. This overlap allows an analysis of cosmic reionization which utilizes a direct combination of CMB and these astrophysical measurements and potentially breaks degeneracies in parameters describing the physics of reionization. For the first time we reconstruct reionization histories by assuming photo-ionization and recombination rates to be free-form and by allowing underlying cosmological parameters to vary with CMB (temperature and polarization anisotropies and lensing) data from Planck 2018 release and a compilation of astrophysical data. We find an excellent agreement between the low-$\ell$ Planck 2018 HFI polarization likelihood and astrophysical data in determining the integrated optical depth. By combining both data, we report for a minimal reconstruction $\tau=0.051^{+0.001+0.002}_{-0.0012-0.002}$ at 68\% and 95\% CL, which, for the errors in the current astrophysical measurements quoted in the literature, is nearly twice better than the projected cosmic variance limited CMB measurements. For the duration of reionization, redshift interval between 10\% and complete ionization, we get $2.9^{+0.12+0.29}_{-0.16-0.26}$ at 68\% and 95\% CL, which improves significantly on the corresponding result obtained by using Planck 2015 data. By a Bayesian analysis of the combined results we do not find evidence beyond monotonic reionization histories, therefore multi-phase reionization scenario is disfavored compared to minimal alternatives.

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.

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.

We explore the constraints on the history of reionization from Planck 2015 Cosmic Microwave Background (CMB) data and we derive the forecasts for future CMB observations. We consider a class of monotonic histories of reionization as parametrized by two additional extra parameters with respect to the average optical depth used in the instantaneous reionization modeling. We investigate the degeneracies between the history of reionization and selected extensions of the standard cosmological model. In particular, we consider the degeneracies with the total mass of the neutrino sector and we discuss the possible correlation between the dark matter annihilation and the duration of reionization in the CMB. We use an extension to poly-reion model that was proposed in Hazra and Smoot, JCAP 1711, 028 (2017). We compare the constraints from Planck 2015 data with the predicted constraints from possible future CMB mission as LiteBIRD, and we also use the proposed CORE-like specifications as an example of what higher resolution can bring in addition. We find that the degeneracy between the average optical depth and the duration of reionization will be substantially removed by both concepts. Degeneracies between the reionization history and either the total neutrino mass and properties of dark matter annihilation will also be improved by future surveys. We find only marginal improvement in the constraints on reionization history for the higher resolution in the case of long duration of reionization.

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 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 constrain the history of reionization using the data from Planck 2015 Cosmic Microwave Background (CMB) temperature and polarization anisotropy observations. We also use prior constraints on the reionization history at redshifts $\sim7-8$ obtained from Lyman-$\alpha$ emission observations. Using the free electron fractions at different redshifts as free parameters, we construct the complete reionization history using polynomials. Our construction provides an extremely flexible framework to search for the history of reionization as a function of redshifts. We present a conservative and an optimistic constraint on reionization that are categorized by the flexibilities of the models and datasets used to constrain them, and we report that CMB data marginally favors extended reionization histories. In both the cases, we find the mean values of optical depth to be larger ($\approx0.09$ and $0.1$) than what we find in standard steplike reionization histories ($0.079\pm0.017$). At the same time we also find that the maximum free electron fraction allowed by the data for redshifts more than 15 is $\sim0.25$ at 95.4\% confidence limit in the case of optimistic constraint.

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 forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

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.

As is well known, the non-Gaussianity parameter $f_{_{\rm NL}}$, which is often used to characterize the amplitude of the scalar bi-spectrum, can be expressed completely in terms of the scalar spectral index $n_{\rm s}$ in the squeezed limit, a relation that is referred to as the consistency condition. This relation, while it is largely discussed in the context of slow roll inflation, is actually expected to hold in any single field model of inflation, irrespective of the dynamics of the underlying model, provided inflation occurs on the attractor at late times. In this work, we explicitly examine the validity of the consistency relation, analytically as well as numerically, away from slow roll. Analytically, we first arrive at the relation in the simple case of power law inflation. We also consider the non-trivial example of the Starobinsky model involving a linear potential with a sudden change in its slope (which leads to a brief period of fast roll), and establish the condition completely analytically. We then numerically examine the validity of the consistency relation in three inflationary models that lead to the following features in the scalar power spectrum due to departures from slow roll: (i) a sharp cut off at large scales, (ii) a burst of oscillations over an intermediate range of scales, and (iii) small, but repeated, modulations extending over a wide range of scales. It is important to note that it is exactly such spectra that have been found to lead to an improved fit to the CMB data, when compared to the more standard power law primordial spectra, by the Planck team. We evaluate the scalar bi-spectrum for an arbitrary triangular configuration of the wavenumbers in these inflationary models and explicitly illustrate that, in the squeezed limit, the consistency condition is indeed satisfied even in situations consisting of strong deviations from slow roll.

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.

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.

Sharp shapes in the inflaton potentials often lead to short departures from slow roll which, in turn, result in deviations from scale invariance in the scalar power spectrum. Typically, in such situations, the scalar power spectrum exhibits a burst of features associated with modes that leave the Hubble radius either immediately before or during the epoch of fast roll. Moreover, one also finds that the power spectrum turns scale invariant at smaller scales corresponding to modes that leave the Hubble radius at later stages, when slow roll has been restored. In other words, the imprints of brief departures from slow roll, arising out of sharp shapes in the inflaton potential, are usually of a finite width in the scalar power spectrum. Intuitively, one may imagine that the scalar bi-spectrum too may exhibit a similar behavior, i.e. a restoration of scale invariance at small scales, when slow roll has been reestablished. However, in the case of the Starobinsky model (viz. the model described by a linear inflaton potential with a sudden change in its slope) involving the canonical scalar field, it has been found that, a rather sharp, though short, departure from slow roll can leave a lasting and significant imprint on the bi-spectrum. The bi-spectrum in this case is found to grow linearly with the wavenumber at small scales, a behavior which is clearly unphysical. In this work, we study the effects of smoothening the discontinuity in the Starobinsky model on the scalar bi-spectrum. Focusing on the equilateral limit, we analytically show that, for smoother potentials, the bi-spectrum indeed turns scale invariant at suitably large wavenumbers. We also confirm the analytical results numerically using our newly developed code BINGO. We conclude with a few comments on certain related points.

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 constrain plausible dark energy models, parametrized by multiple candidate equation of state, using the recently published Cosmic Microwave Background (CMB) temperature anisotropy data from Planck together with the WMAP-9 low-$\ell$ polarization data and data from low redshift surveys. To circumvent the limitations of any particular equation of state towards describing all existing dark energy models, we work with three different equation of state covering a broader class of dark energy models and, hence, provide more robust and generic constraints on the dark energy properties. We show that a clear tension exists between dark energy constraints from CMB and non-CMB observations when one allows for dark energy models having both phantom and non-phantom behavior; while CMB is more favorable to phantom models, the low-z data prefers model with behavior close to a Cosmological Constant. Further, we reconstruct the equation of state of dark energy as a function of redshift using the results from combined CMB and non-CMB data and find that Cosmological Constant lies outside the 1$\sigma$ band for multiple dark energy models allowing phantom behavior. A considerable fine tuning is needed to keep models with strict non-phantom history inside 2$\sigma$ allowed range. This result might motivate one to construct phantom models of dark energy,which is achievable in the presence of higher derivative operators as in string theory. However, disallowing phantom behavior, based only on strong theoretical prior, leads to both CMB and non-CMB datasets agree on the nature of dark energy, with the mean equation of state being very close to the Cosmological Constant. Finally, to illustrate the impact of additional dark energy parameters on other cosmological parameters, we provide the cosmological parameter constraints for different dark energy models.

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.

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

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.

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 present a new and accurate Fortran code, the BI-spectra and Non-Gaussianity Operator (BINGO), for the efficient numerical computation of the scalar bi-spectrum and the non-Gaussianity parameter f_{NL} in single field inflationary models involving the canonical scalar field. The code can calculate all the different contributions to the bi-spectrum and the parameter f_{NL} for an arbitrary triangular configuration of the wavevectors. Focusing firstly on the equilateral limit, we illustrate the accuracy of BINGO by comparing the results from the code with the spectral dependence of the bi-spectrum expected in power law inflation. Then, considering an arbitrary triangular configuration, we contrast the numerical results with the analytical expression available in the slow roll limit, for, say, the case of the conventional quadratic potential. Considering a non-trivial scenario involving deviations from slow roll, we compare the results from the code with the analytical results that have recently been obtained in the case of the Starobinsky model in the equilateral limit. As an immediate application, we utilize BINGO to examine of the power of the non-Gaussianity parameter f_{NL} to discriminate between various inflationary models that admit departures from slow roll and lead to similar features in the scalar power spectrum. We close with a summary and discussion on the implications of the results we obtain.

We explore possibility of using the three dimensional bispectra of the Ly-alpha forest and the redshifted 21-cm signal from the post-reionization epoch to constrain primordial non-Gaussianity. Both these fields map out the large scale distribution of neutral hydrogen and maybe treated as tracers of the underlying dark matter field. We first present the general formalism for the auto and cross bispectrum of two arbitrary three dimensional biased tracers and then apply it to the specific case. We have modeled the 3D Ly-alpha transmitted flux field as a continuous tracer sampled along 1D skewers which corresponds to quasars sight lines. For the post reionization 21-cm signal we have used a linear bias model. We use a Fisher matrix analysis to present the first prediction for bounds on f_{NL} and the other bias parameters using the three dimensional 21-cm bispectrum and other cross bispectra. The bounds on f_{NL} depend on the survey volume, and the various observational noises. We have considered a BOSS like Ly-alpha survey where the average number density of quasars \bar{n} = 10^{-3} Mpc^{-2} and the spectra are measured at a 2-sigma level. For the 21-cm signal we have considered a 4000 hrs observation with a futuristic SKA like radio array. We find that bounds on f_{NL} obtained in our analysis (6 < \Delta f_{NL} < 65) is competitive with CMBR and galaxy surveys. We also find that by combining the various cross-bispectrum estimators it is possible to constrain f_{NL} at a level \Delta f_{NL} ~ 4.7. For the equilateral and orthogonal template we obtain \Delta f_{NL}^{equ}~17 and \Delta f_{NL}^{orth}~13 respectively for the combined estimator. This shall be important in the quest towards understanding the mechanism behind the generation of primordial perturbations.

Features in the primordial scalar power spectrum provide a possible roadway to describe the outliers at the low multipoles in the WMAP data. Apart from the CMB angular power spectrum, these features can also alter the matter power spectrum and, thereby, the formation of the large scale structure. Carrying out a complete numerical analysis, we investigate the effects of primordial features on the formation rates of the halos. We consider a few different inflationary models that lead to features in the scalar power spectrum and an improved fit to the CMB data, and analyze the corresponding imprints on the formation of halos. Performing a Markov Chain Monte Carlo analysis with the WMAP seven year data and the SDSS halo power spectrum from LRG DR7 for the models of our interest, we arrive at the parameter space of the models allowed by the data. We illustrate that, inflationary potentials, such as the quadratic potential with sinusoidal modulations and the axion monodromy model, which generate certain repeated, oscillatory features in the inflationary perturbation spectrum, do not induce a substantial difference in the number density of halos at their best fit values, when compared with, say, a nearly scale invariant spectrum as is generated by the standard quadratic potential. However, we find that the number density and the formation rates of halos change by about 13-22% for halo masses ranging over 10^4-10^14 solar mass, for potential parameters that lie within 2-sigma around the best fit values arrived at from the aforesaid joint constraints. We briefly discuss the implications of our results.

Amongst the multitude of inflationary models currently available, models that lead to features in the primordial scalar spectrum are drawing increasing attention, since certain features have been found to provide a better fit to the CMB data than the conventional, nearly scale invariant, primordial spectrum. In this work, we carry out a complete numerical analysis of two models that lead to oscillations over all scales in the scalar power spectrum. We consider the model described by a quadratic potential which is superposed by a sinusoidal modulation and the recently popular axion monodromy model. Since the oscillations continue even on to arc minute scales, in addition to the WMAP data, we also compare the models with the small scale data from ACT. Though, both the models, broadly, result in oscillations in the spectrum, interestingly, we find that, while the monodromy model leads to a considerably better fit to the data in comparison to the standard power law spectrum, the quadratic potential superposed with a sinusoidal modulation does not improve the fit to a similar extent. We also carry out forecasting of the parameters using simulated Planck data for both the models. We show that the Planck mock data performs better in constraining the model parameters as compared to the presently available CMB datasets.

We investigate the possibility of constraining primordial non-Gaussianity using the 3D bispectrum of Ly-alpha forest. The strength of the quadratic non-Gaussian correction to an otherwise Gaussian primordial gravitational field is assumed to be dictated by a single parameter fnl. We present the first prediction for bounds on fnl using Ly-alpha flux spectra along multiple lines of sight. The 3D Ly-$\alpha$ transmitted flux field is modeled as a biased tracer of the underlying matter distribution sampled along 1D skewers corresponding to quasars sight lines. The precision to which fnl can be constrained depends on the survey volume, pixel noise and aliasing noise (arising from discrete sampling of the density field). We consider various combinations of these factors to predict bounds on fnl. We find that in an idealized situation of full sky survey and negligible Poisson noise one may constrain fnl ~ 23 in the equilateral limit. Assuming a Ly-alpha survey covering large parts of the sky (k_{min} = 8 * 10^{-4} Mpc^{-1}) and with a quasar density of \bar n = 5 * 10^{-3} Mpc^{-2} it is possible to constrain fnl ~ 100 for equilateral configurations. The possibility of measuring fnl at a precision comparable to LSS studies maybe useful for joint constraining of inflationary scenarios using different data sets.

In single field inflationary models, preheating refers to the phase that immediately follows inflation, but precedes the epoch of reheating. During this phase, the inflaton typically oscillates at the bottom of its potential and gradually transfers its energy to radiation. At the same time, the amplitude of the fields coupled to the inflaton may undergo parametric resonance and, as a consequence, explosive particle production can take place. A priori, these phenomena could lead to an amplification of the super-Hubble scale curvature perturbations which, in turn, would modify the standard inflationary predictions. However, remarkably, it has been shown that, although the Mukhanov-Sasaki variable does undergo narrow parametric instability during preheating, the amplitude of the corresponding super-Hubble curvature perturbations remain constant. Therefore, in single field models, metric preheating does not affect the power spectrum of the large scale perturbations. In this article, we investigate the corresponding effect on the scalar bi-spectrum. Using the Maldacena's formalism, we analytically show that, for modes of cosmological interest, the contributions to the scalar bi-spectrum as the curvature perturbations evolve on super-Hubble scales during preheating is completely negligible. Specifically, we illustrate that, certain terms in the third order action governing the curvature perturbations which may naively be expected to contribute significantly are exactly canceled by other contributions to the bi-spectrum. We corroborate selected analytical results by numerical investigations. We conclude with a brief discussion of the results we have obtained.

Certain oscillatory features in the primordial scalar power spectrum are known to provide a better fit to the outliers in the cosmic microwave background data near the multipole moments of $\ell=22$ and 40. These features are usually generated by introducing a step in the popular, quadratic potential describing the canonical scalar field. Such a model will be ruled out, if the tensors remain undetected at a level corresponding to a tensor-to-scalar ratio of, say, $r\simeq 0.1$. In this work, in addition to the popular quadratic potential, we investigate the effects of the step in a small field model and a tachyon model. With possible applications to future datasets (such as PLANCK) in mind, we evaluate the tensor power spectrum exactly, and include its contribution in our analysis. We compare the models with the WMAP (five as well as seven-year), the QUaD and the ACBAR data. As expected, a step at a particular location and of a suitable magnitude and width is found to improve the fit to the outliers (near $\ell=22$ and 40) in all these cases. We point out that, if the tensors prove to be small (say, $r\lesssim 0.01$), the quadratic potential and the tachyon model will cease to be viable, and more attention will need to be paid to examples such as the small field models.