The understanding of stellar structure represents the crossroads of our theories of the nuclear force and the gravitational interaction under the most extreme conditions observably accessible. It provides a powerful probe of the strong field regime of General Relativity, and opens fruitful avenues for the exploration of new gravitational physics. The latter can be captured via modified theories of gravity, which modify the Einstein-Hilbert action of General Relativity and/or some of its principles. These theories typically change the stellar structure equations, thus having a large impact on the astrophysical properties of the corresponding stars and opening a new window to constrain these theories with present and future observations. For relativistic (neutron) stars, the uncertainty on the equation of state of matter at supranuclear densities intertwines with the new parameters of the modified gravity side, providing new phenomenology for the predictions of stellar structure models, such as mass-radius relations, maximum masses, or moment of inertia. For non-relativistic stars (white, brown and red dwarfs), the weakening/strengthening of the gravitational force inside astrophysical bodies may induce changes on the star's mass, radius or luminosity, having an impact, for instance, in the Chandrasekhar's limit for white dwarfs, or in the minimum mass for stable hydrogen burning in brown dwarfs. This work aims to provide a broad overview of the main such results achieved in the recent literature, by combining the results and constraints obtained from the analysis of relativistic and non-relativistic stars in different scenarios. Moreover, we will build a bridge between the efforts of the community working on different theories, formulations, types of stars, theoretical modellings, and observational aspects, highlighting some of the most promising opportunities in the field.

General Relativity yields an analytical prediction of a minimum required mass of roughly $\sim 0.08-0.09 M_{\odot}$ for a star to stably burn sufficient hydrogen to fully compensate photospheric losses and, therefore, to belong to the main sequence. Those objects below this threshold (brown dwarfs) eventually cool down without any chance to stabilize their internal temperature. In this work we consider quadratic Palatini $f(\mathcal{R})$ gravity and show that the corresponding newtonian hydrostatic equilibrium equation contains a new term whose effect is to introduce a weakening/strenghtening of the gravitational interaction inside astrophysical bodies. This fact modifies the General Relativity prediction for this minimum main sequence mass. Through a crude analytical modelling we use this result in order to constraint a combination of the quadratic $f(\mathcal{R})$ gravity parameter and the central density according to astrophysical observations.

Using numerical methods, we investigate the absorption properties of a family of nonsingular solutions {which arise in different metric-affine theories, such as quadratic and Born-Infeld gravity.} These solutions continuously interpolate between Schwarzschild black holes and naked solitons with wormhole topology. The resulting spectrum is characterized by a series of quasibound states excitations, associated with the existence of a stable photonsphere.

We present a novel theory of gravity, namely, an extension of symmetric teleparallel gravity. This is done by introducing a new class of theories where the nonmetricity $Q$ is coupled nonminimally to the matter Lagrangian. This nonminimal coupling entails the nonconservation of the energy-momentum tensor, and consequently the appearance of an extra force. We also present several cosmological applications.

We study the evolution of gravitational waves for non-singular cosmological solutions within the framework of Born-Infeld inspired gravity theories, with special emphasis on the Eddington-inspired Born-Infeld theory. We review the existence of two types of non-singular cosmologies, namely bouncing and asymptotically Minkowski solutions, from a perspective that makes their features more apparent. We study in detail the propagation of gravitational waves near these non-singular solutions and carefully discuss the origin and severity of the instabilities and strong coupling problems that appear. We also investigate the role of the adiabatic sound speed of the matter sector in the regularisation of the gravitational waves evolution. We extend our analysis to more general Born-Infeld inspired theories where analogous solutions are found. As a general conclusion, we obtain that the bouncing solutions are generally more prone to instabilities, while the asymptotically Minkowski solutions can be rendered stable, making them appealing models for the early universe.

