Primordial regular black holes as all the dark matter (I): tr-symmetric metrics
Marco Calzà Davide Pedrotti Sunny Vagnozzi





marco.calza@unitn.it
Department of Physics, University of Trento, Via Sommarive 14, 38123 Povo (TN), Italy
Trento Institute for Fundamental Physics and Applications (TIFPA)-INFN, Via Sommarive 14, 38123 Povo (TN), Italy
M.C. and D.P. contributed equally to this work


davide.pedrotti-1@unitn.it
Department of Physics, University of Trento, Via Sommarive 14, 38123 Povo (TN), Italy
Trento Institute for Fundamental Physics and Applications (TIFPA)-INFN, Via Sommarive 14, 38123 Povo (TN), Italy


sunny.vagnozzi@unitn.it
Department of Physics, University of Trento, Via Sommarive 14, 38123 Povo (TN), Italy
Trento Institute for Fundamental Physics and Applications (TIFPA)-INFN, Via Sommarive 14, 38123 Povo (TN), Italy

§ ABSTRACT 
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 non-singular objects, considering four phenomenological, regular tr-symmetric space-times, featuring either de Sitter or Minkowski cores. We characterize the evaporation of these PBHs and constrain their abundance from γ-ray observations. For three of the metrics, including the well-known Bardeen and Hayward ones, we show 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 up to an order of magnitude, leading to a much larger region of parameter space where PBHs can make up all the DM. A companion paper is instead devoted to non-tr-symmetric metrics, including loop quantum gravity-inspired ones. Our work provides a proof-of-principle for the interface between the DM and singularity problems being a promising arena with a rich phenomenology.



Primordial regular black holes as all the dark matter (I):
 tr-symmetric metrics
    Sunny Vagnozzi
    September 9, 2024
================================================================================





§ 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 <cit.>. 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 <cit.>. On the more theoretical side, continuous gravitational collapse in GR leads to the pathological appearance of curvature singularities <cit.>. The nature of DM and the singularity problem are arguably two among the most important open questions in theoretical physics.

The solution to the former 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 <cit.> (see e.g. Refs. <cit.> 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 <cit.>, turning these objects into extraordinary probes of fundamental physics <cit.>. 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^14 kg≲ M_pbh≲ 10^20 kg <cit.>: 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.

Almost all works on PBHs assume that these are Schwarzschild or Kerr BHs <cit.>. 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”? It is our goal in the present work to systematically address this question, which naturally merges the DM and singularity problems.

We entertain the possibility that PBHs are “regular”, i.e. free of curvature singularities <cit.>, and therefore that DM may be in the form of primordial regular BHs (PRBHs). For concreteness, we consider so-called tr (time-radius)-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 coordinates 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 four 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 <cit.>, Hayward BHs <cit.>, Ghosh-Culetu-Simpson-Visser BHs <cit.>, and a Hayward-like BH which, to the best of our knowledge, is introduced for the first time in this work. These four space-times present a rich phenomenology, including both de Sitter and Minkowski cores, and temperatures which can both decrease and increase with increasing regularizing parameter ℓ. We focus our attention on observational constraints from PRBH evaporation, which set the lower limit 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 primordial regular BHs 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 primordial regular BHs, but rather as a pilot study, pointing towards a direction which has thus far received very little attention and indicating promising directions for further work.

The rest of this paper is then organized as follows. In Sec. <ref> we briefly introduce the regular space-times studied in the rest of the work. Various aspects of our methodology are discussed in Sec. <ref>, with Sec. <ref> devoted to the calculation of the so-called greybody factors, Sec. <ref> to the computation of photon spectra resulting from Hawking evaporation, and Sec. <ref> 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. <ref>. Finally, in Sec. <ref> we draw concluding remarks. A number of more technical aspects concerning the greybody factors computation are discussed in Appendix <ref>. 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 <cit.>: this focuses on non-tr-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 <cit.>.



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

Even 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 holds valid, one can introduce metrics which are free of singularities in the entire space-time, and describe so-called regular BHs (RBHs) <cit.>. 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. <cit.> for an inevitably incomplete selection of examples, as well as Refs. <cit.> for recent reviews on the subject. [Another interesting possibility are gravastars, which are not RBHs in a strict sense <cit.>.] 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 <cit.>.

As 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 either opened or tightened depending on the features of the RBH metric. To the best of our knowledge, we are aware of only four works in this direction <cit.>. Ref. <cit.> studied the thermodynamics of primordial regular BHs, focusing however on the case where they do not evaporate, and therefore did not study constraints on f_pbh. Ref. <cit.> studied the evaporation of a loop quantum gravity-inspired BH, and weakened constraints on f_pbh were reported in a later proceeding (which however does not appear to be widely known). Finally, Ref. <cit.> studied signatures of primordial BHs with magnetic charge, which could be (as is often but necessarily the case) regular. The aim of this pilot study and our companion paper <cit.> 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.

In our work, we shall consider four different non-rotating 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 Boyer-Lindquist coordinates, their line elements take the following general form:


    ds^2 = -f(r)dt^2 + g(r)^-1dr^2 + h(r)dΩ^2 ,


where dΩ^2=dθ^2 +sin^2(θ) dϕ^2 is the metric on the 2-sphere. We also require our space-times to be asymptotically flat, which amounts to the following requirements:


    f(r)  1 ,    g(r)  1 ,    h(r)  r^2 .


In addition, as stated earlier, we require our space-times to be tr-symmetric (the non-tr-symmetric case is covered in a companion paper <cit.>), which imposes the following additional conditions:


    f(r) = g(r) ,    h(r) = r^2 ,


implying that the coordinate r is effectively the areal radius. With the conditions given by Eqs. (<ref>,<ref>) imposed upon Eq. (<ref>), our most general line element therefore takes the following form:


    ds^2 = -f(r) dt^2 + f(r)^-1dr^2 + r^2  ( dθ^2 +sin^2(θ) dϕ^2  )  .


In what follows, we refer to the function f(r) as being the “metric function”. The four different RBH solutions we consider, which we will discuss very shortly in Sections <ref>– <ref>, are characterized by different functional forms of f(r).

An 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. the one evaluated by Wick rotating the metric in the standard way and imposing regularity in the Euclidean period <cit.>. 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. <cit.>). 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:


    T=κ/2π=f'(r)/4π|_r_H ,


