Within the effective field theory approach to cosmic acceleration, the background expansion can be specified separately from the gravitational modifications. We explore the impact of modified gravity in a background different from a cosmological constant plus cold dark matter ($\Lambda$CDM) on the stability and cosmological observables, including covariance between gravity and expansion parameters. In No Slip Gravity the more general background allows more gravitational freedom, including both positive and negative Planck mass running. We examine the effects on cosmic structure growth, as well as showing that a viable positive integrated Sachs-Wolfe effect crosscorrelation easily arises from this modified gravity theory. Using current data we constrain parameters with a Monte Carlo analysis, finding a maximum running $|\alpha_M|\lesssim 0.03$. We provide the modified {\tt hi\_class} code publicly on GitHub, now enabling computation and inclusion of the redshift space distortion observable $f\sigma_8$ as well as the No Slip Gravity modifications.

In general, to avoid a singularity in cosmological models involves the introduction of exotic kind of matter fields, for example, a scalar field with negative energy density. In order to have a bouncing solution in classical General Relativity, violation of the energy conditions is required. In this work, we discuss a case of the bouncing solution in the Brans-Dicke theory with radiative fluid that obeys the energy conditions, and with no ghosts.

Generally the Brans-Dicke theory reduces to General Relativity in the limit $\omega\rightarrow\infty$ if the scalar field goes as $\phi\propto1/\omega$. However, it is also known that there are examples with $\phi\propto1/\sqrt{\omega}$ that does not tend to GR. We discuss another case: a homogeneous and isotropic universe filled with stiff matter. The power of time dependence of these solutions do not depend on $\omega$, and there is no General Relativity limit even though we have $\phi\propto1/\omega$. A perturbative analysis of this exotic case is carried out.

We study the formation of classical singularities in Generalized Brans-Dicke theories that are natural extensions to Brans-Dicke where the kinetic term is modified by a new coupling function $\omega(\varphi)$. We discuss the asymptotic limit $\omega(\varphi)\rightarrow\infty$ and show that the system generically does not approach General Relativity. Given the arbitrariness of $\omega(\varphi)$, one can search for coupling functions chosen specifically to avoid classical singularities. However, we prove that this is not the case. Homogeneous and spherically symmetric collapsing objects form singularities for arbitrary coupling functions. On the other hand, expanding cosmological scenarios are completely free of Big Rip type singularities. In an expanding universe, the scalar field behaves at most as stiff matter, which makes these cosmological solutions asymptotically approach General Relativity.

In the present work we constrain three different profiles of a Lema\^itre-Tolman-Bondi model using supernovae type Ia and baryon acoustic oscillation data. We improve common practice in the literature by carefully calibrating the supernovae in the appropriate inhomogeneous background dynamics. In addition, we address subtle issues in order to propagate the primordial BAO scale to present epoch. The combined analysis of BAO+SNIa offers a stringent test for these models. We use two distinct parameter estimation approaches, namely, the $\chi^2$ and the complete likelihood functional. It has been argued that these two approaches are not equivalent and indeed our analysis shows a specific example of their departure.

In a previous work we obtained a set of necessary conditions for the linear approximation in cosmology. Here we discuss the relations of this approach with the so called covariant perturbations. It is often argued in the literature that one of the main advantages of the covariant approach to describe cosmological perturbations is that the Bardeen formalism is coordinate dependent. In this paper we will reformulate the Bardeen approach in a completely covariant manner. For that, we introduce the notion of pure and mixed tensors, which yields an adequate language to treat both perturbative approaches in a common framework. We then stress that in the referred covariant approach one necessarily introduces an additional hyper-surface choice to the problem. Using our mixed and pure tensors approach, we were able to construct a one-to-one map relating the usual gauge dependence of the Bardeen formalism with the hyper-surface dependence inherent to the covariant approach. Finally, through the use of this map, we define full non-linear tensors that at first order correspond to the three known gauge invariant variables $\Phi$, $\Psi$ and $\Xi$, which are simultaneously foliation and gauge invariant. We then stress that the use of the proposed mixed tensors allows one to construct simultaneously gauge and hyper-surface invariant variables at any order.

In previous works, it was shown that the Lagrangians and Hamiltonians of cosmological linear scalar, vector and tensor perturbations of homogeneous and isotropic space-times with flat spatial sections containing a perfect fluid can be put in a simple form through the implementation of canonical transformations and redefinitions of the lapse function, without ever using the background classical equations of motion. In this paper, we generalize this result to general fluids, which includes entropy perturbations, and to arbitrary spacelike hyper-surfaces through a new method together with the Faddeev-Jackiw procedure for the constraint reduction. A simple second order Hamiltonian involving the Mukhanov-Sasaki variable is obtained, again without ever using the background equations of motion.

A quantum cosmological model with radiation and a dilaton scalar field is analysed. The Wheeler-deWitt equation in the mini-superspace induces a Schr\"odinger equation, which can be solved. An explicit wavepacket is constructed for a particular choice of the ordering factor. A consistent solution is possible only when the scalar field is a phantom field. Moreover, although the wavepacket is time dependent, a Bohmian analysis allows to extract a bouncing behaviour for the scale factor.