Primordial black holes PBHs are usually assumed to be described by the Schwarzschild or Kerr metrics which however feature unwelcome singularities. We study the possibility that PBHs are nonsingular objects considering three phenomenological regular tr timeradialsymmetric spacetimes including the wellknown Bardeen and Hayward ones featuring either de Sitter or Minkowski cores. We characterize the evaporation of these PBHs and constrain their abundance from gammaray observations. For all three metrics we find that constraints on ftextpbh the fraction of dark matter DM in the form of PBHs weaken with respect to the Schwarzschild limits because of modifications to the PBH temperature and greybody factors. This moves the lower edge of the asteroid mass window down by potentially an order of magnitude or more leading to a much larger region of parameter space where PBHs can make up all the DM. A companion paper is devoted to nontextittrsymmetric metrics including loop quantum gravityinspired ones. Our work provides a proofofprinciple for the interface between the DM and singularity problems being a promising arena with a rich phenomenology.
['General Relativity and Quantum Cosmology', 'Astrophysics - Cosmology and Nongalactic Astrophysics', 'High Energy Physics - Phenomenology', 'High Energy Physics - Theory']
Primordial regular black holes as all the dark matter. I. Timeradialsymmetric metrics
{'Primordial regular black holes as all the dark matter. I. Time-radial-symmetric metrics': "Marco Calz'a, 1, 2, ∗ Davide Pedrotti, 1, 2, † and Sunny Vagnozzi 1, 2, ‡  \n1 Department of Physics, University of Trento, Via Sommarive 14, 38123 Povo (TN), Italy § 2 Trento Institute for Fundamental Physics and Applications (TIFPA)-INFN, Via Sommarive 14, 38123 Povo (TN), Italy (Dated: December 9, 2024)  \nPrimordial black holes (PBHs) are usually assumed to be described by the Schwarzschild or Kerr metrics, which however feature unwelcome singularities. We study the possibility that PBHs are non-singular objects, considering three phenomenological, regular tr (time-radial)-symmetric spacetimes (including the well-known Bardeen and Hayward ones), featuring either de Sitter or Minkowski cores. We characterize the evaporation of these PBHs and constrain their abundance from γ -ray observations. For all three metrics we find that constraints on f pbh , the fraction of dark matter (DM) in the form of PBHs, weaken with respect to the Schwarzschild limits, because of modifications to the PBH temperature and greybody factors. This moves the lower edge of the asteroid mass window down by potentially an order of magnitude or more, leading to a much larger region of parameter space where PBHs can make up all the DM. A companion paper is devoted to nontr -symmetric metrics, including loop quantum gravity-inspired ones. Our work provides a proof-ofprinciple for the interface between the DM and singularity problems being a promising arena with a rich phenomenology.", 'I. INTRODUCTION': "The Standard Model of particle physics (SM) and General Relativity (GR) have proven to be extremely successful at describing a huge range of terrestrial, astrophysical, and cosmological observations. However, their successes are limited by a number of shortcomings, potentially (especially in the case of SM) pointing towards the need for new physics which may better describe the matter and gravity sectors. On the more observational/phenomenological side, the SM lacks a candidate for the dark matter (DM) which accounts for ≃ 25% of the energy budget of the Universe [1, 2]. On the more theoretical side, continuous gravitational collapse in GR leads to the pathological appearance of curvature singularities [3, 4]. The nature of DM and the singularity problem are arguably two among the most important open questions in theoretical physics.  \nThe solution to the DM problem could reside in the physics of some of the most peculiar objects in the Universe: black holes (BHs). It has long been realized that primordial BHs (PBHs), hypothetical relics from the primordial Universe formed from the collapse of large density perturbations upon horizon re-entry, are indeed excellent DM candidates [5-104] (see e.g. Refs. [105-114] for reviews): in fact, PBHs are the only viable DM candidate which does not invoke new particles surviving to the present day. Once believed to merely be objects of mathematical speculation, observational effects associated to BHs are now routinely detected [115], turning  \nthese objects into extraordinary probes of fundamental physics [116-166]. As a result, the possibility of PBHs accounting for the entire DM budget is severely constrained by a wide range of considerations and constraints. The only (not entirely undebated) remaining open window of parameter space where PBHs could make up all the DM is the so-called ' asteroid mass window ', roughly for PBH masses 10 17 g ≲ M pbh ≲ 10 23 g [167-182]: lighter PBHs would have evaporated fast enough to either have disappeared by now or overproduced γ -rays in the MeV range, whereas heavier PBHs would have been detected through the microlensing of background stars.  \nAlmost all works on PBHs assume that these are Schwarzschild or Kerr BHs [105-114]. All constraints and considerations on DM potentially being in the form of PBHs are therefore subject to this underlying assumption. The existence of the asteroid mass window, and the extension thereof, is of course no exception. The assumption in question is not at all unreasonable at least from the phenomenological point of view, given that there are at present no signs of tension between astrophysical observations and the Kerr-Newman family of metrics, and more generally the no-hair theorem. Nevertheless, from the theoretical point of view such an assumption might stir some unease, given the appearance of singularities in the Schwarzschild and Kerr metrics. The above considerations naturally lead to the following question: ' what if PBHs are non-singular '? The present work is a pilot study whose goal is to explore this question, which naturally merges the DM and singularity problems.  \nWe entertain the possibility that PBHs are 'regular', i.e. free of curvature singularities [183-187], and therefore that DM may be in the form of primordial regular BHs (PRBHs). For concreteness, we consider so-called tr (time-radial)-symmetric metrics, for which the product of the coefficients for the dt 2 and dr 2 components of the line element in four-dimensional Boyer-Lindquist coordi-  \nnates is equal to -1, and the function which multiplies the angular part of the line element is r 2 , i.e. r is the areal radius. More specifically, we focus on the following three regular, static spherically symmetric space-times, all of which are characterized by an additional regularizing parameter ℓ and recover the Schwarzschild space-time in the ℓ → 0 limit: Bardeen BHs [188], Hayward BHs [189], and Culetu-Ghosh-Simpson-Visser BHs [190-193]. These three space-times present a rich phenomenology, including either de Sitter or Minkowski cores. We focus our attention on observational constraints from PRBH evaporation, which set the lower edge of the asteroid mass window, discussing in detail how the evaporation process is modified with respect to that for Schwarzschild PBHs. We show that, as a result, the phenomenology of PRBHs can be quite different from that of Schwarzschild PBHs, with a larger range of masses where PRBHs could make up the entire DM component, opening up the asteroid mass window by up to an extra decade in mass. Keeping in mind that the metrics in question are phenomenological in nature, our results demonstrate that a common solution to the DM and singularity problems in the form of PRBHs is one which is worth taking seriously, and warrants further investigation, and more generally the interface of these two problems provides a promising arena. We stress that our work should not be intended as a comprehensive analysis of PRBHs, but rather as a pilot study, pointing towards a direction which has thus far received little attention and indicating promising directions for further work.  \nThe rest of this paper is then organized as follows. In Sec. II we briefly introduce the regular space-times studied in the rest of the work. Various aspects of our methodology are discussed in Sec. III, with Sec. III A devoted to the calculation of the so-called greybody factors, Sec. III B to the computation of photon spectra resulting from Hawking evaporation, and Sec. III C to the comparison against observations. The resulting limits on the fraction of DM which may be in the form of primordial regular BHs are then critically discussed in Sec. IV. Finally, in Sec. V we draw concluding remarks. A number of more technical aspects concerning the greybody factors computation are discussed in Appendix A, whereas the time evolution of the primordial regular BHs we consider is studied in Appendix B. Unless otherwise specified, we adopt units where G = c = ℏ = 1. In closing, we note that a related study is being presented in a companion paper [194]: this focuses on nontr -symmetric metrics, which also include loop quantum gravity-inspired metrics, but at the same time complicate the study of the evaporation process. We recommend that the interested reader go through the present work prior to consulting our companion paper [194].", 'II. REGULAR BLACK HOLES': "It is well known, thanks to the Penrose-Hawking singularity theorems, that continuous gravitational collapse in GR sourced by matter contents satisfying reasonable energy conditions leads to the appearance of singularities [3, 4]. These are regions of space-time where curvature invariants, i.e. sets of independent scalars constructed from the Riemann tensor and the metric, diverge (with the archetypal example being the central singularity in the Kerr-Newman family of metrics). These singularities are arguably unsatisfactory as they lead to a potential breakdown in predictivity (see Refs. [195-197] for a different viewpoint). For this reason, they are oftentimes regarded as a manifestation of our lack of knowledge of (new) physics in the high-energy/high-curvature regime. A widespread belief is that quantum gravity effects on these scales would ultimately cure the singularity problem (and potentially lead to observable effects), although this is more of a hope supported only by a few first-principles studies [198-212].  \nEven in the absence of a widely agreed upon theory of quantum gravity, one can still hope to make progress in understanding and taming singularities, while also potentially gaining intuition about the possible features of such a theory, through a more phenomenological approach. Under the assumption that a metric description maintains its validity, one can introduce metrics which are free of singularities in the entire space-time, and describe so-called regular BHs (RBHs) [213]. It is often (albeit not necessarily always) the case that RBH metrics are controlled by an extra parameter, which in what follows we shall refer to as regularizing parameter (and denote by ℓ ), typically recovering the Schwarzschild metric (for non-rotating RBHs) in the limit ℓ → 0. Several RBH metrics have been studied over the past decades, see e.g. Refs. [214-264] for an inevitably incomplete selection of examples, as well as Refs. [183-187] for recent reviews on the subject. 1 While most of these metrics have been introduced on purely phenomenological grounds it is known that possible sources for several RBH metrics lie in theories of non-linear electrodynamics [268-272].  \nAs alluded to earlier, our interest in this work is to explore the possibility that primordial RBHs may play the role of DM. As a proof-of-principle in this sense we will establish constraints on f pbh , the fraction of DM in the form of PRBHs, focusing on the asteroid mass window, whose extent we will show can be further extended. To the best of our knowledge, we are aware of seven works in this little explored direction [273-279]. Ref. [274] studied primordial BHs with de Sitter interiors as DM candidates, but considering the case where the DM is actually constituted of remnants from the evaporation process. Ref. [275] studied the thermodynamics  \nof primordial regular BHs, focusing however on the case where they do not evaporate, and therefore did not study constraints on f pbh . Ref. [276] studied the evaporation of a loop quantum gravity-inspired BH, and weakened constraints on f pbh were reported in a later proceeding [277] (which however does not appear to be widely known)., Ref. [278] studied signatures of primordial BHs with magnetic charge, which could be (as is often but necessarily the case) regular. Finally, Refs. [273, 279] entertained the case where the evaporation times of primordial regular BHs are significantly longer than the standard Schwarzschild case, and did not therefore explicitly compute the evaporation constraints we instead study here. The aim of this pilot study and our companion paper [194] is instead to provide a more comprehensive investigation of primordial regular BHs, considering a more diverse set of metrics and investigating constraints on f pbh in detail.  \nIn our work, we shall consider three different nonrotating RBH metrics, as discussed in more detail in the following subsections. The static, spherically symmetric space-times we investigate are a subset of the Petrov type-D class of metrics. In four-dimensional BoyerLindquist coordinates, their line elements take the following general form:  \nds 2 = -f ( r ) dt 2 + g ( r ) -1 dr 2 + h ( r ) d Ω 2 , (1)  \nwhere d Ω 2 = dθ 2 + sin 2 ( θ ) dϕ 2 is the metric on the 2sphere. We also require our space-times to be asymptotically flat, which amounts to the following requirements:  \nf ( r ) r →∞ ---→ 1 , g ( r ) r →∞ ---→ 1 , h ( r ) r →∞ ---→ r 2 . (2)  \nIn addition, as stated earlier, we require our space-times to be tr -symmetric (the nontr -symmetric case is covered in a companion paper [194]), which imposes the following additional conditions:  \nf ( r ) = g ( r ) , h ( r ) = r 2 , (3)  \nimplying that the coordinate r is effectively the areal radius. With the conditions given by Eqs. (2,3) imposed upon Eq. (1), our most general line element therefore takes the following form:  \nds 2 = -f ( r ) dt 2 + f ( r ) -1 dr 2 + r 2 [ dθ 2 +sin 2 ( θ ) dϕ 2 ] . (4)  \nIn what follows, we refer to the function f ( r ) as being the 'metric function'. The three different RBH solutions we consider, which we will discuss very shortly in Sections II A- II C, are characterized by different functional forms of f ( r ).  \nAn important parameter characterizing the behaviour of BHs is their temperature T . This is particularly crucial when evaluating evaporation constraints on PBHs, given that the temperature controls the strength of the emitted radiation, which in turn can be directly constrained by various observations. We treat the temperature of the RBHs as being the usual Gibbons-Hawking one, i.e.  \nthe one evaluated by Wick rotating the metric in the standard way and imposing regularity in the Euclidean period [280]. The cyclic imaginary time → temperature identification is legitimate if one can formally identify the Euclidean action e -S with the Boltzmann factor e -βH in the partition function, as usually done in finite temperature quantum field theory: in turn, this can be done if one is assuming the standard Boltzmann-Gibbs distribution, but may not be the consistent if other entropies are assumed (see e.g. the recent discussion in Ref. [281]). Since, as we will reiterate later, the RBHs we study are introduced on phenomenological grounds and we remain agnostic as to their theoretical origin (which may in principle be rooted within alternative entropic frameworks), in what follows we assume the Boltzmann-Gibbs distribution, so that the temperature of the RBHs in question is the standard Gibbons-Hawking one, and is given by the following:  \nT = κ 2 π = f ' ( r ) 4 π | r H , (5)  \nwhere the prime indicates a derivative with respect to r , and κ is the BH surface gravity, given by the following:  \nκ = f ' ( r ) 2 | r H . (6)  \nIn Eqs. (5,6), r H is the horizon radius, which is the solution to the following equation:  \ng -1 ( r H ) = f ( r H ) = 0 , (7)  \nwith the first equality following from the choice of focusing on tr -symmetric space-times. In the case of Schwarzschild BHs, where the metric function is f ( r ) = 1 -2 M/r , one recovers the well-known expressions r H = 2 M and T Sch = 1 / 8 πM . However, in more general spacetimes the horizon radius in Eq. (7) is not guaranteed to have a closed form expression, and the same therefore holds for the temperature in Eq. (5). In Fig. 1 we show the evolution of the temperatures (normalized to the temperature of Schwarzschild BHs, T Sch = 1 / 8 πM ) of the three regular space-times we will introduce shortly, as a function of the regularizing parameter ℓ (itself normalized to the event horizon radius r H ). As the Figure clearly shows, for all three space-times the temperature is a monotonically decreasing function of ℓ . 2 We also see  \nFIG. 1. Evolution of the temperatures (normalized to the temperature of Schwarzschild black holes, T Sch = 1 / 8 πM ) as a function of the regularizing parameter ℓ (normalized to the event horizon radius r H ) for the three regular black holes studied in this work: the Bardeen (blue solid curve, Sec. II A), Hayward (red dashed curve, Sec. II B), and Culetu-GhoshSimpson-Visser (green dotted curve, Sec. II C) regular spacetimes. In all cases the temperature is a monotonically decreasing function of ℓ . Note that the range of allowed values of ℓ/r H is different for the three regular black holes.  \n<!-- image -->  \nthat the temperatures of these BHs approach zero in the extremal limit. As has been observed earlier, in this limit the usual constraints on f pbh vanish, since these BHs do not evaporate [275].  \nA final caveat is in order before discussing the RBH metrics we consider. The latter are all regular in the sense of having finite curvature invariants R ≡ g µν R µν , R µν R µν , and K ≡ R µνρσ R µνρσ . However, a more stringent criterion for regularity is that of geodesic completeness, which does not necessarily imply finiteness of curvature invariants and viceversa. We recall that a spacetime is geodesically complete if any geodesic thereon can be extended to arbitrary values of their affine parameter. In other words, all geodesics extend (or can be extended) for all times. A number of 'popular' RBHs have indeed been shown to have finite curvature invariants but to be geodesically incomplete [293]. This includes the well-known Hayward RBH, which is among the ones we shall consider here. Nevertheless, given the significant interest in this metric, the fact that it is widely taken as prototype for RBHs, and our phenomenological goal of going beyond Schwarzschild PBHs, we will take this space-time into consideration, while cautioning the reader about the above issues, and therefore that the Hayward metric (but more generally all the metrics we will study) should be considered nothing more than phenomenological toy models at this stage. Note, in addition, that the stability of RBHs featuring inner horizons is currently a matter of debate in the literature [294-298].", 'A. Bardeen black hole': 'The Bardeen BH is easily one of the best known RBHs, and one of the first ones to ever have been proposed [188]. It is characterized by the following metric function: 3  \nf B ( r ) = 1 -2 Mr 2 ( r 2 + ℓ 2 ) 3 / 2 , (8)  \nwhere, in terms of the BH mass M , the regularizing parameter satisfies ℓ ≤ √ 16 / 27 M ∼ 0 . 77 M in order for the metric to describe a BH and not a horizonless object. Note that the Schwarzschild metric function is recovered in the ℓ → 0 limit. A perhaps physically more motivated choice is to express quantities in units of the horizon radius r H , defined as the largest root of the equation f ( r H ) = 0, in which case the regularizing parameter is subject to the constraint ℓ ≲ 0 . 70 r H (which is the point at which T = 0 in Fig. 1). In order to obtain this limit we have computed the solution to f ( r H ) = 0 fixing M = 1, in order to extrapolate r H ( ℓ ), before analyzing for which real values of the parameter n the equation ℓ = nr H ( ℓ ) admits solutions. We note that the peculiar factor of √ 16 / 27 in the extremality condition appears when one demands that the cubic equation determining the horizon(s) location(s) of the Bardeen RBH, f B ( r ) = 0, admits one real non-zero root: with some algebraic manipulation, one sees that this condition amounts to the requirement that √ 27 ℓ 8 M 4 -16 ℓ 6 M 6 is real, from which the requirement ℓ ≤ √ 16 / 27 M follows.  \nIt is worth noting that the Bardeen BH possesses a de Sitter (dS) core which replaces the central singularity of the Schwarzschild BH. This is evident by noting that, in the limit r → 0, the metric function goes as f B ( r ) ∝ r 2 , exactly as expected for an asymptotically dS space-time. Although originally introduced on phenomenological grounds, it is now known that the Bardeen RBH can emerge from a magnetic monopole source [299], potentially within the context of a specific non-linear electrodynamics theory [300]. Another possible origin for the Bardeen RBH are quantum corrections to the uncertainty principle [301]. Irrespective of its origin, and consistently with the approach pursued for the other spacetimes, we consider this solution as a model-agnostic phenomenological toy model.', 'B. Hayward black hole': 'Another widely known regular space-time is the Hayward RBH [189], characterized by the following metric  \nfunction:  \nf H ( r ) = 1 -2 Mr 2 r 3 +2 Mℓ 2 . (9)  \nIf expressed in terms of BH mass M , the regularizing parameter for the Hayward BH is subject to the same limit as that of the Bardeen BH, i.e. ℓ ≤ √ 16 / 27 , M . On the other hand, if expressed in terms of the more physically motivated horizon radius, the limit is instead ℓ ≲ 0 . 57 r H (again this is the point where T = 0 in Fig. 1). We note that the Schwarzschild metric function is recovered in the ℓ → 0 limit. The factor of √ 16 / 27 in the extremality condition originates from considerations similar to those we made for the Bardeen RBH in Sec. II A: with some algebraic manipulation, one sees that this condition amounts to the requirement that √ 27 ℓ 4 M 2 -16 ℓ 2 M 4 is real, from which the requirement ℓ ≤ √ 16 / 27 M follows.  \nJust as the Bardeen RBH possesses a dS core, so does the Hayward RBH. Indeed, introducing a dS core characterized by a (positive) cosmological constant Λ = 3 /ℓ 2 in order to prevent the central singularity was precisely the original justification for the Hayward BH which, just like its Bardeen counterpart, was proposed on purely phenomenological grounds. Nevertheless, potential theoretical origins for the Hayward BH have been investigated, and range from corrections to the equation of state of matter at high density [302, 303], finite density and finite curvature proposals [304-306], theories of non-linear electrodynamics [307, 308], and more generally as the result of corrections due to quantum gravity [309, 310]. Just as with the Bardeen RBH, we shall here treat the Hayward RBH as a model-agnostic phenomenological toy model for a singularity-free space-time.', 'C. Culetu-Ghosh-Simpson-Visser black hole': "The regular space-times considered so far featured dS cores, which in itself is a very common feature of several RBH metrics. Nevertheless, another interesting phenomenological possibility consists in considering 'hollow' RBHs wherein the central singularity is replaced by an asymptotically Minkowski core, where the associated energy density and pressure asymptote to zero. This is quite unlike the case of the dS core where the energy density asymptotes to a finite value associated to a positive cosmological constant, and the pressure asymptotes to an equal but opposite value. Possible theoretical/mathematical motivations for considering RBHs with Minkowski cores include the fact that the vanishing energy density can significantly simplify the physics in the deep core, whereas the otherwise messy solutions to polynomial equations (which often cannot be written down in closed form) can be traded for arguably more elegant special functions, resulting in the space-time being more tractable. Our physical motivation in considering this class of BHs is instead to broaden the range of physical  \nproperties and phenomenological implications of PRBHs, going beyond the dS core RBHs studied thus far.  \nWith this in mind, we consider a RBH featuring a Minkowski core, independently studied by Culetu [190, 191], Ghosh [192], as well as Simpson and Visser [193]. Although such a RBH does not have any particular name associated to it in the literature, here we refer to it as CGSV BH (from the initials of the four authors above). The space-time is characterized by the following metric function:  \nf CGSV ( r ) = 1 -2 M r exp ( -ℓ r ) . (10)  \nThe horizon radius r H , for which a closed form expression is not available in the Bardeen and Hayward cases, here is given by:  \nr H = -ℓ W ( -ℓ 2 M ) , (11)  \nwhere W denotes the Lambert function. Considering the principal branch W 0 , a real and positive horizon radius is present for:  \nW 0 ( -ℓ 2 M ) ≤ 0 = ⇒ 0 ≤ ℓ < 2 M e , (12)  \nor, alternatively, 0 ≤ ℓ < r H . While the CGSV BH was original introduced purely on phenomenological/mathematical grounds, it was shown in Refs. [311, 312] that such a space-time can emerge within the context of GR coupled to a specific non-linear electrodynamics source. In this case, denoting by g the non-linear electrodynamics coupling constant/charge, the regularizing parameter ℓ is given by ℓ = g 2 / 2 M , with M the BH mass. Nevertheless, as with all the other RBHs considered, here we shall treat the CGSV RBH as a toy model for a regular space-time possessing a Minkowski core.", 'A. Greybody factors': "A set of parameters playing a key role in describing the Hawking radiation spectra emitted from evaporating BHs are the so-called greybody factors (GBFs). These are functions of energy/frequency and angular momentum which govern the deviation of the emitted spectrum from that of a blackbody [313-315]. Although the emitted Hawking radiation at the horizon takes the blackbody form, the potential barrier due to space-time geometry will attenuate the radiation, so that an observer at spatial infinity will measure a spectrum which differs from that of a blackbody by a frequency-dependent function Γ( ω ). GBFs can be characterized by setting up a classical scattering problem around the BH potential barrier, with boundary conditions allowing for incoming wave packets  \nfrom infinity or equivalently, due to the symmetries of the scattering problem, originating from the horizon. The scattering problem is governed by the so-called Teukolsky equation, which is a partial differential equation describing the propagation of perturbations of given spin in the BH background [316].  \nFor the static, spherically and tr -symmetric metrics given by Eq. (4) which we consider, the Teukolsky equation in spherical coordinates is separable. A key role in computing the GBFs is played by the radial Teukolsky equation, which we now report in full generality for the class of metrics in question. Using the Newman-Penrose (NP) formalism, the Teukolsky equation governing the evolution of (massless) perturbations of different spin can be condensed into a single master equation [316]:  \n[ -r 2 f ∂ 2 t + s ( r 2 f ' f -2 r ) ∂ t ] Υ s + [ ( s +1)( r 2 f ' +2 rf ) ∂ r ] Υ s + [ 1 sin θ ∂ θ (sin θ∂ θ ) + 2 is cot θ sin θ ∂ ϕ + 1 sin 2 θ ∂ 2 ϕ -s -s 2 cot 2 θ ] Υ s + [ sr 2 f '' +4 srf ' +2 sf ] Υ s = 0 . (13)  \nHere, Υ s represents a general perturbation of spin s , defined by the NP scalars relative to the respective perturbation. To not make the notation too heavy, we drop the l and m indices labelling the field mode, so Υ s is understood to really mean Υ lm s . We note that Eq. (13) is separable if one makes the following wave ansatz:  \nΥ s = ∑ l,m e -iωt e imϕ S l s ( θ ) R s ( r ) , (14)  \nwhere ω is the perturbation frequency, l is the angular node number, and m is the azimuthal node number.  \nThe functions S l s ( θ ) contribute to defining the socalled spin-weighted spherical harmonics S s l,m ( θ, ϕ ) = ∑ S l s ( θ ) e imϕ , satisfying the following equation [317-320]:  \n( 1 sin θ ∂ θ (sin θ ∂ θ ) + csc 2 θ ∂ 2 ϕ + 2 is cot θ sin θ ∂ ϕ + s -s 2 cot 2 θ + λ s l ) S s l,m = 0 , (15)  \nwhere λ s l ≡ l ( l +1) -s ( s +1) is the separation constant. For the spin 0 case, these functions reduces to the usual spherical harmonics S 0 l,m = Y l,m .  \nAnalogously to the Schwarzschild and Kerr BH cases [316], the decoupled radial Teukolsky equation takes the following form [321, 322]:  \n1 ∆ s ( ∆ s +1 R ' s ) ' + ( ω 2 r 2 f +2 iωsr -isωr 2 f ' f + s (∆ '' -2) -λ s l ) R s = 0 , (16)  \nwhere ∆( r ) ≡ r 2 f ( r ) and ' ≡ ∂ r . We set in purely ingoing boundary conditions, so the asymptotic solutions of Eq. (16) are given by:  \nR s ∼ R in s e -iωr ⋆ r + R out s e iωr ⋆ r 2 s +1 ( r →∞ ) R s ∼ R hor s ∆ -s e -iωr ⋆ ( r → r H ) , (17)  \nwhere r ⋆ is the tortoise coordinate defined by:  \ndr ⋆ dr = 1 f ( r ) . (18)  \nWe note that r ⋆ → r for large values of r , given that the metrics we consider are asymptotically flat.  \nIn general, numerical integration methods are required to compute GBFs for general space-times, and this holds for our tr -symmetric RBHs as well. In our work, we make use of the so-called shooting method (see Appendix A for further details), which has already been successfully applied to these types of calculations in several earlier works [323-330].  \nTo begin with, we rewrite Eq. (16) in terms of the rescaled coordinate x , given by the following:  \nx ≡ r -r H r H , (19)  \nwhere r H is the largest real root of f ( r ) = 0. With this substitution Eq. (16) is rewritten as follows:  \nx 2 ( x +1) 3 f R s +( s +1) x ( x +1) ( 2( x +1) f +( x +1) 2 ˙ f ) ˙ R s + V ( ω, x ) R s = 0 , (20)  \nwhere ˙ ≡ ∂ x , and V ( ω, x ) is given by:  \nV ( ω, x ) = ( ω 2 r 2 H ( x +1) 2 f +2 is ( x +1) ω -isr H ( x +1) 2 ˙ f f ω + s ( 2 f +4( x +1) ˙ f +( x +1) 2 f -2 ) -l ( l +1) + s ( s +1) ) x ( x +1) . (21)  \nIn order to further simplify the problem, we work in units of horizon radius and therefore set r H = 1, so that r = x + 1. In these units, the metric functions of the three RBHs under consideration are given by the following:  \nf B ( x ) = 1 -(1 + ℓ 2 ) 3 / 2 ( x +1) 2 ( ℓ 2 +( x +1) 2 ) 3 2 , f H ( x ) = 1 -( x +1) 2 (1 -ℓ 2 ) ( ( x +1) 3 -ℓ 2 ℓ 2 -1 ) , f CGSV ( x ) = 1 -e ℓ -ℓ x +1 x +1 , (22)  \nfor the Bardeen, Hayward, and CGSV space-times respectively.  \nSetting purely ingoing boundary in proximity of the horizon, the solutions to Eq. (20) can be expressed in the form of a Taylor expansion as follows [323, 324, 331-336]:  \nR s ( x ) = x -s -iω τ ∞ ∑ n =0 a n x n , (23)  \nwhere iω/τ is a function of the field spin and the regularizing parameter, and also depends on the space-time in question. We refer the reader to Appendix A for further details. The a n coefficients can be determined by substituting Eq. (23) in Eq. (20) and iteratively solving the resulting algebraic equations. The near-horizon solution is then used to set the boundary conditions and numerically integrate the radial equation up to large distances, where the general form of the solution is the following:  \nR ( x ) r →∞ ---→ R in s e -iωx x + R out s e iωx x 2 s +1 . (24)  \nThe GBFs can then be computed from the s R lm in ( ω ) coefficient. More specifically, the normalization of the scattering problem is set by requiring a 0 = 1. With this normalization, the GBFs then read:  \nΓ s lm ( ω ) = δ s | s R lm in ( ω ) | -2 , (25)  \nwhere the coefficient δ s is given by:  \nδ s = τ ie iπs (2 ω ) 2 s -1 Γ ( 1 -s -2 iω τ ) Γ ( s -2 iω τ ) , (26)  \nwhere Γ is the Γ function. Using the method discussed above, we compute the GBFs for perturbations of different spin on the backgrounds of the three regular metrics discussed earlier, for different values of the regularizing parameter ℓ . In the specific case s = 1, we have checked that calculating the GBFs up to l = 4 is sufficient for our purposes. The GBFs we calculate are then used to characterize the Hawking evaporation spectra.  \nIn Fig. 2 we show the Γ s =1 l =1 GBFs for the Schwarzschild BH and the three regular space-times we study. For illustrative purposes we have focused on the Γ s =1 l =1 GBFs, since we are interested in photons ( s = 1) and the dominant emission mode is the l = 1 one. We also note that, since we are considering spherically symmetric space-times, the (2 l +1) different m modes are degenerate. We have fixed the regularizing parameter to ℓ = 0 . 3 r H for all three regular space-times, also for illustrative purposes. We see that, for all three regular space-times, the GBFs are slightly higher than their Schwarzschild counterparts (by ≲ 20% at most), and asymptote to the latter for both ω/M ≲ 0 . 3 and ω/M ≳ 0 . 7. While for definiteness we have focused on these specific values of s , l , and ℓ , we have explicitly checked that the features described above are present for other values of the parameters in question as well.  \nFIG. 2. Greybody factors Γ s =1 l =1 as a function of ω/M for Schwarzschild BHs (black curve), as well as the Bardeen (blue solid curve), Hayward (red dashed curve), and Culetu-GhoshSimpson-Visser (green dotted curve) regular space-times. For illustrative purposes we only plot Γ s =1 l =1 , since we are interested in photons ( s = 1) and the dominant emission mode is the l = 1 one. We have fixed the regularizing parameter to ℓ = 0 . 3 r H for all three regular space-times. We see that in all three cases the GBFs are consistently (slightly) higher than their Schwarzschild counterparts. The features shown in this plot do not change sensibly for higher values of l and other values of ℓ .  \n<!-- image -->", 'B. Evaporation spectra': 'We now discuss our computation of the photon spectra resulting from Hawking evaporation of the regular BHs discussed previously. In what follows, we only account for the primary photon spectrum. Nevertheless, we have checked that in the mass region of interest the impact of the secondary component of the spectrum, i.e. that resulting from the decay into photons of other unstable particles which are also produced during the evaporation process, is negligible.  \nThe Hawking radiation rate (number of particles emitted per unit time per unit energy) of a given particle species i with spin s , as a result of Hawking evaporation, is given by the following [337-341]: 4  \nd 2 N i dtdE i = 1 2 π ∑ l,m n i Γ s l,m ( ω ) e ω/T ± 1 , (27)  \nwhere n i is the number of degrees of freedom of the particle in question, ω = E i is the mode frequency (in natural units), Γ s l,m are the GBFs discussed previously, and we have implicitly set k B = 1. Note that  \nFIG. 3. Primary photon spectra resulting from the evaporation of Bardeen black holes of mass 10 16 g for different values of the regularizing parameter ℓ (normalized by the horizon radius r H ): ℓ/r H = 0 . 15 (red dotted curve), 0 . 3 (green dashed curve), and 0 . 45 (magenta dash-dotted curve). The blue solid curve corresponds to the case ℓ/r H = 0, which recovers the Schwarzschild black hole.  \n<!-- image -->  \nFIG. 4. As in Fig. 3, but for Hayward black holes, with identical values of the regularizing parameter ℓ/r H and identical color coding.  \n<!-- image -->  \nthe plus (minus) sign in the denominator is associated to fermions (bosons). Following the methodology discussed in Sec. III A, we calculate the GBFs within all the BH space-times in question for photons ( s = 1), up to l = 4 (note that the angular node number l should not be confused with the regularizing parameter ℓ ). We have checked that adding higher l modes does not appreciably improve the resulting spectra.  \nWe show examples of the resulting evaporation spectra in Figs. 3, 4, and 5. The spectra obviously depend on the mass of the evaporating PRBH, which we have set to M pbh = 10 16 g, as it sits roughly in the middle of  \nFIG. 5. As in Fig. 3, but for Culetu-Ghosh-Simpson-Visser black holes, with identical values of the regularizing parameter ℓ/r H and identical color coding.  \n<!-- image -->  \nthe mass range of interest. Nevertheless, we stress that the features we discuss below do not depend on the chosen mass. The resulting spectra all peak approximately between 5 MeV and 10 MeV.  \nFor the Bardeen, Hayward, and CGSV PRBHs we observe that an increase in the regularizing parameter ℓ leads to a decrease in the intensity of the spectra at all energies . This is na¨ıvely what one could expect by inspecting the temperature evolution shown in Fig. 1, since for all three these PRBHs the temperatures decrease relative to their Schwarzschild counterparts. However, as Eq. (27) shows, the temperature is not the only quantity at play, since the GBFs enter as well. As we have seen in Fig. 2, the GBFs for all three PRBHs are larger than their Schwarzschild counterparts, which would seem to counteract the effect of a lower temperature. A posteriori, however, the dominant effect turns out to be that of a lower temperature. We could actually have expected as much, since the resulting spectra given by Eq. (27) are linear in the GBFs, but depend exponentially on the temperature. Therefore, one might expect that a decrease in the temperature will have a much more dramatic effect than a comparable increase in the GBFs. This is precisely what we observe in Figs. 3, 4, and 5.  \nAside from the amplitude of the spectrum, for the Bardeen and Hayward PRBHs we note that the position of the peak in the spectrum is only mildly affected by the regularizing parameter, an increase in which pushes the peak towards slightly lower energies. On the other hand, for the CGSV BH an increase in the regularizing parameter pushes the peak towards slightly higher energies. We do not exclude that this different behaviour may be related to the type of core being considered, and we defer a more detailed investigation of this point to future work. At any rate, we expect that the behaviour of the spectra discussed above should lead to constraints on f pbh which', 'C. Evaporation constraints': 'The spectra calculated in Sec. III B are then used to set evaporation constraints on f pbh ( M ) ≡ Ω pbh / Ω dm , the fraction of DM in the form of PBHs, where Ω pbh and Ω dm are the PBH and DM density parameters respectively. Specifically, the computed spectra are used to obtain predictions for the flux of photons resulting from Hawking evaporation, which are then directly compared against measurements of the extragalactic photon background across a wide range of energies (see e.g. Ref. [342] for a recent review). Evaporation constraints are the dominant ones in the 10 13 g ≲ M pbh ≲ 10 18 g mass range: the lower limit of the range is set by the requirement that PBHs have not yet evaporated at the time of recombination, whereas the upper limit is defined by measurements of the diffuse extragalactic γ -ray background (EGRB) in the energy range 100 keV ≲ E γ ≲ 5 GeV, given that the intensity of the Hawking radiation flux is inversely proportional to the mass of the evaporating BH. In what follows, we will direct our attention exclusively to PBHs for which M pbh ≳ 10 15 g: these have yet to fully evaporate today and, having formed deep during the radiation domination era, are therefore excellent non-baryonic DM candidates. There is another important reason for focusing on this mass range. As shown in Appendix B, PBHs (either Schwarzschild or the three PRBHs we consider) within this range have lifetimes much longer than the age of the Universe, are far from having fully evaporated at the present time, and have only lost a negligible fraction of their mass from formation until now. When using the symbol M pbh , we are therefore safe in denoting the values of the PBH mass both at formation and now.  \nWe work under the commonly adopted assumption that PBHs are isotropically distributed on sufficiently large scales. Therefore, the flux resulting from their evaporation and reaching us today is given by the redshifted sum of the contributions from all evaporating PBHs in our Universe, and can be used to constrain the average extragalactic distribution of DM in the form of PBHs. Furthermore, we work within the (also commonly adopted) approximation of monochromatic mass distributions (which can be expected if the formation mechanism arises from an amplification of the power spectrum at a very specific scale), although the effect of extended mass distributions is the subject of active research [343357]. Finally, as discussed earlier, we only consider the primary photon contribution, as the secondary component resulting from the decay into photons of other unstable particles is verified a posteriori to be negligible given the mass range of interest. While all these are clearly approximations, albeit widely adopted ones, we are confident that they are appropriate given the aim of our work. Our main goal is in fact to examine how the limits on f pbh change when moving from the Schwarzschild  \nPBH framework to that of the regular metrics presented in Sec. II, altering the asteroid mass window. It is more than reasonable to expect that the shift in constraints relative to the Schwarzschild case, δf pbh , is only weakly affected by the above approximations. In other words we expect these approximations to have similar impacts on the constraints on f pbh relative to Schwarzschild BHs for the different RBHs discussed in Sec. II, therefore leading to negligible effects on the relative shift δf pbh , which is the quantity we are ultimately interested in. At any rate, the adopted approximations also allow for a more direct comparison to several previous works, and we therefore consider them appropriate for our pilot study, while stressing that their impact should definitely be explored in future follow-up works. Finally, note that we are tacitly assuming that PBHs cluster in the galactic halo in the same way as other forms of DM (unless they are extremely large, which is not the case for the mass range of interest).  \nIn what follows, we therefore assume that PBHs all have the same initial mass M pbh . Following Ref. [358] we approximate the number of emitted photons in the logarithmic energy bin ∆ E γ ≃ E γ as being given by ˙ N γ ( E γ ) ≃ E γ ( d ˙ N γ /dE γ ). The emission rate of photons from Hawking evaporation per volume at a cosmological time t is then given by [358]: 5  \ndn γ dt ( E γ , t ) ≃ n pbh ( t ) E γ d 2 N γ dtdE γ ( M pbh , E γ ) , (28)  \nwhere n pbh ( t ) is the PBH number density at time t . By integrating and taking into account the redshift scaling of the photon energy and density we end up with:  \nn γ 0 ( E γ 0 ) = n pbh ( t 0 ) E γ 0 ∫ t 0 t ⋆ dt (1 + z ) d 2 N γ dtdE γ ( M pbh , (1 + z ) E γ 0 ) = n pbh ( t 0 ) E γ 0 ∫ z ⋆ 0 dz H ( z ) d 2 N γ dtdE γ ( M pbh , (1 + z ) E γ 0 ) , (29)  \nwhere t 0 denotes the present time, t ⋆ and z ⋆ are respectively the cosmic time and redshift at recombination, and H ( z ) is the expansion rate. Finally, n γ 0 ( E γ 0 ) is the present number density of photons with energy E γ 0 . The resulting photon flux (more properly, the rate of photons per unit time per unit area per unit solid angle) is then given by:  \nI ( E γ 0 ) ≡ c 4 π n γ 0 ( E γ 0 ) . (30)  \nIt is this quantity which can then be directly compared against observations.  \nWe assume a spatially flat ΛCDM cosmological model in specifying the expansion rate entering into Eq. (29), with the same cosmological parameters as in Ref. [358]. This allows us to robustly cross-check our Schwarzschild constraints on f pbh against those reported in the seminal Ref. [358], although we stress that our constraints are very stable against reasonable changes in the values of the cosmological parameters. Once the cosmological model is fixed, all the relevant quantities in Eq. (29) are known except for the present-day PBH number density, n pbh ( t 0 ), which can be constrained from EGRB observations and is ultimately related to f pbh . More specifically, for any given value of the PBH mass M pbh , through Eqs. (27,29,30) we can compute the unnormalized photon flux I ( E γ 0 ) /n pbh ( t 0 ), and adjust the normalization n pbh ( t 0 ) by comparing against EGRB observations (as we will explain shortly). This procedures gives us an upper limit on n pbh ( t 0 ), which can be translated into an upper limit on f pbh as follows:  \nf pbh ( M pbh ) ≡ Ω pbh Ω dm = n pbh ( t 0 ) M pbh ρ crit , 0 Ω dm , (31)  \nwhere ρ crit , 0 = 3 H 2 0 / 8 πG is the present-day critical density, with H 0 the Hubble constant, and we recall that this procedure is done for various values of M pbh .  \nWe compare our theoretical predictions against various measurements of the EGRB. Specifically, we use observations of the EGRB from the HEAO-1 X-ray telescope in the 3-500 keV range [359], the COMPTEL imaging Compton γ -ray telescope in the 0 . 8-30 MeV range [360], and the EGRET γ -ray telescope [361]. A few comments are in order concerning the adopted datasets. While these are by now a couple of decades old, they basically still represent the state-of-the-art in the energy range of interest. One could entertain other observations, including local galactic measurements of the galactic γ -ray background [362], positron flux [363], 0 . 511 MeV annihilation radiation [364-366], and various other sources. While these galactic observations could lead to potentially stronger limits, they depend strongly on the form of the PBH mass function (assumed to be monochromatic in our study), as well as the clustering properties of these PBHs. On the other hand, our limits on f pbh are effectively testing the average extragalactic distribution of DM. Finally, other measurements of the EGRB are available, e.g. from Fermi-LAT [367], but these are mostly important for energy ranges larger than the ones of interest, and therefore for PBHs lighter than the ones we are considering. Therefore, we believe the choice of datasets (which is the one adopted in several works estimating evaporation limits on PBHs) is appropriate given the objective of our study.  \nTo set upper limits on n pbh ( t 0 ) - and therefore f pbh through Eq. (31) - we adopt the simple method first explained in the seminal Ref. [358], and later adopted in most of the works examining constraints on PBHs from  \nFIG. 6. Photon fluxes resulting from the evaporation of primordial Bardeen black holes with regularizing parameter ℓ = 0 . 3 r H , and masses of M pbh = 10 15 g (brown curve), 2 . 3 × 10 15 g (purple curve), 5 × 10 15 g (red curve), 10 16 g (green curve), 3 × 10 16 g (green curve), and 5 × 10 16 g (blue curve). The triangles indicate 1 σ upper limits on the extragalactic γ -ray background flux as measured by HEAO-1 (maroon), COMPTEL (turquoise), and EGRET (teal). For each value of M pbh , the PBH fraction f pbh has been set to its upper limit, determined as soon as the first datapoint is overshot (the latter is different for different values of M pbh ).  \n<!-- image -->  \nthe EGRB. Specifically, for each value of M pbh , and for given values of the regularizing parameter ℓ , the maximum allowed value of f pbh is determined by requiring that the predicted photon flux does not overshoot any of the ERGB datapoints by more than 1 σ . An example is shown in Fig. 6 for Bardeen PRBHs with regularizing parameter ℓ = 0 . 3 r H : for each of the mass values M pbh represented, the upper limit on f pbh is set as soon as the first datapoint is overshot. As is clear from the Figure, different PBH masses result in different datapoints being overshot. For each of the PRBHs considered, we use this procedure to determine upper limits on f pbh for fixed, representative values of ℓ ( ℓ/r H = 0 . 15, 0 . 3, and 0 . 45), comparing the results to the Schwarzschild case which is recovered when ℓ = 0. 6 We note that the exact origin of the EGRB is currently a matter of debate [368]:  \nFIG. 7. Upper limits on f pbh , the fraction of dark matter in the form of primordial regular Bardeen black holes, as a function of the black hole mass M pbh . The limits are derived for different values of the regularizing parameter ℓ (normalized by the horizon radius r H ), with the shaded regions excluded: ℓ/r H = 0 . 15 (red dotted curve), 0 . 3 (green dashed curve), and 0 . 45 (magenta dash-dotted curve). Note that the blue solid curve corresponds to the case ℓ/r H = 0, which recovers the Schwarzschild black hole, whereas the value of M pbh corresponding to the upper right edge of the f pbh constraints marks the lower edge of the asteroid mass window.  \n<!-- image -->  \nalthough it is commonly believed that distant astrophysical sources such as blazars give a major contribution to the EGRB, there is no complete consensus on the level of this contribution. In this light, our approach of simply requiring that the PRBH evaporation contribution to the EGRB does not exceed any observed datapoint is rather conservative (given that there could in principle be a PBH contribution to the EGRB, should it be conclusively determined that known astrophysical sources are unable to fully account for the latter).', 'IV. RESULTS': "For each of the PRBHs discussed in Sec. II, we now proceed to derive upper limits on f pbh as a function of the PRBH mass M pbh , for different values of ℓ , using the methodology presented in Sec. III C. The results are shown in Figs. 7, 8, and 9 for Bardeen, Hayward, and CGSV PRBHs respectively. For each case, we also plot the constraints on f pbh for the ℓ = 0 case (blue solid curve in all the Figures), which correspond to the standard Schwarzschild PBH scenario widely studied in the  \nFIG. 8. As in Fig. 7, but for primordial regular Hayward black holes, with identical values of the regularizing parameter ℓ/r H and identical color coding.  \n<!-- image -->  \nFIG. 9. As in Fig. 7, but for primordial regular Culetu-GhoshSimpson-Visser black holes, with identical values of the regularizing parameter ℓ/r H and identical color coding.  \n<!-- image -->  \nliterature. As a sanity check, we have verified that our ℓ = 0 constraints exactly recover those of the seminal Ref. [358]. It is worth noting that, for any given value of ℓ , the value of M pbh corresponding to the upper right edge of the f pbh constraints (i.e. the value of M pbh for which the limit reads f pbh < 1) marks the lower edge of the asteroid mass window.  \nFor all PRBHs, we saw earlier that the temperature and photon spectra decrease in intensity with increasing regularizing parameter ℓ (see the discussion in Sec. III B, and Figs. 3-5). As we could have expected, this behaviour leads to overall looser constraints on f pbh (for any given M pbh ) relative to the standard limits reported for Schwarzschild PBHs in the literature. In the case of near-extremal Hayward PRBHs ( ℓ = 0 . 45 r H ) this be-  \nour is somewhat enhanced compared to the nearextremal Bardeen and CGSV PRBHs, with the upper limits on f pbh approximately three orders of magnitude looser than the corresponding Schwarzschild ones: again, this is somewhat unsurprising when comparing Fig. 4 to Figs. 3 and 5. This could also have been expected from Fig. 1, noting that the temperature of Hayward BHs decreases more rapidly with increasing ℓ/r H relative to the Bardeen and especially CGSV ones. Although the temperature is not the only factor at play in determining the resulting evaporation spectra, given that the GBFs also play a key role as per Eq. (27), it is reassuring to see that the temperature behaviour observed in Fig. 1 is qualitatively reflected in the limits on f pbh we derive.  \nAs a result of the shifts discussed above, the lower edge of the asteroid mass window where PBHs could make up the entire DM component is modified for all three metrics considered. We recall that in the Schwarzschild case, the lower edge of the window lies at M pbh ≃ 10 17 g. For the PRBHs we consider, the looser constraints on f pbh result in the asteroid mass window further opening up by up to half a decade in mass or more. The maximum extension of the window is reached for the Hayward PRBH closer to extremality, in which case the lower edge decreases by about an order of magnitude to M pbh ≃ 10 16 g. Overall, we therefore observe that considering PRBHs in place of the standard Schwarzschild ones can relax the resulting constraints on f pbh , further opening up the asteroid mass window. The allowed region for the window lower edge spans over a decade in mass, at least for the PRBHs and range of ℓ considered. We note that, as ℓ moves towards the extremal limit, the temperatures of these PRBHs T approaches zero. This is indeed generally expected from thermodynamical arguments. In this case, which is the one studied in Ref. [275], the PRBHs do not evaporate, and therefore our evaporation constraints do not apply. See also Refs. [279] for a related study appearing after ours.  \nThree comments are in order before concluding. Firstly we note that, for a given PRBH space-time, the curves describing the f pbh ( M ) limits are approximately, but not exactly parallel to the Schwarzschild ones (blue solid curves in Figs. 7-9). The reason is simply that, as ℓ is increased, the datapoint shown in Fig. 6 which is first being overshot and therefore responsible for determining the f pbh limit can potentially change (in part due to the spectrum slightly changing shape).  \nNext, the constraints we have determined on f pbh at a fixed value of ℓ/r H implicitly assume that all PRBHs in the Universe carry the same value of 'hair' parameter ℓ . However, particularly given our agnostic stance with regards to the origin of these space-times, in principle the value of ℓ/r H can vary from PRBH to PRBH. To make an analogy, let us assume for a moment that ReissnerNordstrom BHs are astrophysically relevant. Then, since the electric charge Q is not tied to a universal parameter of the underlying Einstein-Maxwell Lagrangian, there is no reason to expect it to carry the same value across  \nall BHs. In the language of Ref. [282], the regularizing parameter for all three PRBHs considered is a 'specific hair' rather than an 'universal hair' (see Ref. [282] for various examples of BH solutions carrying universal hair), unless one were able to tie ℓ to some fundamental parameter of the underlying theory, which however is not the case in the phenomenological approach we are following. In principle one should therefore account for the (non-monochromatic) ℓ distribution for PRBHs across the Universe to determine constraints on f pbh . We see no obvious way of doing this, while noting that such a procedure would most likely result in upper limits on f pbh lying between the Schwarzschild and extremal cases: this observation suffices for our pilot study, and we defer a more complete investigation to future work.  \nOur final comment concerns the fact that evaporation limits on the PBH abundance are not the only ones at play. Indeed, as recently summarized in Ref. [110], there are essentially four classes of limits, each of which is relevant in a different mass range: evaporation, lensing, dynamical, and accretion constraints. Constraints from the accretion of background gas at early times are relevant in a completely different mass range (10 33 ≲ M pbh / g ≲ 10 40 - see Fig. 7 of Ref. [110] and Fig. 10 of Ref. [369]). Although these have been derived assuming Schwarzschild PBHs, moving to the PRBH picture we have considered will not shift the relevant mass range by the ≳ 18 orders of magnitude required for these constraints to compete with the evaporation ones, unless the physics of gas accreting around RBHs changes drastically with respect to the standard picture, which appears very unlikely. Dynamical constraints, most of which are associated to the destruction of different astronomical objects by the passage of nearby PBHs, are also relevant in a completely different mass range (10 34 ≲ M pbh / g ≲ 10 55 - see Fig. 7 of Ref. [110] and Fig. 10 of Ref. [369]), and considerations completely analogous to those we made for accretion constraints hold. 7  \nOf potentially more relevance to the present work are lensing constraints, which constrain the abundance of PBHs (and more generally massive compact halo objects) with masses M pbh ≳ 10 15 g. Indeed, it is lensing constraints which locate the upper edge of the asteroid mass window where PBHs can make up all the DM. Nevertheless, we expect that these constraints should not change when moving from the Schwarzschild PBH framework to the PRBHs considered in this work. Indeed, with all other quantities being fixed (mass of source, relative  \ndistances, and so on), lensing constraints only depend on the lens mass M , and are unaffected by the metric structure of the lens. Therefore, at fixed mass M , we can assume that the lensing limits on Schwarzschild PBHs hold for our PRBHs as well. Note that, as already pointed out in footnote 2, the parameter M appearing in the RBH metrics can be unambiguously identified with the RBH mass, just as with the parameter M in the Schwarzschild metric. We can therefore conclude that for the PRBHs we are considering it should only be the lower edge of the asteroid mass window which is altered with respect to the Schwarzschild case, but not the upper edge. In other words, space-times for which the lower edge moves towards lower masses (as in the Bardeen, Hayward, and CGSV PRBH cases) genuinely correspond to an enlarged asteroid mass window. Therefore, the window where Bardeen, Hayward, and CGSV PRBHs could account for all the DM is much larger compared to that of Schwarzschild PBHs.  \nOther potentially relevant constraints come from µ -distortions in the Cosmic Microwave Background, and gravitational waves (either a stochastic background due to a population of coalescing PBHs or produced via second-order tensor perturbations generated by the scalar perturbations producing the PBHs, or associated to resolved events). The latter are expected to be relevant in a much higher mass range (again, see Fig. 7 of Ref. [110] and Fig. 10 of Ref. [369]), whereas the former are somewhat dependent on the PBH formation scenario from highσ tails of density fluctuations, and in particular on the shape of the tail. At any rate, while the focus in the present pilot study has been solely on evaporation constraints from the ERGB, revisiting all these other important sources of constraints (including the ones we discussed earlier) is a worthwhile endeavour which we plan to explore in upcoming works.", 'V. CONCLUSIONS': "Over the past decade, primordial black holes have regained tremendous interest as viable dark matter candidates, with the so-called 'asteroid mass window' (10 17 g ≲ M pbh ≲ 10 23 g) where PBHs could potentially account for the entire DM currently still open. Nearly all works on PBHs assume that these are Schwarzschild or Kerr BHs. However, while phenomenologically perfectly valid, such an assumption may stir some unease on the theoretical side, due to the appearance of singularities in these metrics. In our work, we have conducted a pilot study aimed at addressing a question which naturally merges the DM and singularity problems, arguably two among the most important open problems in theoretical physics: ' What if PBHs are non-singular '? Our study of primordial regular BHs (PRBHs) has focused on three so-called tr -symmetric metrics (including the wellknown Bardeen and Hayward space-times), whereas our companion paper [194] considers nontr -symmetric met-  \n, including various metrics inspired from loop quantum gravity.  \nWe show that evaporation constraints on f pbh , the fraction of DM in the form of PRBHs, can be substantially loosened when moving away from the Schwarzschild picture, leading to the asteroid mass window further opening up. For the three PRBHs considered (the Bardeen, Hayward, and CGSV ones) the lower edge of the asteroid mass window is shifted by a decade in mass or more, leading to a larger region of parameter space where PRBHs could account for the entire DM component, which should be the target of the same probes proposed to test the standard window [171, 173, 175177, 179, 180, 377, 378]. The nature of the regular BH core (de Sitter or Minkowski) does not appear to play a significant role in this sense. Overall, we have shown that the phenomenology of primordial regular BHs can be particularly rich, making the associated simultaneous solution to the DM and singularity problems one worthy of further study. We note that the constraints we obtained would become even weaker had we moved closer to the extremal limit for the regularizing parameter ℓ : indeed, in this limit the temperatures of the RBHs we considered approaches zero, implying that there are no constraints from evaporation [275].  \nWe remark that the present work (alongside our companion paper [194]) should be intended as a pilot study, and there are a huge number of interesting follow-up directions. One interesting avenue for further work involves systematically revisiting other sources of constraints which have been studied in the Schwarzschild PBH case, including but not limited to lensing, accretion, and dynamical constraints: while we have argued that these should not alter our considerations on the asteroid mass window, a detailed study which would allow us to extend our constraints over a much larger region of M pbh -f pbh plane is nevertheless in order. In addition, the metrics we have considered are inherently phenomenological in nature, and it would therefore be worth extending our study to non-singular metrics which enjoy a strong theoretical motivation (our companion paper [194] goes partially in this direction), including potentially metrics which are coupled to the cosmological expansion [237, 242, 244, 261, 264]. Moreover, the formation mechanisms for these PRBHs is likely to be much more complex than the corresponding mechanisms for Schwarzschild PBHs: in fact, the space-times considered are not vacuum solutions of GR, which implies that Birkhoff's theorem does not (necessarily) hold, and correspondingly the endpoint of gravitational collapse may not be unique. At this stage it is impossible for us to make definitive statements on the matter without abandoning our agnostic viewpoint and assuming a specific underlying theoretical model, but the issue of PRBH formation (possibly within the context of inflationary models leading to an enhanced spectrum of curvature fluctuations over specific scales) is one which requires further study. Last but definitely not least, if PBHs are truly regular,  \none would hope to ascertain this via signatures complementary to those we have studied: gravitational wave observations, VLBI imaging, motion around BHs, or energy extraction are potentially interesting observables in this sense. For instance, the shadows of Bardeen, Hayward, and CGSV BHs are smaller than their Schwarzschild counterparts by up to 20% [282], whereas the motion of stars around these BHs (take for instance the S2 star orbiting around Sgr A ⋆ ) would be altered compared to the motion around a Schwarzschild BH. These probes are being actively studied within the community as a means of distinguishing Schwarzschild BHs from alternatives thereto, and are likely to provide a promising route towards testing the regular nature of PRBHs, in the event that these are detected in the future and demonstrated to make up the DM, or even just a fraction thereof. We plan to address these and other related points in followup work.", 'ACKNOWLEDGMENTS': "We acknowledge support from the Istituto Nazionale di Fisica Nucleare (INFN) through the Commissione Scientifica Nazionale 4 (CSN4) Iniziativa Specifica 'Quantum Fields in Gravity, Cosmology and Black Holes' (FLAG). M.C. and S.V. acknowledge support from the University of Trento and the Provincia Autonoma di Trento (PAT, Autonomous Province of Trento) through the UniTrento Internal Call for Research 2023 grant 'Searching for Dark Energy off the beaten track' (DARKTRACK, grant agreement no. E63C22000500003). This publication is based upon work from the COST Action CA21136 'Addressing observational tensions in cosmology with systematics and fundamental physics' (CosmoVerse), supported by COST (European Cooperation in Science and Technology).", 'Appendix A: Details on the computation of greybody factors': "Here we provide a few more details on the computation of GBFs. We recall that we expressed the solutions to the radial Teukolsky equation, Eq. (20), in the form of a Taylor expansion as given by Eq. (23). This is also known as a Frobenius series, being a by-product of a method for solving second-order differential equations named after Frobenius. The method applies to equations which take the following form  \nu '' + p ( x ) u ' + q ( x ) u = 0 , (A1)  \nin proximity of its singular points, namely those where p ( x ) and/or q ( x ) diverge. One can notice that Eq. (20) can be rewritten in the form of Eq. (A1), with one of its singular point being at x = 0, i.e. at the event horizon.  \nTo solve the radial Teukolsky equation we therefore proceed as follows:  \n- · We work in units of the event horizon and rewrite Eq. (20) in order to remove the denominators  \nA ( x ) R '' s + B ( x ) R ' s + C ( x ) R s = 0 , (A2)  \nwhere the functions A ( x ), B ( x ), and C ( x ) are given by the following:  \nA ( x ) = f 2 ( x +1) 2 , B ( x ) = ( s +1) f 2 ( x )(2( x +2) + ( x +1) 2 f ' /f ) , C ( x ) = +( x +1) 2 ω 2 +2 is ( x +1) ωf -is ( x +1) 2 ωf ' + sf ( ( x +1) 2 f '' +4( x +1) f ' +2 f -2 ) -l ( l +1) f + s ( s +1) f ,  \n- · The lowest power term around x = 0 of each coefficient can be written in the following form:  \nA ( x ) ∼ x 2 τ 2 , B ( x ) ∼ x ( s +1) τ 2 , C ( x ) ∼ ω 2 -iωsτ ,  \nwhere τ = τ ( ℓ ) depends on the choice of RBH.  \n- · We then build the following characteristic equation:  \nm ( m -1) τ 2 + m ( s +1) τ 2 + ω ( ω -isτ ) = 0 , (A3)  \nwhose solutions are the following:  \nm 1 = -s -iω τ , m 2 = iω τ (A4)  \n- · It is then possible to conclude that Eq. (20) admits solutions near the singular point x = 0 of the form given by Eq. (23).  \nExplicitly, for the three RBHs in question, τ is given by the following:  \nτ B = 1 -2 ℓ 2 ℓ 2 +1 , τ H = (1 -3 ℓ 2 ) , τ CGSV = 1 -ℓ .  \nWe notice that in the Schwarzschild limit ℓ → 0, all of the above reduce to τ = 1 as one could expect.", 'Appendix B: PRBH time evolution': "Throughout the paper, in using M to denote the masses of PRBHs, we never specified whether we were referring to the initial mass, the present mass of the PRBH undergoing evaporation, or another quantity. This question is relevant since the mass of an evaporating BH is of course a monotonically decreasing quantity. As we shall see, this question is intertwined with the question of what  \nis the lifetimes of these PRBHs, and whether they have already evaporated by now: the aim of this Appendix is to address these questions.  \nIn Sec. III B we discussed the spectrum of photons emitted due to Hawking evaporation. In general, however, BHs will emit not only photons, but the whole spectrum of particles of the underlying theory. In fact, no specific coupling is required for this to occur. Enforcing energy conservation leads to the conclusion that the emission comes at the expense of the BH mass. Following the steps of Refs. [339-341, 379, 380], and considering a field of spin s , the associated depletion function is defined as follows:  \nf s = M 2 ∫ ∞ 0 dω d 2 N dtdE ω = M 2 ∑ i,l,m 1 2 π ∫ ∞ 0 dω n i Γ s l,m ( ω ) e ω/T ± 1 ω = M 2 ∑ i,l (2 l +1) 2 π ∫ ∞ 0 dω n i Γ s l ( ω ) e ω/T ± 1 ω, (B1)  \nwhere in the last equality we have exploited the fact that for spherically symmetric space-times the (2 l +1) different m modes are degenerate, and n i is the number of degrees of freedom of the particle in question. It is worth noting that there is no reference to µ , the mass of the field being evaporated by the BH. Given that the emission of particles with masses above the Hawking temperature T is Boltzmann-suppressed, in practice working within the usually adopted approximation where the particles are considered massless for µ < T , and absent at lower temperatures, is sufficiently accurate for our purposes (so we do not need to compute the GBFs for massive fields).  \nAt a given particle temperature, summing over all the degrees of freedom accessible at a given BH temperature one can define the following function:  \nF = N 0 f 0 + N 1 / 2 f 1 / 2 + N 1 f 1 + N 3 / 2 f 3 / 2 + N 2 f 2 , (B2)  \nwhere N s count the degrees of freedom of spin s such that µ < T . This function is related to the BH mass loss rate as a function of time t through the following differential equation:  \ndM dt = -F M 2 , (B3)  \n- [1] A. Arbey and F. Mahmoudi, Prog. Part. Nucl. Phys. 119 , 103865 (2021), arXiv:2104.11488 [hep-ph].\n- [2] M. Cirelli, A. Strumia, and J. Zupan, (2024), arXiv:2406.01705 [hep-ph].\n- [3] R. Penrose, Phys. Rev. Lett. 14 , 57 (1965).\n- [4] S. W. Hawking and R. Penrose, Proc. Roy. Soc. Lond. A 314 , 529 (1970).  \nwhich can easily be integrated numerically with initial condition given by M ( t = 0) = M 0 . The lifetime of a BH is then determined by the time at which M goes to 0.  \nLet us begin by considering the case of evaporating Schwarzschild PBHs. Integrating Eq. (B3) backwards, we compute the initial mass of a Schwarzschild PBH whose lifetime is roughly the age of the Universe ( t U ∼ 13 . 8 Gyr), finding M 0 ∼ 5 × 10 14 g in agreement with Ref. [339]. On the other hand, we find the lifetime of a Schwarzschild PBH of initial mass M 0 = 10 15 g to be ∼ 10 t U . However, of greater relevance to us is the question of how much does its mass vary from M 0 during the first t U , i.e. from formation until today. We find that, during this time, the mass of the PBH is subject to a negligible (sub-% level) decrease. In this light, and given that we have focused on PBHs with masses M > 10 15 g, which are far from having fully evaporated today, for our purposes we can safely approximate the BHs as being quasi-static throughout the lifetime of the Universe, while consequently denoting by M the values of the PBH mass both at formation and today.  \nThe above is valid for Schwarzschild PBHs. However, as we have already seen earlier, the three PRBHs we considered are all colder compared to their Schwarzschild counterparts (at a given mass), and the slight increase in the GBFs is not sufficient to compensate for these decreased temperatures. In other words, the considerations we drew in Sec. III B lead us to conclude that these PRBHs have a longer lifetime compared to their Schwarzschild counterparts at a given mass, expectation that we have checked explicitly. For example, the lifetime of a Bardeen BH of ℓ = 0 . 45 r H and initial mass M 0 = 10 15 g is > 10 3 t U , and during the first t U the mass lost to evaporation is negligible. In short, the PRBHs we considered have lifetimes considerably longer than Schwarzschild PBHs of the same initial mass, and they lose a negligible fraction of mass throughout the age of the Universe: therefore, they are very far from having fully evaporated today, and the quasi-static approximation we have adopted is justified.  \n- [5] G. F. Chapline, Nature 253 , 251 (1975).\n- [6] P. Meszaros, Astron. Astrophys. 38 , 5 (1975).\n- [7] M. Y. Khlopov and A. G. Polnarev, Phys. Lett. B 97 , 383 (1980).\n- [8] M. Y. Khlopov, B. A. Malomed, I. B. Zeldovich, and Y. B. Zeldovich, Mon. Not. Roy. Astron. Soc. 215 , 575 (1985).\n- [9] P. Ivanov, P. Naselsky, and I. Novikov, Phys. Rev. D 50 , 7173 (1994).\n- [10] S. Choudhury and A. Mazumdar, Phys. Lett. B 733 , 270 (2014), arXiv:1307.5119 [astro-ph.CO].\n- [11] K. M. Belotsky, A. D. Dmitriev, E. A. Esipova, V. A. Gani, A. V. Grobov, M. Y. Khlopov, A. A. Kirillov, S. G. Rubin, and I. V. Svadkovsky, Mod. Phys. Lett. A 29 , 1440005 (2014), arXiv:1410.0203 [astro-ph.CO].\n- [12] S. Bird, I. Cholis, J. B. Mu˜noz, Y. Ali-Ha¨ımoud, M. Kamionkowski, E. D. Kovetz, A. Raccanelli, and A. G. Riess, Phys. Rev. Lett. 116 , 201301 (2016), arXiv:1603.00464 [astro-ph.CO].\n- [13] S. Clesse and J. Garc'ıa-Bellido, Phys. Dark Univ. 15 , 142 (2017), arXiv:1603.05234 [astro-ph.CO].\n- [14] V. Poulin, P. D. Serpico, F. Calore, S. Clesse, and K. Kohri, Phys. Rev. D 96 , 083524 (2017), arXiv:1707.04206 [astro-ph.CO].\n- [15] A. Raccanelli, F. Vidotto, and L. Verde, JCAP 08 , 003 (2018), arXiv:1708.02588 [astro-ph.CO].\n- [16] J. L. Bernal, N. Bellomo, A. Raccanelli, and L. Verde, JCAP 10 , 052 (2017), arXiv:1709.07465 [astro-ph.CO].\n- [17] S. Clesse and J. Garc'ıa-Bellido, Phys. Dark Univ. 22 , 137 (2018), arXiv:1711.10458 [astro-ph.CO].\n- [18] K. Kohri and T. Terada, Class. Quant. Grav. 35 , 235017 (2018), arXiv:1802.06785 [astro-ph.CO].\n- [19] L. Liu, Z.-K. Guo, and R.-G. Cai, Phys. Rev. D 99 , 063523 (2019), arXiv:1812.05376 [astro-ph.CO].\n- [20] L. Liu, Z.-K. Guo, and R.-G. Cai, Eur. Phys. J. C 79 , 717 (2019), arXiv:1901.07672 [astro-ph.CO].\n- [21] R. Murgia, G. Scelfo, M. Viel, and A. Raccanelli, Phys. Rev. Lett. 123 , 071102 (2019), arXiv:1903.10509 [astroph.CO].\n- [22] B. Carr, S. Clesse, J. Garc'ıa-Bellido, and F. Kuhnel, Phys. Dark Univ. 31 , 100755 (2021), arXiv:1906.08217 [astro-ph.CO].\n- [23] L. Liu, Z.-K. Guo, R.-G. Cai, and S. P. Kim, Phys. Rev. D 102 , 043508 (2020), arXiv:2001.02984 [astro-ph.CO].\n- [24] M. P. Hertzberg, E. D. Schiappacasse, and T. T. Yanagida, Phys. Rev. D 102 , 023013 (2020), arXiv:2001.07476 [astro-ph.CO].\n- [25] P. D. Serpico, V. Poulin, D. Inman, and K. Kohri, Phys. Rev. Res. 2 , 023204 (2020), arXiv:2002.10771 [astroph.CO].\n- [26] V. De Luca, G. Franciolini, P. Pani, and A. Riotto, JCAP 04 , 052 (2020), arXiv:2003.02778 [astro-ph.CO].\n- [27] V. De Luca, G. Franciolini, P. Pani, and A. Riotto, Phys. Rev. D 102 , 043505 (2020), arXiv:2003.12589 [astro-ph.CO].\n- [28] V. De Luca, G. Franciolini, P. Pani, and A. Riotto, JCAP 06 , 044 (2020), arXiv:2005.05641 [astro-ph.CO].\n- [29] B. Carr, F. Kuhnel, and L. Visinelli, Mon. Not. Roy. Astron. Soc. 501 , 2029 (2021), arXiv:2008.08077 [astroph.CO].\n- [30] S. Bhagwat, V. De Luca, G. Franciolini, P. Pani, and A. Riotto, JCAP 01 , 037 (2021), arXiv:2008.12320 [astro-ph.CO].\n- [31] V. De Luca, V. Desjacques, G. Franciolini, P. Pani, and A. Riotto, Phys. Rev. Lett. 126 , 051101 (2021), arXiv:2009.01728 [astro-ph.CO].\n- [32] K. W. K. Wong, G. Franciolini, V. De Luca, V. Baibhav, E. Berti, P. Pani, and A. Riotto, Phys. Rev. D 103 , 023026 (2021), arXiv:2011.01865 [gr-qc].\n- [33] B. Carr, F. Kuhnel, and L. Visinelli, Mon. Not. Roy. Astron. Soc. 506 , 3648 (2021), arXiv:2011.01930 [astro-  \n- ph.CO].\n- [34] G. Dom'enech, C. Lin, and M. Sasaki, JCAP 04 , 062 (2021), [Erratum: JCAP 11, E01 (2021)], arXiv:2012.08151 [gr-qc].\n- [35] V. De Luca, G. Franciolini, P. Pani, and A. Riotto, JCAP 05 , 003 (2021), arXiv:2102.03809 [astro-ph.CO].\n- [36] A. Arbey, J. Auffinger, P. Sandick, B. Shams Es Haghi, and K. Sinha, Phys. Rev. D 103 , 123549 (2021), arXiv:2104.04051 [astro-ph.CO].\n- [37] G. Franciolini, V. Baibhav, V. De Luca, K. K. Y. Ng, K. W. K. Wong, E. Berti, P. Pani, A. Riotto, and S. Vitale, Phys. Rev. D 105 , 083526 (2022), arXiv:2105.03349 [gr-qc].\n- [38] V. De Luca, G. Franciolini, P. Pani, and A. Riotto, JCAP 11 , 039 (2021), arXiv:2106.13769 [astro-ph.CO].\n- [39] A. Cheek, L. Heurtier, Y. F. Perez-Gonzalez, and J. Turner, Phys. Rev. D 105 , 015022 (2022), arXiv:2107.00013 [hep-ph].\n- [40] A. Cheek, L. Heurtier, Y. F. Perez-Gonzalez, and J. Turner, Phys. Rev. D 105 , 015023 (2022), arXiv:2107.00016 [hep-ph].\n- [41] S. Heydari and K. Karami, Eur. Phys. J. C 82 , 83 (2022), arXiv:2107.10550 [gr-qc].\n- [42] G. Dvali, F. Kuhnel, and M. Zantedeschi, Phys. Rev. D 104 , 123507 (2021), arXiv:2108.09471 [hep-ph].\n- [43] S. Heydari and K. Karami, JCAP 03 , 033 (2022), arXiv:2111.00494 [gr-qc].\n- [44] V. De Luca, G. Franciolini, A. Kehagias, P. Pani, and A. Riotto, Phys. Lett. B 832 , 137265 (2022), arXiv:2112.02534 [astro-ph.CO].\n- [45] L. Liu, X.-Y. Yang, Z.-K. Guo, and R.-G. Cai, JCAP 01 , 006 (2023), arXiv:2112.05473 [astro-ph.CO].\n- [46] A. K. Saha and R. Laha, Phys. Rev. D 105 , 103026 (2022), arXiv:2112.10794 [astro-ph.CO].\n- [47] N. Bhaumik, A. Ghoshal, and M. Lewicki, JHEP 07 , 130 (2022), arXiv:2205.06260 [astro-ph.CO].\n- [48] L. A. Anchordoqui, I. Antoniadis, and D. Lust, Phys. Rev. D 106 , 086001 (2022), arXiv:2206.07071 [hep-th].\n- [49] Y.-F. Cai, X.-H. Ma, M. Sasaki, D.-G. Wang, and Z. Zhou, JCAP 12 , 034 (2022), arXiv:2207.11910 [astroph.CO].\n- [50] M. Oguri, V. Takhistov, and K. Kohri, Phys. Lett. B 847 , 138276 (2023), arXiv:2208.05957 [astro-ph.CO].\n- [51] G. Franciolini, I. Musco, P. Pani, and A. Urbano, Phys. Rev. D 106 , 123526 (2022), arXiv:2209.05959 [astroph.CO].\n- [52] K. Mazde and L. Visinelli, JCAP 01 , 021 (2023), arXiv:2209.14307 [astro-ph.CO].\n- [53] R.-G. Cai, T. Chen, S.-J. Wang, and X.-Y. Yang, JCAP 03 , 043 (2023), arXiv:2210.02078 [astro-ph.CO].\n- [54] L. A. Anchordoqui, I. Antoniadis, and D. Lust, Phys. Lett. B 840 , 137844 (2023), arXiv:2210.02475 [hep-th].\n- [55] L. Liu, Z.-Q. You, Y. Wu, and Z.-C. Chen, Phys. Rev. D 107 , 063035 (2023), arXiv:2210.16094 [astro-ph.CO].\n- [56] C. Fu and S.-J. Wang, JCAP 06 , 012 (2023), arXiv:2211.03523 [astro-ph.CO].\n- [57] S. Choudhury, M. R. Gangopadhyay, and M. Sami, Eur. Phys. J. C 84 , 884 (2024), arXiv:2301.10000 [astroph.CO].\n- [58] T. Papanikolaou, Class. Quant. Grav. 40 , 134001 (2023), arXiv:2301.11439 [gr-qc].\n- [59] J. A. de Freitas Pacheco, E. Kiritsis, M. Lucca, and J. Silk, Phys. Rev. D 107 , 123525 (2023), arXiv:2301.13215 [astro-ph.CO].\n- [60] S. Choudhury, S. Panda, and M. Sami, Phys. Lett. B 845 , 138123 (2023), arXiv:2302.05655 [astro-ph.CO].\n- [61] S. Choudhury, S. Panda, and M. Sami, JCAP 11 , 066 (2023), arXiv:2303.06066 [astro-ph.CO].\n- [62] I. Musco, K. Jedamzik, and S. Young, Phys. Rev. D 109 , 083506 (2024), arXiv:2303.07980 [astro-ph.CO].\n- [63] G.-W. Yuan, L. Lei, Y.-Z. Wang, B. Wang, Y.-Y. Wang, C. Chen, Z.-Q. Shen, Y.-F. Cai, and Y.-Z. Fan, Sci. China Phys. Mech. Astron. 67 , 109512 (2024), arXiv:2303.09391 [astro-ph.CO].\n- [64] S. Choudhury, S. Panda, and M. Sami, JCAP 08 , 078 (2023), arXiv:2304.04065 [astro-ph.CO].\n- [65] A. Ghoshal, Y. Gouttenoire, L. Heurtier, and P. Simakachorn, JHEP 08 , 196 (2023), arXiv:2304.04793 [hep-ph].\n- [66] Y. Cai, M. Zhu, and Y.-S. Piao, Phys. Rev. Lett. 133 , 021001 (2024), arXiv:2305.10933 [gr-qc].\n- [67] S. Choudhury, A. Karde, S. Panda, and M. Sami, JCAP 01 , 012 (2024), arXiv:2306.12334 [astro-ph.CO].\n- [68] H.-L. Huang, Y. Cai, J.-Q. Jiang, J. Zhang, and Y.S. Piao, Res. Astron. Astrophys. 24 , 091001 (2024), arXiv:2306.17577 [gr-qc].\n- [69] S. Choudhury, A. Karde, S. Panda, and M. Sami, Nucl. Phys. B 1007 , 116678 (2024), arXiv:2308.09273 [astroph.CO].\n- [70] G. Bhattacharya, S. Choudhury, K. Dey, S. Ghosh, A. Karde, and N. S. Mishra, Phys. Dark Univ. 46 , 101602 (2024), arXiv:2309.00973 [astro-ph.CO].\n- [71] S. Heydari and K. Karami, JCAP 02 , 047 (2024), arXiv:2309.01239 [astro-ph.CO].\n- [72] S. Heydari and K. Karami, Eur. Phys. J. C 84 , 127 (2024), arXiv:2310.11030 [gr-qc].\n- [73] S. Choudhury, K. Dey, A. Karde, S. Panda, and M. Sami, Phys. Lett. B 856 , 138925 (2024), arXiv:2310.11034 [astro-ph.CO].\n- [74] S. Choudhury, K. Dey, and A. Karde, (2023), arXiv:2311.15065 [astro-ph.CO].\n- [75] A. Ghoshal and A. Strumia, JCAP 07 , 011 (2024), arXiv:2311.16236 [hep-ph].\n- [76] H.-L. Huang, J.-Q. Jiang, and Y.-S. Piao, Phys. Rev. D 109 , 063515 (2024), arXiv:2312.00338 [astro-ph.CO].\n- [77] H.-L. Huang and Y.-S. Piao, Phys. Rev. D 110 , 023501 (2024), arXiv:2312.11982 [astro-ph.CO].\n- [78] L. A. Anchordoqui, I. Antoniadis, and D. Lust, Phys. Rev. D 109 , 095008 (2024), arXiv:2401.09087 [hep-th].\n- [79] S. Choudhury, A. Karde, S. Panda, and M. Sami, JCAP 07 , 034 (2024), arXiv:2401.10925 [astro-ph.CO].\n- [80] V. Thoss, A. Burkert, and K. Kohri, Mon. Not. Roy. Astron. Soc. 532 , 451 (2024), arXiv:2402.17823 [astroph.CO].\n- [81] T. Papanikolaou, X.-C. He, X.-H. Ma, Y.-F. Cai, E. N. Saridakis, and M. Sasaki, Phys. Lett. B 857 , 138997 (2024), arXiv:2403.00660 [astro-ph.CO].\n- [82] S. Choudhury, (2024), 10.1142/S0218271824410074, arXiv:2403.07343 [astro-ph.CO].\n- [83] S. Choudhury, A. Karde, P. Padiyar, and M. Sami, (2024), arXiv:2403.13484 [astro-ph.CO].\n- [84] L. A. Anchordoqui, I. Antoniadis, and D. Lust, Phys. Rev. D 110 , 015004 (2024), arXiv:2403.19604 [hep-th].\n- [85] T. Papanikolaou, S. Banerjee, Y.-F. Cai, S. Capozziello, and E. N. Saridakis, JCAP 06 , 066 (2024), arXiv:2404.03779 [gr-qc].\n- [86] Z. Yin and L. Visinelli, JCAP 10 , 013 (2024), arXiv:2404.10340 [hep-ph].  \n- [87] S. Choudhury, A. Karde, S. Panda, and S. SenGupta, Eur. Phys. J. C 84 , 1149 (2024), arXiv:2405.06882 [astro-ph.CO].\n- [88] S. Heydari and K. Karami, Astrophys. J. 975 , 148 (2024), arXiv:2405.08563 [gr-qc].\n- [89] G. Dvali, J. S. Valbuena-Berm'udez, and M. Zantedeschi, Phys. Rev. D 110 , 056029 (2024), arXiv:2405.13117 [hep-th].\n- [90] A. Boccia, F. Iocco, and L. Visinelli, (2024), arXiv:2405.18493 [astro-ph.CO].\n- [91] H.-L. Huang, J.-Q. Jiang, and Y.-S. Piao, (2024), arXiv:2407.15781 [astro-ph.CO].\n- [92] S. Choudhury, S. Ganguly, S. Panda, S. SenGupta, and P. Tiwari, JCAP 09 , 013 (2024), arXiv:2407.18976 [astro-ph.CO].\n- [93] L. A. Anchordoqui, I. Antoniadis, D. Lust, and K. P. n. Castillo, Phys. Dark Univ. 46 , 101714 (2024), arXiv:2407.21031 [hep-th].\n- [94] C. Yang, S. Wang, M.-L. Zhao, and X. Zhang, JCAP 10 , 083 (2024), arXiv:2408.10897 [astro-ph.HE].\n- [95] A. K. Saha, A. Singh, P. Parashari, and R. Laha, (2024), arXiv:2409.10617 [astro-ph.CO].\n- [96] L. A. Anchordoqui, I. Antoniadis, D. Lust, and K. P. n. Castillo, Phys. Dark Univ. 46 , 101681 (2024), arXiv:2409.12904 [hep-ph].\n- [97] C. Chen, A. Ghoshal, G. Tasinato, and E. Tomberg, (2024), arXiv:2409.12950 [astro-ph.CO].\n- [98] D.-C. Dai and D. Stojkovic, Phys. Dark Univ. 46 , 101662 (2024), arXiv:2409.14321 [gr-qc].\n- [99] H.-L. Huang, Y.-T. Wang, and Y.-S. Piao, (2024), arXiv:2410.05891 [astro-ph.GA].\n- [100] M. Zantedeschi and L. Visinelli, (2024), arXiv:2410.07037 [astro-ph.HE].\n- [101] M. Chianese, A. Boccia, F. Iocco, G. Miele, and N. Saviano, (2024), arXiv:2410.07604 [astro-ph.HE].\n- [102] W. Barker, B. Gladwyn, and S. Zell, (2024), arXiv:2410.11948 [astro-ph.CO].\n- [103] D. Borah and N. Das, (2024), arXiv:2410.16403 [hepph].\n- [104] H.-L. Huang, J.-Q. Jiang, J. He, Y.-T. Wang, and Y.-S. Piao, (2024), arXiv:2410.20663 [astro-ph.GA].\n- [105] M. Y. Khlopov, Res. Astron. Astrophys. 10 , 495 (2010), arXiv:0801.0116 [astro-ph].\n- [106] B. Carr, F. Kuhnel, and M. Sandstad, Phys. Rev. D 94 , 083504 (2016), arXiv:1607.06077 [astro-ph.CO].\n- [107] A. M. Green and B. J. Kavanagh, J. Phys. G 48 , 043001 (2021), arXiv:2007.10722 [astro-ph.CO].\n- [108] B. Carr and F. Kuhnel, Ann. Rev. Nucl. Part. Sci. 70 , 355 (2020), arXiv:2006.02838 [astro-ph.CO].\n- [109] P. Villanueva-Domingo, O. Mena, and S. PalomaresRuiz, Front. Astron. Space Sci. 8 , 87 (2021), arXiv:2103.12087 [astro-ph.CO].\n- [110] B. Carr and F. Kuhnel, SciPost Phys. Lect. Notes 48 , 1 (2022), arXiv:2110.02821 [astro-ph.CO].\n- [111] S. Bird et al. , Phys. Dark Univ. 41 , 101231 (2023), arXiv:2203.08967 [hep-ph].\n- [112] B. Carr, S. Clesse, J. Garcia-Bellido, M. Hawkins, and F. Kuhnel, Phys. Rept. 1054 , 1 (2024), arXiv:2306.03903 [astro-ph.CO].\n- [113] A. Arbey, in 58th Rencontres de Moriond on Very High Energy Phenomena in the Universe (2024) arXiv:2405.08624 [gr-qc].\n- [114] S. Choudhury and M. Sami, (2024), arXiv:2407.17006 [gr-qc].\n- [115] C. Bambi, Annalen Phys. 530 , 1700430 (2018), arXiv:1711.10256 [gr-qc].\n- [116] P. Creminelli and F. Vernizzi, Phys. Rev. Lett. 119 , 251302 (2017), arXiv:1710.05877 [astro-ph.CO].\n- [117] J. Sakstein and B. Jain, Phys. Rev. Lett. 119 , 251303 (2017), arXiv:1710.05893 [astro-ph.CO].\n- [118] J. M. Ezquiaga and M. Zumalac'arregui, Phys. Rev. Lett. 119 , 251304 (2017), arXiv:1710.05901 [astroph.CO].\n- [119] S. Boran, S. Desai, E. O. Kahya, and R. P. Woodard, Phys. Rev. D 97 , 041501 (2018), arXiv:1710.06168 [astro-ph.HE].\n- [120] T. Baker, E. Bellini, P. G. Ferreira, M. Lagos, J. Noller, and I. Sawicki, Phys. Rev. Lett. 119 , 251301 (2017), arXiv:1710.06394 [astro-ph.CO].\n- [121] L. Amendola, M. Kunz, I. D. Saltas, and I. Sawicki, Phys. Rev. Lett. 120 , 131101 (2018), arXiv:1711.04825 [astro-ph.CO].\n- [122] L. Visinelli, N. Bolis, and S. Vagnozzi, Phys. Rev. D 97 , 064039 (2018), arXiv:1711.06628 [gr-qc].\n- [123] M. Crisostomi and K. Koyama, Phys. Rev. D 97 , 021301 (2018), arXiv:1711.06661 [astro-ph.CO].\n- [124] A. Dima and F. Vernizzi, Phys. Rev. D 97 , 101302 (2018), arXiv:1712.04731 [gr-qc].\n- [125] Y.-F. Cai, C. Li, E. N. Saridakis, and L. Xue, Phys. Rev. D 97 , 103513 (2018), arXiv:1801.05827 [gr-qc].\n- [126] A. Casalino, M. Rinaldi, L. Sebastiani, and S. Vagnozzi, Phys. Dark Univ. 22 , 108 (2018), arXiv:1803.02620 [grqc].\n- [127] L. Barack et al. , Class. Quant. Grav. 36 , 143001 (2019), arXiv:1806.05195 [gr-qc].\n- [128] B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett. 123 , 011102 (2019), arXiv:1811.00364 [gr-qc].\n- [129] A. Casalino, M. Rinaldi, L. Sebastiani, and S. Vagnozzi, Class. Quant. Grav. 36 , 017001 (2019), arXiv:1811.06830 [gr-qc].\n- [130] A. Held, R. Gold, and A. Eichhorn, JCAP 06 , 029 (2019), arXiv:1904.07133 [gr-qc].\n- [131] C. Bambi, K. Freese, S. Vagnozzi, and L. Visinelli, Phys. Rev. D 100 , 044057 (2019), arXiv:1904.12983 [grqc].\n- [132] S. Vagnozzi and L. Visinelli, Phys. Rev. D 100 , 024020 (2019), arXiv:1905.12421 [gr-qc].\n- [133] T. Zhu, Q. Wu, M. Jamil, and K. Jusufi, Phys. Rev. D 100 , 044055 (2019), arXiv:1906.05673 [gr-qc].\n- [134] P. V. P. Cunha, C. A. R. Herdeiro, and E. Radu, Universe 5 , 220 (2019), arXiv:1909.08039 [gr-qc].\n- [135] I. Banerjee, S. Chakraborty, and S. SenGupta, Phys. Rev. D 101 , 041301 (2020), arXiv:1909.09385 [gr-qc].\n- [136] I. Banerjee, S. Sau, and S. SenGupta, Phys. Rev. D 101 , 104057 (2020), arXiv:1911.05385 [gr-qc].\n- [137] A. Allahyari, M. Khodadi, S. Vagnozzi, and D. F. Mota, JCAP 02 , 003 (2020), arXiv:1912.08231 [gr-qc].\n- [138] S. Vagnozzi, C. Bambi, and L. Visinelli, Class. Quant. Grav. 37 , 087001 (2020), arXiv:2001.02986 [gr-qc].\n- [139] M. Khodadi, A. Allahyari, S. Vagnozzi, and D. F. Mota, JCAP 09 , 026 (2020), arXiv:2005.05992 [gr-qc].\n- [140] R. Kumar, A. Kumar, and S. G. Ghosh, Astrophys. J. 896 , 89 (2020), arXiv:2006.09869 [gr-qc].\n- [141] M. Khodadi and E. N. Saridakis, Phys. Dark Univ. 32 , 100835 (2021), arXiv:2012.05186 [gr-qc].\n- [142] R. C. Pantig, P. K. Yu, E. T. Rodulfo, and A. Ovgun, Annals Phys. 436 , 168722 (2022), arXiv:2104.04304 [gr-qc]  \n-\n- [143] M. Khodadi, G. Lambiase, and D. F. Mota, JCAP 09 , 028 (2021), arXiv:2107.00834 [gr-qc].\n- [144] R. Roy, S. Vagnozzi, and L. Visinelli, Phys. Rev. D 105 , 083002 (2022), arXiv:2112.06932 [astro-ph.HE].\n- [145] A. Uniyal, R. C. Pantig, and A. Ovgun, Phys. Dark Univ. 40 , 101178 (2023), arXiv:2205.11072 [gr-qc].\n- [146] R. C. Pantig and A. Ovgun, Annals Phys. 448 , 169197 (2023), arXiv:2206.02161 [gr-qc].\n- [147] S. G. Ghosh and M. Afrin, Astrophys. J. 944 , 174 (2023), arXiv:2206.02488 [gr-qc].\n- [148] M. Khodadi and G. Lambiase, Phys. Rev. D 106 , 104050 (2022), arXiv:2206.08601 [gr-qc].\n- [149] R. Kumar Walia, S. G. Ghosh, and S. D. Maharaj, Astrophys. J. 939 , 77 (2022), arXiv:2207.00078 [gr-qc].\n- [150] R. Shaikh, Mon. Not. Roy. Astron. Soc. 523 , 375 (2023), arXiv:2208.01995 [gr-qc].\n- [151] S. D. Odintsov and V. K. Oikonomou, EPL 139 , 59003 (2022), arXiv:2208.07972 [gr-qc].\n- [152] V. K. Oikonomou, P. Tsyba, and O. Razina, Universe 8 , 484 (2022), arXiv:2209.04669 [gr-qc].\n- [153] R. C. Pantig, Chin. J. Phys. 87 , 49 (2024), arXiv:2303.01698 [gr-qc].\n- [154] E. Gonz'alez, K. Jusufi, G. Leon, and E. N. Saridakis, Phys. Dark Univ. 42 , 101304 (2023), arXiv:2305.14305 [astro-ph.CO].\n- [155] S. K. Sahoo, N. Yadav, and I. Banerjee, Phys. Rev. D 109 , 044008 (2024), arXiv:2305.14870 [gr-qc].\n- [156] K. Nozari and S. Saghafi, Eur. Phys. J. C 83 , 588 (2023), arXiv:2305.17237 [gr-qc].\n- [157] A. Uniyal, S. Chakrabarti, M. Fathi, and A. Ovgun, Annals Phys. 462 , 169614 (2024), arXiv:2309.13680 [grqc].\n- [158] A. A. A. Filho, J. A. A. S. Reis, and H. Hassanabadi, JCAP 05 , 029 (2024), arXiv:2309.15778 [gr-qc].\n- [159] M. A. Raza, J. Rayimbaev, F. Sarikulov, M. Zubair, B. Ahmedov, and Z. Stuchlik, Phys. Dark Univ. 44 , 101488 (2024), arXiv:2311.15784 [gr-qc].\n- [160] H. Hoshimov, O. Yunusov, F. Atamurotov, M. Jamil, and A. Abdujabbarov, Phys. Dark Univ. 43 , 101392 (2024), arXiv:2312.10678 [gr-qc].\n- [161] L. Chakhchi, H. El Moumni, and K. Masmar, Phys. Dark Univ. 44 , 101501 (2024), arXiv:2403.09756 [gr-qc].\n- [162] W. Liu, D. Wu, X. Fang, J. Jing, and J. Wang, JCAP 08 , 035 (2024), arXiv:2406.00579 [gr-qc].\n- [163] W. Liu, D. Wu, and J. Wang, (2024), arXiv:2407.07416 [gr-qc].\n- [164] M. Khodadi, S. Vagnozzi, and J. T. Firouzjaee, Sci. Rep. 14 , 26932 (2024), arXiv:2408.03241 [gr-qc].\n- [165] W. Liu, D. Wu, and J. Wang, Phys. Lett. B 858 , 139052 (2024), arXiv:2408.05569 [gr-qc].\n- [166] S. Nojiri and S. D. Odintsov, Phys. Dark Univ. 46 , 101669 (2024), arXiv:2408.05668 [gr-qc].\n- [167] A. Katz, J. Kopp, S. Sibiryakov, and W. Xue, JCAP 12 , 005 (2018), arXiv:1807.11495 [astro-ph.CO].\n- [168] Y. Bai and N. Orlofsky, Phys. Rev. D 99 , 123019 (2019), arXiv:1812.01427 [astro-ph.HE].\n- [169] N. Smyth, S. Profumo, S. English, T. Jeltema, K. McKinnon, and P. Guhathakurta, Phys. Rev. D 101 , 063005 (2020), arXiv:1910.01285 [astro-ph.CO].\n- [170] A. Coogan, L. Morrison, and S. Profumo, Phys. Rev. Lett. 126 , 171101 (2021), arXiv:2010.04797 [astroph.CO].\n- [171] A. Ray, R. Laha, J. B. Mu˜noz, and R. Caputo, Phys. Rev. D 104 , 023516 (2021), arXiv:2102.06714 [astroph.CO].\n- [172] J. Auffinger, Eur. Phys. J. C 82 , 384 (2022), arXiv:2201.01265 [astro-ph.HE].\n- [173] D. Ghosh and A. K. Mishra, Phys. Rev. D 109 , 043537 (2024), arXiv:2208.14279 [hep-ph].\n- [174] A. L. Miller, N. Aggarwal, S. Clesse, and F. De Lillo, Phys. Rev. D 105 , 062008 (2022), arXiv:2110.06188 [grqc].\n- [175] N. P. Branco, R. Z. Ferreira, and J. a. G. Rosa, JCAP 04 , 003 (2023), arXiv:2301.01780 [hep-ph].\n- [176] B. Bertrand, M. Cuadrat-Grzybowski, P. Defraigne, M. Van Camp, and S. Clesse (2023) arXiv:2312.14520 [astro-ph.CO].\n- [177] T. X. Tran, S. R. Geller, B. V. Lehmann, and D. I. Kaiser, Phys. Rev. D 110 , 063533 (2024), arXiv:2312.17217 [astro-ph.CO].\n- [178] M. Gorton and A. M. Green, SciPost Phys. 17 , 032 (2024), arXiv:2403.03839 [astro-ph.CO].\n- [179] J. B. Dent, B. Dutta, and T. Xu, (2024), arXiv:2404.02956 [hep-ph].\n- [180] M. Tamta, N. Raj, and P. Sharma, (2024), arXiv:2405.20365 [astro-ph.HE].\n- [181] P. Tinyakov, (2024), arXiv:2406.03114 [astro-ph.CO].\n- [182] A. Loeb, Res. Notes AAS 8 , 211 (2024), arXiv:2408.10799 [hep-ph].\n- [183] S. Ansoldi, in Conference on Black Holes and Naked Singularities (2008) arXiv:0802.0330 [gr-qc].\n- [184] P. Nicolini, Int. J. Mod. Phys. A 24 , 1229 (2009), arXiv:0807.1939 [hep-th].\n- [185] L. Sebastiani and S. Zerbini, Astronomy 1 , 99 (2022), arXiv:2206.03814 [gr-qc].\n- [186] R. Torres, (2022), arXiv:2208.12713 [gr-qc].\n- [187] C. Lan, H. Yang, Y. Guo, and Y.-G. Miao, Int. J. Theor. Phys. 62 , 202 (2023), arXiv:2303.11696 [gr-qc].\n- [188] J. Bardeen, in Proceedings, 5th International Conference on Gravitation and the theory of relativity: Tbilisi, USSR, December 9-13, 1968 (1968) p. 174.\n- [189] S. A. Hayward, Phys. Rev. Lett. 96 , 031103 (2006), arXiv:gr-qc/0506126.\n- [190] H. Culetu, (2013), arXiv:1305.5964 [gr-qc].\n- [191] H. Culetu, Int. J. Theor. Phys. 54 , 2855 (2015), arXiv:1408.3334 [gr-qc].\n- [192] S. G. Ghosh, Eur. Phys. J. C 75 , 532 (2015), arXiv:1408.5668 [gr-qc].\n- [193] A. Simpson and M. Visser, Universe 6 , 8 (2019), arXiv:1911.01020 [gr-qc].\n- [194] M. Calz'a, D. Pedrotti, and S. Vagnozzi (2024) arXiv:2409.02807 [gr-qc].\n- [195] I. Sachs, M. Schneider, and M. Urban, Phys. Rev. D 104 , 125020 (2021), arXiv:2105.01071 [hep-th].\n- [196] A. Ashtekar, T. De Lorenzo, and M. Schneider, Adv. Theor. Math. Phys. 25 , 1651 (2021), arXiv:2107.08506 [gr-qc].\n- [197] A. Ashtekar, A. del R'ıo, and M. Schneider, Gen. Rel. Grav. 54 , 45 (2022), arXiv:2205.00298 [gr-qc].\n- [198] I. Dymnikova, Gen. Rel. Grav. 24 , 235 (1992).\n- [199] I. Dymnikova and E. Galaktionov, Class. Quant. Grav. 22 , 2331 (2005), arXiv:gr-qc/0409049.\n- [200] A. Ashtekar and M. Bojowald, Class. Quant. Grav. 22 , 3349 (2005), arXiv:gr-qc/0504029.\n- [201] M. V. Bebronne and P. G. Tinyakov, JHEP 04 , 100 (2009), [Erratum: JHEP 06, 018 (2011)],  \n- arXiv:0902.3899 [gr-qc].\n- [202] L. Modesto, J. W. Moffat, and P. Nicolini, Phys. Lett. B 695 , 397 (2011), arXiv:1010.0680 [gr-qc].\n- [203] E. Spallucci and S. Ansoldi, Phys. Lett. B 701 , 471 (2011), arXiv:1101.2760 [hep-th].\n- [204] A. Perez, Class. Quant. Grav. 32 , 084001 (2015), arXiv:1410.7062 [gr-qc].\n- [205] A. Coll'eaux, S. Chinaglia, and S. Zerbini, Int. J. Mod. Phys. D 27 , 1830002 (2018), arXiv:1712.03730 [gr-qc].\n- [206] P. Nicolini, E. Spallucci, and M. F. Wondrak, Phys. Lett. B 797 , 134888 (2019), arXiv:1902.11242 [gr-qc].\n- [207] L. Bosma, B. Knorr, and F. Saueressig, Phys. Rev. Lett. 123 , 101301 (2019), arXiv:1904.04845 [hep-th].\n- [208] K. Jusufi, Annals Phys. 448 , 169191 (2023), arXiv:2208.12979 [gr-qc].\n- [209] G. J. Olmo and D. Rubiera-Garcia, (2022), arXiv:2209.05061 [gr-qc].\n- [210] K. Jusufi, Phys. Dark Univ. 39 , 101156 (2023), arXiv:2212.06760 [gr-qc].\n- [211] A. Ashtekar, J. Olmedo, and P. Singh, (2023), arXiv:2301.01309 [gr-qc].\n- [212] P. Nicolini (2023) arXiv:2306.01480 [gr-qc].\n- [213] C. Bambi, Regular Black Holes , Springer Series in Astrophysics and Cosmology (Springer Singapore, 2023) arXiv:2307.13249 [gr-qc].\n- [214] A. Borde, Phys. Rev. D 55 , 7615 (1997), arXiv:grqc/9612057.\n- [215] E. Ayon-Beato and A. Garcia, Phys. Rev. Lett. 80 , 5056 (1998), arXiv:gr-qc/9911046.\n- [216] E. Ayon-Beato and A. Garcia, Phys. Lett. B 464 , 25 (1999), arXiv:hep-th/9911174.\n- [217] K. A. Bronnikov and J. C. Fabris, Phys. Rev. Lett. 96 , 251101 (2006), arXiv:gr-qc/0511109.\n- [218] W. Berej, J. Matyjasek, D. Tryniecki, and M. Woronowicz, Gen. Rel. Grav. 38 , 885 (2006), arXiv:hepth/0606185.\n- [219] K. A. Bronnikov, R. A. Konoplya, and A. Zhidenko, Phys. Rev. D 86 , 024028 (2012), arXiv:1205.2224 [grqc].\n- [220] M. Rinaldi, Phys. Rev. D 86 , 084048 (2012), arXiv:1208.0103 [gr-qc].\n- [221] Z. Stuchl'ık and J. Schee, Int. J. Mod. Phys. D 24 , 1550020 (2014), arXiv:1501.00015 [astro-ph.HE].\n- [222] J. Schee and Z. Stuchlik, JCAP 06 , 048 (2015), arXiv:1501.00835 [astro-ph.HE].\n- [223] T. Johannsen, Phys. Rev. D 88 , 044002 (2013), arXiv:1501.02809 [gr-qc].\n- [224] R. Myrzakulov, L. Sebastiani, S. Vagnozzi, and S. Zerbini, Class. Quant. Grav. 33 , 125005 (2016), arXiv:1510.02284 [gr-qc].\n- [225] Z.-Y. Fan and X. Wang, Phys. Rev. D 94 , 124027 (2016), arXiv:1610.02636 [gr-qc].\n- [226] L. Sebastiani, S. Vagnozzi, and R. Myrzakulov, Adv. High Energy Phys. 2017 , 3156915 (2017), arXiv:1612.08661 [gr-qc].\n- [227] B. Toshmatov, Z. Stuchl'ık, and B. Ahmedov, Phys. Rev. D 95 , 084037 (2017), arXiv:1704.07300 [gr-qc].\n- [228] S. Chinaglia and S. Zerbini, Gen. Rel. Grav. 49 , 75 (2017), arXiv:1704.08516 [gr-qc].\n- [229] V. P. Frolov, EPJ Web Conf. 168 , 01001 (2018), arXiv:1708.04698 [gr-qc].\n- [230] M. Bertipagani, M. Rinaldi, L. Sebastiani, and S. Zerbini, Phys. Dark Univ. 33 , 100853 (2021), arXiv:2012.15645 [gr-qc].\n- [231] G. G. L. Nashed and E. N. Saridakis, JCAP 05 , 017 (2022), arXiv:2111.06359 [gr-qc].\n- [232] A. Simpson and M. Visser, JCAP 03 , 011 (2022), arXiv:2111.12329 [gr-qc].\n- [233] E. Franzin, S. Liberati, J. Mazza, R. Dey, and S. Chakraborty, Phys. Rev. D 105 , 124051 (2022), arXiv:2201.01650 [gr-qc].\n- [234] L. Chataignier, A. Y. Kamenshchik, A. Tronconi, and G. Venturi, Phys. Rev. D 107 , 023508 (2023), arXiv:2208.02280 [gr-qc].\n- [235] R. Ghosh, M. Rahman, and A. K. Mishra, Eur. Phys. J. C 83 , 91 (2023), arXiv:2209.12291 [gr-qc].\n- [236] M. Khodadi and R. Pourkhodabakhshi, Phys. Rev. D 106 , 084047 (2022), arXiv:2210.06861 [gr-qc].\n- [237] D. Farrah et al. , Astrophys. J. Lett. 944 , L31 (2023), arXiv:2302.07878 [astro-ph.CO].\n- [238] M. Fontana and M. Rinaldi, Phys. Rev. D 108 , 125003 (2023), arXiv:2302.08804 [gr-qc].\n- [239] K. Boshkayev, A. Idrissov, O. Luongo, and M. Muccino, Phys. Rev. D 108 , 044063 (2023), arXiv:2303.03248 [astro-ph.HE].\n- [240] O. Luongo, S. Mancini, and P. Pierosara, Phys. Rev. D 108 , 104059 (2023), arXiv:2304.06593 [gr-qc].\n- [241] O. Luongo and H. Quevedo, Class. Quant. Grav. 41 , 125011 (2024), arXiv:2305.11185 [gr-qc].\n- [242] M. Cadoni, A. P. Sanna, M. Pitzalis, B. Banerjee, R. Murgia, N. Hazra, and M. Branchesi, JCAP 11 , 007 (2023), arXiv:2306.11588 [gr-qc].\n- [243] R. Giamb'o and O. Luongo, Class. Quant. Grav. 41 , 125005 (2024), arXiv:2308.10060 [gr-qc].\n- [244] M. Cadoni, R. Murgia, M. Pitzalis, and A. P. Sanna, JCAP 03 , 026 (2024), arXiv:2309.16444 [gr-qc].\n- [245] O. Luongo, H. Quevedo, and S. N. Sajadi, Gen. Rel. Grav. 56 , 17 (2024), arXiv:2311.13264 [gr-qc].\n- [246] S. N. Sajadi, M. Khodadi, O. Luongo, and H. Quevedo, Phys. Dark Univ. 45 , 101525 (2024), arXiv:2312.16081 [gr-qc].\n- [247] F. Javed, S. Mumtaz, G. Mustafa, F. Atamurotov, and S. G. Ghosh, Chin. J. Phys. 88 , 55 (2024).\n- [248] A. Ditta, T. Xia, R. Ali, G. Mustafa, G. Mustafa, and A. Mahmood, Phys. Dark Univ. 43 , 101418 (2024).\n- [249] A. Al-Badawi, S. Shaymatov, M. Alloqulov, and A. Wang, Commun. Theor. Phys. 76 , 085401 (2024), arXiv:2401.12723 [gr-qc].\n- [250] A. Ovgun, R. C. Pantig, and A. Rinc'on, Annals Phys. 463 , 169625 (2024), arXiv:2402.14190 [gr-qc].\n- [251] D. Corona, R. Giamb'o, and O. Luongo, Int. J. Geom. Meth. Mod. Phys. 21 , 2440019 (2024), arXiv:2402.18997 [gr-qc].\n- [252] P. Bueno, P. A. Cano, and R. A. Hennigar, (2024), arXiv:2403.04827 [gr-qc].\n- [253] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 109 , 104005 (2024), arXiv:2403.07848 [gr-qc].\n- [254] D. Pedrotti and S. Vagnozzi, Phys. Rev. D 110 , 084075 (2024), arXiv:2404.07589 [gr-qc].\n- [255] K. A. Bronnikov, Phys. Rev. D 110 , 024021 (2024), arXiv:2404.14816 [gr-qc].\n- [256] Y. Kurmanov, K. Boshkayev, T. Konysbayev, O. Luongo, N. Saiyp, A. Urazalina, G. Ikhsan, and G. Suliyeva, Phys. Dark Univ. 46 , 101566 (2024), arXiv:2404.15437 [gr-qc].\n- [257] S. V. Bolokhov, K. A. Bronnikov, and M. V. Skvortsova, Grav. Cosmol. 30 , 265 (2024),  \n- arXiv:2405.09124 [gr-qc].\n- [258] A. S. Agrawal, S. Zerbini, and B. Mishra, Phys. Dark Univ. 46 , 101637 (2024), arXiv:2406.01241 [gr-qc].\n- [259] A. Belfiglio, S. M. Chandran, O. Luongo, and S. Mancini, (2024), arXiv:2407.03775 [gr-qc].\n- [260] O. Stashko, Phys. Rev. D 110 , 084016 (2024), arXiv:2407.07892 [gr-qc].\n- [261] V. Faraoni and M. Rinaldi, Phys. Rev. D 110 , 063553 (2024), arXiv:2407.14549 [gr-qc].\n- [262] R. A. Konoplya and O. S. Stashko, (2024), arXiv:2408.02578 [gr-qc].\n- [263] M. Khodadi and J. T. Firouzjaee, Phys. Lett. B 857 , 138986 (2024), arXiv:2408.12873 [gr-qc].\n- [264] M. Calz'a, F. Gianesello, M. Rinaldi, and S. Vagnozzi, (2024), arXiv:2409.01801 [gr-qc].\n- [265] P. O. Mazur and E. Mottola, Universe 9 , 88 (2023), arXiv:gr-qc/0109035.\n- [266] P. O. Mazur and E. Mottola, Proc. Nat. Acad. Sci. 101 , 9545 (2004), arXiv:gr-qc/0407075.\n- [267] E. Mottola, (2023), 10.1007/978-981-99-1596-5 8, arXiv:2302.09690 [gr-qc].\n- [268] K. A. Bronnikov, Int. J. Mod. Phys. D 27 , 1841005 (2018), arXiv:1711.00087 [gr-qc].\n- [269] S. G. Ghosh and R. K. Walia, Annals Phys. 434 , 168619 (2021), arXiv:2109.13031 [gr-qc].\n- [270] K. A. Bronnikov and R. K. Walia, Phys. Rev. D 105 , 044039 (2022), arXiv:2112.13198 [gr-qc].\n- [271] K. A. Bronnikov, (2022), arXiv:2211.00743 [gr-qc].\n- [272] A. Bokulic, E. Franzin, T. Juric, and I. Smolic, Phys. Lett. B 854 , 138750 (2024), arXiv:2311.17151 [gr-qc].\n- [273] D. A. Easson, JHEP 02 , 037 (2003), arXiv:hepth/0210016.\n- [274] I. Dymnikova and M. Khlopov, Int. J. Mod. Phys. D 24 , 1545002 (2015), arXiv:1510.01351 [gr-qc].\n- [275] J. A. d. F. Pacheco, Universe 4 , 62 (2018), arXiv:1805.03053 [gr-qc].\n- [276] A. Arbey and J. Auffinger, Eur. Phys. J. C 81 , 910 (2021), arXiv:2108.02737 [gr-qc].\n- [277] A. Arbey, J. Auffinger, M. Geiller, E. Livine, and F. Sartini, PoS EPS-HEP2021 , 066 (2022).\n- [278] A. Banerjee, L. Singh Bhandari, A. Jain, and A. M. Thalapillil, (2024), arXiv:2406.08728 [astro-ph.CO].\n- [279] P. C. W. Davies, D. A. Easson, and P. B. Levin, (2024), arXiv:2410.21577 [hep-th].\n- [280] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15 , 2738 (1977).\n- [281] H. Lu, S. Di Gennaro, and Y. C. Ong, (2024), arXiv:2407.00484 [gr-qc].\n- [282] S. Vagnozzi et al. , Class. Quant. Grav. 40 , 165007 (2023), arXiv:2205.07787 [gr-qc].\n- [283] D. I. Kazakov and S. N. Solodukhin, Nucl. Phys. B 429 , 153 (1994), arXiv:hep-th/9310150.\n- [284] P. Kocherlakota et al. (Event Horizon Telescope), Phys. Rev. D 103 , 104047 (2021), arXiv:2105.09343 [gr-qc].\n- [285] K. Akiyama et al. (Event Horizon Telescope), Astrophys. J. Lett. 930 , L17 (2022), arXiv:2311.09484 [astroph.HE].\n- [286] K. A. Bronnikov and R. A. Konoplya, (2023), arXiv:2306.11083 [gr-qc].\n- [287] V. Vertogradov and A. Ovgun, Phys. Lett. B 854 , 138758 (2024), arXiv:2404.18536 [gr-qc].\n- [288] V. Vertogradov and A. Ovgun, Phys. Dark Univ. 45 , 101541 (2024), arXiv:2404.04046 [gr-qc].\n- [289] V. Vertogradov, A. Ovgun, and R. C. Pantig, (2024), arXiv:2405.05077 [gr-qc].\n- [290] A. Arbey, J. Auffinger, and J. Silk, Mon. Not. Roy. Astron. Soc. 494 , 1257 (2020), arXiv:1906.04196 [astroph.CO].\n- [291] A. Arbey, J. Auffinger, and J. Silk, Phys. Rev. D 101 , 023010 (2020), arXiv:1906.04750 [astro-ph.CO].\n- [292] A. Arbey, J. Auffinger, and J. Silk, PoS ICHEP2020 , 585 (2021), arXiv:2012.14767 [astro-ph.CO].\n- [293] T. Zhou and L. Modesto, Phys. Rev. D 107 , 044016 (2023), arXiv:2208.02557 [gr-qc].\n- [294] R. Carballo-Rubio, F. Di Filippo, S. Liberati, C. Pacilio, and M. Visser, JHEP 05 , 132 (2021), arXiv:2101.05006 [gr-qc].\n- [295] R. Carballo-Rubio, F. Di Filippo, S. Liberati, C. Pacilio, and M. Visser, JHEP 09 , 118 (2022), arXiv:2205.13556 [gr-qc].\n- [296] A. Bonanno, A.-P. Khosravi, and F. Saueressig, Phys. Rev. D 107 , 024005 (2023), arXiv:2209.10612 [gr-qc].\n- [297] A. Bonanno and F. Saueressig, (2022), arXiv:2211.09192 [gr-qc].\n- [298] R. Carballo-Rubio, F. Di Filippo, S. Liberati, C. Pacilio, and M. Visser, Phys. Rev. D 108 , 128501 (2023), arXiv:2212.07458 [gr-qc].\n- [299] E. Ayon-Beato and A. Garcia, Phys. Lett. B 493 , 149 (2000), arXiv:gr-qc/0009077.\n- [300] E. Ayon-Beato and A. Garcia, Gen. Rel. Grav. 37 , 635 (2005), arXiv:hep-th/0403229.\n- [301] R. V. Maluf and J. C. S. Neves, Int. J. Mod. Phys. D 28 , 1950048 (2018), arXiv:1801.08872 [gr-qc].\n- [302] A. D. Sakharov, Sov. Phys. JETP 22 , 241 (1966).\n- [303] E. B. Gliner, Soviet Journal of Experimental and Theoretical Physics 22 , 378 (1966).\n- [304] M. A. Markov, Soviet Journal of Experimental and Theoretical Physics Letters 36 , 265 (1982).\n- [305] M. A. Markov, Soviet Journal of Experimental and Theoretical Physics Letters 46 , 431 (1987).\n- [306] V. F. Mukhanov and R. H. Brandenberger, Phys. Rev. Lett. 68 , 1969 (1992).\n- [307] A. Kumar, D. V. Singh, and S. G. Ghosh, Annals Phys. 419 , 168214 (2020), arXiv:2003.14016 [gr-qc].\n- [308] S. I. Kruglov, Grav. Cosmol. 27 , 78 (2021), arXiv:2103.14087 [gr-qc].\n- [309] A. Addazi et al. , Prog. Part. Nucl. Phys. 125 , 103948 (2022), arXiv:2111.05659 [hep-ph].\n- [310] R. Alves Batista et al. , (2023), arXiv:2312.00409 [grqc].\n- [311] R. Kumar and S. G. Ghosh, Class. Quant. Grav. 38 , 8 (2021), arXiv:2004.07501 [gr-qc].\n- [312] D. V. Singh, S. G. Ghosh, and S. D. Maharaj, Nucl. Phys. B 981 , 115854 (2022).\n- [313] I. Sakalli and S. Kanzi, Turk. J. Phys. 46 , 51 (2022), arXiv:2205.01771 [hep-th].\n- [314] R. A. Konoplya and A. Zhidenko, JCAP 09 , 068 (2024), arXiv:2406.11694 [gr-qc].\n- [315] R. A. Konoplya and A. Zhidenko, (2024), arXiv:2408.11162 [gr-qc].\n- [316] S. A. Teukolsky, Astrophys. J. 185 , 635 (1973).\n- [317] E. D. Fackerell and R. G. Crossman, Journal of Mathematical Physics 18 , 1849 (1977), https://pubs.aip.org/aip/jmp/articlepdf/18/9/1849/19234633/1849 1 online.pdf.\n- [318] K. G. Suffern, E. D. Fackerell, and C. M. Cosgrove, Journal of Mathematical Physics 24 ,  \n- 1350 (1983), https://pubs.aip.org/aip/jmp/articlepdf/24/5/1350/19310962/1350 1 online.pdf.\n- [319] E. Seidel, Classical and Quantum Gravity 6 , 1057 (1989).\n- [320] E. Berti, V. Cardoso, and M. Casals, Phys. Rev. D 73 , 024013 (2006), [Erratum: Phys.Rev.D 73, 109902 (2006)], arXiv:gr-qc/0511111.\n- [321] C. M. Harris and P. Kanti, JHEP 10 , 014 (2003), arXiv:hep-ph/0309054.\n- [322] A. Arbey, J. Auffinger, M. Geiller, E. R. Livine, and F. Sartini, Phys. Rev. D 103 , 104010 (2021), arXiv:2101.02951 [gr-qc].\n- [323] J. G. Rosa, Phys. Rev. D 95 , 064017 (2017), arXiv:1612.01826 [gr-qc].\n- [324] J. G. Rosa, JHEP 02 , 014 (2013), arXiv:1209.4211 [hepth].\n- [325] M. Calz'a, J. March-Russell, and J. a. G. Rosa, (2021), arXiv:2110.13602 [astro-ph.CO].\n- [326] M. Calz'a, Phys. Rev. D 107 , 044067 (2023), arXiv:2207.10467 [gr-qc].\n- [327] M. Calz'a and J. a. G. Rosa, JHEP 12 , 090 (2022), arXiv:2210.06500 [gr-qc].\n- [328] M. Calz'a, J. a. G. Rosa, and F. Serrano, JHEP 05 , 140 (2024), arXiv:2306.09430 [hep-ph].\n- [329] M. Calz'a and J. a. G. Rosa, JHEP 08 , 012 (2024), arXiv:2311.12930 [gr-qc].\n- [330] M. Calz'a and J. a. G. Rosa, (2023), arXiv:2312.09261 [hep-ph].\n- [331] A. H. Wilson, Proceedings of The Royal Society A: Mathematical, Physical and Engineering Sciences 118 , 617 (1928).\n- [332] W. G. Baber and H. R. Hass'e, Mathematical Proceedings of the Cambridge Philosophical Society 31 , 564-581 (1935).\n- [333] E. W. Leaver, Proc. Roy. Soc. Lond. A 402 , 285 (1985).\n- [334] E. W. Leaver, J. Math. Phys. 27 , 1238 (1986).\n- [335] R. A. Konoplya, D. Ovchinnikov, and B. Ahmedov, Phys. Rev. D 108 , 104054 (2023), arXiv:2307.10801 [grqc].\n- [336] R. A. Konoplya, JCAP 07 , 001 (2023), arXiv:2305.09187 [gr-qc].\n- [337] S. W. Hawking, in 1st Oxford Conference on Quantum Gravity (1975).\n- [338] S. W. Hawking, Commun. Math. Phys. 43 , 199 (1975), [Erratum: Commun.Math.Phys. 46, 206 (1976)].\n- [339] D. N. Page, Phys. Rev. D 13 , 198 (1976).\n- [340] D. N. Page, Phys. Rev. D 14 , 3260 (1976).\n- [341] D. N. Page, Phys. Rev. D 16 , 2402 (1977).\n- [342] J. Auffinger, Prog. Part. Nucl. Phys. 131 , 104040 (2023), arXiv:2206.02672 [astro-ph.CO].\n- [343] F. Kuhnel, C. Rampf, and M. Sandstad, Eur. Phys. J. C 76 , 93 (2016), arXiv:1512.00488 [astro-ph.CO].\n- [344] F. Kuhnel and K. Freese, Phys. Rev. D 95 , 083508 (2017), arXiv:1701.07223 [astro-ph.CO].\n- [345] B. Carr, M. Raidal, T. Tenkanen, V. Vaskonen, and H. Veermae, Phys. Rev. D 96 , 023514 (2017), arXiv:1705.05567 [astro-ph.CO].\n- [346] M. Raidal, V. Vaskonen, and H. Veermae, JCAP 09 , 037 (2017), arXiv:1707.01480 [astro-ph.CO].\n- [347] N. Bellomo, J. L. Bernal, A. Raccanelli, and L. Verde, JCAP 01 , 004 (2018), arXiv:1709.07467 [astro-ph.CO].\n- [348] B. V. Lehmann, S. Profumo, and J. Yant, JCAP 04 , 007 (2018), arXiv:1801.00808 [astro-ph.CO].\n- [349] B. Carr and F. Kuhnel, Phys. Rev. D 99 , 103535 (2019),\n- arXiv:1811.06532 [astro-ph.CO].\n- [350] A. D. Gow, C. T. Byrnes, A. Hall, and J. A. Peacock, JCAP 01 , 031 (2020), arXiv:1911.12685 [astro-ph.CO].\n- [351] V. De Luca, G. Franciolini, and A. Riotto, Phys. Lett. B 807 , 135550 (2020), arXiv:2001.04371 [astro-ph.CO].\n- [352] A. D. Gow, C. T. Byrnes, and A. Hall, Phys. Rev. D 105 , 023503 (2022), arXiv:2009.03204 [astro-ph.CO].\n- [353] A. Ashoorioon, A. Rostami, and J. T. Firouzjaee, Phys. Rev. D 103 , 123512 (2021), arXiv:2012.02817 [astroph.CO].\n- [354] E. Bagui and S. Clesse, Phys. Dark Univ. 38 , 101115 (2022), arXiv:2110.07487 [astro-ph.CO].\n- [355] U. Mukhopadhyay, D. Majumdar, and A. Halder, JCAP 10 , 099 (2022), arXiv:2203.13008 [astro-ph.CO].\n- [356] T. Papanikolaou, JCAP 10 , 089 (2022), arXiv:2207.11041 [astro-ph.CO].\n- [357] Y.-F. Cai, C. Tang, G. Mo, S.-F. Yan, C. Chen, X.-H. Ma, B. Wang, W. Luo, D. A. Easson, and A. Marciano, Sci. China Phys. Mech. Astron. 67 , 259512 (2024), arXiv:2301.09403 [astro-ph.CO].\n- [358] B. J. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, Phys. Rev. D 81 , 104019 (2010), arXiv:0912.5297 [astroph.CO].\n- [359] D. E. Gruber, J. L. Matteson, L. E. Peterson, and G. V. Jung, Astrophys. J. 520 , 124 (1999), arXiv:astroph/9903492.\n- [360] V. Schoenfelder (COMPTEL), Astron. Astrophys. Suppl. Ser. 143 , 145 (2000), arXiv:astro-ph/0002366.\n- [361] A. W. Strong, I. V. Moskalenko, and O. Reimer, Astrophys. J. 613 , 956 (2004), arXiv:astro-ph/0405441.\n- [362] B. J. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, Phys. Rev. D 94 , 044029 (2016), arXiv:1604.05349 [astro-ph.CO].\n- [363] M. Boudaud and M. Cirelli, Phys. Rev. Lett. 122 , 041104 (2019), arXiv:1807.03075 [astro-ph.HE].\n- [364] W. DeRocco and P. W. Graham, Phys. Rev. Lett. 123 , 251102 (2019), arXiv:1906.07740 [astro-ph.CO].  \n- [365] R. Laha, Phys. Rev. Lett. 123 , 251101 (2019), arXiv:1906.09994 [astro-ph.HE].\n- [366] B. Dasgupta, R. Laha, and A. Ray, Phys. Rev. Lett. 125 , 101101 (2020), arXiv:1912.01014 [hep-ph].\n- [367] M. Ackermann et al. (Fermi-LAT), Astrophys. J. 799 , 86 (2015), arXiv:1410.3696 [astro-ph.HE].\n- [368] E. Dwek and F. Krennrich, Astropart. Phys. 43 , 112 (2013), arXiv:1209.4661 [astro-ph.CO].\n- [369] B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, Rept. Prog. Phys. 84 , 116902 (2021), arXiv:2002.12778 [astro-ph.CO].\n- [370] F. Capela, M. Pshirkov, and P. Tinyakov, Phys. Rev. D 87 , 023507 (2013), arXiv:1209.6021 [astro-ph.CO].\n- [371] F. Capela, M. Pshirkov, and P. Tinyakov, Phys. Rev. D 87 , 123524 (2013), arXiv:1301.4984 [astro-ph.CO].\n- [372] P. Pani and A. Loeb, JCAP 06 , 026 (2014), arXiv:1401.3025 [astro-ph.CO].\n- [373] P. W. Graham, S. Rajendran, and J. Varela, Phys. Rev. D 92 , 063007 (2015), arXiv:1505.04444 [hep-ph].\n- [374] R. Ibata, C. Nipoti, A. Sollima, M. Bellazzini, S. Chapman, and E. Dalessandro, Mon. Not. Roy. Astron. Soc. 428 , 3648 (2013), arXiv:1210.7787 [astro-ph.CO].\n- [375] G. Defillon, E. Granet, P. Tinyakov, and M. H. G. Tytgat, Phys. Rev. D 90 , 103522 (2014), arXiv:1409.0469 [gr-qc].\n- [376] P. Montero-Camacho, X. Fang, G. Vasquez, M. Silva, and C. M. Hirata, JCAP 08 , 031 (2019), arXiv:1906.05950 [astro-ph.CO].\n- [377] K. Kainulainen, S. Nurmi, E. D. Schiappacasse, and T. T. Yanagida, Phys. Rev. D 104 , 123033 (2021), arXiv:2108.08717 [astro-ph.HE].\n- [378] D. W. P. Amaral and E. D. Schiappacasse, Phys. Rev. D 110 , 083532 (2024), arXiv:2312.09285 [hep-ph].\n- [379] C. M. Chambers, W. A. Hiscock, and B. Taylor, Phys. Rev. Lett. 78 , 3249 (1997), arXiv:gr-qc/9703018.\n- [380] B. E. Taylor, C. M. Chambers, and W. A. Hiscock, Phys. Rev. D 58 , 044012 (1998), arXiv:gr-qc/9801044."}