Quantum production of black holes at colliders

We investigate black hole production in pp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p\,p$$\end{document} collisions at the Large Hadron Collider by employing the horizon quantum mechanics for models of gravity with extra spatial dimensions. This approach can be applied to processes around the fundamental gravitational scale and naturally yields a suppression below the fundamental gravitational scale and for increasing number of extra dimensions. The results of numerical simulations performed with the black hole event generator BLACKMAX are here reported in order to illustrate the main differences in the numbers of expected black hole events and mass distributions.


Introduction
The possibility to produce black holes (BHs) at particle colliders is directly related to the question whether the fundamental gravitational scale is somewhere in the few TeV range, as was suggested in scenarios with extra spatial dimension, like the ADD model [1,2] and the RS model [3,4] (see also Refs. [5,6] for a comprehensive review). Above the gravitational scale, it is generally expected that BHs can be created and finding signatures of their decays would be evidence in favour of these extra-dimensional models [5,7]. During the recent years, it was proposed that high energy particle colliders could turn out to be huge BH factories [8,9], and there have actually been many searches to observe the production and decay of semiclassical and quantum BHs at a e-mail: nicusorarsene@spacescience.ro b e-mail: casadio@bo.infn.it c e-mail: octavian.micu@spacescience.ro the LHC. 1 The ATLAS collaboration looked for events with jet+leptons [12,13] or dimuon [14] in the final state of pp collisions at a centre-of-mass energy of √ s = 8 TeV. At the same time, the CMS collaboration was searching for energetic multi-particle final states, as well as for resonances and quantum black holes using the dijet mass spectra at √ s = 8 TeV [15,16]. These searches and their results are very important for the community, especially in the context of the existing extra-dimensional models. They represent direct comparisons between an experiment and the theoretical predictions for new physics at these energies and can be used to constrain the parameters of the models [17]. For example, the CMS collaboration [15] excluded the production of quantum/semiclassical BHs with masses below 4.3 to 6.2 TeV (depending on the models) with 95 % confidence level, while ATLAS results indicate the threshold mass of the quantum BH to be larger than 5.3 TeV [13]. However, this exclusion limits are strongly dependent on the BH production cross section and different decay modes.
Here, we analyse the BH production by employing the modified cross section in the ADD model [1,2] obtained in Ref. [18] from the Horizon Quantum Mechanics (HQM) of localised sources [19][20][21][22][23][24][25]. In fact, this approach was specifically devised to yield the probability that a particle is a BH, and is therefore perfectly suited to address this issue. To perform the analysis, we then adapt the BLACKMAX code [26], one of the most powerful and widely used BH event generators, which includes different scenarios like tension/tensionless rotating/non-rotating BHs. 2 The results of our findings will be presented in Sect. 3. Before that, we familiarise the reader with the HQM and provide some use- 1 A BH is considered semiclassical if it decays via Hawking radiation, whereas it is generically called quantum if the decay is not thermal, including the case of a stable remnant [10,11]. 2 Examples of other available BH event generators are CHARYB-DIS2 [27], QBH [28] and CATFISH [29]. ful references for a more in-depth study of the formalism in Sect. 2.

Horizon quantum mechanics
The HQM for static sources [19][20][21][22][23][24] was extended to higher dimensions in Refs. [18,25], which can be naturally applied to BHs in the ADD scenario. Let us start by considering the wave-function for a localised massive particle as given by a spherically symmetric Gaussian wave-packet of width in D spatial dimensions [representing a source in a (D + 1)dimensional space-time], whose form in momentum space is where =h/ = m D D / , m D is the fundamental gravitational mass and D =h/m D represents the corresponding length scale. We study the simplest case and assume that when a BH forms, it will be described by the (D + 1)dimensional Schwarzschild metric, where the classical horizon radius is given by and G D = D−2 D /m D represents the fundamental gravitational constant in this ADD scenario.
We now consider the mass-shell relation in flat space, where m is the rest mass of the source, and express E in terms of the above horizon radius (4), r H = R D (E). After using these results in Eq. (2), the normalised horizon wave-function reads [18] where the normalisation N is obtained from the Schrödinger scalar product in D spatial dimensions and the step function ensures that the gravitational radius r H ≥ R D (m), since m is the minimum energy eigenvalue contributing to the packet. We can now calculate the probability for the particle to be a BH, where P < (r < r H ) represents the probability density for the particle to be inside its own gravitational radius r H and is the product of two factors: the probability for the particle to be located inside a D-ball of radius r H and the probability density that the horizon radius equals r H . In this particular case, the BH probability depends on the Gaussian width , particle mass m and number of spatial dimensions D. We can further assume = λ m = m D D /m is the Compton length of the source, which represents the minimum uncertainty in its size, so that = m and the probability only depends on m and the number of dimensions D [18], where (a, b) is the upper incomplete gamma function and γ (a, b) the lower incomplete gamma function. The above expression can be computed numerically and is displayed in Fig. 1 for D = 5, 7 and 9.
There are a few important observations regarding this result. First of all, like in D = 3, there is no sharp threshold for BH formation, but the BH probability drops very fast for m < m D (or, equivalently, > D ). Moreover, for any given mass, say m m D , the probability P BH (m; D) decreases for increasing values of D. In the next section, we will focus on expressing these differences in a more quantitative way.

