We present an up to cubic curvature correction to General Relativity with the following features: (i) its vacuum spectrum solely consists of a graviton and is ghost-free, (ii) it possesses well-behaved black hole solutions which coincide with those of Einsteinian cubic gravity, (iii) its cosmology is well-posed as an initial value problem and, most importantly, (iv) it entails a geometric mechanism triggering an inflationary period in the early universe (driven by radiation) with a graceful exit to a late-time cosmology arbitrarily close to $\Lambda$CDM.

We present a new parameterization for the equation of state (EoS) $\omega_X=P_X/\rho_X$, which can reproduce a $f(R)$-like evolution with a precision between $[0.5\%-0.8\%]$ over the numerical solutions. Also, our proposal can render a variety of popular $f(R)$ models that are considered as viable candidates for the cosmic late time acceleration. By using observational data from baryonic acoustic oscillations, supernovae and cosmic chronometers we investigate the constraints on the new EoS parameters. This proposal set a EoS formulation which can be used in an efficient way and makes a good candidate to be implemented in a variety of surveys in order to test the $f(R)$ generic behaviour.

We analyze the evolution of a Friedmann-Robertson-Walker spacetime within the framework of $f(R)$ metric gravity using an exponential model. We show that $f(R)$ gravity may lead to a vanishing effective cosmological constant in the far future (i.e. $R\rightarrow 0$) and yet produce a transient accelerated expansion at present time with a potentially viable cosmological history. This is in contrast with several $f(R)$ models which, while viable, produce in general a non-vanishing effective cosmological constant asymptotically in time ($R\rightarrow 4\Lambda_{\rm eff}$). We also show that relativistic {stars in asymptotically flat spacetimes can be supported within this framework without encountering any singularity, notably in the Ricci scalar $R$.

We derive a class of non-static inhomogeneous dust solutions in f(R) gravity described by the Lemaitre-Tolman-Bondi (LTB) metric. The field equations are fully integrated for all parameter subcases and compared with analogous subcases of LTB dust solutions of GR. Since the solutions do not admit regular symmetry centres, we have two possibilities: (i) a spherical dust cloud with angle deficit acting as the source of a vacuum Schwarzschild-like solution associated with a global monopole, or (ii) fully regular dust wormholes without angle deficit, whose rest frames are homeomorphic to the Schwarzschild-Kruskal manifold or to a 3d torus. The compatibility between the LTB metric and generic f(R) ansatzes furnishes an "inverse procedure" to generate LTB solutions whose sources are found from the f(R) geometry. While the resulting fluids may have an elusive physical interpretation, they can be used as exact non--perturbative toy models in theoretical and cosmological applications of f(R) theories.

We discuss with a rather critical eye the current situation of black hole (BH) solutions in $f(R)$ gravity and shed light about its geometrical and physical significance. We also argue about the meaning, existence or lack thereof of a Birkhoff's theorem in this kind of modified gravity. We focus then on the analysis and quest of $non-trivial$ (i.e. hairy) $asymptotically\,\,flat$ (AF) BH solutions in static and spherically symmetric (SSS) spacetimes in vacuum having the property that the Ricci scalar does $not$ vanish identically in the domain of outer communication. To do so, we provide and enforce the $regularity\,\,conditions$ at the horizon in order to prevent the presence of singular solutions there. Specifically, we consider several classes of $f(R)$ models like those proposed recently for explaining the accelerated expansion in the universe and which have been thoroughly tested in several physical scenarios. Finally, we report analytical and numerical evidence about the $absence$ of $geometric\,\,hair$ in AFSSSBH solutions in those $f(R)$ models. First, we submit the models to the available no-hair theorems, and in the cases where the theorems apply, the absence of hair is demonstrated analytically. In the cases where the theorems do not apply, we resort to a numerical analysis due to the complexity of the non-linear differential equations. Within that aim, a code to solve the equations numerically was built and tested using well know exact solutions. In a future investigation we plan to analyze the problem of hair in De Sitter and Anti-De Sitter backgrounds.

We perform the two$-$point diagnostic for the $Om(z)$ function proposed by Sahni ${\it et al}$ in 2014 for the Starobinsky and Hu & Sawicki models in $f(R)$ gravity. We show that the observed values of the $Omh^2$ function can be explained in $f(R)$ models while in LCDM the $Omh^2$ funticon is expected to be a redshift independent number. We perform the analysis for some particular values of $\Omega_m^0$ founding a cumulative probability ($P(\chi^2 \leq \chi^2_{{\it model}})$) $P \sim 0.16$ or $\sim0.09$ for the better cases versus a cumulative probability of $P \sim 0.98$ in the $\Lambda$CDM scenario. We also show that these models present a characteristic signature around the interval between $z\sim 2$ and $z\sim 4$, that could be confronted with future observations using the same test.

A general formulation to compute habitable zones (HZ) around binary stars is presented. A HZ in this context must satisfy two separate conditions: a radiative one and one of dynamical stability. For the case of single stars, the usual concept of circumstellar habitable zone is based on the radiative condition only, as the dynamical stability condition is taken for granted. For the radiative condition, we extend the simple formulation of the circumstellar habitable zone for single stars, to the case of eccentric stellar binary systems, where two sources of luminosity at different orbital phases contribute to the irradiance of their planetary circumstellar and circumbinary regions. Our approach considers binaries with eccentric orbits and guarantees that orbits in the computed habitable zone remain within it at all orbital phases. For the dynamical stability condition, we use the approach of invariant loops developed by Pichardo et al. 2005 to find regions of stable, non-intersecting orbits, which is a robust method to find stable regions in binary stars, as it is based in the existence of integrals of motion. We apply the combined criteria to calculate HZ for 64 binary stars in the solar neighborhood with known orbital parameters, including some with discovered planets. Formulae and interpolating tables are provided, so the reader can compute the boundaries of the HZ for an arbitrary binary system, using the stellar flux limits they prefer. Together with the formulae provided for stable zones, these allow the computation of both regions of stability and habitability around any binary stellar system. We found 56% of the cases we consider can satisfy both restrictions, this is a very important constriction to binary systems. Nevertheless, we conclude that these systems where a dynamical and radiative safe zone exists, must be considered strong candidates in the search for habitable planets

We provide further numerical evidence which shows that R^n models in f(R) metric gravity whether produces a late time acceleration in the Universe or a matter domination era (usually a transient one) but not both. Our results confirm the findings of Amendola et al. (2007), but using a different approach that avoids the mapping to scalar-tensor theories of gravity, and therefore, dispense us from any discussion or debate about frames (Einstein vs Jordan) which are endemic in this subject. This class of models has been used extensively in the literature as an alternative to the dark energy, but should be considered ruled out for being inconsistent with observations. Finally, we discuss a caveat in the analysis by Faraoni (2011), which was used to further constrain these models by using a chameleon mechanism.

Using an approach that treats the Ricci scalar itself as a degree of freedom, we analyze the cosmological evolution within an f(R) model that has been proposed recently (exponential gravity) and that can be viable for explaining the accelerated expansion and other features of the Universe. This approach differs from the usual scalar-tensor method and, among other things, it spares us from dealing with unnecessary discussions about frames. It also leads to a simple system of equations which is particularly suited for a numerical analysis.

Using the results of Pichardo \et (2005,2008), we determine regions of dynamical stability where planets (or discs in general) could survive in stable orbits around binary stellar systems. We produce this study for 161 binary stars in the Solar neighborhood with known orbital parameters. Additionally, we constructed numerically the discs (invariant loops) around five binary systems with known orbital parameters and with confirmed planets: HIP 10138, HIP 4954, HIP 67275, HIP 116727 and Kepler 16, as a test to the approximation of Pichardo et al. (2005,2008). In each single case, the reported position of the planets lay within our calculated stability regions. This study intends to provide a guide in the search for planets around binary systems with well know orbital parameters, since our method defines precise limits for the stable regions, where discs may have established and planets formed.

We consider a class of metric f(R) modified gravity theories, analyze them in the context of a Friedmann-Robertson-Walker cosmology and confront the results with some of the known constraints imposed by observations. In particular, we focus in correctly reproducing the matter and effective cosmological constant eras, the age of the Universe, and supernovae data. Our analysis differs in many respects from previous studies. First, we avoid any transformation to a scalar-tensor theory in order to be exempted of any potential pathologies (e.g. multivalued scalar potentials) and also to evade any unnecessary discussion regarding frames (i.e. Einstein vs Jordan). Second, based on a robust approach, we recast the cosmology equations as an initial value problem subject to a modified Hamiltonian constraint. Third, we solve the equations numerically where the Ricci scalar itself is one of the variables, and use the constraint equation to monitor the accuracy of the solutions. We compute the "equation of state" (EOS) associated with the modifications of gravity using several inequivalent definitions that have been proposed in the past and analyze it in detail. We argue that one of these definitions has the best features. In particular, we present the EOS around the so called "phantom divide" boundary and compare it with previous findings.

We consider metric f(R) theories of gravity without mapping them to their scalar-tensor counterpart, but using the Ricci scalar itself as an "extra" degree of freedom. This approach avoids then the introduction of a scalar-field potential that might be ill defined (not single valued). In order to explicitly show the usefulness of this method, we focus on static and spherically symmetric spacetimes and deal with the recent controversy about the existence of extended relativistic objects in certain class of f(R) models.