Abstract: 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 Tolman-Oppenheimer-Volkoff equations of stellar's hydrostatic equilibrium, 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 of different types of stars. For relativistic stars, such as neutron stars, the uncertainty on the equation of state of matter at supranuclear densities intertwines with the new parameters coming from the modified gravity side, providing a whole new phenomenology for the typical predictions of stellar structure models, such as mass-radius relations, maximum masses, or moment of inertia. For non-relativistic stars, such as white, brown and red dwarfs, the weakening/strengthening of the gravitational force inside astrophysical bodies via the modified Newtonian (Poisson) equation may induce changes on the star's mass, radius, central density 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 high-mass brown dwarfs. This work aims to provide a broad overview of the main such results achieved in the recent literature for many such modified theories of gravity, 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 modelings, and observational aspects, highlighting some of the most promising opportunities in the field.

Abstract: We show that the effective brane-world and the loop quantum cosmology background expansion histories can be reproduced from a modified gravity perspective in terms of an f (R) gravity action plus a g(R) term non-minimally coupled with the matter Lagrangian. The reconstruction algorithm that we provide depends on a free function of the matter density that must be specified in each case and allows to obtain analytical solutions always. In the simplest cases, the function f (R) is quadratic in the Ricci scalar, R, whereas g(R) is linear. Our approach is compared with recent results in the literature. We show that working in the Palatini formalism there is no need to impose any constraint that keeps the equations second order, which is a key requirement for the successful implementation of the reconstruction algorithm.

Abstract: We revisit the gauge symmetry related to integrable projective transformations in metric-affine formalism, identifying the gauge field of the Weyl (conformal) symmetry as a dynamical component of the affine connection. In particular, we show how to include the local scaling symmetry as a gauge symmetry of a large class of geometric gravity theories, introducing a compensator dilaton field that naturally gives rise to a Stuckelberg sector where a spontaneous breaking mechanism of the conformal symmetry is at work to generate a mass scale for the gauge field. For Ricci-based gravities that include, among others, General Relativity, f(R) and f(R, R μnu R μnu) theories and the EiBI model, we prove that the on-shell gauge vector associated to the scaling symmetry can be identified with the torsion vector, thus recovering and generalizing conformal invariant theories in the Riemann-Cartan formalism, already present in the literature.

Abstract: We explore equilibrium solutions of spherically symmetric boson stars in the Palatini formulation of f (R) gravity. We account for the modifications introduced in the gravitational sector by using a recently established correspondence between modified gravity with scalar matter and general relativity with modified scalar matter. We focus on the quadratic theory f (R) = R + xi R-2 and compare its solutions with those found in general relativity, exploring both positive and negative values of the coupling parameter xi. As matter source, a complex, massive scalar field with and without self-interaction terms is considered. Our results show that the existence curves of boson stars in Palatini f (R) gravity are fairly similar to those found in general relativity. Major differences are observed for negative values of the coupling parameter which results in a repulsive gravitational component for high enough scalar field density distributions. Adding self-interactions makes the degeneracy between f (R) and general relativity even more pronounced, leaving very little room for observational discrimination between the two theories.

Abstract: 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 similar 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.