Non-gravitational interaction between dark matter and dark energy has been considered in a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) universe. The interaction rate is assumed to be linear in the energy densities of dark matter and dark energy and it is also proportional to the Hubble rate of the FLRW universe. This kind of interaction model leads to an autonomous linear dynamical system, and depending on the coupling parameters, could be solved analytically by calculating the exponential of the matrix, defining the system. We show here that such interaction rate has a very deep connection with some well known cosmological theories. We then investigate the theoretical bounds on the coupling parameters of the interaction rate in order that the energy densities of the dark sector remain positive throughout the evolution of the universe and asymptotically converge to zero at very late times. Our analyses also point out that such linear interacting model may encounter with finite time future singularities depending on the coupling parameters as well as the dark energy state parameter.

We review the well-known Peebles-Vilenkin (PV) quintessential inflation model and discuss its possible improvements in agreement with the recent observations. The improved PV model depends only on two parameters: the inflaton mass $m$, and another smaller mass $M$; where the latter has to be chosen in order to undertake that, at present time, the dark energy density of the universe is approximately about 70\% of the total energy budget of the universe. The value of the inflaton mass $m$ is calculated using the observational value of the power spectrum of the scalar perturbations, and the value of mass $M$, which depends on the reheating temperature, is calculated by solving the corresponding dynamical system whose initial conditions are taken at the matter-radiation equality and are obtained from three observational data: the red shift at the matter-radiation equality, the ratio of the matter energy density to the critical one at the present time and the current value of the Hubble parameter.

The PLANCK collaboration has determined values for the spectral parameters of the CMB radiation, namely the spectral index $n_s$, its running $\alpha_s$, the running of the running $\beta_s$, using a growing body of measurements of CMB anisotropies by the Planck satellite and other missions. These values do not follow the hierarchy of sizes predicted by single field, slow roll inflationary theory, and are thus difficult to fit for such inflation models. In this work we present first a study of 49 single field, slow roll inflationary potentials in which we assess the likelyhood of these models fitting the spectral parameters to their currently most accurate determination given by the PLANCK collaboration. We check numerically with a MATLAB program the spectral parameters that each model can yield for a very broad, comprehensive list of possible parameter and field values. The comparison of spectral parameter values supported by the models with their determinations by the PLANCK collaboration leads to the conclusion that the data provided by PLANCK2015 TT+lowP and PLANCK2015 TT,TE,EE+lowP taking into account the running of the running disfavours 40 of the 49 models with confidence level at least 92.8\%. Next, we discuss the reliability of the current computations of these spectral parameters. We identify a bias in the method of determination of the spectral parameters by least residue parameter fitting (using MCMC or any other scheme) currently used to reconstruct the power spectrum of scalar perturbations. This bias can explain the observed contradiction between theory and observations. Its removal is computationally costly, but necessary in order to compare the forecasts of single field, slow roll theories with observations.

It is shown that teleparallel $F({\mathcal T})$ theories of gravity combined with Loop Quantum Cosmology support a Matter Bounce Scenario which is an alternative to the inflation scenario in the Big Bang paradigm. It is checked thatthese bouncing models provide theoretical data that fits well with the current observational data, allowing the viability of the Matter Bounce Scenario.

The CMB map provided by the Planck project constrains the value of the ratio of tensor-to-scalar perturbations, namely $r$, to be smaller than $0.11$ (95% CL). This bound rules out the simplest models of inflation. However, recent data from BICEP2 is in strong tension with this constrain, as it finds a value $r=0.20^{+0.07}_{-0.05}$ with $r=0$ disfavored at $7.0 \sigma$, which allows these simplest inflationary models to survive. The remarkable fact is that, even though the BICEP2 experiment was conceived to search for evidence of inflation, its experimental data matches correctly theoretical results coming from the matter bounce scenario (the alternative model to the inflationary paradigm). More precisely, most bouncing cosmologies do not pass Planck's constrains due to the smallness of the value of the tensor/scalar ratio $r\leq 0.11$, but with new BICEP2 data some of them fit well with experimental data. This is the case with the matter bounce scenario in the teleparallel version of Loop Quantum Cosmology.