where the prime indicates a derivative with respect to r, and κ is the BH surface gravity, given by the following:


    κ=f'(r)/2|_r_H .


In Eqs. (<ref>,<ref>), r_H is the horizon radius, which is the solution to the following equation:


    g^-1(r_H)=f(r_H)=0 ,


with 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-2M/r, one recovers the well-known expressions r_H=2M and T_Sch=1/8π M. However, in more general space-times the horizon radius in Eq. (<ref>) is not guaranteed to have a closed form expression, and the same therefore holds for the temperature in Eq. (<ref>). In Fig. <ref> we show the evolution of the temperatures (normalized to the temperature of Schwarzschild BHs, T_Sch=1/8π M) of the four regular space-times we will introduce shortly, as a function of the regularizing parameter ℓ (normalized by the event horizon radius r_H). As the Figure clearly shows, for three of these space-times (Bardeen, Hayward, and Ghosh-Culetu-Simpson-Visser) the temperature is a monotonically decreasing function of ℓ, whereas for the fourth (Hayward-like), after a small initial dip, the temperature grows monotonically with growing ℓ. As we will see, these different behaviour turn out to have important phenomenological consequences.



A 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. A number of “popular” RBHs have indeed been shown to have finite curvature invariants but to be geodesically incomplete <cit.>. This includes the well-known Bardeen and Hayward RBHs, which are among the ones we shall consider here. However, given the significant interest in these metrics, the fact that they are widely taken as prototypes for RBHs, and our phenomenological goal of going beyond Schwarzschild PBHs, we will take these space-times into consideration, while cautioning the reader about the above issues, and therefore that these metrics 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 <cit.>.



 §.§ 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 <cit.>. It is characterized by the following metric function: [We note that in all the metrics considered, the parameter M appearing in the metric function can always be unambiguously identified with the BH mass (either the Komar, ADM, Misner-Sharp-Hernandez, or Brown-York mass). This is important as the later constraints on f_pbh as a function of M_pbh will identify the latter with M.]


    f_B(r)=1-2Mr^2/(r^2+ℓ^2)^3/2 ,


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