Cross section p p → BH
We will now focus on the implications for BH searches at the LHC. As stated in the Introduction, we performed the numer- continuous black lines represent the modified case (9) and the events ratio gives the corresponding suppression factor. Each curve is an average over 10 4 simulations of tensionless non-rotating BHs in BLACK-MAX 2.02.0 ical simulations using BLACKMAX 2.02.0 and considering tensionless non-rotating BHs. In the standard configuration, BLACKMAX employs the BH production cross section where b D = 2 1 + (D − 1) 2 /4 1 2−D and R D is the horizon radius (4). Among other parameters, one can set the values of the fundamental gravitational scale m D and the minimum BH mass m min . In fact, it is important to remark that no threshold of BH production is fixed in the standard scenarios, although one expects that BHs do not form below a certain mass because of quantum fluctuations, and one can at best constrain m min from the data. Typically, we shall consider m min m D in order to ensure that no BH is produced with a mass below m D in the standard case.
In the HQM picture of Refs. [18,25], the effective BH production cross section is instead given by where P BH (E) = P BH (m = E; D) is the probability (7) for a particle with energy E to be a BH, while σ BH is still given by Eq. (8). Note that there is now no need for imposing a minimum BH mass, since P BH acts as a proper quantum regulator. In order to implement this improved cross section, we considered the fundamental gravitational mass scale to have the same value as in the standard case, and set a minimum BH mass m min = 0.2 m D for computational convenience. 3 Finally, we added a subroutine to BLACKMAX which weighs the standard BH mass distributions with the probability P BH (E). We first illustrate the typical differences between the simulations that employ the standard cross section (8) and the HQM cross section (9) in Figs. 2 and 3. Later on, we will investigate more general cases and include comparisons to the current bounds on m D and m min by the ATLAS and CMS collaborations. The blue dashed lines are obtained by employing the standard BH production cross section (8) for p p collisions at √ s = 8 TeV, with m D = 1 TeV, and setting the minimum BH mass to m min = m D (in Fig. 2) or m min = 2 m D (in Fig. 3). The continuous black lines in the same plots represent the analogous BH mass distribution derived from the modified production cross section (9). First of all, in agreement with Fig. 1 and the HQM approach [18,[21][22][23], BHs with masses below the fundamental scale of gravity are now possible and no sharp threshold effect like the one forced in the standard case exists. For m min = m D = 1 TeV, the HQM cross section (9) leads to a significant suppression of BH production, whereas for m min = 2 m D = 2 TeV the situation is reversed. Besides investigating how the differential production cross section varies with the value of the resulting BH mass, we can also compare the total cross sections in the two cases. This comparison is given by the events ratio, the ratio between the HQM total production cross section relative to the standard case, which we also display in the plots. In particular, the event ratio is smaller than one (thus signalling a suppression) for m min = m D = 1 TeV, but larger than one (indicating an enhancement) for m min = 2 m D = 2 TeV. We also notice that this ratio is always smaller for larger D, again in agreement with Fig. 1. This can be viewed as a check which makes us confident that our numerical simulations are accurate.
In light of these preliminary results, it appeared interesting to study how the events ratio depends on the value of the fundamental gravitational scale for different numbers of spatial dimensions. This analysis is presented in Fig. 4 for m D = m min . We see that, in this case, regardless of the value of m D , this ratio is always smaller for larger number of extra-dimensions (at the same value of the gravity scale). Another feature we notice is that, for all numbers of spatial dimensions, from around m D = 2 TeV for √ s = 8 TeV (respectively m D = 4 TeV for √ s = 13 TeV) the events ratio starts to increase with m D , eventually crossing unity from below. Even though it seems that more BHs are produced in the HQM scenario than in the standard case for higher values of m D , we have to remember that the total cross section decreases with D and the number of expected BH events remains very small. This dumping of the total cross sections for increasing m D is exemplified in Fig. 5 for D = 6 (but the same behaviour holds in all cases): for smaller values of m D , the standard production is larger than the HQM expectation, and the two cross at a relatively large value of m D (which depends on D). As one can see, even where σ p+ p→BH predicted by the HQM is larger, the actual values of the total  Figure 5 shows, for instance, that when √ s = 8 TeV, the cross sections are of the order of 10 4 pb for m D = 1 TeV, but they reduce to about 1 pb for m D = 4 TeV.
So far we mostly analysed cases with m min = m D in the standard scenario and compared with the HQM predictions for the same m D . The tables in Appendix A present the dependence of the events ratio on m min and m D , taken to be independent parameters for the standard scenario, for the same m D in the HQM case. The thick black lines in the tables show where the events ratio crosses over one: below the lines, the HQM predicts less BH events, whereas the standard simulations predict less such events above the lines. It is in particular interesting to compare the number of events predicted by the HQM with the number of events one expects to see in the standard case for the current lower bounds imposed by the LHC collaborations on m D and on the minimum mass that BHs can have [12][13][14][15]17]. In Fig. 6, we show a comparative analysis between the BH production cross sections (which ultimately translate into the number of events expected to be produced at the LHC), and the distributions of BH masses for these two cases. In the plots, we assumed the strongest lower bounds on m D and m min available [15,17], and as usual compared with the HQM predictions for the same m D . It immediately appears that the HQM predicts more BH events. Upon examining Fig. 6 further, one also notices that the BH mass distributions differ: the HQM predicts that most BHs are produced with masses smaller than the values expected in the standard scenario. This is no surprise, given that we assumed the same value for m D in the two scenarios, the lower bounds imposed by the LHC groups on m min are stronger (higher values) than the bounds imposed on m D , and that BHs with mass below m D are possible in the HQM. Since the HQM yields no minimum BH mass, the above comparison is not completely significant, and one should actually constrain only m D using experimental data in this scenario. We then determined the value of m D in the HQM formalism for which the number of BH events is expected to be the same as in the standard case at the current bounds. This means we set the BLACKMAX parameters for the standard case equal to the LHC bounds, and changed the value of m D in the HQM simulations. The results are presented in Fig. 7, from which one can see the events ratio is roughly equal to one for m 5 6.8 TeV, m 7 6.4 TeV, m 9 5.3 TeV. In all cases, the HQM lower bounds on m D therefore appear stronger than in the standard case. The number of BH events expected at the LHC is about the same as in the corresponding standard scenario, but with a very different distribution of masses (most of the BHs are produced with masses lower than the minimum BH mass in the standard case). We thus caution the reader that in deriving these HQM lower bounds on m D , we neglected the impact that the different HQM distributions of BH masses may have on the likelihood for the LHC collaborations to detect them. In fact, before the LHC collaborations reached the current limits, they also scanned the parameter space below those values.
With that disclaimer in mind, we can still consider the HQM bounds on m D and simulate the BH production at √ s = 13 TeV. Figure 8 shows the expected distribution of BH masses in this case, with cross sections of the order of a few times 10 −2 pb.

