BoseEinstein graviton condensate in a Schwarzschild black hole
Abstract
We analyze in detail a previous proposal by Dvali and Gómez that black holes could be treated as consisting of a BoseEinstein condensate of gravitons. In order to do so we extend the EinsteinHilbert action with a chemical potentiallike term, thus placing ourselves in a grandcanonical ensemble. The form and characteristics of this chemical potentiallike piece are discussed in some detail. We argue that the resulting equations of motion derived from the action could be interpreted as the GrossPitaevskii equation describing a graviton BoseEinstein condensate trapped by the black hole gravitational field. After this, we proceed to expand the ensuring equations of motion up to second order around the classical Schwarzschild metric so that some nonlinear terms in the metric fluctuation are kept. Next we search for solutions and, modulo some very plausible assumptions, we find out that the condensate vanishes outside the horizon but is nonzero in its interior. Inspired by a linearized approximation around the horizon we are able to find an exact solution for the meanfield wave function describing the graviton BoseEinstein condensate in the black hole interior. After this, we can rederive some of the relations involving the number of gravitons and the black hole characteristics along the lines suggested by Dvali and Gómez.
May 2017
ICCUB16032
1 Introduction
In an interesting saga of papers, Dvali, Gómez and others[1, 2] have put forward an intriguing suggestion: black holes (BH) could perhaps be understood as BoseEinstein condensates of gravitons. If correct, this would reveal a deep quantum nature of such fascinating objects and could lead to an alternative understanding of some of the most striking features of BH; for instance, Hawking radiation [3] could be understood as being due to leakage from the condensate. Besides, this picture brings new ideas about the Bekenstein entropy [4], the absence of hair [5], as well as the quantum nature of information storage and the possible information loss in BHs [6].
The main point of these works is that the physics of BH can be understood in this picture in terms of a single number , the number of (offshell) gravitons contained in the BoseEinstein condensate (BEC). These condensed gravitons have a wavelength , being the Planck length; they have a characteristic interaction strength and the leakage leads to a Hawking temperature of order , equal to the inverse of . The mass of the BH is and its Schwarzschild radius therefore is given by , thus agreeing with the Compton wavelength of the quantum gravitons , in accordance with the uncertainty principle that dictates in the ground state of the quantum system. Therefore, up to various factors of the Planck mass everything is governed by , the number of intervening gravitons. Modulo some assumptions, all these results stem from the basic relation (unless otherwise stated we work in units where )
(1) 
that relates the number of gravitons , the mass of the gravitating object and its gravitational radius . For a Schwarzschild BH .
We found these results intriguing and we set up to try and understand them in a different language and, if possible, attempt to be more quantitative. On the process we have separated slightly from the original line of thought of the authors. The original approach of [1, 2] is not geometric at all. There is no mention of horizon or metric. Here we will adopt a more conservative approach. We will assume the preexistence of a classical gravitational field created by an unspecified source that generates the Schwarzschild metric. As it is known, nothing can classically escape from a BH so if we wish to interpret this in potential terms (which of course is not correct but serves us for the purpose of creating a picture of the phenomenon) it would correspond to a confining potential. On and above this classical potential one can envisage a number of quantum fields being trapped. For sure there is a gravity quantum field present inside the horizon; hence gravitons. Other quanta may get trapped by the BH potential as well, but a possible graviton BEC seems particularly challenging to treat and this is the purpose of this work.
Continuing with our semiclassical analogy, these trapped gravitons as well as other quanta present have had a long time to thermalize in a (eternal) Schwarzschild BH and it is therefore natural to expect that after cooling they can form a BoseEinstein condensate. Of course these ‘gravitons’ are in no way freely propagating transverse gravitons. They are necessarily offshell () and have some sort of effective mass. It may help to get a picture of the phenomenon to think of them as quasiparticles.
It is well known in BEC theory that Bose–Einstein condensation of a spatially homogeneous gas with attractive interactions is precluded by a conventional phase transition into either a liquid or solid. When confined to a trap; however, such a condensate can form, if its occupation number is low. Repulsive forces act to stabilize the condensate against collapse. Knowing this, one would immediately think that gravitons do not have repulsive interactions, at least naively, and that therefore a BEC is impossible to sustain, particularly as we expect to be very large. To this objection one could reply in a twofold way. First, it is up to the equations to determine whether such a condensate is possible or not (we will see below that indeed it is, at least in the case of vanishing angular momentum). Secondly one might also answer that, in theories of emergent gravity, the ultimate nature of gravitons may be some type of fermionic degrees of freedom (such as e.g. in the model suggested in [7]). Then, repulsion is assured at some scale and fundamental collapse prevented.
We will try to identify and propose a consistent set of equations describing a BEC constructed on top of the classical field created by a BH^{1}^{1}1A classical field is not the same as a quantum condensate, although the latter may trigger the former.. We will succeed and we will see that remarkably enough the characteristics of the resulting BEC are uniquely described in terms of the Schwarzschild radius of the BH and the value of a dimensionless parameter, interpreted as a chemical potential. A condensate appears not only to be possible but actually intimately related to the classical field that sustains it and determines its characteristics.
It would therefore be tempting to go one step beyond and reverse the order of the logical implication and eventually attempt to derive the classical field as a sort of mean field potential à la HartreeFock. However we will stop short of doing so here. In any case, even without embarking in the discussion just mentioned, we are aware of the speculative nature of the present study.
2 Building up a condensate over a Schwarzschild background
In what follows we shall adhere to the following notational conventions: the Einstein tensor is constructed with the metric in the usual way. We will denote by the background metric that in our case will invariably be the Schwarzschild metric. Perturbations above this background metric will be denoted by , so . We will use the Minkowskian metric convention .
To initiate our program we should, first of all, identify an equation (or set of equations) that could provide a suitable description of a BEC in the present context. In other words, we have to find the appropriate generalization of the GrossPitaevskii equation. The graviton condensate has necessarily to be described by a tensor field that within our philosophy has to be connected necessarily with a perturbation of the classical metric.
In order to keep things as simple as possible we will attempt to describe only condensates with quantum states having . This will translate to spherically symmetric perturbations of the gravity field only
(2) 
The Einstein tensor derived from the previous metric will be expanded up to second order in to retain the leading nonlinearities (selfinteractions of the desired condensate).
It is well known that Birkhoff’s theorem [8, 9] guarantees the uniqueness of the solution of Einstein’s equations in vacuum with the properties of having spherical symmetry and being static. As every dimensionful quantity can be expressed as a function of the Schwarzschild radius , the difference between a given Schwarzschild metric and any perturbed solution built upon it must necessarily correspond to a change in the mass to . The fact that spherical symmetric perturbations are related to shifts in mass, indicates a correspondence between any perturbation (i.e. each ‘graviton’) and a certain amount of energy that is reflected in a change of the BH mass.
2.1 Einstein equations as GrossPitaevskii equations
The familiar GrossPitaevskii equation [10] employed to describe BoseEinstein condensates is a nonlinear Schrödinger equation; i.e. an equation of motion that contains selfinteractions (hence the nonlinearity), a confining potential for the atoms or particles constituting the condensate, and a chemical potential, the conjugate thermodynamic variable of the number of particles or atoms contained in the condensate.
Among all these ingredients, the Einstein equations perturbed around a BH metric contain most of them. They are already nonlinear and while there is no an explicit confining potential (as befits a relativistic theory) they do confine particles, at least classically, because if the selected background corresponds to a Schwarzschild BH, the strong gravitational field classically traps particles inside the horizon. In addition, they are the only known consistent equations describing a ranktwo tensor. Finally, the background (i.e. the BH metric) is a solution of Einstein equations; therefore they are necessarily part of the answer. However, there is one ingredient still missing, namely the equivalent of the chemical potential. Therefore we have to extend the formulation of perturbations around a classical BH solution to the grandcanonical ensemble by adding to the appropriate action a chemical potential term.
As it is well known since the early days of quantum field theory [11] there are no conserved currents or continuity equations for fundamental fields that are chargeless (such as a real KleinGordon field). Therefore there is no way of defining a number operator for freely propagating gravitons or photons. Nonetheless, in the picture of [1] and [2] the situation is different. If a BEC is present with ‘gravitons’ acquiring all the same momenta and being weakly interacting for a macroscopic BH (recall ) the total energy stored in the condensate should be and this quantity would be a conserved one. Therefore, lest change, would be conserved. The previous reasoning shows very clearly that the ‘gravitons’ contemplated in the present scenario, if realized, have nothing to do with freely propagating gravitons.
The energy contained in a given volume occupied by a noninteracting scalar field is
(3) 
by analogy, in the present case
(4) 
where we assume that the energy per graviton is constant and approximately given by , and
(5) 
with . While there is no formally conserved current, the above quantity can be interpreted as a ‘graviton’ number density, and the integral of the graviton density (5) in the interior of the BH has to be interpreted as the number of constituents of the BEC.
The above considerations can now be phrased in a Lagrangian language. The chemical potential term in the action should be related to the graviton density of the condensate inside a differential volume element . In order to respect the basic symmetry of General Relativity (GR), the simplest form of introducing such a term is by means of
(6) 
This term does of course resemble a mass term for the spin2 excitation and indeed it is some sort of effective mass in practice as the ‘gravitons’ in the BEC are quasinoninteracting. However we will see that GR eventually requires for the quantity to actually be position dependent, i.e. , and transform as a scalar.
Given the resemblance of the term that we interpret as chemical potential to a mass term, it is legitimate to ask why not use the FierzPauli form of the mass term, that is known to be able to provide ghostfree propagation of gravitons (see e.g. [12] for a detailed discussion of massive gravity). The answer is simple: our ‘gravitons’ correspond to timeindependent solutions; they do not propagate and the issue of ghosts is totally irrelevant in the present discussion. We do not advocate adding (6) as a fundamental ingredient of gravity, but only as a means of describing in an effective and thermodynamical way the formation of a BEC condensate —exactly what the GrossPitaevskii equation is meant to implement. From a more pragmatic point of view, modifications of (6) in the sense of making this term look like a FierzPauli mass term do not really change the results at the qualitative level, but we emphasize that it is not really necessary to worry about the ghost issue at this point.
Although the additional term may look odd, it is invariant under the gauge group of diffeomorphisms. Part of this statement is shown in [13]: under an infinitesimal displacement in the coordinates of the form the full metric changes as . The same gauge transformation rules the background metric. As the perturbation is defined by , under the same perturbation of the coordinate system, the same rule of covariant transformation is obtained. Together with the fact that the chemical potential behaves as a scalar under a general coordinate transformation, ensure automatically the general covariance of the theory. However, while the formalism is diff invariant, it is not background independent. The separation leads to an action that depends on the choice of the background metric, which is the one of the BH. Likewise also will depend on the background metric. We will not postulate any particular dependence of on . In fact, at this point it would be conceivable that the only consistent solution implies . This will not be the case however and will be determined in the subsequent from the requirement of selfconsistency of the proposal. Therefore an appropriate action for the field is
(7) 
Indices are raised and lowered using the full metric . A completely covariant expansion in powers of can be performed up to the desired order of accuracy. By construction these equations will be nonlinear, but we will keep only the leading nonlinearities to maintain the formalism simple and analytically tractable. It is worth noting again that the action (7) is as shown reparametrization invariant, but it is not background independent. The fact that fluctuations take place above a BH background –in whatever coordinates one chooses to describe it– does matter.
The action principle yields the two equations of motion for the field , namely
(8) 
where is the Kronecker delta and . It is important to keep in mind the following: we are working here in the grand canonical ensemble; this implies that the magnitude is an external field and does not vary in the action. In particular, for these equations of motion is independent of , so . This should be the way of introducing the chemical potential. Otherwise, were not an external field, it would be necessary to take it into account when performing variations to derive the equations of motion. An equation of motion would be obtained for and this equation would nullify automatically the perturbation as well as the chemical potential itself, i.e. and . Therefore, the external field is introduced in the theory as some kind of Lagrange multiplier. As Lagrange multipliers have constraint equations, the chemical potential of the theory may not be arbitrary at all. It must satisfy binding conditions with . The main difference is that these constraints are implicit in a GR theory. The restriction of is through the general covariance conditions for the action. In other words, the diffeomorphisms covariance implies the covariant conservation of the Einstein tensor, and this in itself entails the same for the chemical potential term
(9) 
The covariant derivative is defined using the full metric . This differential equation of motion is valid up to every order in perturbation theory.
Hereinafter, we will proceed to find acceptable solutions of these equations and interpret them. We will separate the problem into two regions: outside and inside the BH horizon. We certainly expect that the graviton condensate will disappear quickly in the outside region; imposing this as a boundary condition for we will see that in fact the condensate is identically zero on this side of the horizon. On the contrary, a unique nontrivial solution will be found in the interior of the BH.
3 Outside the horizon
In this section it will be shown that even after the inclusion of there is not other normalizable solution outside the black hole horizon than the trivial solution for the perturbation. If we place ourselves far away from the BH, in the limit, one is able to keep only dominant terms in the Einstein equations. In this perturbative treatment the chemical potential, far away from the condensate, is expected to behave as a perturbative small parameter. Therefore, on the LHS side of the equations of motion (8), possible terms of the form will be considered as and neglected. We will make use of another equation that, although it does not contribute with new information, will make explicit how the components of the perturbation behaves between each other; namely, any of the two angular components of the Einstein tensor . In this limit, as we want a vanishing perturbation at infinity, the following ansatz is imposed in the faraway region:
(10) 
Before proceeding, an additional consideration is needed: since for us represents a localized BEC, in analogy with the wave function of a confined particle, the perturbation must be a squareintegrable function. This is
(11) 
There are two reasons that lead us to the previous requirement. Let us review first why one has to request squareintegrability. Note that the solutions of (8) should not be understood as a quantum field, but rather as solutions of the GrossPitaevskii equation. Here we adopt Bogoliubov theory [14] and its interpretation of the GP wave function. This is the commonly accepted interpretation and essentially it boils down to the fact that in the large occupation limit () the creation and annihilation operators of the ground state can be approximately treated as commuting cnumbers, and hence the many body quantum problem is described by a classical function called the ‘macroscopic wave function’ or simply the ‘order parameter’. Because the modulus square of this quantity is proportional to , it is therefore proportional to . Hence the order parameter itself is proportional to . If the potential is not uniform, may of course depend on the coordinates. From this the need to require squareintegrability follows. This interpretation is in the present case also supported by the dimensionality of .
Why then this precise condition? The leading role of the geometry has been to provide, in an indirect manner, a “confining potential” that traps gravitons inside the horizon. Therefore, from this point of view when giving a physical interpretation to the chemical potential (see section 4) one implicitly assumes that the geometry of the spacetime is flat and that geometry acts via the external potential and the chemical one. In the present situation, it turns out that coincides with the one corresponding to a 3 dimensional spatial flat metric in spherical coordinates; therefore is reparametrization invariant in 3 dimensions. This interpretation is equivalent to defining a wave function normalized by the temporal component of the metric tensor, . Then, any magnitude computed by means of an integration over a 3spatial volume is effectively written in a fully covariant way as , since the 4dimensional metric can be decomposed as , with if the metric is diagonal. Therefore, the volume element is coordinate dependent as well as the square modulus of the macroscopic wave function , but the combination of both is not; and we end up finally with the corresponding volume element to the previously discussed flat spacetime.
Whatever the case, in this perturbative limit the components of the Schwarzschild metric accomplish ; then the integral in Eq. (11) should fulfill
(12) 
For this to happen, if a power law ansatz at infinity is imposed, the exponents must obey . Then the equations corresponding to , and respectively become
(13) 
The requirement of the solution being squareintegrable exclude that and could be zero. This implies automatically the neglection of the terms proportional to the Schwarzschild radius. Even more, with the previous considerations, (13) reduce to the following asymptotic identities
(14) 
(15) 
From the angular equation (15), it is mandatory for both terms to contribute at infinity. If this is not the case, this would nullify (or ), and then, via both equations in (14), fix (or ) and ; i.e. . The competition of both terms is possible if and only if . This condition modifies the angular equation (15) as
(16) 
Then, from the second equation in (14) the following relationship can be read
(17) 
This automatically leads to a null chemical potential as the only possible solution for this region with the proposed ansatz.
We have also explored the possibility of exponentially vanishing chemical potential when with analogous conclusions. If we change the ansatz and impose an exponential decreasing solution for the perturbations at infinity,
(18) 
the linear order dominates in front of the higher ones. As decreasing solutions are expected, the parameters , and must be not null. If one use this ansatz in the Einstein equations partially evaluated at infinity, one finds that in the angular equation there is no way for both leading terms to compete between each other
(19) 
In conclusion, depending on the fact whether if is bigger or not than , or . For the first case, if , the leading terms of the temporal equation
(20) 
fixes . Being both, and null, the chemical potential disappears from the theory in this region. For the second case, if , the same equation nullifies (as no null solutions for are expected). No perturbation makes the chemical potential senseless. To sum up, this ansatz implies as well that the only solution at infinity is a null perturbation, , and the disappearance of the chemical potential from the theory, . Likewise it can be seen that a much faster Gaussian decay is also excluded with analogous calculations.
Before concluding this section we will give additionally a numerical argument that confirms the vanishing of the condensate and the chemical potential in the outer region when the perturbations are null at infinity. The following change of variables allows rewriting the relevant equations (8) in terms of dimensionless quantities
(21) 
This quantities are universal as and are the only physical parameters. This redefinition maps the exterior part into a compact interval: is transformed into and into . It is well known that Schwarzschild spacetime is asymptotically flat; therefore, at infinity () it is expected for the perturbation to vanish (). Nevertheless equations are difficult to deal with because they are singular at , i.e. it is not possible to isolate the higher order derivatives.
Let us see how we can proceed. As mentioned before, we have to settle for a point (but ) to set up boundary conditions for the numerical integration procedure to start. Likewise setting at any arbitrary value of outside the event horizon, immediately triggers divergences at the first step of the routine. Therefore we take a “small” initial value for the perturbation in a point near infinity (), and then decrease this initial value; this is .
With the set of solutions for the components of obtained for each of the initial condition imposed, it is possible to compute the quantity defined in Eq. (11); the integral extends over the exterior of the BH (i.e. runs between 0 and 1). In Figure 1 the behaviour of the integral of outside the BH is presented when the initial condition for the integration approaches zero at a fixed point far away from the condensate (near infinity or ). The numerical analysis results very stable and the initial condition can be reduced as many orders of magnitude as desired. The graph is presented in a log–log scale, hence the linear behaviour. The points are fitted by a linear function with a slope approaching as more points with lower initial condition are added. This implies that the integral converges quadratically to zero as the initial condition for nullifies. This seems to confirm —also numerically— that there is no condensate, i.e. and no dimensionless chemical potential in the outer region of the BH.
4 Inside the horizon
In order to study the behaviour of our equations when the horizon is crossed, it is convenient to keep working in the and variables introduced in the previous section and redefine the perturbations as
(22) 
With the new definitions the Einstein equations become more compact. At the same order as before, the temporal, radial and angular components of the Einstein tensor read as
(23) 
4.1 Linearization near the event horizon
In order to get a feeling for possible solutions to these equations we consider their linearized approximation. Only terms linear in , and their derivatives are kept in the region . No a priori assumption for the dimensionless chemical potential is made, therefore is a priori considered as a linear contribution in the perturbation. That is
(24) 
Let us now make the following ansatz for the new perturbations
(25) 
while we leave the chemical potential as a free function. The three equations take respectively the following form
(26) 
It is worth noting that if , all terms contribute as we get closer to the event horizon and behaves as a constant. In this situation, we obtain three equations for the coefficients
(27) 
The system of equations is algebraic for the variable (that as we have seen should behave as a constant as ); therefore, it is possible to eliminate by combining the temporal and radial equation. In the present ansatz this leads to
(28) 
Together with the angular equation, this determines the solutions up to a single constant. There are two possible solutions for this system