We present a novel theory of gravity by considering an extension of symmetric teleparallel gravity. This is done by introducing, in the framework of the metric-affine formalism, a new class of theories where the nonmetricity $Q$ is nonminimally coupled to the matter Lagrangian. More specifically, we consider a Lagrangian of the form $L \sim f_1(Q) + f_2(Q) L_M$, where $f_1$ and $f_2$ are generic functions of $Q$, and $L_M$ is the matter Lagrangian. This nonminimal coupling entails the nonconservation of the energy-momentum tensor, and consequently the appearance of an extra force. The motivation is to verify whether the subtle improvement of the geometrical formulation, when implemented in the matter sector, would allow more universally consistent and viable realisations of the nonminimal curvature-matter coupling theories. Furthermore, we consider several cosmological applications by presenting the evolution equations and imposing specific functional forms of the functions $f_1(Q)$ and $f_2(Q)$, such as power-law and exponential dependencies of the nonminimal couplings. Cosmological solutions are considered in two general classes of models, and found to feature accelerating expansion at late times.

General Relativity has shown an outstanding observational success in the scales where it has been directly tested. However, modifications have been intensively explored in the regimes where it seems either incomplete or signals its own limit of validity. In particular, the breakdown of unitarity near the Planck scale strongly suggests that General Relativity needs to be modified at high energies and quantum gravity effects are expected to be important. This is related to the existence of spacetime singularities when the solutions of General Relativity are extrapolated to regimes where curvatures are large. In this sense, Born-Infeld inspired modifications of gravity have shown an extraordinary ability to regularise the gravitational dynamics, leading to non-singular cosmologies and regular black hole spacetimes in a very robust manner and without resorting to quantum gravity effects. This has boosted the interest in these theories in applications to stellar structure, compact objects, inflationary scenarios, cosmological singularities, and black hole and wormhole physics, among others. We review the motivations, various formulations, and main results achieved within these theories, including their observational viability, and provide an overview of current open problems and future research opportunities.

We study a general class of gravitational theories formulated in the Palatini approach and derive the equations governing the evolution of tensor perturbations. In the absence of torsion, the connection can be solved as the Christoffel symbols of an auxiliary metric which is non-trivially related to the space-time metric. We then consider background solutions corresponding to a perfect fluid and show that the tensor perturbations equations (including anisotropic stresses) for the auxiliary metric around such a background take an Einstein-like form. This facilitates the study in a homogeneous and isotropic cosmological scenario where we explicitly establish the relation between the auxiliary metric and the space-time metric tensor perturbations. As a general result, we show that both tensor perturbations coincide in the absence of anisotropic stresses.

In the framework of Born-Infeld inspired gravity theories, which deviates from General Relativity (GR) in the high curvature regime, we discuss the viability of Cosmic Inflation without scalar fields. For energy densities higher than the new mass scale of the theory, a gravitating dust component is shown to generically induce an accelerated expansion of the Universe. Within such a simple scenario, inflation gracefully exits when the GR regime is recovered, but the Universe would remain matter dominated. In order to implement a reheating era after inflation, we then consider inflation to be driven by a mixture of unstable dust species decaying into radiation. Because the speed of sound gravitates within the Born-Infeld model under consideration, our scenario ends up being predictive on various open questions of the inflationary paradigm. The total number of e-folds of acceleration is given by the lifetime of the unstable dust components and is related to the duration of reheating. As a result, inflation does not last much longer than the number of e-folds of deceleration allowing a small spatial curvature and large scale deviations to isotropy to be observable today. Energy densities are self-regulated as inflation can only start for a total energy density less than a threshold value, again related to the species' lifetime. Above this threshold, the Universe may bounce thereby avoiding a singularity. Another distinctive feature is that the accelerated expansion is of the superinflationary kind, namely the first Hubble flow function is negative. We show however that the tensor modes are never excited and the tensor-to-scalar ratio is always vanishing, independently of the energy scale of inflation.

Recently, the phenomenology of f(R) gravity has been scrutinized motivated by the possibility to account for the self-accelerated cosmic expansion without invoking dark energy sources. Besides, this kind of modified gravity is capable of addressing the dynamics of several self-gravitating systems alternatively to the presence of dark matter. It has been established that both metric and Palatini versions of these theories have interesting features but also manifest severe and different downsides. A hybrid combination of theories, containing elements from both these two formalisms, turns out to be also very successful accounting for the observed phenomenology and is able to avoid some drawbacks of the original approaches. This article reviews the formulation of this hybrid metric-Palatini approach and its main achievements in passing the local tests and in applications to astrophysical and cosmological scenarios, where it provides a unified approach to the problems of dark energy and dark matter.

