#### Abstract

We present a graviatom with de Sitter interior as a new candidate to atomic dark matter generically related to a vacuum dark energy through its de Sitter vacuum interior. It is a gravitationally bound quantum system consisting of a nucleus represented by a regular primordial black hole (RPBH), its remnant or gravitational vacuum soliton G-lump, and a charged particle. We estimate probability of formation of RPBHs and G-lumps in the early Universe and evaluate energy spectrum and electromagnetic radiation of graviatom which can in principle bear information about a fundamental symmetry scale responsible for de Sitter interior and serve as its observational signatures.

#### 1. Introduction

Nonluminous atomic dark matter includes a wide range of candidates starting from the historically first mirror dark matter [1–3]. Presented here is graviatom, a gravitationally bound quantum system consisting of a regular black hole or gravitational soliton G-lump and a captured charged particle, which can be formally classified as atomic dark matter with an additional intrinsic dark feature presented by dark energy interior of a certain fundamental scale of de Sitter vacuum. Electromagnetic radiation of a graviatom bearing information on this scale fits in the ultrahigh gamma range and can be regarded as its independent observational signature.

The idea of an atom with a black hole as a nucleus goes back to 1971 when Hawking put forward the idea that primordial charged black holes could capture free charged particle forming neutral and nonrelativistic ultraheavy *black hole atoms* [4]. The minimal mass of a primordial black hole that would not have evaporated entirely by now is g [5]. The term *gravitational atom* was introduced in [6] for gravitationally bound neutral black hole and a charged particle. Similar quantum system called a *hypothesized gravitational atom* was considered in [7] by solving the Schrödinger equation with the Newtonian gravitational potential of a point-like mass modified by the presence of extra compactified dimensions and a captured neutron [7].

Quantum levels in the field of a black hole have been studied in detail starting from 1987 [8–10]. The gravitational analogue of the hydrogen atom orbitals was considered in [11] where the spectrum of normalizable fermion bound states was calculated in the Schwarzschild background. Perturbative modification of electronic orbitals around charged black holes was proposed in [12]. It is clear that a particle undergoing transitions between quantum levels should radiate. However, this question has been much less studied in the literature than the question of quantum spectrum of a black hole.

Quantum radiation of a charged nonrelativistic particle in the Schwarzschild field was first studied in 1990 in [13] with taking into account the DeWitt conservative self-force whose origin is interaction of an electrostatic field of a charged particle with a gravitational field of a source [14]. It was shown that in this case the DeWitt self-force always exceeds the nonconservative radiation damping force; the energy spectrum and intensity of the electric dipole radiation were calculated for the case of black hole with masses ranging from 10^{−5} g to 10^{14} g [13] (see also [15]). The conditions of existence of a gravitationally bound quantum system with a neutral black hole as a nucleus were formulated in [16, 17] where it was called *graviatom*.

In this paper we study graviatoms with regular black holes with de Sitter interior [18], their remnants [19, 20], and vacuum gravitational solitons G-lumps [19, 21] as nuclei.

Primordial black hole remnants left after the Hawking evaporation have been considered as a source of dark matter for more than two decades [22–29] (for a review see [30–32]). They could also be a source of additional information about the primordial power spectrum [29, 33, 34] and about their effect on big bang nucleosynthesis [29, 35]. However, the question whether some remnants leave after the Hawking evaporation of a singular black hole still remains open. The generalized uncertainty principle requires existence of a black hole remnant as a Planck size black hole [25]. On the other hand, no evident symmetry or quantum number was found which would prevent a complete evaporation [36]. Another question in the case of a singular black hole is how to evaporate a singularity? The complete evaporation would involve a serious and unclear change in spacetime structure which would have to evolve from a singularity to a maximally symmetric Minkowski (or de Sitter) space [37].

These problems do not arise in the case of a regular black hole with de Sitter centre instead of a singularity which leaves behind a thermodynamically stable double-horizon remnant [20] and does not involve a dramatic change in spacetime structure [37].

