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.

We discuss the coupling of the electromagnetic field with a curved and torsioned Lyra manifold using the Duffin-Kemmer-Petiau theory. We will show how to obtain the equations of motion and energy-momentum and spin density tensors by means of the Schwinger Variational Principle.

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

We develop a model for spiral galaxies based on a nonlinear realization of the Newtonian dynamics starting from the momentum and mass conservations in the phase space. The radial solution exhibits a rotation curve in qualitative accordance with the observational data.

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.