It 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 <cit.>, potentially within the context of a specific non-linear electrodynamics theory <cit.>. Another possible origin for the Bardeen RBH are quantum corrections to the uncertainty principle <cit.>. Irrespective of its origin, and consistently with the approach pursued for the other space-times, we consider this solution as a model-agnostic phenomenological toy model.



 §.§ Hayward black hole



Another widely known RBH space-time is the Hayward RBH <cit.>, characterized by the following metric function:


    f_H(r)=1-2Mr^2/r^3+2Mℓ^2 .


If 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. We note that the Schwarzschild metric is function is recovered in the ℓ→ 0 limit.

Just 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 introduced 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 <cit.>, finite density and finite curvature proposals <cit.>, theories of non-linear electrodynamics <cit.>, and more generally as the result of corrections due to quantum gravity <cit.>. 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.



 §.§ Ghosh-Culetu-Simpson-Visser black hole



The regular space-times considered so far feature dS cores, which in itself is a very common feature of several RBH metrics. Nevertheless, another interesting phenomenological possibility consists in considering “hollow” RBHs in which the central singularity is replaced by an asymptotically Minkowski core, where the associated energy density and associated 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 polynomical 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 properties and phenomenological implications of PRBHs, going beyond the dS core RBHs studied thus far.

With this in mind, we consider a RBH featuring a Minkowski core, independently studied by Ghosh <cit.>, Culetu <cit.>, as well as Simpson and Visser <cit.>. Although such a RBH does not have any particular name associated to it in the literature, here we shall conform to the name introduced in Ref. <cit.>, referring to it as GCSV BH (from the initials of the four authors above). This space-time is characterized by the following metric function:


    f_GCSV(r)= 1-2M/rexp ( -ℓ/r )  .


The horizon radius r_H, for which a closed form expression is not available in the Bardeen and Hayward cases, is given by:


    r_H=-ℓ/W  ( -ℓ/2 M )  ,


where W denotes the Lambert function. Considering the principal branch W_0, a real and positive horizon radius is present for:


    W_0  ( -ℓ/2 M ) ≤ 0  0 ≤ℓ < 2M/e ,


or, alternatively, 0 ≤ℓ < r_H. While the GCSV BH was original introduced purely on phenomenological/mathematical grounds, it was shown in Refs. <cit.> 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/2M, with M the BH mass. Nevertheless, as with all the other RBHs considered, here we shall treat the GCSV RBH as a toy model for a regular space-time possessing a Minkowski core.



 §.§ Hayward-like black hole



For all of the three RBHs discussed earlier, the temperature decreases with increasing regularizing parameter, while converging to the Schwarzschild value T_Sch=1/8π M as ℓ→ 0, see Fig. <ref>. It is therefore not unreasonable to expect that evaporation constraints on f_pbh may be weakened if the DM consists of these PRHBs (this will of course be checked explicitly later). We have verified that such a behaviour of temperature versus regularizing parameter holds for most PRBHs of phenomenological interest considered in the literature. [This is actually not unrelated to the fact that most RBHs, and more generally hairy BH solutions, feature shadows whose size decreases relative to the size of the Schwarzschild BH shadow, i.e. 3√(3)M, see Ref. <cit.> for further more detailed discussions (see also Refs. <cit.>). In fact, the equations governing BH temperature and shadow size bear some resemblance, so an increase/decrease in one is expected to lead to an increase/decrease in the other. Ref. <cit.> also presents a few space-times whose shadow radii increase with increasing hair parameter, and for these we would expect the temperatures to increase as well. However, for a few of these space-times we explicitly tested (which, we stress, are phenomenological in nature), the Teukolsky equation used to calculate the greybody factors turned out to be numerically very complicated, which is why we opted for proposing the equally phenomenological Hayward-like BH discussed here.] For phenomenological reasons, especially in light of our expectation (verified a posteriori) that such a behaviour would lead to weaker bounds on f_pbh and thereby the asteroid window further opening, it could therefore be desirable to also contemplate a case where the BH temperature actually increases with increasing regularizing parameter.

With the previous considerations in mind, here we introduce a new RBH solution metric which somewhat resembles the Hayward BH. This Hayward-like RBH is characterized by the following metric function:


    f_H-l(r)=1 - 2 M r^2/r^3+ℓ ( 1-ℓ r  )  .