The idea of de Sitter interior goes back to the heuristic hypothesis that de Sitter vacuum could be a final state in a collapse [38], to the idea of a self-regulatory behavior of space-time geometry at achieving the Planckian densities [39], and to the idea of the existence of the limiting curvature of the Planck scale [40]. Second-order phase transition during evaporation of a black hole with de Sitter interior [19, 20] suggests possibility for a phase transition to de Sitter vacuum in the origin [41].

Arguments in favour of a regular black hole are provided by a loop quantum gravity [42, 43]. The “renormalization group improving” approach based on the running Newton constant, applied to the Schwarzschild spacetime, predicts an appearance of a smooth de Sitter core replacing a singularity [44]. The noncommutative geometry approach (for a review see [45]) applied to the Schwarzschild black hole leads to a regular de Sitter core at short distances from the origin [46]. Appearance of de Sitter core was found also for a cosmological noncommutative black hole of positive mass with the Gaussian density profile [47].

The Einstein equations admit the class of regular spherically symmetric solutions asymptotically de Sitter as with , of a certain fundamental scale and at infinity [21, 48, 49]. A source term for this class is specified by and provides a unified model-independent description of dark ingredients in the Universe by a vacuum dark fluid [50, 51]. In this approach a vacuum dark energy is described by a time evolving and spatially inhomogeneous cosmological term [48], while dark matter is represented by compact objects generically related to vacuum dark energy through de Sitter vacuum interior. They include regular black holes [18, 52, 53], their remnants [20], and gravitationally bound vacuum structures without black hole horizon (dark particles or dark stars, dependently on a mass) [19, 21], called G-lumps, since they hold themselves together by their own gravity [21]. For this class a static metric is described by the line element: In the asymptotically flat case to be studied in this paper, the metric function reads

A metric function is asymptotically Schwarzschild at , where is the characteristic length of de Sitter-Schwarzschild spacetime [18, 19], , and is characteristic scale for the interior de Sitter vacuum related to its density by . (Characteristic length scales essential for graviatom are summarized in the Appendix.) The scale comes as the radius of the direct matching of the de Sitter and Schwarzschild metrics, ; that is, [39]. It is also the scale at which the scalar curvature changes sign [19]. The mass is related to interior de Sitter vacuum and to breaking of spacetime symmetry from the de Sitter group at [21, 54, 55].

Regular black hole with de Sitter centre has two horizons [19, 21] which coalesce in the course of evaporation at a certain mass corresponding to a thermodynamically stable double-horizon remnant [19–21, 46, 56, 57]. For , the metric (1)-(2) describe G-lumps. The typical behaviour of the metric function [21] is shown in Figure 1.

Primordial black holes can be formed by various mechanisms (for a review see [35, 58]). Most general possibility involves primordial density inhomogeneities forming overdense regions which can stop expanding and collapse [4, 59]. The particular origin of inhomogeneities can be quantum fluctuations arising in various inflationary scenarios [35]. PBH formation would be enhanced if inflation is followed by a dust-like stage dominated by nonrelativistic particles [58]. This point is of special interest in the context of graviatoms. As we will see below, the case of graviatom with RPBH remnant or G-lump and captured GUT particles, − GeV, which are considered as superheavy dark matter, is the most promising [60–64].

In the literature two mechanisms of producing such particles have been considered. One is gravitational production towards the end of inflation [60–62, 65, 66]. The other is decay of the inflationary de Sitter vacuum during symmetry-breaking phase transitions at the GUT scale: statistical mechanics approach with model-independent description of cosmological background shows that appearance of GUT particles from the GUT scale de Sitter vacuum looks like evaporation of Bose-condensate: particles satisfy the equation of state of an ideal quantum degenerate Bose gas with the Gibbons-Hawking temperature related to the de Sitter Horizon [67–71]. GUT particles could be captured by RPBH remnants or G-lumps and survive to the present epoch as constituents of graviatoms.

