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.

Einstein-Cartan theory is an extension of the standard formulation of General Relativity characterized by a non-vanishing torsion. The latter is sourced by the matter fields via the spin tensor, and its effects are expected to be important at very high spin densities. In this work we analyze in detail the physics of Einstein-Cartan theory with Dirac and Maxwell fields minimally coupled to the spacetime torsion. This breaks the $U(1)$ gauge symmetry, which is suggested by the possibility of a torsion-induced phase transition in the early Universe. The resulting Dirac-like and Maxwell-like equations are non-linear with self-interactions as well as having fermion-boson non-minimal couplings. We discuss several cosmological aspects of this theory, including bounces, acceleration phases and matter-antimatter asymmetry in the torsion era, as well as late-time effects such as the generation of an effective cosmological constant, dark energy, and future bounces within cyclic solutions.

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.

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 explore the dynamics and evolution of the Universe at early and late times, focusing on both dark energy and extended gravity models and their astrophysical and cosmological consequences. Modified theories of gravity not only provide an alternative explanation for the recent expansion history of the universe, but they also offer a paradigm fundamentally distinct from the simplest dark energy models of cosmic acceleration. In this review, we perform a detailed theoretical and phenomenological analysis of different modified gravity models and investigate their consistency. We also consider the cosmological implications of well motivated physical models of the early universe with a particular emphasis on inflation and topological defects. Astrophysical and cosmological tests over a wide range of scales, from the solar system to the observable horizon, severely restrict the allowed models of the Universe. Here, we review several observational probes -- including gravitational lensing, galaxy clusters, cosmic microwave background temperature and polarization, supernova and baryon acoustic oscillations measurements -- and their relevance in constraining our cosmological description of the Universe.

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.

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