Primordial Black Holes from Spatially Varying Cosmological Constant Induced by Field Fluctuations in Extra Dimensions

The origin and evolution of supermassive black holes (SMBHs) in our universe have sparked controversy. In this study, we explore the hypothesis that some of these black holes may have seeded from the direct collapse of dark energy domains with density significantly higher than the surrounding regions. The mechanism of the origin of such domains relies on the inflationary evolution of a scalar field acting in D dimensions, which is associated with the cosmological constant in our four-dimensional spacetime manifold. Inner space quantum fluctuations of the field during inflation are responsible for the spatial variations of the dark energy density in our space. This finding holds particular significance, especially considering recent evidence from pulsar timing array observations, which supports the existence of a stochastic gravitational wave background consisting of SMBH mergers.


Introduction
Primordial black holes (PBHs) have been extensively studied over the decades [1][2][3][4][5][6][7][8][9][10] and offer a scenario with the potential to leave distinct imprints on cosmic history.Depending on the ratio of their abundance relative to the overall dark matter (DM), f PBH = Ω PBH /Ω DM , the range of possible PBH masses M PBH spans a wide spectrum, including PBHs of small masses [9] which have undergone scrutiny through various observations (for comprehensive reviews, refer to [8,10]).Additionally, since PBHs formed during the early stages of the Universe, they have the capacity to develop bound binaries via multiple mechanisms [8,11,12].As these binaries become close, they emit gravitational waves (GWs) continuously until a final dramatic burst occurs at the point of their ultimate merger.Notably, for black holes of stellar mass, such mergers have already been detected by ground-based interferometers [13,14].Moreover, it is plausible that several of these observed mergers might be attributed to the coalescence of PBHs [15][16][17][18][19][20].
PBHs with masses exceeding 10 2 M ⊙ hold particular significance due to their impact on the growth of massive objects during the evolution of the early Universe.Notably, it is wellestablished that galactic nuclei host supermassive black holes (SMBHs) with masses surpassing 10 6 M ⊙ [21][22][23].It has been theorized that PBHs could be their progenitors, achieving such masses through processes like merging, accretion [24][25][26][27][28][29], or the direct collapse of primordial fluctuations [30,31].In the latter scenario, SPBHs are constrained to constitute less than O(0.1%) of dark matter (DM).As they have been present since the dawn of the matterdominated era, they can serve as cosmic seeds, enhancing galaxy formation [32,33].Furthermore, different observations have provided evidence for the existence of intermediate-mass black holes (IMBHs) with masses ranging from 10 4 M ⊙ to 10 6 M ⊙ [34].Additionally, a subdominant fraction of dark matter may consist of immensely massive PBHs, exceeding 10 12 M ⊙ [35], capable of traversing the intergalactic medium.
SMBHs may be responsible for the generation of early galaxies reported by JWST [36].They can bind in binary systems which leads to late time merging and radiation of gravitational waves in the nHz frequency range that are detectable by pulsar timing array (PTA) experiments [37][38][39][40][41][42][43].The results from the PTA observations have been extensively analyzed and interpreted in numerous studies, including recent ones such as  and earlier analyses such as [114][115][116][117][118][119][120][121][122][123][124][125][126], which are related to the previously published NANOGrav signal evidence [127].These interpretations and effects may independently explain the PTA gravitational wave signal.They can also be considered in combination with the modeling of gravitational waves originating from supermassive black hole binaries (SMBHBs).
In this paper, we propose and validate a mechanism for the formation of PBHs based on the generation of specially varying cosmological constants, which may be generic for theories with compact extra dimensions.Considering extra dimensions allows us to examine fluctuations of fields within the internal space during inflation in addition to fluctuations of ordinary scalar fields.However, the fate of these field fluctuations differs significantly from those associated with conventional four-dimensional scalar fields.While conventional field fluctuations transform their energy density into radiation during the FRW stage through decay into other particle-like species, the energy density associated with scalar fields within the internal space remains unchanged, effectively stored within the scalar field itself, manifesting as a cosmological constant.Fluctuations of the scalar field within the internal space manifest as spatial variations of the local Λ term.Domains containing an extraordinarily high cosmological constant may collapse into PBHs.After the end of inflation, the horizon expands and the particle energy density decreases, approaching its present-day value.Simultaneously, the energy density associated with the Λ term, being dependent on the Hubble parameter, also decreases over time, converging to its present-day value, which equals the dark energy density.There must exist a moment in time when both the energy density of matter and the energy density associated with the Λ term are equal.Evidently, since the value of Λ varies across space coordinates, this equality primarily arises within the densest regions, implying the existence of a density contrast close to unity in domains with high values of Λ.Once such a domain becomes encompassed by the cosmological horizon, it may evolve into a PBH.We assert that PBHs formed through the proposed mechanism are cosmologically feasible candidates for seeding SMBHs and explaining the observed IMBHs.
The flexible metrics characterizing extra dimensions constitute a continuous set of static classical solutions derived from the generalized Einstein igat [128,129], and they share fixed Lagrangian parameters.This approach, distinct from the brane world model, renders the extra dimensions invisible due to their small size.However, unlike Kaluza-Klein geometries, these dimensions exhibit inhomogeneity.The concept of such geometries was initially introduced in [130] with further discussion in [131].Subsequent research, as presented in [132], applied a top-down approach to elucidate observed physical laws.It demonstrated that incorporating quantum corrections to initial parameters established at high energies eases the renormalization procedure.
Investigating the evolution of extra field distributions leading to a static state is a crucial endeavor.This aspect has been explored in previous works such as [133].The outcomes reveal that the resultant metric and field distribution are contingent upon both model parameters and initial conditions.Notably, the extra-dimensional stationary field distributions evolve in tandem with the energy density across distinct volumes below the horizon, which are replicated during inflation.This particular aspect forms the focal point of our investigation.Furthermore, our current investigation is grounded in nonlinear f (R) gravity, as extensively discussed in reviews such as [134,135].This framework holds significant potential for diverse cosmological implications, with one notably remarkable consequence being the emergence of dark matter [136].Several viable f (R) models in 4D space aligning with observational constraints have been proposed in works such as [128,129,[137][138][139].
This paper is structured as follows: In Section 2, we provide a concise overview of the mathematical setup employed in the extra-dimensional framework under consideration.Section 3 is dedicated to exploring the distinct behaviors of fluctuations in our space compared to those in extra dimensions.Section 4 is focused on deriving the conditions essential for the formation of PBHs and estimating their mass spectrum.The conclusions of our study are summarized in Section 5.

Static Field Distribution in Internal Space
The primary objective of this section is to revisit the foundational concepts of extra-dimensional frameworks that give rise to a continuum set of static metric distributions.This issue has been elaborated in our previous papers [133], and we refer the reader to them for details.Consider f (R) gravity with a minimally coupled scalar field ζ in a D = 4 + n-dimensional where g D ≡ Detg M N ; M, N = 1, D; X A = (x µ , y a ); the coordinate set x µ , µ = 1, 2, 3, 4 describes the four-dimensional space M 4 , and the set y a , a = 5, 6, ..., n describes the n-dimensional manifold M n , which is assumed to be a closed manifold without boundary; f (R) is a function of the D-dimensional Ricci scalar R; and m D is the D-dimensional Planck mass.Below, we will work in the units m D = 1.Note that the main results of this work hold even for the simplest form of the potential V .
The metric is postulated to have the form Such a metric ansatz has been extensively studied within the realm of linear gravity [140][141][142][143], particularly in addressing the hierarchy problem [133,[144][145][146].Our approach is based on the concept of compact extra dimensions.A preliminary investigation suggests that their scale could be as small as 10 −28 cm or even smaller.This implies that the extra dimensions remain invisible to our instruments, and our rulers and clocks do not measure intervals of space and time at a specific value of u.Instead, all metric functions, such as the function e γ(u) , should be averaged over the extra space.The way to achieve this is discussed in [133] and briefly presented below.
The equations of motion, see [128,129], represented by possess a continuum set of solutions just as the differential equations do.We choose those solutions that exhibit homogeneity in the spatial coordinates x and inhomogeneity in the internal coordinates y.We consider only those solutions that refer to the compact extra space.
This means that the metric function r(u) must have two zeros.This condition is fulfilled at the coordinates u min and u max , i.e., r(u min ) = r(u max ) = 0, which is the result of numerical calculations.These coordinate values depend on additional conditions which are different in different space domains.
After integration over the extra-dimensional coordinates, the action (1) reduces to the ef- The term c ef f represents the cosmological constant Λ: assuming the scalar function ζ is homogeneous in a 3-dimensional space under horizon.This value varies in different space regions due to the fluctuations at the inflationary stage.We are interested in those space domains where the effective parameter Λ is considerably large as compared to values in the surrounding space.Here, g 4 is the determinant of the 4D metric: The effective parameters are expressed as follows is the volume of n − 1-dim sphere.
The parameter u is expressed in D-dimensional Planck units.
The right-hand side of Equation ( 8) is expressed in units where m D = 1.This relation allows us to articulate the D-dimensional Planck mass in terms of the four-dimensional Planck mass m P .In this context, we assume that the functions γ(u), r(u), ζ(u), R(u) constitute a specific solution to the system (4), with details available in [130] for a specific value of the Hubble parameter H. Figure 1 illustrates some examples of static distributions.This approximation remains valid during the inflationary period and at the present time, particularly when the Hubble parameter remains nearly constant.
Our comprehension of the specific value of the energy density, denoted as ρ Λ = Λm 2 P /(8π) = −c ef f m 2 P /(16π), is quite limited.Observational constraints provide an upper limit of approximately 10 −123 m 2 P at the present time.Understanding this density during inflation is even more uncertain, with the sanity bound being ρ Λ (H ≃ 10 −6 m P ) ≪ H 2 ≃ 10 −12 m 2 P , implying its negligible impact on the inflation rate.Post-inflation, considering the variation in the Hubble parameter becomes crucial, and establishing a connection between ρ Λ and this parameter remains elusive.Multiple factors contribute to the complexity of this issue, including quantum corrections, the influence of other fields, and the effects of averaging after the horizon crossing.Furthermore, obtaining an accurate solution to the dynamic equations during the reheating stage appears challenging.The subsequent section is dedicated to a detailed discussion of these aspects.
Clearly, both the extra dimensions and the scalar field experience fluctuations in the Ddimensional space.Quantum fluctuations during inflation are expected to induce significant deviations from their initial values (distributions), especially during the later stages of inflation when scales much smaller than those relevant for CMB observations exit the inflating Hubble patches.Specifically, fluctuations in the parameter c ef f during inflation can lead to spatial variations in the cosmological constant.The cosmological effects of these variations are the primary focus of investigation in the subsequent sections of the paper.These fluctuations may be substantial, giving rise to domains where the density of dark energy is significantly higher compared to the surrounding regions.Such domains could persist for an extended duration, provided gravity exerts a strong influence.

Inflationary Field Dynamics in Extra Space
The evolution of the Universe is significantly influenced by field fluctuations during inflation.
Following the completion of inflation, the rapid decrease of the Hubble parameter induces vigorous damped fluctuations of the field, ultimately converging asymptotically to one of the minima of its potential.These inhomogeneities, influenced by gravitational effects, give rise to a large-scale structure after the conclusion of the radiation-dominated stage.
We consider the presence of extra spatial dimensions, allowing for fluctuations in fields within this inner space.The destiny of these field fluctuations differs fundamentally from those inherent to the usual, four-dimensional scalar fields described above.This distinction arises from the fact that stationary distributions of fields constitute a set of measured continuum, as established in the early study [130].
Similar to the situation with usual fields, during the FRW epoch, the asymptotic distribution of the fields in the inner space undergoes variations in causally disconnected regions due to random fluctuations during inflation.However, a significant distinction arises: while the energy density of usual fields is transformed into radiation at the FRW stage through the decay of the fields into other particle-like species, the energy density associated with the scalar fields in the inner space remains in its initial form, being stored in the scalar field.Therefore, the evolution of the energy density of the fields in the inner space is still governed by the Hubble parameter, resulting in a slower decrease compared to the energy density stored in particle-like species generated from the decay of typical scalar fields, which could exist during the inflationary epoch.To distinguish between these two kinds of energy densities, we use the notation ρ for the energy stored in the ordinary fields fluctuating in the observable three-dimensional space, which is eventually converted into particle-like species, and ρ Λ for the energy remaining stored in the scalar field exhibiting inhomogeneities in the inner space.By choosing the model parameters such that ρ ≫ ρ Λ , we ensure that the impact of the field distribution in the inner space on the expansion rate, as well as on the rate of horizon growth, can be safely neglected.
After the end of inflation, the horizon expands and the particle energy density ρ decreases, approaching its present-day value ρ(t 0 ).Simultaneously, ρ Λ , being dependent on the Hubble parameter, also decreases over time, converging to the present-day value ρ Λ (t 0 ), which equals the dark energy density.Since ρ Λ (t 0 ) > ρ(t 0 ), there must exist a moment in time, denoted by t * , when both densities are equal, ρ(t * ) = ρ Λ (t * ).Evidently, since the value ρ Λ (t * ) varies across space coordinates, this equality primarily arises within the densest regions, implying a density contrast δρ/ρ ≃ 1.

Formation of PBHs Induced by Inhomogeneous Cosmological Constants
As indicated in the preceding section, our setup involves the total energy density, which comprises the 4D energy density ρ, represented either by a scalar field or particle-like species created after the conversion of this field into radiation at the reheating stage, and ρ Λ , associated with the energy density emerging from the scalar field distribution in the inner space.During the inflationary stage and for some period afterward, the Universe was dominated by the 4D energy density, such that ρ ≫ ρ Λ , while at the present time, ρ ≲ ρ Λ obs , where the current observable value of the dark energy density is given by ρ Λ obs ∼ 10 −123 m 4 Pl .Therefore, at some moment t * during the evolution of the Universe, both types of energy density become equal, resulting in If the equality condition (9) occurs within a causally connected domain, it implies that the density contrast, expressed in this particular case as Let us consider a scenario in which fluctuations of scalar fields in the internal space during inflation lead to the formation of a domain with size R(t end ), determined at the end of inflationary epoch t end , filled with a cosmological constant Λ 1 that exceeds its observable average value Λ obs .After the inflationary period, during the FRW epoch, the domain of size R(t end ) undergoes simple conformal stretching due to the expansion of the Universe where a(t) is the scale factor.It is evident that at a time t 1 ≳ t * , ensuring δρ ρ(t * ) > 1, the domain reaches a radius of R(t 1 ) as described by Equation (11), acquiring the mass M 1 .Subsequently, it becomes encompassed by a Hubble radius H −1 (t 1 ) = H 1 , thereby becoming detached from the cosmological expansion and initiating collapse.Within about a Hubble time, it will convert into a black hole (BH) of mass M PBH = ξM 1 .Below, we assume that almost the entirety of the energy contained in the domain is deposited into the BH, so that ξ ≃ 1.
Since the interior of such a domain can exert repulsive gravity due to its substantial energy dominance within the encompassing Hubble horizon, it may maintain a negative pressure, particularly if the density contrast δρ ρ(t 1 ) exceeds a threshold of around 10.In such a scenario, the domain enters the Hubble radius at t 1 ≫ t * and starts expanding faster than the background, eventually reaching the inflationary vacuum and potentially developing a wormhole to a baby universe.Such a wormhole would appear as a BH in the FRW Universe.In this paper, we focus on the regime with the most plausible collapse rather than expansion, where t 1 ≈ t * , and thus, the size of the domain filled with Λ 1 is close to the Hubble radius at the moment when local dominance of the cosmological term occurs, i.e., when the condition δρ ρ(t * ) ≈ 1 is reached.If the cosmological constant Λ 1 substantially exceeds its universe-averaged value Λ obs , then the component ρ Λ obs (t * ) can be neglected in (10).Therefore, the conditions for reaching a density contrast (10) exceeding unity can be described as where G stands for the Newtonian constant, H 1 denotes the Hubble rate at t 1 , and the Λ term is measured in units of the Planck mass squire.
Thus, an overdense object with a size given by is formed, with its mass as measured by a distant observer being expressed as This mass is determined by the localized value of the cosmological constant Λ 1 within a specific domain, which surpasses its universe-average value outside the domain.The validity of Equation ( 14) in the case of an overdense domain emerging due to fluctuations in extra dimensions is rigorously proven in Appendix A. Thus, it appears that the Schwarzschild radius of the above object, given by exceeds the size of the Λ 1 overdense domain given by Equation (13).Hence, assuming that the spherical shape of the domain is not significantly disturbed, it will be converted into a BH.To account for the population of unevaporated black holes, it is instructive to express the mass of such Λ-term-induced PBHs (ΛPBHs) in units of solar masses: If we consider that the growth of supermassive black holes (SMBHs) observed today originated from seed black holes, then this process must have commenced in the early Universe, approximately 3 million years after the Big Bang, with seeds heavier than 10 2 M ⊙ .Additionally, there is evidence of the existence of intermediate-mass black holes (IMBHs) with masses ranging from 10 2 M ⊙ up to about 10 6 M ⊙ .Thus, attributing the seeding objects with masses from 10 2 M ⊙ to 10 6 M ⊙ to ΛPBHs implies that they appeared as a result of the collapse of domains containing Λ-terms spanning the range In the subsequent analysis, we examine the comparability of the spectrum of the population of ΛPBHs with the constraints on the abundance of PBHs within the considered mass range.A domain of radius R ≈ H −1 inf , filled with Λ ̸ = Λ obs , that emerges at the time moment t Λ , during inflation with a total duration t inf , when the Universe is yet to inflate over ∆N Λ = H inf (t inf − t Λ ) = N inf − N Λ e-folds, undergoes stretching during expansion as The number of domains created in a comoving volume dV within an e-fold interval dN Λ is determined by where Γ Λ represents the formation rate of domains with Λ per Hubble time-space volume H −4 inf .By expressing N Λ from (18), we can derive the number distribution of domains with respect to their physical radius R as Therefore, the number density in the physical inflationary volume dV inf = e 3N Λ dV is In the context of the setup discussed in this section, where domains of high-density contrast and obeying the condition ( 12) are considered, the distribution ( 21) covers a range of scales from , where N Λ 1 represents the number of e-folds when the probability of the appearance of at least one domain with Λ 1 becomes significant.This probability becomes notable over the course of the progression of inflation, which lasts for a sufficient number of N inf e-folds necessary to address the horizon and flatness problems.
It is worth noting that if inflation were to occur above the TeV scale, the comoving Hubble scale at the end of inflation would be less than one astronomical unit.Consequently, a causally connected patch could encompass our entire observable Universe today, which has a size of about 30 Gpc, if there were more than 40 e-folds of inflation.Similarly, if inflation occurred at the GUT scale (≃10 16 GeV), then it would require more than 60 e-folds.The upper bound on the value of the Hubble scale during slow-roll inflation provided by Planck [147] is The mass distribution of black holes formed during the collapse of domains with values of Λ deviating from Λ obs is determined by the size distribution (21), scaled with respect to the expansion of the Universe (11).This distribution can be expressed as at the equality time t eq = 51kyr.A convenient characteristic of this distribution, which facilitates comparison of the PBH yield with constraints on their abundance in different mass ranges (see, for instance, Figure 18 in [10]), is the mass density of PBHs per logarithmic mass interval, expressed in units of the total density of the Universe: where ρ eq = m 2 Pl /(6πt 2 eq ) represents the matter density at the time of equality.Using (23), we can obtain where ρ Λ is the energy density contained in the domain filled with Λ term which reads as where we recall that Λ is expressed in units of m 2 Pl .Thus, (24) can be expanded as For those values of Λ within the domains of inhomogeneities, as indicated in (17), the rate Γ Λ can be approximated as (a detailed derivation is provided in Appendix B): where Substituting this into expression (27) and using relations ( 16), (28), and (29), we finally obtain Comparing distribution (30) with the model used in [10] to quote the constraints on the density fraction β deposited in PBHs at the moment of their formation, we arrive at the following condition Condition ( 32) is useful for assessing the consistency of ΛPBH formation with cosmological constraints on the abundance of PBHs across different mass ranges.By analyzing the combined constraints on β for a monochromatic mass function, as presented in Figure 18 of [10], we can verify the consistency of considering ΛPBHs as candidates for seeding SMBHs and IMBHs, taking into account the CMB constraints on the inflation scale.For the seeding masses M ΛPBH ≈ 10 2 M ⊙ , the abundance is constrained to the level of β ≈ 10 −14 [10], which is saturated at the inflation energy scale H inf ≃ 10 −5 m Pl .At this level of precision, this can be considered as the saturation point, ensuring that it does not exceed the CMB Planck limit (22).Similar estimates of the inflation scale can be obtained for IMBHs in the mass range 10 2 M ⊙ ≤ M ΛPBH ≤ 5 × 10 5 M ⊙ , where β ≈ 10 −15 [10].For IMBH with masses M ΛPBH ≈ 10 6 M ⊙ , the energy scale H inf ≃ 10 −7 m Pl saturates the constraint β ≈ 10 −7 [10].Similarly, the constraint β ≈ 3 × 10 −7 [10] imposed for the mass scale M ΛPBH ≈ 10 10 M ⊙ , which is typically relevant for currently observed biggest SMBHs, is saturated at the inflation scale H inf ≃ 10 −7 m Pl .Therefore, it can be concluded that ΛPBHs are cosmologically consistent for serving as seeds for SMBHs as well as explaining the observed IMBHs.
Additionally, ΛPBHs with masses as low as M ΛPBH ≈ 10 −2 M ⊙ remain compatible with the constraint β ≃ 10 −11 [10] imposed by the upper CMB inflation scale limit (22).This scenario corresponds to Λ 1 ≃ 10 −76 m 2 Pl , where the formation mechanism of ΛPBHs would generate about 10 times the mass of the Jupiter PBHs.Such PBHs could potentially account for a component of DM.

Concluding Remarks
It is theorized that the large-scale structure of the Universe was shaped by quantum fluctuations of scalar fields and/or metrics during inflation.These fluctuations, scaled exponentially with conserved amplitude, gave rise to primordial inhomogeneities, culminating in the formation of the cosmic web that represents the Universe's structure.In theories involving extra dimensions, fluctuations of fields within these dimensions can also be considered.However, the fate of these field fluctuations differs significantly from those associated with conventional four-dimensional scalar fields.
While conventional field fluctuations transform their energy density into radiation during the FRW stage through decay into other particle-like species, the energy density associated with scalar fields within the internal space remains unchanged, effectively stored within the scalar field itself, manifesting as a cosmological constant.Fluctuations of the scalar field within the internal space manifest as spatial variations of the local Λ term.Domains containing an extraordinarily high cosmological constant may collapse into PBHs.
Upon investigating the mass distribution of such ΛPBHs, we find that it may satisfy existing cosmological constraints on the abundance of PBHs without conflicting with the upper bound on the inflation energy scale inferred from CMB measurements, within the mass range from 10 −2 M ⊙ up to 10 10 M ⊙ .Of particular interest is the possibility of associating ΛPBHs with masses of 10 2 M ⊙ with seeds or supermassive black holes (SMBHs) and associating those with masses spanning the range from 10 2 M ⊙ to 10 6 M ⊙ with intermediate-mass black holes (IMBHs).
The lightest ΛPBHs of masses 10 −2 M ⊙ can potentially contribute to the dark matter budget of the Universe.
The inevitable clustering of PBHs formed by the connected mechanism, driven by inflationary dynamics, leads to the formation of a Swiss cheese-like special structure of domains filled with high values of the Λ term.This clustering may impact the characteristics of the observable spectrum of gravitational waves in the nanohertz frequency band, which are believed to be a signal from SMBHBs [148].

Acknowledgments
The A Justification of Formula ( 14) for Mass Measurement by a Distant Observer Here, we explore the conditions under which Equation ( 14) can be reliably used.We consider a three-dimensional space of volume Λ 1 > Λ 0 .Our objective is to estimate the mass of such a region as observed by a distant observer.For our estimation, we make several assumptions: the field distribution varies slowly, allowing us to neglect its time dependence; we operate far below the inflationary scale, implying that the term aR 2 is negligible; the domain with the higher energy density has approximately spherical geometry; and the amount of ordinary matter is negligible.
We set the four-dimensional effective action Equation ( 5) as and assume that the D dimensional metric depends on the radial coordinate v, leading to the interval in the form which evolves the action (33) into the following expression where the term aR 2 is neglected.Now, the situation is essentially simplified, allowing for the analytical evaluation of the mass using the nontrivial equations of the theory B 2 Subtracting the first two equations yields so that Equations ( 36) and ( 37) are reduced to a single equation with the solution We assume that the value of Λ remains constant within the sphere of radius l * , denoted as Λ(l * ) = Λ 1 .Therefore, inside the sphere, ( 41) is converted into the expression Since, in our setup, there is no point-like mass producing the singularity in the center, C 1 = 0, and hence, we obtain the well-known de Sitter metric.For a distant observer in the Minkowski space we assume Λ(v) = 0 everywhere outside the sphere of radius l * so that one can derive Thus, the mass of the domain filled with Λ 1 reads 2M ∞ = Λ 1 l 3 * m 2 P /3 and finally can be expressed as providing the definition Λ 1 = 8πρ Λ /m 2 P .Note that taking into account ( 39), (42), and (43), Equation ( 38) becomes an identity.

B Probability of Domain Formation with Specific Energy Density during Inflation
To estimate the volume fraction dP (Λ) filled with a specific value of the energy density ρ Λ = Λm 2 P /8π, we need to relate the field fluctuations in the extra dimensions during the inflationary stage to the parameter Λ.Following the approach outlined in [149], we consider the extra metric as the background one, allowing the parameter c eff to vary together with the scalar field ζ.The scalar field fluctuates intensively during inflation.Therefore, our first approximation is to assume a pure de Sitter metric, such that the function ζ(x, y) is governed by the equation This equation, one of Equation set (4), has an asymptote of ζ 1 (u) in a chosen 3D volume under the horizon.The surrounding 3D space is characterized by another static configuration ζ 0 (u) ̸ = ζ 1 (u).
The parameter value c eff [ζ 0 ] expressed by the last equation in Equation (8), is assumed to be small in order to avoid disrupting the dynamics of the inflationary stage.Additionally, it is assumed to tend to a post-inflationary value of 10 −123 m 2 P .Therefore, the quantity is responsible for the excess energy density.
The scalar field action (1) can be reduced to the standard form as follows from ( 5) and (51).The excitation amplitudes ϕ (a) act as independent free fields with an initial amplitude equal to ϕ 0 (u) [149].The probability of finding a set of functions ϕ (a)   dP ({ϕ a }) ≃ a dϕ (a) • q a /π exp − a q a ϕ (a)2 , q a = 4π 2 3 According to (52), The final expression for the probability is where the pre-exponent follows from the normalization.The number of domains with specific Λ is related to the probability as dn = e 3N dP (Λ).(56)
work of S.G.R. was funded by the Ministry of Science and Higher Education of the Russian Federation, Project "New Phenomena in Particle Physics and the Early Universe" FSWU-2023-0073, and the Kazan Federal University Strategic Academic Leadership Program.The work of A.A.P. was funded by the development program of the Volga Region Mathematical Center (agreement No. 075-02-2023-944).