The main goal of this paper is to study typical features of graviatoms radiation which could result in their observational signatures as dark matter candidates.

In Section 2 we estimate the probability of formation of a RPBH and G-lump in a quantum collapse of a primordial fluctuation. In Section 3 we introduce the Schrödinger equation for graviatom and study typical behaviour of its potential. Section 4 is devoted to quantum spectra and mechanisms of graviatom radiation. In Section 5 we summarize the results.

#### 2. Formation of the Regular Primordial Black Holes and G-Lumps with de Sitter Interior

In the classical approach geometry of a collapsing body in the comoving coordinates can be presented by the FRW spacetime with the positive curvature [72]. In quantum cosmology the wave function is defined on a superspace which is the space of all 3-dimensional geometries and matter field configuration , . It satisfies the Wheeler-DeWitt equation . In the case of the FRW cosmology with the only dynamical variable, the scale factor , the Wheeler-DeWitt equation reduces to the simple case called minisuperspace model. For the Friedmann closed model it reads [73] where In general case a collapsing quantum fluctuation can be described as a superposition of partial energy densities of possible matter contributions [74]: The parameter is related to the parameter in the equation of state , as . Contributions to density (5) include an inflationary vacuum (), ensembles of domain walls () and strings (), nonrelativistic matter (dust with ), ultrarelativistic gas (), perfect gas (), and ultrastiff matter ().

The Wheeler-DeWitt equation (3) reduces to where the potential expressed in terms of the dimensionless variable has the form [75] Energy of a fluctuation in (6) is Coefficients are related to the coefficients in (5) by .

To model a collapse to a compact object with the de Sitter interior, we modify the potential (7) by introducing a cutoff related to characteristic scale for an interior de Sitter vacuum in the terms of the negative powers of corresponding to high density contributions to (5), in accordance with the underlying hypothesis [38–40] that transition to interior de Sitter vacuum may occur at achieving a certain high density during contraction. Modified potential has the form

The potential (9) tends to minus infinity for and . Near , the potential (9) has a maximum , and the process of a quantum collapse involves tunnelling through the potential barrier.

For big values of the variable , in the limit , the Schrödinger equation (6) reduces to which has the solution [76]: with the asymptotics [75]: Choosing integration constant proper for a collapse, we get the WKB wave function corresponding to deflation. Indeed, in the limits and , the action , where , , and [75]. The classical solutions are (inflation) and (deflation).

For very close to , , the Schrödinger equation takes the form Its solution, satisfying the boundary condition , is given by

After tunnelling through the barrier, a collapsing object appears in the potential well. In this region, for the Schrödinger equation (6) reduces to [75] whose solution, satisfying the boundary condition , reads where is the degenerate hypergeometric function, , , , , and .

The energy spectrum is given by

The discrete spectrum does exist only if the following condition is satisfied:

This allows one to estimate the probability of tunnelling through the barrier and range of masses of a collapsing object. The tunnelling probability is determined by the WKB penetration factor: where , , and .

The factor is estimated as , which gives .

The Schrödinger equation (6) describes the planckeon on the energy level given by (17); . Its total energy includes the rest mass, , so the mass of an object is given by . Taking into account and , we estimate the probability of the tunnelling for the compact object with the mass : It seems for the first sight that smaller objects collapse with bigger probability. However, this estimate does not involve the scale for de Sitter interior . General constraint on the mass involving is obtained from (18). It gives where . Ultimately we get

As we see, possibilities of formation of compact objects with the de Sitter interior depend on the scales of the interior vacuum and of the inflationary vacuum . The scale of inflationary vacuum can be adopted as the GUT scale GeV. In the frame of the hypothesis of arising of interior de Sitter vacuum due to the phase transition at the GUT scale [41], constraint on the mass is g. In this case only regular primordial black holes can be formed in a collapse of primordial fluctuations. In the frame of hypothesis of self-regulation of geometry near the Planck scale [39] or existence of limiting curvature of the Planck scale [40], mass range for collapsing objects admits G-lumps and RPBH including those with masses sufficiently small to evaporate and produce remnants to the end of inflation. Objects with small masses are produced with the bigger probability (20).