Conclusions
We investigated the implications of the HQM on the BH production cross sections at the LHC in the context of the ADD models with extra spatial dimensions. We used BLACKMAX to perform numerical simulations allowing us to compare both the production cross sections and the resulting BH mass distributions. The events ratio was used to express quantitatively the differences in the BH production cross sections in the two cases. We find that this ratio is always smaller for larger number of extra-dimensions, and that in each case it eventually increases with m D until it becomes larger than one. A wide range of cases are presented in Fig. 4 and the tables in Appendix A. When looking at the distribution of BH masses, in particular, we find that the HQM predicts most BHs are produced with masses smaller than the values expected in the standard scenario.
We also compared the cross section predicted by the HQM with the standard one for the current lower bounds imposed by the LHC collaborations on the fundamental gravity scale m D and m min . We found that in the HQM case more BHs are produced. We thus determined new lower bounds on m D , by finding the value of the fundamental gravity scale at which the number of BH events is expected to be the same as in the standard case at the current bounds: m 5 6.8 TeV, m 7 6.4 TeV, m 9 5.3 TeV. Finally, we calculated the BH production cross sections for these new lower bounds at √ s = 13 TeV. It will also be interesting to investigate the implications of the HQM for BH remnants at the LHC [10,30], or other modified decay channels [31], but we leave this analysis for future work.