To the best of our knowledge, this Hayward-like RBH has never been proposed before. Just as the Hayward BH, it is evident that this space-time possesses a dS core, with an effective cosmological constant given by Λ=6M/ℓ. We explicitly verify the finiteness of the curvature invariants, which take the following very simple forms:


    R   =   g^μνR_μν=24 M/ℓ ,    R_μνR^μν=144 M^2/ℓ^2 ,
     K   =    R_μνσρR^μνσρ=96 M^2/ℓ^2


We later explicitly verify that the intensity of Hawking radiation of the Hayward-like RBH increases with increasing regularizing parameter, therefore tightening constraints on f_pbh, and further closing the asteroid mass window. This behavior is similar to that of Kerr BHs, which display an enhancement of the primary photon emission spectrum with increasing spin parameter <cit.>. A similar behaviour has also been observed in the context of the Kazakov-Solodukhin BH <cit.> in Ref. <cit.>, consistently with the expectation laid out in footnote 3 (note, in fact, that the shadow of the Kazakov-Solodukhin BH increases with increasing regularizing parameter, see e.g. Fig. 1 of Ref. <cit.>, Fig. 18 of Ref. <cit.>, and Fig. 9 of Ref. <cit.>). Finally, we have checked that the allowed limit for the regularizing parameter in terms of the horizon radius is 0≤ℓ≤ r_H.



§ METHODOLOGY





 §.§ 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 <cit.>. 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 from 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 <cit.>.