One more possibility for a wide range of masses is related to production of compact objects with de Sitter interior during later stages of the early Universe. The standard model of particle physics predicts a phase transition at the QCD scale of 100–200 MeV (see [77] and references therein) which can lead to a second inflationary stage with duration of about 7–10 e-foldings [78–80]. In this case the formation constraint (21) admits any mass.

We conclude that regular primordial black holes, their remnants, and G-lumps can arise in the early Universe during first and second inflationary stages, so that they can capture available particles and form graviatoms.

#### 3. The Schrödinger Equation for Graviatom and Conditions of Its Existence

The problem of a motion of a nonrelativistic charged particle in a curved spherically symmetric spacetime reduces to the Schrödinger equation [81]: where is the mass of a particle. The potential includes contribution of gravity and of the DeWitt conservative self-force acting on a charged particle with the mass and the charge . It is given by [14] and directed outward the field centre. Here is a classical electromagnetic radius of a particle. In a weak gravitational field the conservative self-force always exceeds the nonconservative radiation damping force , where is the particle velocity [13]. In this approximation the gravitational part of the potential is given by . For a metric of the class (1) the gravitational potential reads .

This approach corresponds to the Pauli approximation to the Dirac equation [16, 17]. As we will see below, characteristic orbits of particles in the field of graviatom are in most cases located far from the characteristic Schwarzschild-de Sitter radius which justifies application of the Pauli approximation for analysis of basic properties of graviatom. In what follows we will consider the range of frequencies relevant for this approximation.

With taking into account the DeWitt self-force we have for the potential

The behaviour of the potential (24) is determined by typical behaviour of the metric function : in the asymptotically flat case it has not more than two zeros, one maximum at where and one minimum [21]. This follows directly from the Einstein equation for the metric (1) which reads

The pressures are given by For the solutions satisfying the weak energy condition (nonnegative density for an observer on a timelike curve), [21]. In the extremum of the metric function, , and the type of extremum is determined by the sign of the tangential pressure . In the asymptotically flat spacetime with de Sitter centre there is one scale of a vacuum energy, at , and at this point has a maximum. The next extremum is a minimum [21], in the region where is already positive (since and there). In an extremum of the potential (24) the condition gives It follows that has an extremum in the region where and the metric function is growing. In particular, the extremum of is always outside the double horizon on which . The second derivative is given by, with taking into account (25) and (27),

According to (27), at the extremum of , and also there; hence, the second derivative of the potential is always positive in its extremum and the potential can thus have only one minimum. For , the potential goes to infinity. For , the metric function goes asymptotically to the Schwarzschild metric, so that the potential for goes to and therefore as . The only minimum, , is thus always negative.

The radial wave function in (22) satisfies the equation The potential can be presented as

Equation (30) describes a captured particle in the effective potential:

In the Schwarzschild limit, In this limit we see clearly that in the case of a graviatom the fine structure constant is replaced with the effective gravitational fine structure constant [13]: The corresponding Bohr radius for the graviatom with the potential (33) is given by [16, 17] where is the Compton wavelength of a captured particle.

The condition (the Bohr orbit in the R-region outside the horizon) gives Equation (36) gives the constraint on the mass of the graviatom nucleus: where is the mass of the electron.

The main condition for the existence of a graviatom is the geometrical condition: Here is the characteristic size of a graviatom, is the size of its nucleus, and of a captured particle. For a hydrogen-like graviatom, it reads .

For the metric function (2), the black hole horizon, , satisfies since, for any density profile satisfying the weak energy condition , is monotonically growing, and . The restriction following from will also hold for to guarantee , so we can take . The condition (38) in this case, , with taking into account (35), gives which follows in the constraint:

Fulfillment of the conditions (39) guarantees fulfillment of the constraint (36).

In the case of G-lump, its characteristic size is of the order of [19, 21], and the condition of the existence of the hydrogen-like graviatom reads Introducing the dimensionless quantities and , we get the inequality . With the new variable , it reduces to . This is the 4th order polynomial; complex conjugate roots come in pairs and make this polynomial positive, so that it can be negative only if real roots exist. The condition of their existence, , results in the following constraint: This is the important constraint; only when (41) is satisfied, the geometric condition (40) holds. In this case the above polynomial has three negative and one positive root. The only possibility to satisfy (40) results in the second constraint: With these two constraints we estimate

For further analysis we adopt some slowly falling density profile which should vanish for large faster than to guarantee finiteness of the mass [21]. If we choose , we obtain The requirement for the Bohr orbit (35) to be outside the horizon, , results in where is the double horizon radius. It is calculated from , which gives The first inequality in (45) gives, with taking into account (31), the constraint where is the mass of the electron and is its Compton wavelength.

For the GUT scale interior GeV, this gives for the mass of the captured particle GeV and for the Planck scale interior .

The second inequality in (45) gives, with taking into account (46), the constraint It follows that the more the mass contained under the double horizon is, the better graviatom with a regular black hole can be approximated by the Schwarzschild potential (33). It suggests also that more quickly falling density profile would be more favorable for the existence of the Bohr orbits in the R-region outside the horizon.

For detailed estimates we adopt the density profile [18]: which describes vacuum polarization effects leading to de Sitter interior [39] in the semiclassical model for vacuum polarization in the gravitational field [19]. The metric function is given by A black hole exists for [19]. The boundary value corresponds to the mass g. With using (50) we get for the horizon

The most essential difference between a graviatom and the hydrogen atom is that, in the graviatom, due to the presence of the DeWitt force and modification of the Coulomb-like part of the potential, the effective potential has the minimum. This leads to an appearance of oscillatory spectrum and radiation from oscillatory levels.

With the density profile (49) the effective potential (32) takes the form The condition for a minimum of the function (52), , gives the relation where and Equation (53) can be written as

For we have and . For , , oscillatory levels are determined by the effective classical electromagnetic radius , and the second term in (54) would provide quantum corrections to oscillatory levels.

For we have and

This is the most promising case since quantum radiation which we present in the next section carries information about de Sitter interior of the graviatom nucleus.

In the general case the position of the minimum is determined by the interpolation formula:

The hydrogen-like radiation dominates when the effect of the DeWitt contribution is less that of the Coulomb-like term in (33), . Radiation from the oscillatory levels dominates in the opposite case, .

#### 4. Graviatom Radiation

##### 4.1. Hydrogen-Like Spectrum and Radiation

The constraint (38) leads, with taking (35) into account, to .

For the density profile (49) ; as a result for the graviatom with a regular black hole and its remnant, and spectrum of radiation can be approximated with the high accuracy by that for the Schwarzschild potential.

In this case the energy spectrum is given by [16, 17] where , , . Here and are the principal and orbital quantum numbers, respectively. This is general formula including both the Bohr levels and oscillatory levels.

The radiation from transitions between the Bohr levels dominates when . The condition of its domination, gives the constraint g.

In the case of a black hole, the requirement () results in the estimate:

In the limit (60) the energy spectrum reduces to

In this case we have the dipole electromagnetic radiation of graviatom (dominating for the Bohr levels) similar to that for a hydrogen-like atom. Frequencies for the transitions and are given by and intensities are related by [16, 17].

