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.

We will briefly discuss the necessary conditions for stability of polytropies in $f(\hat R)$ Palatini gravity and the differences with the General Relativity ones.

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 study static symmetric solutions in the context of a gravitational theory based on a action-dependent Lagrangian. Such theory has been designed as a setup to implement dissipative effects into the traditional principle of least action. Dissipation appears therefore from the first principles and has a purely geometric origin. An interesting feature of this theory is the existence of a coupling four-vector $\lambda_{\mu}$, which in an expanding background is related to cosmological bulk viscosity. General Relativity is recovered with a vanishing $\lambda_{\mu}$. We analyse the existence of equilibrium solutions of static configurations aiming to describe astrophysical objects. We find out that the existence of static spherically symmetric configurations occurs only in the particular scenario with vanishing $\lambda_t$, $\lambda_r$ and $\lambda_{\phi}$ components i.e, $\lambda_{\mu}=\{0,0,\lambda_{\theta},0\}$. Thus, the component $\lambda_{\theta}$ is the unique available parameter of the theory in the astrophysical context. This result severely constrains the existence of this sort of gravitational theories. We proceed then verifying the impact of $\lambda_{\theta}$ on the stability and the mass-radius configurations for a reasonable equation of state for the cold dense matter inside compact stars. We further investigate the relativistic spherical collapse in order to track the structure of geometrical singularities appearing in the theory.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Two exact solutions for $n=0$ and $n=1$ of the Palatini-modified Lane-Emden equation are found. We have employed these solutions to describe a Palatini-Newtonian neutron star and compared the result with the pure Newtonian counterpart. It turned out that for the negative parameter of the Starobinsky model the star is heavier and larger.

We have derived a modified Lane-Emden equation for the Starobinsky model in Palatini gravity which is numerically solvable. Comparing the results to the ones provided by General Relativity we observe a significant difference depending on the theory parameter for the $M-R$ relations.

We study static, spherically symmetric equilibrium configurations in extended theories of gravity (ETG) following the notation introduced by Capozziello et {\it al}. We calculate the differential equations for the stellar structure in such theories in a very generic form i.e., the Tolman-Oppenheimer-Volkoff generalization for any ETG is introduced. Stability analysis is also investigated with special focus on the particular example of scalar-tensor gravity.

We investigate further (cf. arXiv:1512.01199, JCAP01 (2016) 040) Starobinsky cosmological model $R+\gamma R^2$ in the Palatini formalism with Chaplygin gas and baryonic matter as a source. For this aim we use dynamical system theory. The dynamics is reduced to the 2D sewn dynamical system of a Newtonian type (a piecewise-smooth dynamical system). We classify all evolutional paths in the model as well as trajectories in the phase space. We demonstrate the presence of a degenerate freeze singularity (glued freeze type singularities) for the positive $\gamma$. In this case it is a generic feature of early evolution of the universe. We point out that a degenerate type III of singularity can be considered as an endogenous model of inflation between the matter dominating epoch and the dark energy phase. We also investigate cosmological models with negative $\gamma$. It is demonstrated that $\gamma$ equal zero is a bifurcation parameter and dynamics qualitatively changes in comparison to positive $\gamma$. Instead of the big bang the sudden bounce singularity of a finite scale factor appears and there is a generic class of bouncing solutions sewn along the line $a=a_{\text{sing}}$. And we argue that the presence of sudden singularities in an evolutional scenario of the Universe falsifies the negative $\gamma$ in the Palatini cosmology. Only very small values of $\Omega_{\gamma}$ parameter are admissible if we requires that agreements physics with the $\Lambda$CDM model. From the statistical analysis of astronomical observations, we deduce that the case of negative values of $\Omega_\gamma$ can be rejected even if it may fit better to the data.

We present a simple generalisation of the $\Lambda$CDM model which on the one hand reaches very good agreement with the present day experimental data and provides an internal inflationary mechanism on the other hand. It is based on Palatini modified gravity with quadratic Starobinsky term and generalized Chaplygin gas as a matter source providing, besides a current accelerated expansion, the epoch of endogenous inflation driven by type III freeze singularity. It follows from our statistical analysis that astronomical data favors negative value of the parameter coupling quadratic term into Einstein-Hilbert Lagrangian and as a consequence the bounce instead of initial Big-Bang singularity is preferred.

Lie symmetries are discussed for the Wheeler-De Witt equation in Bianchi Class A cosmologies. In particular, we consider General Relativity, minimally coupled scalar field gravity and Hybrid Gravity as paradigmatic examples of the approach. Several invariant solutions are determined and classified according to the form of the scalar field potential. The approach gives rise to a suitable method to select classical solutions and it is based on the first principle of the existence of symmetries.

We test a free {\it ad hoc} parametrization of the Tolman-Oppenheimer-Volkoff (TOV) equation. We do not have in mind any specific extended theory of gravity (ETG) but each new parameter introduced has a physical interpretation. Our aim is fully pedagogical rather than a proposal for a new ETG. Given a realistic neutron star equation of state we map the contributions of each new parameter into a shift in trajectories of the mass-radius diagram. This exercise allows us to make the correspondence between each TOV sector with possible modifications of gravity and clarifies how neutron star observations are helpful for distinguishing theories.

Symmetries play a crucial role in physics and, in particular, the Noether symmetries are a useful tool both to select models motivated at a fundamental level, and to find exact solutions for specific Lagrangians. In this work, we consider the application of point symmetries in the recently proposed metric-Palatini Hybrid Gravity in order to select the $f({\cal R})$ functional form and to find analytical solutions for the field equations and for the related Wheeler-DeWitt (WDW) equation. We show that, in order to find out integrable $f({\cal R})$ models, conformal transformations in the Lagrangians are extremely useful. In this context, we explore two conformal transformations of the forms $d\tau=N(a) dt$ and $d\tau=N(\phi) dt$. For the former conformal transformation, we found two cases of $f({\cal R})$ functions where the field equations admit Noether symmetries. In the second case, the Lagrangian reduces to a Brans-Dicke-like theory with a general coupling function. For each case, it is possible to transform the field equations by using normal coordinates to simplify the dynamical system and to obtain exact solutions. Furthermore, we perform quantization and derive the WDW equation for the minisuperspace model. The Lie point symmetries for the WDW equation are determined and used to find invariant solutions.