For the static, spherically and tr-symmetric metrics given by Eq. (<ref>) 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 <cit.>, the Teukolsky equation governing the evolution of (massless) perturbations of different spin can be condensed into a single master equation <cit.>:


    [ - r^2/f∂_t^2 + s  ( r^2 f'/f -2 r  ) ∂_t  ] Υ_s 
       +  [ (s+1) (r^2 f' + 2 r f) ∂_r  ] Υ_s 
       + [ 1/sinθ∂_θ (sinθ∂_θ) + 2 i s θ/sinθ∂_ϕ. 
       . + 1/sin^2θ∂^2_ϕ -s -s^2 ^2θ ] Υ_s 
       + [ s r^2 f” + 4 s r f' + 2 s f  ] Υ_s=0 .


Here, Υ_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. (<ref>) is separable if one makes the following wave ansatz:


    Υ_s= ∑_l,m e^-i ω t  e^i m ϕ S^l_s(θ) R_s(r) ,


where ω is the perturbation frequency, l is the angular node number, and m is the azimuthal node number.

The functions S^l_s(θ) contribute to defining the so-called spin-weighted spherical harmonics S^s_l,m(θ, ϕ)=∑ S^l_s(θ) e^imϕ, satisfying the following equation <cit.>:


    ( 1/sinθ∂_θ(sinθ ∂_θ)+^2θ ∂_ϕ^2 . 
       . + 2isθ/sinθ∂_ϕ+s-s^2^2θ+λ_l^s  ) S_l,m^s=0 ,


where λ_l^s≡ l(l+1)-s(s+1) is the separation constant. For the spin 0 case, these functions reduces to the usual spherical harmonics S_l,m^0=Y_l,m.

Analogously to the Schwarzschild and Kerr BH cases <cit.>, the decoupled radial Teukolsky equation takes the following form <cit.>:


    1/Δ^s(Δ^s+1R'_s)' 
       +(ω^2r^2/f+2iω sr-isω r^2f'/f+s(Δ”-2)-λ_l^s)R_s=0 ,


where Δ(r)≡ r^2f(r) and ' ≡∂_r. We set in purely ingoing boundary conditions, so the asymptotic solutions of Eq. (<ref>) are given by:


    R_s ∼ R^in_s e^-iω r^⋆/r+ R^out_s e^iω r^⋆/r^2s+1   (r→∞) 
        R_s ∼ R^hor_s Δ^-s e^-i ω r^⋆           (r → r_H) ,


where r^⋆ is the tortoise coordinate defined by:


    dr^⋆/dr=1/f(r) .


We note that r^⋆→ r for large values of r, given that the metrics we consider are asymptotically flat.

In 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 <ref> for further details), which has already been successfully applied to these types of calculations in several earlier works <cit.>.

To begin with, we rewrite Eq. (<ref>) in terms of the rescaled coordinate x, given by the following:


    x ≡r-r_H/r_H ,


where r_H is the largest real root of f(r)=0. With this substitution Eq. (<ref>) is rewritten as follows:


    x^2(x+1)^3 f R̈_̈s̈
       + (s+1) x(x+1)  ( 2(x+1)f+(x+1)^2 ḟ ) Ṙ_̇ṡ
       + V(ω,x)R_s=0 ,


where ≡∂_x, and V(ω,x) is given by:


    V(ω,x) =  ( ω^2 r_H^2 (x+1)^2/f + 2 i s (x+1) ω. 
        - i s r_H (x+1)^2 ḟ/fω + s  ( 2 f + 4 (x+1) ḟ.
       . + (x+1)^2 f̈  -2  ) -l(l+1) + s(s+1) ) x (x+1) .


In 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 four RBHs under consideration are given by the following:

    f_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_GCSV(x)=1- e^ℓ - ℓ/x+1/x+1 ,
        f_H-l(x)=1-(1+ℓ-ℓ^2)(x+1)^2/ ( (x+1)^3 + ℓ - ℓ^2 (x+1)  )  ,

for the Bardeen, Hayward, GCSV, and Hayward-like space-times respectively.

Setting purely ingoing boundary in proximity of the horizon, the solutions to Eq. (<ref>) can be expressed in the form of a Taylor expansion as follows <cit.>:


    R_s(x)= x^-s- i ω/τ∑_n=0^∞ a_n x^n ,


where 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 <ref> for further details. 
The a_n coefficients can be determined by substituting Eq. (<ref>) in Eq. (<ref>) 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:


    R(x)  R^in_s e^-i ω x/x+R^out_s e^i ω x/x^2s+1 .


The GBFs can then be computed from the _s R^l m_in(ω) coefficient. More specifically, the normalization of the scattering problem is set by requiring a_0=1. With this normalization, the GBFs then read:


    Γ^s_l m(ω)=δ_s | _s R^l m_in(ω)|^-2 ,


where the coefficient δ_s is given by:


    δ_s = τ i e^i π s (2 ω)^2s-1Γ ( 1-s- 2 i ω/τ ) /Γ ( s-2 i ω/τ )


Using the method discussed above, we compute the GBFs for perturbations of different spin on the backgrounds of the four RBH space-times 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, as we will discuss shortly.



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

The 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 <cit.>: [This expression implicitly assumes that the particles emitted by the BH are not coupled to the regularizing parameter ℓ, an assumption which is reasonable.]


    d^2N_i/dtdE_i=1/2π∑_l,mn_iΓ^s_l,m(ω)/ e^ω/T± 1 ,


where 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 the plus (minus) sign in the denominator is associated to fermions (bosons). Following the methodology discussed in Sec. <ref>, 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.









We show examples of the resulting evaporation spectra in Figs. <ref>, <ref>, <ref>, and <ref>. The spectra obviously depend on the mass of the evaporating PBH, which we have set to M_pbh=10^13 kg, as it sits roughly in the middle of the 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. For the Bardeen, Hayward, and GCSV RBHs we observe that an increase in the regularizing parameter ℓ leads to a decrease in the intensity of the spectra at all energies. These behaviours are consistent with the temperature evolution shown in Fig. <ref>, although we stress that studying the temperature alone is not sufficient to draw these conclusions as the GBFs also play a key role in determining the shape and intensity of the resulting spectra, as Eq. (<ref>) makes very clear. For the Bardeen and Hayward RBHs 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 GCSV BH an increase in the regularizing parameter pushes the peak towards 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 are loosened for these classes of primordial RBHs.

Contrarily to what we observe for these three RBHs, we see that the intensity of the spectra for Hayward-like RBHs substantially increase with increasing ℓ for sufficiently large energies E ≳ O(MeV), whereas for lower energies they very slightly decrease. This behaviour agrees qualitatively with the expectations which motivated us to construct the Hayward-like RBH in first place, and with the temperature evolution shown in Fig. <ref>. In addition, we observe that an increase in the regularizing parameter leads to a much richer structure in the spectra: the observed multi-peak structure is the result of the contribution of different l modes in the emission being disentangled more clearly. A similar behaviour is actually observed in the  Kerr BHs evaporating spectra which, by virtue of their rotation, lead to a more complex spectrum where the contribution of high-l modes, while still subdominant with respect to the l=1 one, emerges more clearly (although we have explicitly checked that adding l>4 modes has virtually no effect on our results). In the case of Hayward-like RBHs, we remain agnostic as to the physical reason why such a structure is observed, and leave a more detailed study to follow-up work.



 §.§ Evaporation constraints



The spectra calculated in Sec. <ref> 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. <cit.> for a recent review). Evaporation constraints are the dominant ones in the 10^10 kg≲ M_pbh≲ 10^15 kg 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^12 kg: these have yet to fully evaporate today and, having formed deep during the radiation domination era, are therefore excellent non-baryonic DM candidates.

We 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. We also 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 <cit.>. 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 to examine how the limits on f_pbh change when moving from the Schwarzschild PBH framework to that of the regular metrics presented in Sec. <ref>, potentially opening or closing 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 contraints on f_pbh relative to the Schwarzschild BHs and the RBHs discussed in Sec. <ref>, therefore leading to negligible effects on the shift δ f_pbh, in which we are ultimately interested. In any case, such approximations also allow for a more direct comparison to several previous works and therefore we consider them appropriate for our pilot study, but their impact should definitely be explored in 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).

In what follows, we therefore assume that PBHs all have the same initial mass M_pbh. Following Ref. <cit.> we approximate the number of emitted photons in the logarithmic energy bin Δ E_γ≃ E_γ as being given by Ṅ_γ(E_γ) ≃ E_γ(dṄ_γ/dE_γ). The emission rate of photons from Hawking evaporation per volume at a cosmological time t is then given by <cit.>:[In the original paper M_pbh is assumed to be a function of time, due to the evaporation process. However, for the PBH masses considered in the present work (M > 10^12kg) we can safely assume M to be roughly constant during the evaporation process <cit.>.]


    dn_γ/dt(E_γ,t) ≃ n_pbh(t)E_γd^2N_γ/dt dE_γ(M_pbh,E_γ) ,


where 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:


    n_γ 0(E_γ 0)
       = n_pbh(t_0) E_γ 0∫^t_0_t_⋆ dt(1 + z) d^2N_γ/dt dE_γ(M_pbh,(1+z)E_γ 0)
        = n_pbh(t_0) E_γ 0∫^z_⋆_0dz/H(z)d^2N_γ/dt dE_γ(M_pbh,(1+z)E_γ 0)


where 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:


    I(E_γ 0) ≡c/4πn_γ 0(E_γ 0) .


It is this quantity which can then be directly compared against observations.

We assume a spatially flat ΛCDM cosmological model in specifying the expansion rate entering into Eq. (<ref>), with the same cosmological parameters as in Ref. <cit.>. This allows us to robustly cross-check our Schwarzschild constraints on f_pbh against those reported in the seminal Ref. <cit.>, 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. (<ref>) 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. (<ref>,<ref>,<ref>) we can compute the unnormalized photon flux I(E_γ 0)/n_pbh(E_γ 0), and adjust the normalization n_pbh(E_γ 0) by comparing against EGRB observations (as we will explain shortly). This procedures gives us an upper limit on n_pbh(E_γ 0), which can be translated into an upper limit on f_pbh as follows:


    f_pbh(M_pbh) ≡Ω_pbh/Ω_dm = n_pbh(t_0)M_pbh/ρ_crit,0Ω_dm ,


where ρ_crit,0=3H_0^2/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.

We 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 <cit.>, the COMPTEL imaging Compton γ-ray telescope in the 0.8-30 MeV range <cit.>, and the EGRET γ-ray telescope <cit.>. 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 <cit.>, positron flux <cit.>, 0.511 MeV annihilation radiation <cit.>, and various other sources. While these galactic observations could give 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 <cit.>, 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.

To set upper limits on n_pbh(t_0) – and therefore f_pbh through Eq. (<ref>) – we adopt the simple method first explained in the seminal Ref. <cit.>, and later adopted in most of the works examining constraints on PBHs from the 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. <ref> for a Bardeen PRBH 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 for all BHs, as well as 0.6 and 0.75 only for the Hayward-like BH), comparing the results to the Schwarzschild case which is recovered when ℓ=0. [Clearly, from the statistical point of view, more precise analyses are possible. For instance, one could construct a metric to be minimized (χ^2 or similar), or adopt a fully-fledged Bayesian approach exploring the joint M-f_pbh-ℓ posterior. Nevertheless, we believe our approach is sufficient for the purposes of our work, for several reasons. First and foremost, as with the approximations discussed earlier, adopting this method allows for a more direct comparison to several previous works. Furthermore, for most of these older datasets, often only the datapoints shown in Fig. <ref> are available, with no further available details on aspects which would be required to properly build a χ^2 or likelihood (e.g. correlation between the datapoints, instrumental details, and so on). Finally, we expect that the relative shift in f_pbh limits with respect to their Schwarzschild counterparts will be largely unaffected by the adopted methodology. For all these reasons, and especially being ours a pilot study, we believe the adopted methodology is appropriate for the purposes of our work.] We note that the exact origin of the EGRB is currently a matter of debate <cit.>: although it is 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 cannot fully account for the latter).





§ RESULTS











For each of the PRBHs discussed in Sec. <ref>, 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. <ref>. The results are shown in Figs. <ref>,  <ref>,  <ref>, and <ref> for the Bardeen, Hayward, GCSV, and Hayward-like BHs 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 literature. As a sanity check, we have verified that our ℓ=0 constraints exactly recover those of the seminal Ref. <cit.>. 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 (modified – either enlarged or contracted) asteroid mass window.

We begin by discussing the Bardeen, Hayward, and GCSV PRBHs, for which we saw earlier that the temperature and photon spectra decrease in intensity with increasing regularizing parameter ℓ (see the discussion in Sec. <ref>, and Figs. <ref>–<ref>). 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.45r_H) this behaviour is somewhat enhanced compared to the near-extremal Bardeen and GCSV 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. <ref> to Figs. <ref> and <ref>. This could also have been expected from Fig. <ref>, noting that the temperature of Hayward BHs decreases more rapidly with increasing ℓ/r_H relative the Bardeen and especially GCSV ones. Although the temperature is not the only factor at play in determining the resulting evaporation spectrum, given that the GBFs also play a key role as per Eq. (<ref>), it is reassuring that the temperature behaviour observed in Fig. <ref> is qualitatively reflected in the limits on f_pbh we derive.

For the Hayward-like BH, we observe exactly the opposite trend, with overall tighter constraints on f_pbh relative to the Schwarzschild ones. Again, such a behaviour is in line with expectations given the behaviour of the temperature observed in Fig. <ref>, and of the spectra as shown in Fig. <ref>. The relative shift in f_pbh constraints for a given value of M_pbh and ℓ/r_H is somewhat less dramatic than what we observed for the Bardeen, Hayward, and GCSV BHs, with shifts of less than two orders of magnitude.

As 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 four metrics considered. We recall that in the Schwarzschild case, the lower edge of the window lies at M_pbh≃ 10^14 kg. For the Bardeen, Hayward, and GCSV PRBHs, the looser constraints on f_pbh result in the asteroid mass window further opening up by approximately 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^13 kg. The opposite behaviour is of course observed for Hayward-like PRBHs, in which case the asteroid mass window further closes down, although by a lesser extent with respect to the previous cases. For instance, for ℓ/r_H=0.75, the lower edge increases to M_pbh≃ 2× 10^14 kg. Overall, we therefore observe that considering primordial regular BHs in place of the standard Schwarzschild ones can move the resulting constraints on f_pbh in either direction, further opening up or closing down the asteroid mass window, with the allowed region for the window lower edge spanning over a decade in mass, at least for the PRBHs considered.

Three 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. <ref>–<ref>). The reason is simply that, as ℓ is increased, the datapoint shown in Fig. <ref> 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, especially in the Hayward-like case where we saw that an increase in ℓ lead to a richer peak structure due to the effect of different l modes).