For the case of a black hole there is the natural stability condition [16, 17]: and where is the graviatom lifetime, is the black hole lifetime, and is the particle lifetime. The condition of indestructibility, , where is the destructive energy (due to tidal forces, the Hawking radiation and “ionization”-eventual removing of a charged captured particle) and is the binding energy, was thoroughly studied in [16, 17]. For the case of the hydrogen-like graviatom with the Schwarzschild potential (), conditions of stability and indestructibility constrain the gravitational fine structure constant within the narrow range [16, 17] . For this range of the coupling the applied Pauli approximation holds for the frequencies of the hydrogen-like radiation (63).

In the case of G-lump the estimate (43) restricts the correction to the potential for the density profile (49), However, we cannot approximate the radiation spectrum by that for the Schwarzschild case with this accuracy since the energy is proportional to and thus to .

In the case of G-lump we should have to solve a separate problem with the potential (32). Preliminary consideration can be made by replacing the mass in (62) with . The rough estimate with gives the correction of the order of . This estimate is relevant also for a remnant since its characteristic size, , is close to . This allows us to conclude that radiation from the Bohr levels in the case of G-lump and remnant will depend on the density of the interior de Sitter vacuum and would result in observational signatures for them as dark matter candidates.

##### 4.2. Oscillatory Spectrum and Radiation

In the case when the DeWitt force is dominating, and we deal with an oscillatory spectrum and radiation.

Near the minimum, the potential reduces to

Introducing in (30) with the effective potential (67) the new variables we obtain the equation of a harmonic oscillator in one dimension: The energy spectrum is given by The wave function being a solution to (69) has the form [82]: and satisfies the normalization condition: Matrix elements of the coordinate are The intensity of an electric dipole radiation has the form [83]: where is the frequency of a transition and are the matrix elements of the charge coordinate in the transition from the state to the state .

Using (73), we obtain, for and , the expression

For the graviatom with the de Sitter interior in the limit (66), there are two possibilities. For , the correction to the Coulomb part potential in (32) satisfies by virtue of . The graviatom can be described with the high accuracy by the Schwarzschild potential (33) which does not provide information about internal structure of the graviatom nucleus. The minimum of the potential (33) is located at . The series in gives the oscillatory spectrum: where .

For a black hole as a nucleus, the condition for the minimum outside of horizon, , leads to the constrain g. The condition (66) gives the restriction which results in for particles captured by a black hole. Only atomic nuclei can be captured with [13]. However, the captured nuclei with are quickly destroyed by the Hawking radiation from the black hole horizon. The condition for existence of such a graviatom, where is a binding energy, is not satisfied for [16, 17].

In the case of a G-lump or remnant as a graviatom nucleus, the condition (66) of capture of a particle on a level near the minimum leads to the constraint which excludes a nucleus with . This is not a big disaster since a spectrum in this case does not carry information about an internal structure. From this point of view, the most interesting is the case (56) when and thus frequency depend on the interior vacuum scale .

The second possibility is specified by For the case of the density profile (49) the minimum of the effective potential (52) in the limit (66) is given by (56), and the oscillatory spectrum the intensity of the dipole radiation and the frequency depend on and hence on .

The geometrical condition requires . It gives us the upper limit for the intensity
where and erg s^{−1}. In the case of G-lump with cm we have erg s^{−1}, and for a remnant and near extremal black hole is bigger. Formula (81) is valid for particles with the mass which follows from the restriction of applicability of nonrelativistic quantum mechanics.

The energy of radiated quanta exceeds the average energy of the quantum radiation from the Schwarzschild black hole and the average energy from the de Sitter horizon of the interior scale .

The ratio of the intensity of radiated quanta in our case to the Hawking radiation of the Schwarzschild black hole is estimated as The ratio of our intensity to intensity of the Hawking radiation of the de Sitter horizon is

##### 4.3. Availability of Information on Graviatom Interior

