Compact Objects in Alternative Gravities

We address neutron stars and black holes in alternative gravities, after recalling their basic properties in General Relativity. Among the plethora of interesting alternative gravities we here focus on an interesting set of scalar-tensor theories. We discuss the phenomenon of spontaneous scalarization, that is matter induced for neutron stars and curvature induced for black holes. Along with other relevant physical properties, we address the quasi-normal modes of these compact objects. In particular, we consider \textit{universal relations} of neutron stars to largely reduce the dependence on the equation of state, and we briefly address the shadow of black holes.


Introduction
Alternative gravities have moved into the focus since many years for a variety of reasons, as discussed in numerous contributions at this conference, where cosmological considerations are the focal point. The present contribution is meant as a supplement to these considerations, putting emphasis on compact objects in alternative gravities (see e.g. [1][2][3]). With the advent of gravitational wave observations by LIGO/VIRGO, in combination with multi-messenger astronomy [4][5][6] and first observations by the EHT consortium [7], the study of black holes and neutron stars has entered a new golden age of gravitational physics. The new data and the enhanced future data from upgraded and new observatories will allow to put (stronger) constraints on numerous alternative gravities, that will complement the available solar system data and the data from pulsar observations (see e.g. [8][9][10]).
In this contribution we will first address neutron stars and then black holes. In both cases we will set the stage by recalling properties of neutron stars and black holes in General Relativity (GR). Subsequently we will discuss corresponding properties for a selection of alternative gravities, where we focus on metric theories with an additional scalar degree of freedom. In the case of neutron stars we will put emphasis on universal relations, that allow to minimize the effects of the persisting uncertainty of the equation of state of nuclear matter, like the I-Love-Q relations and universal relations for quasi-normal modes (QNMs) (see e.g. [11,12]). The theories we will consider include scalar-tensor theories (STTs) [13][14][15] with spontaneous scalarization [16], R 2 gravity [17][18][19][20][21], or Chameleon screening [22]. For black holes we will address besides their basic properties their shadow and also their QNMs. We will discuss scalarization in Einstein-scalar-Gauss-Bonnet (EsGB) theories [23][24][25], where the scalar could be a dilaton or some other field, depending on the coupling function. In the latter case curvatureinduced spontaneously scalarized black holes can arise [26][27][28], where for rapid rotation also spin-induced scalarization can be present [29].

Neutron Stars in General Relativity
Neutron stars are highly compact objects as demonstrated in Fig. 1, where their massradius relation is shown for various equations of state. Despite strong efforts both from the observational and the theoretical side, there are still large uncertainties concerning the proper description of nuclear matter under the extreme conditions of strong gravity as found in neutron stars (see, e.g., the reviews [30][31][32]), although the observations of pulsars with masses of about 2 solar masses have already ruled out numerous equations of state [33][34][35]. The strong dependence of neutron star properties on the equation of state may seem a serious obstacle to employ neutron star data to extract constraints for many of the alternative gravities. However, here the universal relations can come to the rescue [11,12]. These represent relations between neutron star properties, that are to a large extent independent of the neutron star's internal structure. The most impressive relations are those between the moment of inertia I, the Love number, i.e., the tidal deformability λ, and the quadrupole moment Q [36].
Like the mass and the radius, also the moment of inertia, the tidal deformability and the quadrupole moment show a strong dependence on the equation of state (see e.g. [37]). However, when instead the dimensionless quantitiesĪ = I/M 3 ,λ = λ/M 5 andQ = QM/J 2 (with spin J) are considered, the dependence on the equation of state almost disappears [11,36]. The dimensionless quantities satisfy simple relations of the type the so-called I-Love-Q relations, that are satisfied to 1% or better [11,36]. Moreover, besides the lowest multipole moments also the higher multipole moments satisfy such relations, although the deviations increase somewhat [11,36]. (See for instance [38,39] for some analytical results on such universal relations.) Whereas the multipole moments of the neutron stars represent important characteristics of the stationary solutions and their tidal reaction to an external companion, further relevant characteristics of neutron stars are obtained from asteroseismology, i.e., the study of their spectra. The neutron star matter possesses a large variety of modes, that may be excited by external interactions. Of foremost importance are of course gravitational waves, that will be emitted in the collision of neutrons stars. In such processes the inspiral phase is followed by the merger and the ringdown, providing ample data for crucial analysis [6].
The spectrum of neutron stars can be investigated by calculating their QNMs (see e.g. [2,40,41]), that should be discernible by future gravitational wave detectors in the ringdown phase after collisions. Theoretically QNMs are obtained by employing perturbation theory. So far, most studies of QNMs of neutron stars have considered static spherically symmetric background solutions, making use of the simplifications arising due to symmetry. In particular, the perturbation equations then decouple for modes with odd (axial) and even (polar) parity, and they depend only on the angular number l (and not the azimuthal number m).
The perturbations concern the metric g µν = g and the nuclear matter with pressure p, density ρ and velocity u, and the index (0) denotes the background values.
The resulting system of perturbation equations then corresponds to an eigenvalue problem with a complex eigenvalue ω, whose real part ω R and imaginary part ω I describe the frequency and the decay rate of the gravitational wave, respectively. As noticed early on [40,42], also the QNMs of neutron stars exhibit universal relations. For better illustration, we show a set of l = 2 polar QNMs in Fig. 2 [43]. These are the lowest mode, called the fundamental f-mode, the first pressure p-mode and the first space-time w-mode. As seen in Fig. 2(a), where we exhibit the frequency ω R versus the compactness M/R, these modes typically exhibit a considerable dependence on the equation of state. Fig. 2(b) then demonstrates for the w-mode that by scaling the frequency with the radius a universal relation results. Universal relations arise also for the axial modes.