Next, 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 Reissner-Nordström 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 all BHs. In the language of Ref. <cit.>, the regularizing parameter for all four RBHs considered is a “specific hair” rather than an “universal hair” (see Ref. <cit.> 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.

Our 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. <cit.>, 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^30≲ M_pbh/kg≲ 10^37 – see Fig. 7 of Ref. <cit.> and Fig. 10 of Ref. <cit.>). 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^31≲ M_pbh/kg≲ 10^52 – see Fig. 7 of Ref. <cit.> and Fig. 10 of Ref. <cit.>), and considerations completely analogous to those we made for accretion constraints hold. [Potential exceptions to this mass range for dynamical constraints are those from capture of PBHs by white dwarfs or neutron stars at the centres of globular clusters <cit.>, or from supernovae explosions resulting from transit of a PBH through a white dwarf <cit.>. However, these limits are highly disputed because of uncertainties in the dark matter density in globular clusters <cit.>, or based on the results of hydrodynamical simulations <cit.>. For these reason, we will not consider the previously mentioned limits in our discussion.]

Of 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^12 kg. 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 distances, 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 is only 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 GCSV PRBH cases) genuinely corresponds to an enlarged asteroid mass window, and conversely for the Hayward-like case where the lower edge moves towards higher masses. Therefore, the window where Bardeen, Hayward, and GCSV PRBHs could account for all the DM is larger compared to that of Schwarzschild PBHs.

