We develop an action principle to construct the dynamics that give rise to a minimal generalization of Einstein's equations, where the speed of light ($c$), the gravitational constant ($G$) and the cosmological constant ($\Lambda$) are allowed to vary. Our construction preserves general covariance of the theory, which yields a general dynamical constraint on $c$, $G$ and $\Lambda$. This action is general and can be applied to describe different cosmological solutions. We apply this formulation to the initial condition puzzles of the early universe and show that it generates a dynamical mechanism to obtain the homogeneous and flat universe we observe today. We rewrite the conditions necessary to solve the horizon and flatness problems in this framework, which does not necessarily lead to an accelerated expansion as in inflation. Then, we show how the dynamics of the scalar field that represents $c$ or $G$ (and $\Lambda$) can be used to solve the problems of the early universe cosmology by means of different ways to c-inflate the horizon in the early universe. By taking $\Lambda = 0$, we show that the dynamics of the scalar field representing $c$ can be described once a potential is given.

We propose a cosmological scenario involving a scalar field, $\varphi$, that is a source of Dark Matter as well as of Dark Energy. Besides $\varphi$, the Lagrangian of the field theory envisaged in our scenario contains a second field $\chi$, for simplicity assumed to be a scalar, too. For fixed values of $\chi$, the potential term decays exponentially at large positive values of $\varphi$. While $\varphi$ is not coupled to Standard Model fields, $\chi$ is assumed to be coupled to them, and the Green functions of $\chi$ depend on the cosmological redshift in the expanding universe. We assume that the term in the Lagrangian coupling $\chi$ to $\varphi$ is such that, at redshifts $z$ larger than some critical redshift $z_c$, $\varphi$ is trapped near $\varphi=0$, and oscillations of $\varphi$ about $\varphi=0$ describing massive scalar particles give rise to Dark Matter. At redshifts below $z_c$, the field $\varphi$ is no longer trapped near the origin and starts to `roll' towards large field values. A homogenous component of $\varphi$ emerges that acts as Dark Energy. Within over-dense regions, such as galaxies and galaxy clusters, the redshifting of $\chi$ stops, and $\varphi$ therefore remains trapped near $\varphi=0$ as long as $z_c$ is smaller than the redshift when structures on galactic scales decouple from the Hubble flow. Thus, at the present time, $\varphi$ describes both Dark Energy and Dark Matter.

An extension of the Starobinsky model is proposed. Besides the usual Starobinsky Lagrangian, a term proportional to the derivative of the scalar curvature, $\nabla_{\mu}R\nabla^{\mu}R$, is considered. The analyzis is done in the Einstein frame with the introduction of a scalar field and a vector field. We show that inflation is attainable in our model, allowing for a graceful exit. We also build the cosmological perturbations and obtain the leading-order curvature power spectrum, scalar and tensor tilts and tensor-to-scalar ratio. The tensor and curvature power spectrums are compared to the most recent observations from BICEP2/Keck collaboration. We verify that the scalar-to-tensor rate $r$ can be expected to be up to three times the values predicted by Starobinsky model.

We study cosmological perturbations in a model of unified dark matter and dark energy with a sharp transition in the late-time universe. The dark sector is described by a dark fluid which evolves from an early stage at redshifts $z > z_C$ when it behaves as cold dark matter (CDM) to a late time dark energy (DE) phase ($z < z_C$) when the equation of state parameter is $w = -1 + \epsilon$, with a constant $\epsilon$ which must be in the range $0 < \epsilon < 2/3$. We show that fluctuations in the dark energy phase suffer from an exponential instability, the mode functions growing both as a function of comoving momentum $k$ and of conformal time $\eta$. In order that this exponential instability does not lead to distortions of the energy density power spectrum on scales for which we have good observational results, the redshift $z_C$ of transition between the two phases is constrained to be so close to zero that the model is unable to explain the supernova data.

The equivalence between theories depending on the derivatives of $R$, i.e. $f\left( R,\nabla R,...,\nabla^{n}R\right) $, and scalar-multi-tensorial theories is verified. The analysis is done in both metric and Palatini formalisms. It is shown that $f\left( R,\nabla R,...,\nabla^{n}R\right) $ theories are equivalent to scalar-multi-tensorial ones resembling Brans-Dicke theories with kinetic terms $\omega_{0}=0$ and $\omega_{0}= - \frac{3}{2}$ for metric and Palatini formalisms respectively. This result is analogous to what happens for $f(R)$ theories. It is worthy emphasizing that the scalar-multi-tensorial theories obtained here differ from Brans-Dicke ones due to the presence of multiple tensorial fields absent in the last. Furthermore, sufficient conditions are established for $f\left( R,\nabla R,...,\nabla^{n}R\right) $ theories to be written as scalar-multi-tensorial theories. Finally, some examples are studied and the comparison of $f\left( R,\nabla R,...,\nabla^{n}R\right) $ theories to $f\left( R,\Box R,...\Box^{n}R\right) $ theories is performed.

We investigate the influence of massive photons on the evolution of the expanding universe. Two particular models for generalized electrodynamics are considered, namely de Broglie-Proca and Bopp-Podolsky electrodynamics. We obtain the equation of state (EOS) $P=P(\varepsilon)$ for each case using dispersion relations derived from both theories. The EOS are inputted into the Friedmann equations of a homogeneous and isotropic space-time to determine the cosmic scale factor $a(t)$. It is shown that the photon non-null mass does not significantly alter the result $a\propto t^{1/2}$ valid for a massless photon gas; this is true either in de Broglie-Proca's case (where the photon mass $m$ is extremely small) or in Bopp-Podolsky theory (for which $m$ is extremely large).

The propagation of non-linear electromagnetic waves is carefully analyzed on a curved spacetime created by static spherically symmetric mass and charge distribution. We compute how non-linear electrodynamics affects the geodesic deviation and the redshift of photons propagating near this massive charged object. In the first order approximation, the effects of electromagnetic self-interaction can be distinguished from the usual Reissner-Nordstr\"om terms. In the particular case of Euler-Heisenberg effective Lagrangian, we find that these self-interaction effects might be important near extremal compact charged objects.

We propose a phenomenological unified model for dark matter and dark energy based on an equation of state parameter $w$ that scales with the $\arctan$ of the redshift. The free parameters of the model are three constants: $\Omega_{b0}$, $\alpha$ and $\beta$. Parameter $\alpha$ dictates the transition rate between the matter dominated era and the accelerated expansion period. The ratio $\beta / \alpha$ gives the redshift of the equivalence between both regimes. Cosmological parameters are fixed by observational data from Primordial Nucleosynthesis (PN), Supernovae of the type Ia (SNIa), Gamma-Ray Bursts (GRB) and Baryon Acoustic Oscillations (BAO). The calibration of the 138 GRBs events is performed using the 580 SNIa of the Union2.1 data set and a new set of 79 high-redshift GRBs is obtained. The various sets of data are used in different combinations to constraint the parameters through statistical analysis. The unified model is compared to the $\Lambda$CDM model and their differences are emphasized.

This paper analyses the cosmological consequences of a modified theory of gravity whose action integral is built from a linear combination of the Ricci scalar $R$ and a quadratic term in the covariant derivative of $R$. The resulting Friedmann equations are of the fifth-order in the Hubble function. These equations are solved numerically for a flat space section geometry and pressureless matter. The cosmological parameters of the higher-order model are fit using SN Ia data and X-ray gas mass fraction in galaxy clusters. The best-fit present-day $t_{0}$ values for the deceleration parameter, jerk and snap are given. The coupling constant $\beta$\ of the model is not univocally determined by the data fit, but partially constrained by it. Density parameter $\Omega_{m0}$ is also determined and shows weak correlation with the other parameters. The model allows for two possible future scenarios: there may be either a premature Big Rip or a Rebouncing event depending on the set of values in the space of parameters. The analysis towards the past performed with the best-fit parameters shows that the model is not able to accommodate a matter-dominated stage required to the formation of structure.

We construct a phenomenological theory of gravitation based on a second order gauge formulation for the Lorentz group. The model presents a long-range modification for the gravitational field leading to a cosmological model provided with an accelerated expansion at recent times. We estimate the model parameters using observational data and verify that our estimative for the age of the Universe is of the same magnitude than the one predicted by the standard model. The transition from the decelerated expansion regime to the accelerated one occurs recently (at $\sim9.3\;Gyr$).

Early universe equations of state including realistic interactions between constituents are built up. Under certain reasonable assumptions, these equations are able to generate an inflationary regime prior to the nucleosynthesis period. The resulting accelerated expansion is intense enough to solve the flatness and horizon problems. In the cases of curvature parameter \kappa equal to 0 or +1, the model is able to avoid the initial singularity and offers a natural explanation for why the universe is in expansion.

This work is an application of the second order gauge theory for the Lorentz group, where a description of the gravitational interaction is obtained which includes derivatives of the curvature. We analyze the form of the second field strenght, $G=\partial F +fAF$, in terms of geometrical variables. All possible independent Lagrangians constructed with quadratic contractions of $F$ and quadratic contractions of $G$ are analyzed. The equations of motion for a particular Lagrangian, which is analogous to Podolsky's term of his Generalized Electrodynamics, are calculated. The static isotropic solution in the linear approximation was found, exhibiting the regular Newtonian behaviour at short distances as well as a meso-large distance modification.

The high precision attained by cosmological data in the last few years has increased the interest in exact solutions. Analytic expressions for solutions in the Standard Model are presented here for all combinations of $\Lambda = 0$, $\Lambda \ne 0$, $\kappa = 0$ and $\kappa \ne 0$, in the presence and absence of radiation and nonrelativistic matter. The most complete case (here called the $\Lambda \gamma CDM$ Model) has $\Lambda \ne 0, \kappa \ne 0$, and supposes the presence of radiation and dust. It exhibits clearly the recent onset of acceleration. The treatment includes particular models of interest such as the $\Lambda$CDM Model (which includes the cosmological constant plus cold dark matter as source constituents).