The most promising case is when the frequency of oscillatory transitions (80) depends on the scale of the interior de Sitter vacuum. It could be an observational signature of a graviatom with de Sitter interior as a dark matter candidate. In this case we can estimate frequency as Photons with frequencies of order GeV were first observed in ultra-high-energy cosmic rays in 1962 [84]. Current experiments allow detection of photons up to GeV (see [85] and references therein). Present observational possibilities prefer thus graviatoms with the GUT scale interior although those with the Planck scale interior can exist in principle, and probabilities of their production in a collapse are bigger, but their typical frequency, GeV, is far from the today observational range.

The most promising case can be realized when two conditions are satisfied: the condition (66) of dominance of the oscillatory levels and the condition (77) defining appropriate location of the minimum of the effective potential.

To tell something about a graviatom nucleus, we can take roughly and apply the requirement that the minimum of potential must be located in the R-region outside the event horizon. The requirement leads for the density profile (49) to . This means that the most favorite situation can take place for G-lump.

Precise estimate gives , , for a double-horizon remnant with the density profile (49). The minimum of the potential is outside of the double horizon. The Hawking temperature is zero, so that it does not disturb the oscillatory radiation of the graviatom. The case of almost extremal black hole is also possible, until , in such a case temperature is quickly going to zero, and oscillatory spectrum cannot be destroyed by the Hawking radiation.

Our predictions concerning captured particles are restricted here by the applied approach. The basic requirement of applicability of the Pauli approximation gives immediately in (84) the constraint on the mass of a captured particle GeV for the GUT scale of the interior de Sitter vacuum, so we can trust the above conclusions made in the frame of nonrelativistic quantum mechanics but cannot exclude that a more precise approach would open additional possibilities.

To evaluate which particle can be a captured, we apply the constraint . The condition of the proper location of the minimum (77) requires . This gives cm for the GUT scale of the interior de Sitter vacuum. What a particle it can be? From cm, we get which means GeV.

This range admits GUT particles with masses which could be captured by RPBH at the end of inflation. However, here the problem arises. Formation constraint restricts masses of RPBH with GUT scale interior by g. For such a black hole, even with very quickly falling density profile (49), proper for the frequency (84), minimum (56) is located deeply inside T-region, so a particle would not be captured on an orbit but would be swallowed by a black hole. G-lumps and remnants which do not encounter this problem can be formed at the first inflationary stage only with Planck scale interior. They can capture GUT particles, but their radiation would be then in the range unavailable for near future observations.

Question arises where regular black hole remnants with GUT scale interior produced in evaporation, as well as remnants and G-lumps produced at the second inflationary stage, can capture particles with masses of the order of the GUT mass.

Speaking about RPBH produced during first inflation, one can imagine a situation when it would capture a GUT particle on some higher oscillatory level. Simple estimate with the spectrum (78) shows that it is possible, but lifetime of such a level should have to be comparable to the evaporation time [35]. In our case it gives s.

Such a case could be realized if some metastable level(s) with appropriate quantum number(s) would exist in the oscillatory spectrum, which cannot be excluded *a priori*.

Good news is that the level with energy GeV would not be destroyed by the Hawking radiation. Spherically symmetric regular black hole with de Sitter interior evolves during evaporation to the double-horizon remnant with zero temperature, while its maximal temperature is given by GeV [19, 37, 41].

Another channel is related to leptoquarks which can survive in galactic halos [86]. The leptoquarks arising at the GUT epoch are decayed into quarks and leptons. The quarks form nucleons in three minutes after the Big Bang. However, a part of leptoquarks survives and is accumulated in the galactic halos contributing to a dark matter [86]. The leptoquarks in the galactic halos can be captured by any of the considered here objects with de Sitter interior: by near extreme primordial black holes, by remnants, and by G-lumps.

#### 5. Summary

General constraint on formation of a compact object with the de Sitter interior results in the following options: at the first inflationary stage at the GUT scale, regular black holes with masses g can be formed with the GUT scale interior and any objects with the Planck scale interior including G-lumps. There are no constraints on formation of RPBH and G-lumps during the second inflationary stage at the QCD scale.