(29)

(30)
In any case, at least one of the perturbations is divergent over the event horizon. Nonetheless, the first solution appears to be integrable, while this is not the case for the second one. In case i) while if ii) is taken as solution (i.e. no chemical potential at all).^{2}^{2}2Solution ii) represents however a volumepreserving fluctuation at the linear order.
However, because these solutions do not vanish when , justified doubts can be cast on the relevance of the linearized equations. Let us examine this point taking into consideration only the solution i), which is integrable according to the considerations of the previous section.
In order to see if this solution is modified when nonlinearities are switched on, we substitute it back in the equations system where selfinteractions matters, and see if the solution survives or how would it get modified. In view of the foregoing, near the full temporal, radial and angular equations read
(31) 
(32) 
(33) 
respectively. Replacing the possible solution, —this eliminates any derivative of — and into these equations, we obtain the following relations
(34) 
(35) 
(36) 
Therefore, quite surprisingly, the linear solution is still an exact solution of the nonlinear quadratic differential equations. Thus, we conclude that
(37) 
are solutions of the second order system of equations.
In addition, this exercise gives an interesting result: , where is so far arbitrary, also in sign. Note that at the linear level we got and the fact that the quadratic equation gives an correction to this result is consistent as we are implicitly assuming that , i.e. . The constant itself is arbitrary and is not determined by the structure of the equations. Changes in appear as an overall factor in the solution.
So far we have seen that the solution given in (37) satisfies not only the linearized approximation but also the full quadratic equations. This conclusion is reinforced after performing a numerical integration of the basic equations (8) expanded up to . We found that the numerical study reproduces the general features of the analytical study: as well as the function related to the dimensionless chemical potential , turns out to be constant in the interior of the BH. No other solutions are found.
The equations look simpler if expressed in terms of the functions but to understand what this solution means is better to undo the redefinition of the components of the wave function (i.e. ) by means of (22) for the solution (37). It is of interest to raise one index of the components of the perturbation with the inverse full metric . The function so obtained happens to be constant throughout the interior of the BH
(38) 
The solution for the perturbation turns out to be proportional to the corresponding metric element where each of the two belongs; i.e. and with constant. The angular degrees of freedom remain unchanged, .
4.2 Exact solution
Inspired by the previous analysis, we reformulate our main equations in terms of the metric fluctuation with mixed indices (one covariant, one contravariant). Let us write
(39) 
The full metric would become
(40) 
and seems to impose an upper limit^{3}^{3}3The upper bound for the wave function will become clear when computing the number of constituents of the condensate; such value for the wave function corresponds to the limit . for the constant values and . The exact equations of motion for the theory simplify enormously and reduce to the following ones
(41)  
These Eqs. are valid also for the external solution; however, in the case of the external sector, we know that Minkowski metric must be recovered far from the sources. This condition implies that as explained in Section 3.
Among the two possible solutions of the latter algebraical system of equations, only one is compatible with an acceptable limit for small perturbations, namely (we expect ). All things considered, the two equations in (41) are linearly dependent, hence the resulting equation of motion for a constant becomes
(42) 
Within our philosophy (8) is understood as a GrossPitaevskii equation for a condensate wave function . It is a nonlinear Schrödingerlike equation that produces a unique solution for the condensate and the chemical potential described by Eq. (42). Because the solution is constant, the ‘kinetic’ term drops and one is left with a purely algebraic, meanfieldlike, equation^{4}^{4}4One usually thinks of the GP equation as a nonlinear Schrödinger equation, hence with second order derivatives. Actually this is not always so; for uniform gas of interacting atoms the GP equation is simply , where is the interacting (repulsive) constant and , needless to say, is the chemical potential..
The previous equation defines a dimensionless chemical potential that behaves as a negative constant and, for small perturbations, is related to the meanfield solution by
(43) 
Notice that we are retrieving the solution found in the perturbative analysis when plays the role of the integration constant .
Before moving on, it is mandatory to make a comment on the covariant conservation of our equations. At the end of section 2.1 we have pointed out that the diffeomorphism invariance entails a differential equation for the chemical potential, namely (9). The covariant conservation of the Einstein tensor implies automatically the same for the LHS of our equations of motion
(44) 
The latter equation is a set of 4 equations, , but only one is non trivial; this equation is the radial one, . When the perturbations are equal and constant, i.e. , the general covariance condition yields the following differential equation for the chemical potential
(45) 
The integration is direct and the only degree of freedom is a boundary condition when integrating the differential equation that governs the chemical potential
(46) 
This value should coincide with the value of our negative and constant dimensionless chemical potential in (43); this is .
The qualitative result is the existence of a normalizable solution that can be interpreted as the collective wave function of a graviton condensate. A unique relation is obtained between this (constant) wave function and a (also constant) dimensionless chemical potential.
5 Connection with previous proposals.
It is immediate to see that the solution found for the wave function is of finite norm. Taking into account that the perturbation is null from the event’s horizon onwards, the endpoint on the integration limit can be fixed at the Schwarzschild radius. This way, the integral is
(47) 
and states that the integral of the square modulus of the wave function has a constant value. The volume element for the full metric has been used. The fact that this magnitude is constant automatically ensures a constant behaviour for the probability density of the wave function defined in (5), as . Retrieving the missing constants, we can relate the latter quantity in (47) to the integral of the density . Then, we are able to compute the total number of gravitons of the condensate. From the arguments in section 2
(48) 
Here again the upper limit for the wave function enters explicitly; if , then and the metric becomes singular. At this point we should attempt to make contact with the results of [1, 2]. Under the maximum packaging condition our previous relation agrees nicely with their proposal. The rest of relations of their work can be basically derived from this.
Possibly our more striking results are that the dimensionless chemical potential stays constant and nonzero throughout the interior of the BH, and that so does the quantity previously defined and entirely determined by the value of the dimensionless chemical potential . Therefore, it is totally natural to interpret as the variable conjugate to , the number of gravitons.
As seen above the dimensionless chemical potential has a rather peculiar behaviour. As is a constant function, then and it is not null over the event horizon. Outside it appears to be exactly zero. Let us now for a moment forget about the geometrical interpretation of BH physics and let us treat the problem as a collective many body phenomenon. It is clear why gravitons are trapped behind the horizon: the jump of the chemical potential at would prevent the ‘particles’ inside to reaching infinity. From this point of view it is quite natural to have a lower chemical potential inside the horizon than outside (where is obviously zero) as otherwise the configuration would be thermodynamically unstable. In the present solution particles (‘gravitons’ in our case) cannot escape.
However this is not completely true as the picture itself suggests that one of the modes can scape at a time without paying any energy penalty if the maximum packaging condition is verified. Let us do a semiclassical calculation inspired by this picture; using :
(49) 
To estimate (which is negative) we can use geometrical arguments to determine the flux. If we assume that for a given value of only one mode can get out (as hypothesized above) and that propagation takes place at the speed of light, elementary considerations^{5}^{5}5To determine the rate of variation of we have to multiply the surface () times the flux; i.e. the density of the mode times the velocity, assumed to be in our units. Since the density of the mode is constant in the interior, it is just . lead to
(50) 
This agrees with the results of [2] —for instance Eq. (35)— and yields . Within this picture, several questions concerning the longstanding issue of loss of information may arise as the outcoming state looks thermal [15] but apparently is not; or at least not totally so. However we shall refrain of dwelling on this any further at this point.
The main result from the previous rather detailed analysis is that the BH is able to sustain a graviton BEC and surely similar BECs made of other quanta[16]. But is that condensate really present? Our results do not answer that question, but if we reflect on the case where the limit is taken, without disturbing the BH geometry (i.e. keeping constant) this requires taking in a way that the ratio is fixed. Then one gets . This value appears to be universal and independent of any hypothesis. The metric is 100% Schwarzschild everywhere. We conclude that a BH produces necessarily a (trapping) nonzero chemical potential when the physical system is expressed in terms of the grandcanonical ensemble. Outside the BH, ().
Another way of reaching this conclusion is by taking a closer look to our exact equations in the previous section. If one sets then necessarily and one gets the classical Schwarzschild solution everywhere, but the reverse is not true. One can have but this does not imply . Let us emphasize that these results go beyond the second order perturbative expansion used in parts of this article.
As gravitons cross the horizon and are trapped by the BH classical gravitational field they eventually thermalize and form a BEC. The eventual energy surplus generated in this process is used to increase the mass and therefore the Schwarzschild radius of the BH. is now nonzero; it is directly proportional to , the number of gravitons, and the dimensionless chemical potential departs from the value , presumably increasing it in modulus. As soon as the metric inside the BH is not anymore Schwarzschild (but continues to be Schwarzschild outside). Note that the metric is destabilized and becomes singular for , so surely this is an upper limit where .
Yet another possible interpretation, that we disfavour, could be the following. Each BH has associated a given constant value of , hence of . Then Eqs. (48) would imply that after the emission of each graviton, the value of is readjusted. The problem with this interpretation is that it would require a new dimensionless magnitude () to characterize a Schwarzschild BH; something that most BH practitioners would probably find hard to accept.
6 Conclusions and outlook
We shall now conclude. The purpose of this work was to have some insights in the reformulation of a quantum theory of black holes in the language of condensed matter physics. The key point of the theory is to identify the black hole with a BoseEinstein condensate of gravitons.
We have conjectured the set of equations that play the role of the GrossPitaevskii nonlinear equation; they are derived from the EinsteinHilbert Lagrangian after adding a chemical potentiallike term. We have used a number of different techniques to analyze these equations when the perturbation (i.e. the tentative condensate) has spherical symmetry. The equations appear at first sight rather intractable, but by doing a perturbative analysis around the BH Schwarzschild metric at quadratic order (i.e. including the leading nonlinearities) we found that the chemical potential necessarily vanishes in the exterior of the BH. On the contrary, in the interior we have found two sets of solutions, one of them has to be discarded as producing a nonnormalizable result. The other one leads to a nonzero chemical potential in the interior of the BH that behaves as . Therefore, there is a finite jump on the chemical potential at the BH horizon. Surprisingly —or maybe not so— this solution modifies the coefficient of the and terms in the Schwarzschild solution, but not its functional form. Of course if there is no chemical potential at all, the modification vanishes, in accordance with well known theorems. However, if the former is non zero, the modification affecting the metric is also necessarily nonzero.
The perturbative analysis triggers a unique physical solution for the nonperturbative (exact) theory. This solution is characterized by a constant density of the wave function for the condensate. From the existence and knowledge of this solution, an unambiguous relation between the number of gravitons and the geometric properties of the BH is obtained. Hence, we find an expression for the Schwarzschild radius that involves an a priori independent and tuneable parameter, the dimensionless chemical potential (related to the meanfield wave function of the condensate ). We find this somewhat strange as this would be a new black hole parameter. Therefore we favour the universal value as discussed in the text. From this expression for the Schwarzschild radius, most relations obtained in [1, 2] can be rederived.
As should be obvious to the reader who has followed our discussion, our approach is somewhat different from the one developed in the initial papers by Dvali, Gómez and coworkers. We assume from the start the existence of a classical geometry background that acts as confining potential for the condensate. The fact that the functional form of the metric perturbation induced by the condensate is exactly the same as the original background, of course gives a lot of credence to the possibility of deriving the latter from the former in a sort of selfconsistent derivation. We have not explored this possibility in detail yet.
It is quite plausible that one could entertain the presence of condensates of other quantum fields inside the BH horizon (why only gravitons?). While we do not expect much of a conceptual difference, it would be very interesting to see the similarities and differences with the case of quantum gravitons.
Acknowledgments
We acknowledge the financial support of projects FPA201346570C21P (MINECO) and 2014SGR104 (Generalitat de Catalunya). This research was supported in part though the E.U. EPLANET exchange program FP7PEOPLE2009IRSES, project number 246806. The work of J.A. is partially supported by grants Fondecyt 1150390 and CONICYTPIAACT14177 (Government of Chile).
References
 [1] G. Dvali and C. Gomez, Black Hole’s 1/N Hair. Phys. Lett. 719, 419 (2013), arXiv:1203.6575; G. Dvali and C. Gomez, LandauGinzburg Limit of Black Hole’s Quantum Portrait: Self Similarity and Critical Exponent. Phys. Lett. B 716, 240 (2012), arXiv:1203.3372; G. Dvali, D. Flassig, C. Gomez, A. Pritzel, N. Wintergerst, Scrambling in the Black Hole Portrait. Phys. Rev. D 88, 124041 (2013), arXiv:1307.3458; G. Dvali and C. Gomez, Black Holes as Critical Point of Quantum Phase Transition. Eur. Phys. J. C 74, 2752 (2014), arXiv:1207.4059.
 [2] G. Dvali and C. Gomez, Black Hole’s Quantum N–Portrait. Fortsch. Phys. 61, 742 (2013), arXiv:1112.3359.
 [3] S.W. Hawking, Particle Creation by Black Holes. Commun. Math. Phys. 43, 199 (1975).
 [4] J.D. Bekenstein, Black Holes and Entropy. Phys. Rev. D 7, 2333 (1973); J.D. Bekenstein, Generalized second law of thermodynamics in blackhole physics. Phys. Rev. D 9, 3292 (1974); S.W. Hawking, Black holes and thermodynamics. Phys. Rev. D 13, 191 (1976).
 [5] J.D. Bekenstein, Phys. Rev. Lett. 28, 452 (1972); J.D. Bekenstein, Phys. Rev. D 5, 1239 (1972); J.D. Bekenstein, Phys. Rev. D 5, 2403 (1972); C. Teitelboim, Phys. Rev. D 5 (1972) 294; J. Hartle, Phys. Rev. D 3, 2938 (1971).
 [6] S.W. Hawking, The Unpredictability of Quantum Gravity. Commun. Math. Phys. 87, 395 (1982). For a review on information loss paradox, see, J. Preskill, Do Black Holes Destroy Information? International Symposium on Black Holes, Membranes, Wormholes, and Superstrings, Houston Advanced Research Center (1992), arXiv:hepth/9209058.
 [7] J. Alfaro, D. Espriu, D. Puigdomènech, Spontaneous generation of geometry in four dimensions. Phys. Rev. D 86 025015 (2012), arXiv:1201.4697; J. Alfaro, D. Espriu, D. Puigdomènech, The emergence of geometry: a twodimensional toy model. Phys. Rev. D 82, 045018 (2010), arXiv:1004.3664.
 [8] T. Regge and J.A. Wheeler, Stability of a Schwarzschild Singularity. Phys. Rew. 108, 1063 (1957);
 [9] G.D. Birkhoff and R. Langer, Relativity and Modern Physics (Harvard Univ. Press, 1923).
 [10] F. Dalfovo, S. Giorgini, L.P. Pitaevskii and S. Stringari, Theory of BoseEinstein condensation in trapped gases. Rev. Mod. Phys. 71, 463 (1999).
 [11] W. Pauli and V. Weisskopf, Helv. Phys. Acta 7, 709 (1934).
 [12] K. Hinterbichler, Rev. Mod. Phys. 84, 671 (2012), arXiv:1105.3735.
 [13] S. Weinberg, Cosmology (Oxford Univ. Press, 2008).
 [14] N. N. Bogoliubov, J. Phys. (USSR), 11, (1947).
 [15] S.W. Hawking, Breakdown of predictability in gravitational collapse. Phys. Rev. D 14 2460 (1976).
 [16] S. Das and R.K. Bhaduri, Class. Quant. Grav. 32, 105003 (2015), arXiv:1411.0753.