We consider two types of modifications of Born-Infeld gravity in the Palatini formulation and explore their dynamics in the early universe. One of these families considers $f(R)$ corrections to the Born-Infeld Lagrangian, which can be seen as modifications of the dynamics produced by the quantum effects of matter, while the other consists on different powers of the elementary building block of the Born-Infeld Lagrangian, which we denote by $|\hat\Omega|$. We find that the two types of nonsingular solutions that arise in the original Born-Infeld theory are also present in these extensions, being bouncing solutions a stable and robust branch. Singular solutions with a period of approximate de Sitter inflation are found even in universes dominated by radiation.

We generalize the ultraviolet sector of gravitation via a Born-Infeld action using lessons from massive gravity. The theory contains all of the elementary symmetric polynomials and is treated in the Palatini formalism. We show how the connection can be solved algebraically to be the Levi-Civita connection of an effective metric. The non-linearity of the algebraic equations yields several branches, one of which always reduces to General Relativity at low curvatures. We explore in detail a {\it minimal} version of the theory, for which we study solutions in the presence of a perfect fluid with special attention to the cosmological evolution. In vacuum we recover Ricci-flat solutions, but also an additional physical solution corresponding to an Einstein space. The existence of two physical branches remains for non-vacuum solutions and, in addition, the branch that connects to the Einstein space in vacuum is not very sensitive to the specific value of the energy density. For the branch that connects to the General Relativity limit we generically find three behaviours for the Hubble function depending on the equation of state of the fluid, namely: either there is a maximum value for the energy density that connects continuously with vacuum, or the energy density can be arbitrarily large but the Hubble function saturates and remains constant at high energy densities, or the energy density is unbounded and the Hubble function grows faster than in General Relativity. The second case is particularly interesting because it could offer an interesting inflationary epoch even in the presence of a dust component. Finally, we discuss the possibility of avoiding certain types of singularities within the minimal model.

Motivated by the properties of matter quantum fields in curved space-times, we work out a gravity theory that combines the Born-Infeld gravity Lagrangian with an $f(R)$ piece. To avoid ghost-like instabilities, the theory is formulated within the Palatini approach. This construction provides more freedom to address a number of important questions such as the dynamics of the early universe and the cosmic accelerated expansion, among others. In particular, we consider the effect that adding an $f(R)=a R^2$ term has on the early-time cosmology. We find that bouncing solutions are robust against these modifications of the Lagrangian whereas the solutions with {\it loitering} behavior of the original Born-Infeld theory are very sensitive to the $R^2$ term. In fact, these solutions are modified in such a way that a plateau in the $H^2$ function may arise yielding a period of (approximately) de Sitter inflationary expansion. This inflationary behavior may be found even in a radiation dominated universe.

