In a recent work, we had constructed a model consisting of two fields---a canonical scalar field and a non-canonical ghost field---that had sourced a symmetric matter bounce scenario. The model had involved only one parameter, viz. the scale associated with the bounce. For a suitable value of the parameter, the model had led to strictly scale invariant power spectra with a COBE normalized scalar amplitude and a rather small tensor-to-scalar ratio. In this work, we extend the model to achieve near-matter bounces, which contain a second parameter apart from the bounce scale. As the new model does not seem to permit analytical evaluation of the scalar modes near the bounce, with the aid of techniques which we had used in our earlier work, we compute the scalar and the tensor power spectra numerically. For appropriate values of the additional parameter, we find that the model produces red spectra with a scalar spectral tilt and a small tensor-to-scalar ratio which are consistent with the recent observations of the anisotropies in the cosmic microwave background by Planck.

We revisit the evolutions of scalar perturbations in a non-singular Galileon bounce. It is known that the second order differential equation governing the perturbations is numerically unstable at a point called $\gamma$-crossing. This instability is usually circumvented using certain gauge choices. We show that the perturbations can be evolved across this point by solving the first order differential equations governing suitable gauge invariant quantities without any instabilities. We demonstrate this method in a matter bounce scenario described by the Galileon action.

Certain features in the primordial scalar power spectrum are known to improve the fit to the cosmological data. We examine whether bouncing scenarios can remain viable if future data confirm the presence of such features. In inflation, the fact that the trajectory is an attractor permits the generation of features. However, bouncing scenarios often require fine tuned initial conditions, and it is only the ekpyrotic models that allow attractors. We demonstrate, for the first time, that ekpyrotic scenarios can generate specific features that have been considered in the context of inflation.

Matter bounces refer to scenarios wherein the universe contracts at early times as in a matter dominated epoch until the scale factor reaches a minimum, after which it starts expanding. While such scenarios are known to lead to scale invariant spectra of primordial perturbations after the bounce, the challenge has been to construct completely symmetric bounces that lead to a tensor-to-scalar ratio which is small enough to be consistent with the recent cosmological data. In this work, we construct a model involving two scalar fields (a canonical field and a non-canonical ghost field) to drive the symmetric matter bounce and study the evolution of the scalar perturbations in the model. If we consider the scale associated with the bounce to be of the order of the Planck scale, the model is completely described in terms of only one parameter, viz the value of the scale factor at the bounce. We evolve the scalar perturbations numerically across the bounce and evaluate the scalar power spectra after the bounce. We show that, while the scalar and tensor perturbation spectra are scale invariant over scales of cosmological interest, the tensor-to-scalar ratio proves to be much smaller than the current upper bound from the observations of the cosmic microwave background anisotropies by the Planck mission. We also support our numerical analysis with analytical arguments.

We evaluate the dimensionless non-Gaussianity parameter $h_{_{\rm NL}}$, that characterizes the amplitude of the tensor bispectrum, numerically for a class of two field inflationary models such as double inflation, hybrid inflation and aligned natural inflation. We compare the numerical results with the slow roll results which can be obtained analytically. In the context of double inflation, we also investigate the effects on $h_{_{\rm NL}}$ due to curved trajectories in the field space. We explicitly examine the validity of the consistency relation governing the tensor bispectrum in the squeezed limit. Lastly, we discuss the contribution to $h_{_{\rm NL}}$ due to the epoch of preheating in two field models.