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.

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.

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.

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$).

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.