Spectra and radiation of graviatoms are studied in the Pauli approximation to the Dirac equation. The gravitational coupling is much bigger than the fine structure constant . The energy levels in a graviatom are determined by , and relativistic correction which would allow for a particle spin is of the order of . The condition of the particle localization dimension exceeding the Compton wavelength is valid all over the paper. This means that the problem is solvable in the framework of nonrelativistic quantum mechanics, although it gives the restriction on the mass of a captured particle, .

The electromagnetic radiation of charged particles in the field of a regular black hole exceeds the Hawking radiation from Schwarzschild black hole of the same mass. The essential difference is that a graviatom emits spectral lines which can bear information on its de Sitter interior.

Hydrogen-like spectrum and radiation of graviatom with RPBH as a nucleus is described with the high accuracy by those for the Schwarzschild black hole. Requirements of stability of graviatom with respect to tidal forces and Hawking radiation constraint the gravitational coupling within the narrow range . Dominating is the dipole radiation which does not carry information about de Sitter interior of graviatom. In the case of G-lump and RPBH remnant, approximation of the potential by the corrected Schwarzschild potential suggests that hydrogen-like radiation will depend on the density of the interior de Sitter vacuum and result in observational signatures for G-lump and remnant as dark matter candidates. This case needs further detailed investigation.

The oscillatory radiation can depend essentially on the interior vacuum scale and gets in the range of ultrahigh cosmic ray energies.

Typical features of graviatoms radiation can result in observational signatures for G-lumps and RPBH remnants as dark matter candidates, as well as provide information about their interior de Sitter vacuum.

Let us note that regular primordial black holes, their remnants, and G-lumps can also serve as a very sensitive universal probe for a scale of inhomogeneities in the very early Universe being in this context the most elusive among dark matter candidates [87].

#### Appendix

#### Characteristic Length Scales of Graviatoms

Graviatoms have rich internal structure, characterized by several essential length scales which we summarize below. Characteristic size of a graviatom satisfies the geometrical condition where is the size of the nucleus and is the size of a captured particle.

The characteristic length of any de Sitter-Schwarzschild spacetime comes as the radius of the direct matching of the de Sitter and Schwarzschild metrics, . It depends on both gravitational radius and characteristic scale of the interior de Sitter vacuum related to its density by . For the GUT scale *≃*10^{15} GeV interior, cm. For the Planck scale, . Characteristic scale gives rough estimate for the size of a G-lump. More precise estimate involves critical gravitational radius , related to critical mass , actually mass of the remnant which gives the boundary between black holes () and G-lumps (). It is given by where the numerical coefficient depends on the model. For the density profile (49), and . In this case, for a G-lump, and for a black hole.

The black hole horizon satisfies always . Double horizon gives roughly a size of RBH remnant. Masses of remnants () range from for the Planck scale interior, through 10^{3} g for the GUT scale interior, up to 10^{32} g ( cm) for the QCD scale interior. Regular remnants are stable, and their population includes all remnants of RPBH evaporated till now.

The gravitational analogue of the Bohr radius is given by , where is the Compton wavelength of a captured particle. The hydrogen-like radiation dominates in the case when , where is electromagnetic radius of a captured particle.

In the case , the oscillatory radiation is dominating. Minimum of the graviatom potential , responsible for its existence, is always outside the double horizon. For the case , oscillatory levels are determined by , with quantum corrections depending on ; in this case, . In the most promising for the observational manifestations case, when , location of the minimum depends directly on and thus on the interior vacuum density, for the density profile (49). For a remnant with the density profile (49), , , , and .

#### Acknowledgments

This work was supported by the Polish Ministry of Science and Education for the research project “Globally Regular Configurations in General Relativity Including Classical and Quantum Cosmological Models, Black Holes and Particle-Like Structures (solitons)” in the frame of the “Polish-Russian Agreement for collaboration in the Field of Science and Technology” and by the Polish National Science Center through the Grant 5828/B/H03/2011/40.