We write down the Lagrangian bias expansion in general relativity up to 4th order in terms of operators describing the curvature of an early-time hypersurface for comoving observers. They can be easily expanded in synchronous or comoving gauges. This is necessary for the computation of the one-loop halo bispectrum, where relativistic effects can be degenerate with a primordial non-Gaussian signal. Since the bispectrum couples scales, an accurate prediction of the squeezed limit behavior needs to be both non-linear and relativistic. We then evolve the Lagrangian bias operators in time in comoving gauge, obtaining non-local operators analogous to what is known in the Newtonian limit. Finally, we show how to renormalize the bias expansion at an arbitrary time and find that this is crucial in order to cancel unphysical $1/k^2$ divergences in the large-scale power spectrum and bispectrum that could be mistaken for a contamination to the non-Gaussian signal.

The large scale structure bispectrum in the squeezed limit couples large with small scales. Since relativity is important at large scales and non-linear loop corrections are important at small scales, the proper calculation of the observed bispectrum in this limit requires a non-linear relativistic calculation. We compute the matter bispectrum in general relativity in the weak field approximation. The calculation is as involved as existing second-order results. We find several differences with the Newtonian calculation such as the non-cancellation of IR divergences, the need to renormalize the background, and the fact that initial conditions must be set at second order in perturbation theory. For the bispectrum, we find relativistic corrections to be as large as the newtonian result in the squeezed limit. In that limit relativistic one-loop contributions, which we compute for the first time, can be as large as tree level results and have the same $1/k^2$ dependence as a primordial local non-Gaussianity signal where $k$ is the momentum approaching zero. Moreover, we find the time dependence of the relavistic corrections to the bispectrum to be the same as that of a primordial non-Gaussianity signal.

We compare supernovae and Baryon Acoustic Oscillations data to the predictions of a cosmological model of interacting dark matter and dark energy. This theoretical model can be derived from the effective field theory of Einstein-Cartan gravity with two scaling exponents $\delta_G$ and $\delta_{\Lambda}$, related to the interaction between dark matter and dark energy. We perform a $\chi^2$ fit to the data to compare and contrast it with the standard $\Lambda$CDM model. We then explore the range of parameter of the model which gives a better $\chi^2$ than the standard cosmological model. All those results lead to tight constraints on the scaling exponents of the model. Our conclusion is that this class of models, provides a decent alternative to the $\Lambda$CDM model.

Loop quantum cosmology tries to capture the main ideas of loop quantum gravity and to apply them to the Universe as a whole. Two main approaches within this framework have been considered to date for the study of cosmological perturbations: the dressed metric approach and the deformed algebra approach. They both have advantages and drawbacks. In this article, we accurately compare their predictions. In particular, we compute the associated primordial tensor power spectra. We show -- numerically and analytically -- that the large scale behavior is similar for both approaches and compatible with the usual prediction of general relativity. The small scale behavior is, the other way round, drastically different. Most importantly, we show that in a range of wavenumbers explicitly calculated, both approaches do agree on predictions that, in addition, differ from standard general relativity and do not depend on unknown parameters. These features of the power spectrum at intermediate scales might constitute a universal loop quantum cosmology prediction that can hopefully lead to observational tests and constraints. We also present a complete analytical study of the background evolution for the bouncing universe that can be used for other purposes.