Neutron Stars beyond General Relativity
Let us now consider the properties of neutron stars in some alternative gravities and begin with a set of STTs, where in addition to the tensor a gravitational scalar degree of freedom is introduced [13][14][15]. When the STT action is formulated in the physical Jordan frame, the nuclear matter is minimally coupled to the Jordan frame metricg µν . Subsequent transformation to the Einstein frame then leads to a coupling of the scalar to the nuclear matter since the Jordan frame metric and the Einstein frame metric g µν are related bỹ where ϕ denotes the scalar field in the Einstein frame, and the coupling function A(ϕ) depends on the respective STT. Whereas the old Brans-Dicke theory is basically ruled out by observations in the solar system by now, other STTs are still interesting candidate theories, that may be constrained by observational data from compact objects. Therefore theoretical studies of the consequences of such theories for neutron star properties and gravitational wave emission are of importance, in particular, since the existence of universal relations may help largely to pinpoint effects of these alternative gravities, when the universal relations vary distinctly with the coupling constants of these theories.
We first address STTs that allow for matter induced spontaneous scalarization of neutron stars [16]. To understand this phenomenon, we inspect the coupled set of field equations in the Einstein frame with stress energy tensor T µν , trace T, and α = d(ln A(ϕ))/dϕ, and the scalar equation has been rewritten in terms of an effective mass m 2 eff . For spontaneous scalarization to occur, the coupling function A(ϕ) has to be chosen such that for ϕ = 0 also α = 0. In that case the GR equations are satisfied, when the scalar field vanishes, and the GR neutron star solutions remain solutions of the STT. However, in addition to the GR solutions also scalarized neutron star solutions can arise, when the scalar field equation develops a tachyonic instability. For instance, the choice [16] allows for a tachyonic instability, m 2 eff = −4πGβ 0 T < 0, when both the coupling constant β 0 and the trace T = 3p − ρ are negative, and form a sufficiently strong source for the scalar field.
Typically, neutron stars possess T < 0, therefore the coupling constant β 0 should be sufficiently negative to achieve spontaneous scalarization. We note, however, for later reference, that the pressure-dominated core of neutron stars with very high densities may lead to regions with T > 0, and therefore spontaneous scalarization is also possible for such stars if the coupling constant is sufficiently positive, β 0 > 0 [44,45]. Returning to β 0 < 0, we exhibit in Fig. 3 the effect of matter-induced scalarization for several equations of state [46]. Fig. 3(a) shows the mass-radius relation for several equations of state for two values of β 0 together with the GR relation. Moving along a family of stars the GR neutron stars develop an unstable mode when scalarization sets. In contrast, the scalarized neutron stars are stable and represent the preferred configurations. Scalarization is the stronger the more negative the coupling constant β 0 becomes. Fig. 3(b) illustrates a kind of universal relation that arises when the value of the . β cr denotes the onset of scalarization [46]. scalar field at the center ϕ(0) is considered versus the gravitational potential at the center (as contained in g tt (0)).
In such STTs the presence of the scalar degree of freedom would lead to an appreciable amount of dipolar radiation and thus to a more rapid orbital decay, while in GR only quadrupole (and higher multipole) radiation can be emitted. Therefore observational data from binary pulsars have been used to put constraints on the parameters of these theories (see e.g. [47,48]), leading recently even to the exclusion of the possibility of spontaneous scalarization for the above discussed class of massless STTs [49]. This conclusion can be avoided, however, when the scalar field has a larger mass, since the dipolar radiation is effectively shut off when the Compton wavelength of the scalar field is small in comparison to the orbital separation of a binary system [50]. The effects of scalar field masses and of scalar self-interaction have, for instance, been studied in [51][52][53], where also universal relations were obtained, that were distinct from those of GR.
We next consider neutron stars in a particular f (R) theory, Mathematically f (R) theories are equivalent to STTs, and can be transformed to the Einstein frame. The theory (10) then involves a coupling function A(ϕ) and a potential V(ϕ), where the scalar field mass m ϕ is related to the coupling constant a by m ϕ = 1/ √ 6a [20,54]. Then a coupling constant of a = 1 corresponds to a mass of m ϕ = 1.08 neV, while a = 10000 to m ϕ = 0.0108 neV. Such a mass range represents an interesting observational window [50] inviting theoretical studies [20,[54][55][56][57][58][59].
Concerning universal relations, the I-Q relations have been obtained for a = 10 4 not only for the slowly rotating case but also for rapidly rotating neutron stars, where the additional dependence on the dimensionless angular momentum J/M 2 leads to a universal relation surface [56] as in GR [60,61]. However, for such a value of the scalar mass, the GR relation and the R 2 relation differ distinctly [56]. Of interest are also the universal relation for QNMs. We exhibit examples for axial l = 2 modes in Fig. 4, where the scaled frequency ω R R and the scaled damping time M/τ are exhibited versus the compactness M/R for a = 10 and a = 10 5 and compared with the corresponding GR modes [57].  Axial modes do not involve perturbations of the matter and of the scalar field [57]. This is in contrast to polar modes, which exhibit a much richer spectrum [58,59]. Moreover, the additional scalar degree of freedom allows for ϕ modes, and the l = 0 normal modes of GR become propagating QNMs in R 2 gravity. The detection of gravitational monopole or dipole radiation (i.e., l < 2) would be a clear sign of the new gravitational degree of freedom.
n=0 φ-mode  It also compares with the massless case and with GR. In GR the F mode is a normal mode, that becomes a zero mode at the neutron star maximum mass, turning into an unstable mode beyond the maximum mass. In contrast, in R 2 gravity the additional degree of freedom leads for the stable neutron stars to QNMs, that possess a very small imaginary part of the eigenvalue ω, making them ultra long-lived [58].
The F mode is shown in Fig. 5(b) as a function of the scalar Compton wavelength for a given neutron star central density. Also shown are the excited H modes and the scalar ϕ mode. Interestingly, the scale of the frequency is set by (multiples of) the size of the star as long as the Compton wavelength is small, whereas for larger Compton wavelength the scale is set by the scalar field mass, as indicated by the curve 2πω R = m ϕ [58]. Universal relations for the polar modes can be obtained as well, as preliminary results including the l = 1 and l = 2 modes suggest.
As a final example of alternative gravities we here briefly address a STT with Chameleon screening [22]. We therefore consider neutron stars in a STT with Einstein frame coupling function A(ϕ) and potential V(ϕ) where α 0 is the coupling constant, and Λ sets the chameleon energy scale [62,63]. In Chameleon screening the scalar field obtains an effective mass whose value is given by the surrounding matter density whereT denotes the trace of the stress energy tensor of the matter in the Jordan frame.
Since the scalar field satisfies a Klein-Gordan equation for an effective potential the scalar field will assume the minimum value of the effective potential, thus depending strongly on the surrounding matter. The Chameleon field is then short-ranged in regions of high density and long-ranged on cosmological scales. For such a theory to remain viable, however, also compact objects like neutron stars should be well described and in accordance with observations. As shown in [62,63] such a Chameleon theory leads to stable neutron stars. This remains the case even when the trace of the stress energy tensor changes sign and the neutron stars develop a pressure dominated core with T > 0, and therefore a partial unscreening of the scalar field arises in the interior [62,63]. Predictions concerning the observability of scalar radiation from neutron stars in such a Chameleon theory must rely, however, on extrapolation, since reliable calculations could be performed only for Λ = 175 GeV,ρ ∞ = 6.5 × 10 +10 g/cm 3 (whereρ ∞ is the density outside the star), with the physical values of these quantities differing by many orders of magnitude, Λ = 2.4 meV,ρ ∞ = 1.67 × 10 −20 g/cm 3 . Extrapolation then predicts a peak in the mHz band in the sensitivity range of LISA, but far below its sensitivity [63].

Black Holes in General Relativity
The Kerr black hole of GR is considered of utmost astrophysical relevance [64,65]. It represents a fascinating family of asymptotically flat vacuum solutions that are regular on and outside their event horizon, and that are uniquely specified by their mass M and angular momentum J (see e.g. [66]). Kerr black holes carry no hair. Their higher mass moments M l and current moments S l are simply given in terms of the mass M 0 = M and the angular momentum in which odd mass moments and even current moments vanish. The Kerr paradigm supposes, that indeed all astrophysical black holes are well described by Kerr black holes [69]. Testing the no-hair hypothesis and the Kerr paradigm with current and future observations is clearly a central issue [70,71]. EHT observations of the shadow of the supermassive black hole at the center of the elliptic galaxy M87 are in good agreement with theoretical predictions obtained with full-fledged general relativistic magneto-hydrodynamics (GRMHD) calculations [7]. Of course, the shadow of a Kerr black hole itself and the dependence of its shape on the angular momentum has been known since long [72].
On the other hand, gravitational wave observations from merger events seen by LIGO/VIRGO employ standard GR templates based on Kerr black holes in their analyses [4,5]. Whereas the merger itself requires numerical relativity, the ringdown involves the QNMs of Kerr black holes, that are well-known (see e.g. [73,74]). An interesting properties of these Kerr QNMs is their isospectrality, i.e., the eigenvalues of polar and axial perturbations are identical.
Besides the Kerr black holes and their static limit, the Schwarzschild black holes (as well as the charged Kerr-Newman and Reissner-Nordström black holes), GR features also asymptotically flat black holes with hair (see, e.g., [75][76][77]). While the hairy black holes with Standard model fields [75,76] appear to be astrophysically irrelevant, this could be different for Kerr black holes with scalar hair [77] or Proca hair [78]. While the scalar and Proca fields needed would be new types of particles beyond the Standard Model, the corresponding rotating black holes represent interesting alternatives to the Kerr black holes and so far viable astrophysical candidates.

Black Holes beyond General Relativity
The STTs dicussed above, that possess matter-induced spontaneously scalarized neutron stars, do not support scalarized black holes. The no-hair theorems of GR forbid such solutions, and only Schwarzschild and Kerr black holes prevail. These remain also solutions of f (R) theories, as evident in the Einstein frame. The presence of the additional degree of freedom will, however, modify their QNM spectra and thus possibly implicate observable effects [79].
Here we would like to address another type of alternative gravity that includes a scalar degree of freedom and leads to black holes with scalar hair. In particular, we focus on a set of higher curvature theories, Einstein-scalar-Gauss-Bonnet (EsGB) theories [23][24][25], where the Einstein action is supplemented by the Gauss-Bonnet (GB) invariant coupled to a scalar field ϕ with a coupling function f (ϕ), and also a scalar kinetic term and, possibly, a scalar potential. The action leads to second order field equations without ghosts and therefore represents a certain type of Horndeski theory [80]. Depending on the coupling function f (ϕ), different types of scalarization arise, as seen from the field equations. Let us therefore take a closer look at the scalar field equation GR black hole solutions can not remain solutions, when d f dϕ = 0, since the scalar field would have a non-vanishing source term, and the resulting finite scalar field would contribute in the Einstein equations. In this case only hairy black holes arise. When, however, d f dϕ = 0 for ϕ = 0, then the GR black holes remain solutions of the theory, but additional scalarized solutions may arise. When the scalar field is considered to be a dilaton, as motivated by string theory, the coupling function f (ϕ) = α exp (−ϕ)/4 is of the first type, and only hairy black holes arise [81][82][83][84][85][86]. Fig. 6(a) shows the scaled event horizon radius r H / √ α versus the scaled mass M/ √ α of the family of static dilatonic black hole solutions, comparing them with the Schwarzschild solutions. Clearly, at fixed coupling constant α, large dilatonic black holes approach the Schwarzschild black holes. Smaller dilatonic black holes deviate significantly, however, and cease to exist beyond a minimal mass. The reason is of purely theoretical nature. The horizon expansion of the dilaton involves a square root, 1 − 6α 2 e 2φ H /r 4 H , whose radicand vanishes when the minimal mass is reached [81].
When these black holes are set into rotation, the dilatonic generalizations of the Kerr black holes are obtained [84]. Their domain of existence is exhibited in Fig. 6(b), where the scaled horizon area a H = A H /16πM 2 is shown versus the scaled angular momentum j = J/M 2 . The left boundary corresponds to the static dilatonic black holes, the upper boundary to the Kerr black holes, and the lower boundary to the critical solutions (where a radicand vanishes). Interestingly, these dilatonic black holes can slightly violate the Kerr bound, J/M 2 = 1. As expected from the static case, the area of the rotating dilatonic black holes is smaller than the area of Kerr black holes with the same angular momentum. The shadow is only slightly smaller than the shadow of comparable Kerr black holes [87]. The entropy of dilatonic black holes is, however, larger since the entropy receives an extra contribution from the GB term [88].
QNMs of static dilatonic GB black holes have been considered widely [89][90][91][92][93][94][95]. The scalar field allows for l = 0 and l = 1 radiation, and also for scalar-led higher l modes in addition to the grav-led modes, that reduce to the GR modes for vanishing GB coupling. As an example we exhibit in Fig. 7 the QNM eigenvalue for the polar l = 2 scalar-led and grav-led modes versus the strength of the GB coupling constant. The latter is scaled by the mass, ζ = α/M 2 . The frequency ω R and the decay rate ω I are scaled by the respective Schwarzschild values [92]. Close to the critical solution the deviations from Schwarzschild become more pronounced. The presence of the scalar field breaks isospectrality of the modes. Recent developments concern the perturbative inclusion of rotation for QNMs [96], and time evolution of the gravitational wave signal [97,98]. We now turn to the second type of coupling function f (ϕ), and require that the coupling function vanishes for vanishing scalar field. The resulting EsGB theory then retains the Schwarzschild and Kerr black holes as solutions. When the coupling is such, that the scalar field equation gets from the coupling function a term that is linear in ϕ, the resulting Klein-Gordan type equation contains an effective mass term. This equation is very reminiscent of Eq. (8) considered in the case of matter-induced spontaneously scalarized neutron stars. Indeed as realized in [26][27][28] this type of coupling function can give rise to curvature-induced spontaneously scalarized black holes, with the matter source term now replaced by the GB term, i.e., a higher curvature term.
We exhibit the domain of existence of static spontaneously scalarized black holes in Fig. 8 for two coupling functions, for exponential coupling f (ϕ) = λ 2 12 1 − e −6ϕ 2 in Fig. 8(a) and quadratic coupling f (ϕ) = λ 2 2 ϕ 2 in Fig. 8(b). The symmetry ϕ → −ϕ of the coupling functions leads to symmetric sets of EsGB solutions. When the curvature source term for the scalar field becomes sufficiently strong, a tachyonic instability sets in, and the fundamental branch (n = 0) of static scalarized black holes bifurcates from the Schwarzschild solution, whose GB term is given by R 2 GB = 48M 2 /r 6 . With increasing curvature, and thus smaller horizon radii and masses, a countable family of radially excited (n > 0) EsGB black holes branches off from the Schwarzschild solution. Inspection of the figures shows, that only the fundamental black hole branch of the exponential coupling features a concave set of solutions. All other sets of solutions are convex. This distinction has serious consequences for the stability of the solutions. As demonstrated in [26,99], only the n = 0 branch is stable against radial perturbations. All other branches possess radially unstable modes, that branch off from the unstable Schwarzschild modes. We recall, that the Schwarzschild solution itself has a zero mode at each bifurcation point, that subsequently turn into unstable modes.
When following the fundamental branch to smaller values of the mass for fixed coupling constant, a point r S1 is reached where the radial perturbation equations lose hyperbolicity, i.e., the employed formalism breaks down. To investigate, whether the solutions on the fundamental branch are mode stable between the bifurcation point from Schwarzschild r B and r S1 , the polar and axial modes should be free of instabilities. As shown in [100,101] this is indeed the case, except that the axial perturbation equations lose hyperbolicity slightly earlier at r S2 . The lowest polar and axial l = 2 grav-led modes are shown in Fig. 9 and compared to the Schwarzschild modes. For the scalarized black holes isospectrality is clearly broken. We note, that stability can also be achieved for a quadratic coupling function, when a potential for the scalar is included [102].
We finally address the rotating spontaneously scalarized EsGB solutions and their properties [29,[103][104][105][106][107][108]. To this end we inspect the GB term of the Kerr solutions (a = J/M), since in the case of rotation the Kerr black holes should be considered in the source term of the scalar equation. Obviously, the static spontaneously scalarized EsGB black holes can be set into rotation, and thus there are rotating EsGB black holes for positive coupling constant [103,104]. However, as already noted in the case of neutron stars, there is a second type of scalarized black holes, that arise for negative coupling constant [29,[105][106][107][108]. In this case, the GB invariant should be negative in some regions of space with high curvature, which can be achieved for sufficiently high spins. In fact, J/M 2 > 0.5 is needed for this spin-induced spontaneous scalarization to occur [29,105]. As expected from these considerations the domain of existence of rotating EsGB black holes with positive coupling decreases, when the spin increases [103,104]. This effect is most pronounced for the exponential coupling, where the domain of existence of static black holes is much larger than for the case of the quadratic coupling, shown in Fig. 10(a). Moreover, these rotating solutions arise largely from stable static solutions. Thus they are expected to retain stability. Calculation of their entropy in fact also points to their stability, since their entropy is always larger than the entropy of comparable Kerr black holes [103]. This is in contrast to the case of the quadratic coupling where the entropy is smaller [104]. Of course, in this case already the static solutions are unstable. Since in the exponential case, the horizon area of the scalarized black holes can be considerably smaller than the horizon area of comparable Kerr black holes, this opens the possibility to test their shadow by observations and thus put bounds on the coupling constant [103].
The domain of existence of spin-induced EsGB black holes has been determined in [107] for the exponential coupling and [108] for the quadratic coupling. Interestingly, in this case both coupling functions lead to similar results. In particular, in both cases the horizon area of spin-induced black holes is smaller than the horizon area of comparable Kerr black holes, whereas their entropy is larger. This indicates that the rapidly spinning EsGB black holes with spin-induced spontaneous scalarization might be stable even in the case of the quadratic coupling. However, a stability analysis of the rotating EsGB black holes is still missing.
We end this discussion by mentioning that recent studies have applied techniques of numerical relativity to EsGB theories with spontaneous scalarization of black holes. These studies address, for instance, dynamical scalarization and descalarization in black hole mergers or the dynamical formation of scalarized black hole in the collapse of burned out stars [109][110][111][112][113][114].

Conclusions
New observations on the ground and in space are giving us unprecedented insight on the Universe. This concerns both the evolution of the Universe as well as the celestial bodies inside it and their evolution and dynamics. Many of the associated theoretical studies have focused in recent years on alternative gravities to overcome the limitations and difficulties encountered. Indeed, alternative gravities have open up many promising directions in the ongoing quest (see e.g. [1][2][3]).
Here we have addressed neutron stars and black holes in a set of alternative gravities, that involve an additional scalar degree of freedom. The presence of such a new degree of freedom then allows for a variety of interesting phenomena. Foremost the spectrum of compact objects is considerably widened, since beside the gravitational quadrupole (and higher mode) radiation, also monopole and dipole radiation is allowed, while the l ≥ 2 gravitational-led modes are supplemented by scalar-led modes. Viable alternative gravities must respect the associated bounds set by current and future observations. Also, depending on the theory, neutron stars can, for instance, undergo matter-induced spontaneous scalarization, that is recurrent in black holes in the form of curvature-induced spontaneous scalarization. The resulting compact objects may then exhibit genuine properties distinguishing them from their GR brethren. Neutron stars may then possess deviating universal relations, while black holes may exhibit distinct shadows, etc., to be compared with observations, instrumental to obtain bounds.
Of considerable future relevance should be the communication between the groups working on alternative gravities as applied to cosmological issues and to compact objects. In particular, cosmologically attractive theories should be scrutinized with respect to their predictions for compact objects. Here first steps have been undertaken, and the future will certainly see more of these promising lines of research, also encouraged by the organizers of ALTECOSMOFUN'21 by the schedule of the meeting.
Author Contributions: All authors contributed substantially.
Funding: This research was funded in part by DFG RTG 1620 Models of Gravity and FCT project PTDC/FIS-AST/3041/2020.