We show that the integral method used to prove the no-hair theorem for Black Holes (BH's) in spherically symmetric and static spacetimes within the framework of general relativity with matter composed by a complex-valued scalar-field does not lead to a straightforward conclusion about the absence of hair in the stationary and rotating (axisymmetric) scenario. We argue that such a failure can be used to justify in a simple and heuristic way the existence of non-trivial boson clouds or hair found numerically by Herdeiro and Radu [1,2] and analytically by Hod in the test field limit [3-5]. This is due to the presence of a contribution that is negative when rotation exists which allows for an integral to vanish even when a non-trivial boson hair is present. The presence of such a negative contribution that depends on the rotation properties of the BH is perfectly correlated with the eigenvalue problem associated with the boson-field equation. Conversely, when the rotation is absent the integral turns to be composed only by non negative (i.e. positive semidefinite) terms and thus the only way it can vanish is when the hair is completely absent. This analysis poses serious challenges and obstructions towards the elaboration of no-hair theorems for more general spacetimes endowed with a BH region even when including matter fields that obey the energy conditions. Thus rotating boson stars, if collapsed, may lead indeed to a new type of rotating BH, like the ones found in [1,2]. In order to achieve this analysis we solve numerically the eigenvalue problem for the boson field in the Kerr-BH background by imposing rigorous regularity conditions at the BH horizon for the non-extremal case ($0<a<M$) which include the near extremal one in the limit $M\rightarrow a$, as well as the small BH limit $M\rightarrow a\rightarrow 0$.

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

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.

Using a 3+1 decomposition of spacetime, we derive a new formula to compute the gravitational light shifts as measured by two observers which are normal to the spacelike hypersurfaces defining the foliation. This formula is quite general and is also independent of the existence of Killing fields. Known examples are considered to illustrate the usefulness of the formula. In particular, we focus on the Sachs-Wolfe effect that arises in a perturbed Friedman-Robertson-Walker cosmology.

Quantum fields in compact stars can be amplified due to a semiclassical instability. This generic feature of scalar fields coupled to curvature may affect the birth and the equilibrium structure of relativistic stars. We point out that the semiclassical instability has a classical counterpart, which occurs exactly in the same region of the parameter space. For negative values of the coupling parameter the instability is equivalent to the well-known "spontaneous scalarization" effect: the plausible end-state of the instability is a static, asymptotically flat equilibrium configuration with nonzero expectation value for the quantum fields, which is compatible with experiments in the weak-field regime and energetically favored over stellar solutions in general relativity. For positive values of the coupling parameter the new configurations are energetically disfavored, and the end-point of the instability remains an open and interesting issue. The vacuum instability may provide a natural mechanism to produce spontaneous scalarization, leading to new experimental opportunities to probe the nature of vacuum energy via astrophysical observations of compact stars.

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.

Two first order strongly hyperbolic formulations of scalar-tensor theories of gravity (STT) allowing nonminimal couplings (Jordan frame) are presented along the lines of the 3+1 decomposition of spacetime. One is based on the Bona-Masso formulation while the other one employs a conformal decomposition similar to that of Baumgarte-Shapiro-Shibata-Nakamura. A modified Bona-Masso slicing condition adapted to the scalar-tensor theory is proposed for the analysis. This study confirms that the scalar-tensor theory posses a well posed Cauchy problem even when formulated in the Jordan frame.

The 3+1 formulation of scalar-tensor theories of gravity (STT) is obtained in the physical (Jordan) frame departing from the 4+0 covariant field equations. Contrary to the common belief (folklore), the new system of ADM-like equations shows that the Cauchy problem of STT is well formulated (in the sense that the whole system of evolution equations is of first order in the time-derivative). This is the first step towards a full first order (in time and space) formulation from which a subsequent hyperbolicity analysis (a well-posedness determination) can be performed. Several gauge (lapse and shift) conditions are considered and implemented for STT. In particular, a generalization of the harmonic gauge for STT allows us to prove the well posedness of the STT using a second order analysis which is very similar to the one used in general relativity. Some spacetimes of astrophysical and cosmological interest are considered as specific applications. Several appendices complement the ideas of the main part of the paper.

It was recently shown by Whinnett & Torres \cite{WT} that the phenomenon of spontaneous scalarization (SC) in compact objects (polytropes) was accompanied also by a {\it spontaneous violation of the weak energy condition} (WEC). Notably, by the encounter of negative-energy densities at several star places as measured by a static observer. Here we argue that the negativeness of such densities is not generic of scalar tensor theories of gravity (STT). We support this conclusion by numerical results within a class of STT and by using three realistic models of dense matter. However, we show that the ``angular parts'' of the additional conditions for the WEC to hold $\rho_{\rm eff} + T^{i ({\rm eff})}_{i} \geq 0$ tend to be ``slightly violated'' at the outskirts of the star.

Using the 3+1 formalism of general relativity we obtain the equations governing the dynamics of spherically symmetric spacetimes with arbitrary sources. We then specialize for the case of perfect fluids accompanied by a flow of interacting massless or massive particles (e.g. neutrinos) which are described in terms of relativistic transport theory. We focus in three types of coordinates: 1) isotropic gauge and maximal slicing, 2) radial gauge and polar slicing, and 3) isotropic gauge and polar slicing.

We analyze the observational constraints on the effective Brans-Dicke parameter and on the temporal variation of the effective gravitational constant within the context of the oscillating G model, a cosmological model based on a massive scalar field non-minimally coupled to gravity. We show that these local constraints cannot be satisfied simultaneously once the values of the free parameters entering the model become fixed by the global attributes of our Universe. In particular, we show that the lower observational bound for the effective Brans-Dicke parameter and the upper bound of the variation of the effective gravitational constant lead to a specific value of the oscillation amplitude which lies well below the value required to explain the periodicity of 128 Mpc h^{-1} in the galaxy distribution observed in the pencil beam surveys.

We discuss scenarios in which the galactic dark matter in spiral galaxies is described by a long range coherent field which settles in a stationary configuration that might account for the features of the galactic rotation curves. The simplest possibility is to consider scalar fields, so we discuss in particular, two mechanisms that would account for the settlement of the scalar field in a non-trivial configuration in the absence of a direct coupling of the field with ordinary matter: topological defects, and spontaneous scalarization.

We analyze the possibility that the spatial periodicity of 128 h^{-1}Mpc in the galaxy count number observed recently in deep pencil-beam surveys could be explained by the mere existence of structure with the appropriate scale in the galaxy distribution in the Universe. We simulate a universe where the distribution of galaxies has an intrinsic length scale and then investigate the probability of observing in a given direction a spatial periodicity similar to that observed in the pencil-beam surveys of the galactic polar regions. A statistical analysis shows that this probability is of the order of 10^{-8}, a value which excludes the afore mentioned explanation. This result contrasts with the estimates based on Voronoi foam models which yield probabilities of the order 10^{-1} - 10^{-2}. The conclusion that emerges from these studies put together with the present one, is that although they can be obtained in some models, the presence of an intrinsic scale in the structure is not sufficient to achieve reasonable probabilities for observation of the detected periodicities.

We perform a numerical analysis of the gravitational field of a global monopole coupled nonminimally to gravity, and find that, for some given nonminimal couplings (in constrast with the minimal coupling case), there is an attractive region where bound orbits exist. We exhibit the behavior of the frequency shifts that would be associated with `rotation curves' of stars in circular orbits in the spacetimes of such global monopoles.