Other 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. <cit.> and Fig. 10 of Ref. <cit.>), 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 discussed earlier) is a worthwhile endeavour which we plan to explore in upcoming works.



§ 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^14 kg≲ M_pbh≲ 10^20 kg) 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 four so-called tr-symmetric metrics (including the well-known Bardeen and Hayward space-times), whereas our companion paper <cit.> considers non-tr-symmetric metrics, including various metrics inspired from loop quantum gravity.

We show that evaporation constraints on f_pbh, the fraction of DM in the form of PRBHs, can be substantially altered in either direction when moving away from the Schwarzschild picture, leading to the asteroid mass window further opening up or closing down depending on the direction of the shifts in f_pbh limits. For three of the PRBHs we considered (the Bardeen, Hayward, and GCSV space-times) the lower edge of the asteroid mass window is shifted by nearly a decade in mass, 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 <cit.>. The opposite trend takes place with the Hayward-like BH we constructed, and we have argued that part of this different behaviour can be traced back to the evolution (increase or decrease) of the PRBH temperature as the regularizing parameter is increased. On the other hand, 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 remark that the present work (alongside our companion paper <cit.>) 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 <cit.> goes partially in this direction), including potentially metrics which are coupled to the cosmological expansion. Last but definitely not least, if PBHs are truly regular, one would hope to ascertain this via signatures complementary to those we have studied (for instance through gravitational wave signatures, VLBI imaging, particle motion, or energy extraction). We plan to address these and other related points in follow-up work.


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