The dynamical generation of wormholes within an extension of General Relativity (GR) containing (Planck's scale-suppressed) Ricci-squared terms is considered. The theory is formulated assuming the metric and connection to be independent (Palatini formalism) and is probed using a charged null fluid as a matter source. This has the following effect: starting from Minkowski space, when the flux is active the metric becomes a charged Vaidya-type one, and once the flux is switched off the metric settles down into a static configuration such that far from the Planck scale the geometry is virtually indistinguishable from that of the standard Reissner-Nordstr\"om solution of GR. However, the innermost region undergoes significant changes, as the GR singularity is generically replaced by a wormhole structure. Such a structure becomes completely regular for a certain charge-to-mass ratio. Moreover, the nontrivial topology of the wormhole allows to define a charge in terms of lines of force trapped in the topology such that the density of lines flowing across the wormhole throat becomes a universal constant. To the light of our results we comment on the physical significance of curvature divergences in this theory and the topology change issue, which support the view that space-time could have a foam-like microstructure pervaded by wormholes generated by quantum gravitational effects.

It has recently been suggested that Einstein-Rosen (ER) bridges can be interpreted as maximally entangled states of two black holes that form a complex Einstein-Podolsky-Rosen (EPR) pair. This relationship has been dubbed as the ER = EPR correlation. In this work, we consider the latter conjecture in the context of quadratic Palatini theory. An important result, which stems from the underlying assumptions about the geometry on which the theory is constructed, is the fact that all the charged solutions of the quadratic Palatini theory possess a wormhole structure. Our results show that spacetime may have a foam-like microstructure with wormholes generated by fluctuations of the quantum vacuum. This involves the spontaneous creation/annihilation of entangled particle-antiparticle pairs, existing in a maximally entangled state connected by a non-traversable wormhole. Since the particles are produced from the vacuum and therefore exist in a singlet state, they are necessarily entangled with one another. This gives further support to the ER=EPR claim.

We give precise details to support that observations of gravitational lensing at scales of individual, groups and clusters of galaxies can be understood in terms of non-Newtonian gravitational interactions with a relativistic structure compatible with the Einstein Equivalence Principle. This result is derived on very general grounds without knowing the underlying structure of the gravitational field equations. As such, any developed gravitational theory built to deal with these astrophysical scales needs to reproduce the obtained results of this article.

The well-formulation and the well-posedness of the Cauchy problem is discussed for {\it hybrid metric-Palatini gravity}, a recently proposed modified gravitational theory consisting of adding to the Einstein-Hilbert Lagrangian an $f(R)$ term constructed {\it \`{a} la} Palatini. The theory can be recast as a scalar-tensor one predicting the existence of a light long-range scalar field that evades the local Solar System tests and is able to modify galactic and cosmological dynamics, leading to the late-time cosmic acceleration. In this work, adopting generalized harmonic coordinates, we show that the initial value problem can always be {\it well-formulated} and, furthermore, can be {\it well-posed} depending on the adopted matter sources.

Generally, the dynamics of test particles around galaxies, as well as the corresponding mass deficit, is explained by postulating the existence of a hypothetical dark matter. In fact, the behavior of the rotation curves shows the existence of a constant velocity region, near the baryonic matter distribution, followed by a quick decay at large distances. In this work, we consider the possibility that the behavior of the rotational velocities of test particles gravitating around galaxies can be explained within the framework of the recently proposed hybrid metric-Palatini gravitational theory. The latter is constructed by modifying the metric Einstein-Hilbert action with an $f(R)$ term in the Palatini formalism. It was shown that the theory unifies local constraints and the late-time cosmic acceleration, even if the scalar field is very light. In the intermediate galactic scale, we show explicitly that in the hybrid metric-Palatini model the tangential velocity can be explicitly obtained as a function of the scalar field of the equivalent scalar-tensor description. The possibility of constraining the form of the scalar field and the parameters of the model by using the stellar velocity dispersions is also analyzed. Furthermore, the Doppler velocity shifts are also obtained in terms of the scalar field. All the physical and geometrical quantities and the numerical parameters in the hybrid metric-Palatini model can be expressed in terms of observable/measurable parameters, such as the tangential velocity, the baryonic mass of the galaxy, the Doppler frequency shifts, and the stellar dispersion velocity, respectively. Therefore, the obtained results open the possibility of testing the hybrid metric-Palatini gravitational models at the galactic or extra-galactic scale by using direct astronomical and astrophysical observations.

The non-equivalence between the metric and Palatini formalisms of $f(R)$ gravity is an intriguing feature of these theories. However, in the recently proposed hybrid metric-Palatini gravity, consisting of the superposition of the metric Einstein-Hilbert Lagrangian with an $f(\cal R)$ term constructed \`{a} la Palatini, the "true" gravitational field is described by the interpolation of these two non-equivalent approaches. The theory predicts the existence of a light long-range scalar field, which passes the local constraints and affects the galactic and cosmological dynamics. Thus, the theory opens new possibilities for a unified approach, in the same theoretical framework, to the problems of dark energy and dark matter, without distinguishing a priori matter and geometric sources, but taking their dynamics into account under the same standard.

We find that the end state of black hole evaporation could be represented by non-singular and without event horizon stable solitonic remnants with masses of the order the Planck scale and up to 16 units of charge. Though these objects are locally indistinguishable from spherically symmetric, massive electric (or magnetic) charges, they turn out to be sourceless geons containing a wormhole generated by the electromagnetic field. Our results are obtained by interpreting semiclassical corrections to Einstein's theory in the first-order (Palatini) formalism, which yields second-order equations and avoids the instabilities of the usual (metric) formulation of quadratic gravity. We also discuss the potential relevance of these solutions for primordial black holes and the dark matter problem.

A new class of modified theories of gravity, consisting of the superposition of the metric Einstein-Hilbert Lagrangian with an $f(\cal R)$ term constructed \`{a} la Palatini was proposed recently. The dynamically equivalent scalar-tensor representation of the model was also formulated, and it was shown that even if the scalar field is very light, the theory passes the Solar System observational constraints. Therefore the model predicts the existence of a long-range scalar field, modifying the cosmological and galactic dynamics. An explicit model that passes the local tests and leads to cosmic acceleration was also obtained. In the present work, it is shown that the theory can be also formulated in terms of the quantity $X\equiv \kappa^2 T+R$, where T and R are the traces of the stress-energy and Ricci tensors, respectively. The variable X represents the deviation with respect to the field equation trace of general relativity. The cosmological applications of this hybrid metric-Palatini gravitational theory are also explored, and cosmological solutions coming from the scalar-tensor representation of f(X)-gravity are presented. Criteria to obtain cosmic acceleration are discussed and the field equations are analyzed as a dynamical system. Several classes of dynamical cosmological solutions, depending on the functional form of the effective scalar field potential, describing both accelerating and decelerating Universes are explicitly obtained. Furthermore, the cosmological perturbation equations are derived and applied to uncover the nature of the propagating scalar degree of freedom and the signatures these models predict in the large-scale structure.

Recently, a modified theory of gravity was presented, which consists of the superposition of the metric Einstein-Hilbert Lagrangian with an $f(\cal R)$ term constructed \`{a} la Palatini. The theory possesses extremely interesting features such as predicting the existence of a long-range scalar field, that explains the late-time cosmic acceleration and passes the local tests, even in the presence of a light scalar field. In this brief report, we consider the possibility that wormholes are supported by this hybrid metric-Palatini gravitational theory. We present here the general conditions for wormhole solutions according to the null energy conditions at the throat and find specific examples. In the first solution, we specify the redshift function, the scalar field and choose the potential that simplifies the modified Klein-Gordon equation. This solution is not asymptotically flat and needs to be matched to a vacuum solution. In the second example, by adequately specifying the metric functions and choosing the scalar field, we find an asymptotically flat spacetime.

We present a novel approach to modified theories of gravity that consists of adding to the Einstein-Hilbert Lagrangian an f(R) term constructed a la Palatini. Using the respective dynamically equivalent scalar-tensor representation, we show that the theory can pass the Solar System observational constraints even if the scalar field is very light. This implies the existence of a long-range scalar field, which is able to modify the cosmological and galactic dynamics, but leaves the Solar System unaffected. We also verify the absence of instabilities in perturbations and provide explicit models which are consistent with local tests and lead to the late-time cosmic acceleration.

We briefly review the need to perform renormalization of inflationary perturbations to properly work out the physical power spectra. We also summarize the basis of (momentum-space) renormalization in curved spacetime and address several misconceptions found in recent literature on this subject.

We point out that if quantum field renormalization is taken into account the predictions of slow-roll inflation for both the scalar and tensorial power spectra change significantly for wavelengths that today are at observable scales

We study the generation of primordial perturbations in a (single-field) slow-roll inflationary universe. In momentum space, these (Gaussian) perturbations are characterized by a zero mean and a non-zero variance $\Delta^2(k, t)$. However, in position space the variance diverges in the ultraviolet. The requirement of a finite variance in position space forces one to regularize $\Delta^2(k, t)$. This can (and should) be achieved by proper renormalization in an expanding universe in a unique way. This affects the predicted scalar and tensorial power spectra (evaluated when the modes acquire classical properties) for wavelengths that today are at observable scales. As a consequence, the imprint of slow-roll inflation on the CMB anisotropies is significantly altered. We find a non-trivial change in the consistency condition that relates the tensor-to-scalar ratio $r$ to the spectral indices. For instance, an exact scale-invariant tensorial power spectrum, $n_t=0$, is now compatible with a non-zero ratio $r\approx 0.12\pm0.06$, which is forbidden by the standard prediction ($r=-8n_t$). The influence of relic gravitational waves on the CMB may soon come within the range of planned measurements, offering a non-trivial test of the new predictions.

We point out that if quantum field renormalization is taken into account, and the counterterms are evaluated at the Hubble-radius crossing time or few e-foldings after it, the predictions of slow-roll inflation for both the scalar and tensorial power spectrum change significantly. This leads to a change in the consistency condition that relates the tensor-to-scalar amplitude ratio with spectral indices. A reexamination of the potentials $\bf{\phi^2, \phi^4}$, shows that both are compatible with five-year WMAP data. Only when the counterterms are evaluated at much larger times beyond the end of inflation one recovers the standard predictions. The alternative predictions presented here may soon come within the range of measurement of near-future experiments.

We find that the amplitude of quantum fluctuations of the invariant de Sitter vacuum coincides exactly with that of the vacuum of a comoving observer for a massless scalar (inflaton) field. We propose redefining the actual physical power spectrum as the difference between the amplitudes of the above vacua. An inertial particle detector continues to observe the Gibbons-Hawking temperature. However, although the resulting power spectrum is still scale-free, its amplitude can be drastically reduced since now, instead of the Hubble's scale at the inflationary period, it is determined by the square of the mass of the inflaton fluctuation field.

We discuss two aspects of f(R) theories of gravity in metric formalism. We first study the reasons to introduce a scalar-tensor representation for these theories and the behavior of this representation in the limit to General Relativity, f(R)--> R. We find that the scalar-tensor representation is well behaved even in this limit. Then we work out the exact equations for spherically symmetric sources using the original f(R) representation, solve the linearized equations, and compare our results with recent calculations of the literature. We observe that the linearized solutions are strongly affected by the cosmic evolution, which makes very unlikely that the cosmic speedup be due to f(R) models with correcting terms relevant at low curvatures.

We study modified theories of gravity of the f(R) type in Palatini formalism. For a generic f(R) lagrangian, we show that the metric can be solved as the product of a scalar function times a rank-two tensor (or auxiliary metric). The scalar function is sensitive to the local energy-momentum density. The auxiliary metric satisfies a set of equations very similar to Einstein's equations and, for weak sources, it can be approximated by the Minkowski metric. According to this, the metric coupled to the matter strongly departs from the Minkowskian one in the neighbourhood of any microscopic physical system. As a consequence, new gravitationally-induced interactions arise and lead to observable effects at microscopic and macroscopic scales. In particular, test body trajectories experience self-accelerations which depend on the internal structure and composition of the body. These facts make very unlikely the viability of Palatini f(R) models designed to change the late-time cosmic evolution.

We compute the complete post-Newtonian limit of the Palatini form of f(R) gravities using a scalar-tensor representation. By comparing the predictions of these theories with laboratory and solar system experiments, we find a set of inequalities that any lagrangian f(R) must satisfy. The constraints imposed by those inequalities allow us to find explicit functions that bound from above and from below the possible nonlinear terms of the lagrangian. We conclude that the lagrangian f(R) must be almost linear in R and that nonlinear corrections that grow at low curvatures are incompatible with observations. This result shows that modifications of gravity at very low cosmic densities cannot be responsible for the observed cosmic speed-up.

We compute the complete post-Newtonian limit of the metric form of f(R) gravities using a scalar-tensor representation. By comparing the predictions of these theories with laboratory and solar system experiments, we find a set of inequalities that any lagrangian f(R) must satisfy. The constraints imposed by those inequalities allow us to find explicit bounds to the possible nonlinear terms of the lagrangian. We conclude that the lagrangian f(R) must be almost linear in R and that corrections that grow at low curvatures are incompatible with observations. This result shows that modifications of gravity at very low cosmic densities cannot be responsible for the observed cosmic speed-up.

In this work we show that the gravity lagrangian f(R) at relatively low curvatures in both metric and Palatini formalisms is a bounded function that can only depart from the linearity within the limits defined by well known functions. We obtain those functions by analysing a set of inequalities that any f(R) theory must satisfy in order to be compatible with laboratory and solar system observational constraints. This result implies that the recently suggested f(R) gravity theories with nonlinear terms that dominate at low curvatures are incompatible with observations and, therefore, cannot represent a valid mechanism to justify the cosmic speed-up.