§ DETAILS ON THE COMPUTATION OF GBFS



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. (<ref>), in the form of a Taylor expansion as given by Eq. (<ref>). 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

    u” +p(x) u' + q(x) u=0 ,


in proximity of its singular points, namely those where p(x) and q(x) diverge. One can notice that Eq. (<ref>) can be rewritten in the form of Eq. (<ref>), with one of its singular point being at x=0, i.e. at the event horizon.

To solve the radial Teukolsky equation we therefore proceed as follows:


    
  * We work in units of the event horizon and rewrite Eq. (<ref>) in order to remove the denominators
        
    A(x)R”_s + B(x)R'_s + C(x) R_s=0,

    where the functions A(x), B(x), and C(x) are given by the following:
        
    A(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 i s (x+1) ω f -i s (x+1)^2 ω  f' 
                    +s f (    (x+1)^2 f”+4 (x+1) f'+2 f-2      )
                    -l (l+1) f +s (s+1) f   ,

    
  * The lowest power term around x=0 of each coefficient can be written in the following form:
    
    A(x)∼ x^2 τ^2   ,
         B(x)∼ x (s+1) τ^2    ,
         C(x) ∼ω^2 - i ω s τ  ,

        where τ=τ(ℓ) depends on the choice of RBH.
    
  * We then build the following characteristic equation:
        
    m(m-1) τ^2  + m (s+1) τ^2 + ω (ω - i s τ)=0 ,

    whose solutions are the following:
    
    
    m_1=-s - i ω/τ   ,     m_2=i ω/τ

    

  * It is then possible to conclude that Eq. (<ref>) admits solutions near the singular point x=0 of the form given by Eq. (<ref>).


Explicitly, for the four RBHs in question, τ is given by the following:

        
    τ_B=1-2 ℓ^2/ℓ^2+1  ,
       τ_H= (1-3  ℓ^2)  ,
       τ_GCSV=1-ℓ  ,
       τ_H-l= (ℓ-1)^2/1-ℓ(ℓ-1)  .


We notice that in the Schwarzschild limit ℓ→ 0, all of the above reduce to τ=1 as one could expect.



