in an interesting saga of papers , dvali , gmez and others@xcite have put forward an intriguing suggestion : black holes ( bh ) could be understood as bose - einstein condensates of gravitons .
if this view is 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 @xcite could be understood as being due to evaporation ( leakage as the authors put it ) of the condensate . besides , this picture brings new ideas about the bekenstein entropy @xcite , the absence of hair @xcite , as well as the quantum nature of information storage and the possible information loss in bhs @xcite .
the main point of these works is that the physics of bh can be understood in this picture in terms of a single number @xmath0 , the number of gravitons contained in the bose - einstein condensate ( bec ) .
these condensed gravitons have a wave length @xmath1 , @xmath2 being the planck length ; they have a characteristic interaction strength @xmath3 ( reminiscent of meson interactions in large @xmath4 gauge theories ) and the leakage leads to a hawking temperature of order @xmath5 , equal to the inverse of @xmath6 .
the mass of the bh is @xmath7 and its schwarzschild radius therefore is given by @xmath8 , thus agreeing with the compton wavelength of the quantum gravitons @xmath6 , in accordance with the uncertainty principle that dictates @xmath9 in the ground state of the quantum system .
therefore , up to various factors of the planck mass everything is governed by @xmath0 , the number of intervening gravitons .
modulo some assumptions , all these results stem from the basic relation ( unless otherwise stated we work in units where @xmath10 ) @xmath11 that relates the number of gravitons @xmath0 , the mass of the gravitating object @xmath12 and its gravitational radius @xmath13 . for a schwarzschild bh @xmath14 .
we found these results intriguing and we set up to try to 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 @xcite 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 pre - existence 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 ; hence gravitons ( longitudinal gravitons that is ) . continuing with our semi - classical analogy ,
these gravitons have had a long time to thermalize and it is therefore expected that they can form a bose - einstein condensate .
of course , one may argue at once that gravitons do not have repulsive interactions , at least naively , and that therefore a bec is impossible to sustain .
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 @xcite ) .
then , repulsion is assured at some scale .
therefore we will try to identify a consistent set of equations describing a bec constructed on top of this classical field .
we will see that remarkably enough the characteristics of the 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 that determines its characteristics .
it is therefore 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 hartree - fock .
however we will stop short of doing so here .
in the subsequent we present the results of our investigations .
in what follows we shall adhere to the following notational conventions : @xmath15 will be the einstein tensor : @xmath16 , where the ricci tensor and scalar curvature are constructed with the metric @xmath17 in the usual way .
we will denote by @xmath18 the background metric that in our case it will invariably be the schwarzschild metric .
perturbations above this background metric will be denoted by @xmath19 , so @xmath20 .
we will use the minkowskian metric convention @xmath21 .
we will leave for later discussion whether the indices of @xmath19 have to be raised or lowered with the background metric @xmath18 or the full one @xmath17 . likewise for the corresponding volume element ( @xmath22 versus @xmath23 ) . as we will see the results depend to some extent on the prescription that is chosen .    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 gross - pitaevskii equation to this case .
the graviton condensate has necessarily to be described by a tensor field that in our philosophy has to be considered as a perturbation of the classical metric .
only spherically symmetric perturbations will be considered to keep things as simple as possible .
as we will see , every dimensionful quantity can be expressed as a function of the schwarzschild radius @xmath24 .
let us consider a spherically symmetric perturbation on the schwarzschild metric ; that is we want to consider a quantum state a condensate having @xmath25 .
then @xmath26 we will now expand the perturbed einstein tensor up to second order in @xmath27 .
this will allow us to retain the leading non - linearities .
four of the components of the einstein tensor are non - null but only three are linearly independent ; both angular components coincide when using the form @xmath28 for the einstein tensor .    taking into account
that @xmath29 , we have @xmath30 and analogous equations for the radial and angular components ( they are given in appendix [ einec ] ) .    it is well known that birkhoff s theorem @xcite guarantees that any spherically symmetric solution of the vacuum field equations must be static and asymptotically flat .
so this means that the solution of the equations @xmath31 with @xmath18 being the schwarzschild s metric , must necessarily be of the schwarzschild type again . in conclusion
, a perturbation of a schwarzschild metric in vacuum must end in another schwarzschild metric but taking into account a change in the parameters responsible of inducing the perturbation . stating this in other words , a spherically symmetric gravitational field should be produced by some massive object at the origin characterized by the schwarzschild radius , hence a perturbation in itself it
can only be produced by an alteration of this radius and therefore a change in its mass .
indeed , imposing the following ansatz for the components of the perturbation @xmath32 and substituting @xmath33 in ( [ gtt schw ] ) is possible to obtain the coefficients for each term of the radial expansion .
we obtain @xmath34 & + \frac{14a^5 + 25r_sa^4 + 20r_s^2a^3 + 10r_s^3a^2 + 5r_s^4a}{r^5}+ ... \,,\end{split}\ ] ] where @xmath35 .
it is obvious that @xmath36 is a solution as implies no perturbation of the metric at all .
but if @xmath37 the rest of the coefficients of the expansion are fixed once @xmath38 is given .
the same happens for @xmath39 using the radial equation ( [ grr schw ] ) given in appendix [ einec ] .
@xmath40 taking into account the accuracy of the expansion in @xmath27 the components of the metric become @xmath41 @xmath42 being the new mass associated to the bh and @xmath43 the perturbed schwarzschild radius .
now we will use similar techniques in the ` quantum ' case and as we will see the results are dramatically different .
however it is interesting to note that any perturbation @xmath27
( i.e. each ` graviton ' ) has associated a certain amount of energy that is reflected in a change of the bh mass .      perturbing the metric as @xmath44 ,
a completely covariant perturbative expansion for @xmath27 is obtained . by construction
these equations are non - linear ( we keep only the leading non - linearities , to keep the formalism simple and analytically tractable as we will see later ) .
we note that the usual gross - pitaevskii equation @xcite employed to describe bose - einstein condensates is a non - linear schrdinger equation ; i.e. an equation of motion that contains self - interactions ( hence the non - linearity ) , a confining potential for the atoms or particles constituting the condensate , and a chemical potential that is conjugate to the number of particles or atoms contained in the condensate .    among all these ingredients ,
the perturbed einstein equations already contain most of them .
they are already non - linear and while there is no confining potential explicitly included ( 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 .
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 grand - canonical ensemble by adding to the appropriate action a chemical potential term such as @xmath45 that is conjugate the quantity @xmath46 , which should be related to the graviton density of the condensate inside a differential volume element @xmath47 . taking into account that @xmath48 is not directly interacting with the background , but a field that propagates over the background
, it may seem reasonable to take the volume element as the one corresponding to the background . still , the new term is completely covariant as can be seen in appendix [ scovariant ] .
irrespective of this chemical potential term an appropriate action has to be chosen for the field @xmath19 .
it is only natural to select the einstein - hilbert action expanded around the classical schwarzschild background @xmath18 to describe fluctuations of the metric in a reparametrization invariant way .
one can then expand the einstein - hilbert lagrangian @xmath49 , with @xmath50 in powers of @xmath27 up to the desired order ( where @xmath51 is the planck mass ) . measure associated to the full
metric has to be taken ( i.e. @xmath52 ) however this is not the only possibility to obtain diff - invariant equations ; for instance @xmath53 also could yield a covariant set of equations . ]
an oddity of eq .
( [ fullaction ] ) is that , as defined , @xmath54 has dimensions of [ energy]@xmath55 .
this is because the field @xmath27 does not have the canonical dimension ; the field that has canonical dimensions is @xmath56 . after scaling the two powers of @xmath57 ,
the chemical potential would have dimensions of [ energy]@xmath58 .
still , these are not quite the correct dimensions for a chemical potential
. as we shall see below the ` natural ' quantity to be identified with a chemical potential will in fact be the dimensionless combination @xmath59 that will be conjugate to @xmath0 , the number of gravitons in the bec condensate with @xmath60 .
note that we assume that the chemical potential is @xmath61-dependent .
the actual equations of motion derived from the previous action may take slightly different forms depending on the choice of volume elements .
for the time being , let us take the simplest possibility @xmath62 where indices are assumed to be raised with the background metric only .
actually eq .
( [ eeq ] ) can be derived from the action we have chosen only if we neglect terms of the form @xmath63 .
we will return to this point later .
it is important to keep in mind the following : in the chemical potential part @xmath54 is considered in the grand canonical ensemble ; this implies that this magnitude is an external field and does not vary in the action ; in particular , for these equations of motion , it is independent of @xmath27 , so @xmath64 . under these considerations
this term resembles the chemical potential term of the gross - pitaevskii equation
.    this should be the way of introducing the chemical potential .
otherwise , were @xmath65 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 @xmath54 and this equation would nullify automatically the perturbation as well as the chemical potential itself , i.e. @xmath66 and @xmath67 .
therefore , the external field @xmath54 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 @xmath27 .
the main difference is that this constrains are implicit in a gr theory .
the restriction of @xmath54 is through the general covariance conditions for the action .
in other words , the diffeomorphism covariance implies the covariant conservation of the einstein tensor , and this in itself entails the same for the chemical potential term , @xmath68 this differential equation of motion for the latter potential 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 possible solutions into two regions : outside and inside the bh horizon .
we certainly expect that the graviton condensate will disappear quickly in the outside region as @xmath69 and , imposing this as a boundary condition , we will see that in fact the condensate is identically zero on this side of the horizon .
on the contrary , a unique non - trivial solution will be found in the interior of the bh .
in this section it will be shown that even after the inclusion of @xmath70 there is not other normalizable solution outside the black hole horizon than the trivial solution for the perturbation previously discussed .      in the limit @xmath71
, one needs only to keep the terms that are dominant in the einstein equations ( [ gtt schw ] ) , ( [ grr schw ] ) and ( [ gthth schw ] ) .
then the equations ( [ eeq ] ) associated to @xmath72 , @xmath73 and @xmath74 in this limit are , respectively will be neglected as we are interested in the simplest approximation to a gross - pitaevskii equation .
the form of these terms depends on which convention is chosen for the volume element and to raise / lower indices . ]
\frac{r_s}{r^3}h_{tt } -\frac{1}{r}{h_{tt } } ' -\frac{1}{r^2}h_{rr } + \frac{r_s}{r^3}{h_{tt}}^2 -\frac{1}{r}h_{tt}{h_{tt } } ' -\frac{r_s}{r^3}h_{rr}h_{tt } + \frac{1}{r}h_{rr}{h_{tt } } ' + \frac{1}{r^2}{h_{rr}}^2 = \mu\,h_{rr } \\[4ex ] \frac{2r_s}{r^3}h_{tt } + \frac{2}{r}{h_{tt } } ' + 2\,{h_{tt } } '' + \frac{2r_s}{r^3}h_{rr } + \frac{2}{r}{h_{rr } } ' + \frac{2r_s}{r^3}{h_{tt}}^2 -\frac{4r_s}{r^3}{h_{rr}}^2 + \frac{2}{r}h_{tt}{h_{tt } } ' + 2\,h_{tt}{h_{tt } } '' \\ + { h_{tt}}'^2 -\frac{2r_s}{r^3}h_{rr}h_{tt } -\frac{2}{r}h_{rr}{h_{tt } } ' -2h_{rr}{h_{tt } } '' + \frac{r_s}{r^2}{h_{rr}}'h_{tt } -{h_{rr}}'{h_{tt } } ' -\frac{4}{r}h_{rr}{h_{rr } } ' = 0   \,.\end{split}\ ] ]    as we want a vanishing perturbation at infinity , the following ansatz is imposed in the faraway region for its components and for the chemical potential : @xmath76 before proceeding , an additional consideration is needed ; namely , the perturbation must be a square - integrable function , i.e. @xmath77 as in the schwarzschild metric @xmath78 , the integral should fulfill @xmath79 for this to happen , if a power law ansatz at infinity is imposed , the exponents must obey @xmath80 . then the three equations in ( [ g schw infty ] ) are @xmath81 \left(\frac{r_s}{r } + \alpha\right)\frac{a}{r^{\alpha+2}}-\frac{b}{r^{\beta+2 } } + \left(\frac{r_s}{r}+\alpha\right)\frac{a^2}{r^{2\alpha+2 } } -\left(\frac{r_s}{r } + \alpha\right)\frac{a\,b}{r^{\alpha+\beta+2 } } + { \frac{b^2}{r^{2\beta+2 } } } = \frac{b\,c}{r^{\beta+\rho } }   \\[4ex ] \left(\frac{r_s}{r}+\alpha^2\right)\frac{2\,a}{r^{\alpha+2 } } + \left(\frac{r_s}{r}-\beta\right)\frac{2\,b}{r^{\beta+2 } } + \left(\frac{2r_s}{r}+3\,\alpha^2\right)\frac{a^2}{r^{2\alpha+2 } } -\left({\frac{r_s}{r}}-\beta\right)\frac{4\,b^2}{r^{2\beta+2 } } \\ -\left[\left(2+\beta\right)\frac{r_s}{r}+2\,\alpha^2+\alpha\,\beta\right]\frac{a\,b}{r^{\alpha+\beta+2 } } = 0 \,.\end{split}\ ] ] the requirement of the solution being square - integrable exclude that @xmath82 and @xmath83 could be zero .
this implies automatically the neglection of the terms proportional to the schwarzschild radius .
even more , as this numbers @xmath82 and @xmath83 are positive and larger than @xmath84 , the following asymptotic identities result @xmath85 @xmath86 from the angular equation ( [ gthth schw powerlaw ansatz o(1 ) ] ) , is mandatory for both terms to contribute at infinity .
if this is not the case , this would nullify @xmath38 ( or @xmath87 ) , and then , via both equations in ( [ gtt schw powerlaw ansatz o(1 ) ] ) , fix @xmath88 ( or @xmath36 ) and @xmath89 , i.e @xmath90 .
the competition of both terms of the angular equation is possible if and only if @xmath91 .
this condition modifies the angular equation ( [ gthth schw powerlaw ansatz o(1 ) ] ) as @xmath92 then from the second equation in ( [ gtt schw powerlaw ansatz o(1 ) ] ) , the following relationship can be read @xmath93 this automatically leads to a null chemical potential as the only possible solution for this region with the proposed ansatz .
this is due to the fact that if @xmath88 , as @xmath87 is proportional to @xmath38 , then @xmath38 is also zero .    in appendix [ ansatzinfty ]
we explore the possibility of exponentially vanishing chemical potential when @xmath94 with analogous conclusions .
likewise it can be seen that a much faster gaussian decay is also excluded ( not reported here ) .
finally , let us turn to the numerical integration of the equations in the outer part of the bh .
on the one hand , the following change of variables @xmath95 in eqs .
( [ eeq ] ) , maps the exterior part into a compact interval ; so as @xmath96 is transformed into the finite boundary @xmath97 ( and @xmath98 into @xmath99 ) .
the numerical integration will run between @xmath100 $ ] . on the other hand
, this redefinition allows rewriting the relevant equations in terms of dimensionless quantities , so they are universal as @xmath101 and @xmath70 are the only physical parameters .    in spite of these simplifications ,
the equations are difficult to deal with because they are singular , i.e. it is not possible to isolate the higher order derivative at @xmath97 or @xmath99 .
let us see how we can proceed .
it is well known that schwarzschild space - time is asymptotically flat , therefore , at infinity ( @xmath97 ) it is expected for the perturbation to vanish ( @xmath66 ) . as mentioned before ,
the equations are singular at @xmath97 and we have to settle for a point @xmath102 but @xmath103 to set up boundary conditions for the numerical integration procedure to start . likewise setting @xmath104 at an arbitrary value for @xmath105 , even if @xmath106 immediately triggers divergences at the first step of the routine .
therefore we take a  small " initial value for the perturbation in a point near infinite ( @xmath97 ) and then , decrease this initial value until the numerical routine becomes unstable . with the set of solutions for the components of @xmath27 obtained for each of the initial condition imposed , it is possible to compute the quantity @xmath107 , defined in eq .
( [ h2integral ] ) , and the integral extends over the exterior of the bh ( i.e. @xmath105 runs between 0 and 1 ) .    in figure [ nout ] the behaviour of the integral of @xmath108 outside the bh is presented when the initial condition for the integration approaches zero at a fixed point near infinity , @xmath109 , i.e. @xmath110 .
it is observed that this integral is invariably small and actually vanishes as the initial condition for @xmath33 nullifies .     outside the event horizon .
the closer the limit to an asymptotically flat space - time is ( i.e. decreasing the initial condition near infinity ) , the smaller this integral is . ]
this seems to confirm also numerically that the only possible solution in the outer region is the one with @xmath66 and therefore @xmath111 throughout the exterior of the bh .
in order to study the behaviour of our equations when the horizon is crossed , it is convenient to keep working in the @xmath112 variable and redefine the perturbations as @xmath113 in the new variables , the horizon corresponds to @xmath99 while @xmath114 is mapped onto @xmath115 . in the new perturbation the equations
become more compact . at the same order as before ,
the temporal , radial and angular equations read @xmath116 ( z-1)\gamma_{rr } + z{\gamma_{tt } } ' + z\gamma_{tt}\,{\gamma_{tt } } ' -z(z-1)^2\gamma_{rr}\,{\gamma_{tt } } ' -(z-1)^3{\gamma_{rr}}^2 = -x(z-1)\gamma_{rr } \\[3ex ] 2z(z-2)(z-1)\gamma_{rr } + z(z-2)(z-1)^2{\gamma_{rr } } ' -z(5z-2){\gamma_{tt } } ' -2z^2(z-1){\gamma_{tt } } '' \\
-2z(z-2)(z-1)^4\gamma_{rr}{\gamma_{rr } } ' + z(7z-2)(z-1)^2\gamma_{rr}{\gamma_{tt } } ' -z(5z-2)\gamma_{tt}\,{\gamma_{tt } } '   \\   + z^2(z-1)^3{\gamma_{rr}}'\,{\gamma_{tt } } ' -z^2(z-1){{\gamma_{tt}}'}^2 + 2z^2(z-1)^3\gamma_{rr}\,{\gamma_{tt } } ''   \\ -2z^2(z-1)\gamma_{tt}\,{\gamma_{tt } } '' = 0 \,.\end{split}\ ] ]      in order to get a feeling for possible solutions to these equations we consider their linearized approximation .
only terms linear in @xmath117 , @xmath118 and their derivatives are kept in the region @xmath119 .
that is @xmath120 ( z-1)\gamma_{rr}+{\gamma_{tt}}'= -x(z-1)\gamma_{rr } \\[3ex ] 2(z-1)\gamma_{rr}+(z-1)^2{\gamma_{rr}}'+3{\gamma_{tt}}'+2(z-1){\gamma_{tt}}''=0 \ , .
\end{split}\ ] ] let us now make the following ansatz @xmath121 the three equations take the following form @xmath122 b(z-1)^{\beta+1 } + a\alpha(z-1)^{\alpha-1}=-x b(z-1)^{\beta+1 } \\[2ex ] b(2+\beta)(z-1)^{\beta+1 } + a\alpha(1 + 2\alpha)(z-1)^{\alpha-1}=0 \,,\end{split}\ ] ] respectively .
it is worth noting that if @xmath123 , i.e. @xmath124 , all the terms have a contribution as we get closer to the event horizon and @xmath125 behaves as a constant . in this situation
, we obtain three equations for the coefficients @xmath126=0 \,.\ ] ] the system of equations is algebraic for the variable @xmath125 ( that as we have seen should behave as a constant as @xmath119 ) ; therefore , it is possible to eliminate @xmath125 by combining the temporal and radial equation . in the present ansatz
this leads to @xmath127 together with the angular equation , this determines the solutions up to a constant .
there are two possible solutions for this system    *   *     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 ) @xmath128 while if ii ) is taken as solution @xmath111 ( i.e. no chemical potential at all ) .
however these solutions do not vanish when @xmath119 ( as one could perhaps have hoped ) , thus casting doubts on the relevance of the linearized equations .
let us examine this point and let us only take into consideration the solution i ) , which is integrable . to make sure that this solution is little modified for @xmath119 when non - linearities are switched on , we substitute the solution back in the non - linear equation and verify that the most singular terms are still satisfied .
calculations are presented in appendix [ solinfulleq ] .
in addition , this exercise gives one interesting result : @xmath129 , where @xmath38 is so far arbitrary , also in sign .
note that at the linear level we got @xmath128 and the fact that the quadratic equation gives an @xmath130 correction to this result is consistent as we are implicitly assuming that @xmath131 , i.e. @xmath132 .      next we consider the region @xmath133 .
we will use another dimensionless variable for this purpose : @xmath134 . in these variables
the by now familiar three equations read in their linearized version @xmath135 @xmath136 let us now employ the ansatz @xmath137 while we leave @xmath138 as a free function .
this function can be eliminated by combining @xmath72 and @xmath73 , leading to one equation that together with @xmath139 leads to several constraints involving the parameters of the ansatz . of all possibilities , the one that gives the same weight to all terms in the equations is the one corresponding to @xmath140 , @xmath141 .
then , close to the origin it would appear as if @xmath142 it remains to be seen if this solution matches with the one obtained in the inner vicinity of the horizon .
so far we have seen that the solution i ) given in eq .
( [ behaviourw=1 ] ) satisfies not only the linearized approximation but also the full quadratic equations .
could it happen that the previous analysis misses some solution ? to try to answer this question we have performed a numerical integration of the basic equations ( [ eeq ] ) .    the integration will be done in the variable @xmath143 .
it is impossible to start integrating from either @xmath144 ( the bh center ) or @xmath145 ( the bh horizon ) exactly because as can be seen from eqs .
( [ gtt schw ] ) , ( [ grr schw ] ) or ( [ gthth schw ] ) the system of equations is singular at @xmath114 and @xmath98 ; i.e. one can not determine in either point the value of the higher order derivatives because precisely at the origin and at the horizon they are multiplied by factors that vanish or diverge .    in order to proceed we start the integration procedure at a small distance @xmath146 from the origin and following ( [ behaviourw=0 ] ) set @xmath147
we could have equally started the numerical integration routine from a point close to @xmath145 , i.e. in the vicinity ( the inner vicinity that is ) of the horizon .
the boundary conditions according to ( [ behaviourw=1 ] ) would then be @xmath148 the respective remaining degrees of freedom are the constants @xmath38 , @xmath149 and @xmath150 , and obviously they must be chosen in such a way that the solution near @xmath145 ( @xmath99 ) and the one near @xmath144 ( @xmath115 ) match .
it turns out that there are unique values of @xmath38 and @xmath151 that make both solutions match .
namely , @xmath152 .
even more , when the previous boundary conditions are assumed for the numerical integration , @xmath117 as well as the function related to the chemical potential @xmath125 , turns out to be constant in the interior of the bh .
the constant @xmath38 itself is arbitrary and is not determined by the structure of the equations .
changes in @xmath38 appear as an overall factor in the solution . for a certain value of @xmath38
the curves presented in figure [ gammasandxin ] are obtained .
the numerical solution reproduces the general features of the analytical study .
is represented by a black line , @xmath118 by a dashed line and the dimensionless chemical potential @xmath125 by a dotted line .
the results depend on an arbitrary constant as discussed in the text .
for the plot we have taken @xmath153 that implies @xmath154 and @xmath155 inside the bh . note that the numerical integration does not extend all the way to @xmath144 and @xmath145 because these are singular points in the differential equation . ]
it is most interesting to plot the results of the numerical integration for @xmath156 and @xmath157 ( i.e. considering @xmath158 and then raising one index , taking into account that in the current approach @xmath159 ) . in figure [ hcte ] , a plot of the numerical solution reveals that to a very good approximation throughout the interior of the bh horizon @xmath160     ( represented by a gray solid line ) and @xmath161 ( by a black dashed line ) , as obtained by numerical integration of the @xmath162 equations , are constant and equal . ]    actually raising an index with @xmath163 or @xmath164 makes no qualitative difference at all because the solution for the perturbation @xmath19 turns out to be proportional to the initial background metric element itself , as can be seen from the results obtained in the previous section for @xmath117 and @xmath118 ( even though these results in principle take into account the leading behaviour when @xmath119 ) . using the full metric to raise the index changes the numerical value ; yet , the qualitative conclusions do not change .      inspired by the previous numerical solution we reformulate our main equations derived at @xmath162 and write them in terms of the metric fluctuation with mixed indices ( one covariant , one contravariant ) .
let us write @xmath165 the full metric would become @xmath166 \,.\ ] ] if indices are lowered and raised with the full metric itself . otherwise , if the background metric is used to raise and lower indices , one obtains @xmath167 \,.\ ] ] and the corresponding volume elements are @xmath168 or @xmath169 , respectively .
if one assumes constant values for @xmath170 and @xmath171 , as the solution in the previous subsection suggests , the corresponding @xmath162 equations simplify enormously and reduce to the three following ones @xmath172 and @xmath173 therefore we are retrieving the solution found in the previous subsection : @xmath174 plays the role of the integration constant @xmath38 .    now let us see that the problem actually admits an exact solution if we assume that the solution is of the form @xmath175 with @xmath174 being a constant .
in fact , there are several possible solutions depending on how one chooses to treat the separation between background and fluctuation metrics .
if one chooses to consistently raise and lower indices with the background metric only , one has to plug eq .
( [ metricbackground ] ) into the einstein tensor and gets the following closed algebraic equation : @xmath176 that agrees with the @xmath162 perturbative solution found above . note that @xmath111 implies @xmath177 but not the other way round .
however , in deriving the equations of motion we have actually neglected a possible contribution of the corresponding volume measures , assuming that eq .
( [ eeq ] ) holds exactly .
if we decide to raise and lower indices with the full metric , it is natural to keep the volume element as the one given by @xmath22 , @xmath17 being the full metric .
the lhs of the equation of motion now reads @xmath178 ( this is obtained after substituting eq .
( [ metricfull ] ) into the einstein tensor ) .
the rhs of the equation also gets modified . on the one hand @xmath23
is replaced by @xmath22 , as mentioned , that depends on @xmath19 . on the other hand , because the indices of @xmath19 are raised with the full metric ( that depends on @xmath19 itself ) , deriving @xmath108 w.r.t .
@xmath19 is not totally immediate .
all things considered , the resulting equation of motion for a constant @xmath174 is now @xmath179 that can be interpreted as a mean field - like gross - pitaevskii equation for the condensate wave function @xmath174 . using the fact that @xmath180 , we ( re)obtain that @xmath125 is a constant and @xmath181 the previous results are valid to all orders in the expansion .
note that @xmath182 .    before moving on
, it is mandatory to make a comment on the covariant conservation of our equations . at the end of section 2.2
we have pointed out that the diffeomorphism invariance entails a differential equation for the chemical potential , namely ( [ diffeqchempot ] ) . in appendix [ diffecforchempot ] a short discussion of this equation
is presented .
we see that the actual relation between @xmath174 and @xmath125 depends from the precise way the fluctuations around the background metric are assumed to interact with the metric itself . in any case
the qualitative result that there is a normalizable solution that can be interpreted as the collective wave function of a graviton condensate remains .
a unique relation is obtained between this ( constant ) wave function and a dimensionless chemical potential .
in order to present a physical interpretation for the obtained results , undoing the redefinition of the components of the wave function ( i.e. @xmath183 ) is mandatory .
the left graph in figure [ hxh2forh ] presents the curves for the two covariant components of the perturbation .
the dimensionful chemical potential is plotted also for comparison . as @xmath184 is a constant function , @xmath185 and it is _ not null _ over the event horizon .
of course it is simpler to represent just @xmath174 or @xmath108 , which are just constants .
( black solid ) , @xmath33 ( dashed ) line and the dimensionful chemical potential @xmath70 ( dotted line ) . on the right plot ,
@xmath186 is shown . as already emphasized ,
a whole family of solutions are obtained , parametrized by a real arbitrary constant .
one can simply trade this constant for the value of @xmath174 .
]    it is immediate to see that the solution found is of finite norm .
this can be computed numerically using the numerical data into the equation ( [ h2integral ] ) .
taking into account that the perturbation @xmath27 is null from the event s horizon onwards , the endpoint on the integration limit can be chosen as the schwarzschild radius .
this way , the integral is @xmath187\simeq h^2 \ , r_s^3 \,.\ ] ] then the integral over the interior of the bh of @xmath108 is given by @xmath188 .
this quantity should , in a way that we will discuss below , be related to the total number of gravitons of the condensate .
@xmath108 would be related to its density in the bh interior that turns out to be constant .
note that @xmath19 is dimensionless and that in order to have a properly normalized kinetic term we have to divide the einstein - hilbert action by @xmath189 .    at this point
we should attempt to make contact with the results of @xcite .
possibly our more striking results are that the dimensionless chemical potential @xmath190 stays constant and non - zero throughout the interior of the bh and that so does the quantity @xmath108 previously defined .
it is totally natural to interpret @xmath125 as the variable conjugate to @xmath0 , the number of gravitons .
how could we determine the value of @xmath0 from our solution ?
even after using a properly normalized @xmath191 , the dimensions of ( [ integralh2 ] ) are not appropriate to deliver to us the value of @xmath0 .
as it is well known since the early days of quantum field theory @xcite there are no conserved currents or continuity equations for fundamental fields that are chargeless ( such as a real klein - gordon field ) .
this also applies to @xmath192 .
however , the quantity @xmath193 with @xmath6 being the graviton wave length , has the right ingredients to play the role of probability density in the present context . note that in the bec @xmath194 . if we assume this , then @xmath195 that agrees nicely with @xcite under the maximum packaging condition @xmath194 . recall that @xmath196 is a constant that is entirely determined by the value of the dimensionless chemical potential @xmath125 . the rest of relations can be basically derived from this .    indeed within the present philosophy we assign to each graviton in the bec an energy @xmath197 . while there is no formally conserved current ,
the total energy of the system will be given , by analogy with the case of a massless neutral scalar field , @xmath198 by @xmath199 if the energy @xmath200 is a constant and given by @xmath201 , then eq . ( [ h2integral ] )
has to be interpreted as the number of gravitons in the bec , i.e. the integral of the graviton density in the interior of the bh .
as seen above the dimensionless chemical potential @xmath125 has a rather peculiar behaviour .
it has a constant value inside the bh and it appears to be exactly zero outside .
this behaviour is summarized in figure [ muxvsw ] .
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 .
from this figure it is clear why gravitons are trapped behind the horizon : the jump of the chemical potential at @xmath98 would prevent the ` particles ' inside to reaching infinity .
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 @xmath202 : @xmath203 to estimate @xmath204 ( which is negative ) we can use geometrical arguments to determine the flux .
if we assume that for a given value of @xmath101 only one mode can get out ( as hypothesized above ) and that propagation takes place at the speed of light , elementary considerations we have to multiply the surface ( @xmath205 ) times the flux ; i.e. the density of the mode times the velocity , assumed to be @xmath206 in our units .
since the density of the mode is constant in the interior , it is just @xmath207 . ]
lead to @xmath208 this agrees with the results of @xcite for instance eq .
( 35) and yields @xmath209 .
the approximations made in the discussion could only modify the coefficient by a numerical factor of @xmath210 .    within this picture , several questions concerning
the long - standing issue of loss of information may arise as the outcoming state looks thermal @xcite but apparently is not ; or at least not totally so .
however we shall refrain of dwelling on this any further at this point .     in solid line and the dimensionless @xmath125 in dashed line , across the horizon of the bh ( dotted line ) .
in this case the jump at the horizon amounts to @xmath211 . ]
the main result from the previous rather detailed analysis is that the bh is able to sustain a graviton bec ( and surely similar bec s made of other quanta ) .
but is that condensate really present ?
our results do not answer that question , but if we reflect on the case where the limit @xmath212 is taken , without disturbing the bh geometry ( i.e. keeping @xmath101 constant ) this requires taking @xmath213 in a way that the ratio @xmath214 is fixed
. then one gets @xmath128 .
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 ) non - zero chemical potential when the physical system is expressed in terms of the grand - canonical ensemble . outside the bh , @xmath111 .
another way of reaching this conclusion is by taking a closer look to our exact equations in the previous section . if one sets @xmath111 then necessarily @xmath177 and one gets the classical schwarzschild solution everywhere , but the reverse is not true .
one can have @xmath177 but this does not imply @xmath111 .
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 .
@xmath174 is now non - zero ; it is directly proportional to @xmath215 , the number of gravitons , and the chemical potential departs from the value @xmath128 , presumably increasing it .
as soon as @xmath216 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 @xmath217 , so surely this is an upper limit .    yet
another possible interpretation , that we disfavour , could be the following . each bh has associated a given constant value of @xmath125 , hence of @xmath174 . then eqs . ( [ mistery ] ) would imply that after the emission of each graviton , the value of @xmath101 is readjusted .
the problem with this interpretation is that it would require a new dimensionless magnitude ( @xmath125 ) to characterize a schwarzschild bh ; something that most bh practitioners would probably find hard to accept . in either case
the metric in the bh interior is different from the schwarzschild one .
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 bose - einstein condensate of gravitons .
we have conjectured the set of equations that play the role of the gross - pitaevskii non - linear equation ; they are derived from the einstein - hilbert lagrangian after adding a chemical potential - like 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 non - linearities ) we found that the chemical potential necessarily vanishes in the exterior of the bh . on the contrary , in the interior
we have found two set of solutions , one of them has to be discarded as producing a non - normalizable result .
the other one leads to a non - zero chemical potential in the interior of the bh that behaves as @xmath218 .
there is a finite jump of the chemical potential at the bh horizon .
surprisingly or maybe not so this solution modifies the coefficient of the @xmath219 and @xmath220 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 is also necessarily non - zero . from the existence and knowledge of this solution , most relations obtained in @xcite can be rederived .
we have seen that there is some ambiguity in the way the effect of the background metric is taken into account and probably more work is needed to produce a unique fully convincing answer .
these ambiguities get reflected in the relation between the number of gravitons and the schwarzschild radius , where the value of the ( constant ) condensate enters .
the relation between the number of gravitons and the geometric properties of the bh involve an a priori independent and tuneable parameter , the dimensionless chemical potential @xmath125 ( related to the mean - field wave function of the condensate @xmath174 ) .
we find this somewhat strange as this would be a new black hole parameter .
therefore we favour the universal value @xmath128 as discussed in the text .
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 , gmez 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 self - consistent 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 .
we acknowledge the financial support of projects fpa2013 - 46570-c2 - 1-p ( mineco ) and 2014sgr104 ( generalitat de catalunya ) . this research was supported in part though the e.u .
eplanet exchange program fp7-people-2009-irses , project number 246806 .
the work of j.a . is partially supported by grants fondecyt 1150390 and conicyt - pia - act14177 ( government of chile ) .
99 g. dvali and c. gomez , black hole s 1/n hair .
lett . _ * 719 * , 419 ( 2013 ) ;    g. dvali and c. gomez , landau - ginzburg limit of black hole s quantum portrait : self similarity and critical exponent .
b * 716 * , 240 ( 2012 ) ;    g. dvali , d. flassig , c. gomez , a. pritzel , n. wintergerst , scrambling in the black hole portrait .
d * 88 * , 124041 ( 2013 ) ;    g. dvali and c. gomez , black holes as critical point of quantum phase transition .
phys . j. c * 74 * , 2752 ( 2014 ) .
g. dvali and c. gomez , black hole s
quantum n - portrait .
fortsch . phys . *
61 * , 742 ( 2013 ) .
j. alfaro , d. espriu , d. puigdomnech , spontaneous generation of geometry in four dimensions .
d _ * 86 * 025015 ( 2012 ) ;    j. alfaro , d. espriu , d. puigdomnech , the emergence of geometry : a two - dimensional toy model . _
* 82 * , 045018 ( 2010 ) .
s.w . hawking , particle creation by black holes .
43 * , 199 ( 1975 ) .
bekenstein , black holes and entropy .
d * 7 * , 2333 ( 1973 ) ;    j.d .
bekenstein , generalized second law of thermodynamics in black - hole physics .
d * 9 * , 3292 ( 1974 ) ;    s.w . hawking , black holes and thermodynamics .
d * 13 * , 191 ( 1976 ) .
bekenstein , phys .
* 28 * , 452 ( 1972 ) ;    j.d .
bekenstein , phys .
d * 5 * , 1239 ( 1972 ) ;    j.d .
bekenstein , phys .
d * 5 * , 2403 ( 1972 ) ;    c. teitelboim , phys .
d * 5 * ( 1972 ) 294 ;    j. hartle , phys .
d * 3 * , 2938 ( 1971 ) .
hawking , the unpredictability of quantum gravity .
phys . * 87 * , 395 ( 1982 ) .
for an excellent 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 : hep - th/9209058 .
birkhoff and r. langer , _ relativity and modern physics _ ( harvard univ . press , 1923 ) .
f. dalfovo , s. giorgini , l.p .
pitaevskii and s. stringari , theory of bose - einstein condensation in trapped gases .
phys . _ * 71 * , 463 ( 1999 ) .
w. pauli and v. weisskopf , helv .
acta * 7 * , 709 ( 1934 ) .    s.w . hawking , breakdown of predictability in gravitational collapse .
d * 14 * 2460 ( 1976 ) .
the two remaining components of the einstein tensor are @xmath221 = 0 \end{split}\ ] ] and @xmath222=0 \,.\end{split}\ ] ]    here @xmath174 stands for any of the angular coordinates : @xmath223 or @xmath224 .
it is worth to note that if a rescaling on the schwarzschild radius @xmath225 is done , fixing it to zero , then we end up with minkowski space - time in the three components of @xmath226 : ( [ gtt schw ] ) , ( [ grr schw ] ) and ( [ gthth schw ] ) .
in the following part , the invariance under the gauge group of diffeomorphisms of the chemical potential term in the action ( [ fullaction ] ) will be shown .
let s start by pointing out that both metrics , @xmath227 and @xmath228 , as well as the perturbation , transform covariantly .
in fact , the full metric @xmath44 , under an infinitesimal displacement d in the coordinates of the form @xmath229 x^\mu=\xi^\mu$ ] , changes as @xmath230    the same gauge transformation rules the background metric under the general coordinate transformation mentioned above .
this is , @xmath231    as the perturbation is defined by @xmath232 , under the same perturbation of the coordinate system , from the subtraction of ( [ fullmetrictransf ] ) and ( [ backgroundmetrictransf ] ) , the following equality is obtained @xmath233    the latter identity implies @xmath27 has also a lie derivative transformation . the fact that the perturbation transforms in the same way than the metric tensor , and the chemical potential @xmath54 behaves as a scalar under a general coordinate transformation ,
ensure automatically the general covariance of the theory .
let s change the ansatz and impose an exponential decreasing solution for the perturbations at infinity : @xmath234    as before for this type of ansatz , the linear order dominates in front of the higher ones . inserting this ansatz into a linearization of the einstein equations partially evaluated at infinity ( [ g schw infty ] ) @xmath235    @xmath236    @xmath237    as decreasing solutions
are expected , the parameters @xmath82 , @xmath83 and @xmath238 must be not null .
this allows one to neglect all the terms besides the principal ones inside each parenthesis . in the angular equation
there is no way for both terms to compete between each other @xmath239    in conclusion , depending on the fact whether if @xmath82 is bigger or not than @xmath83 , @xmath36 or @xmath88 .
for the first case , if @xmath36 , the temporal equation @xmath240    fixes @xmath88 .
being both , @xmath38 and @xmath87 null , the chemical potential disappears form the theory in this region . for the second case , if @xmath88 , the same equation @xmath72 nullifies @xmath38 ( as no null solutions for @xmath241 are expected )
however , no perturbation makes the chemical potential senseless . to sum up
, this ansatz implies that the only solution at infinity is a null perturbation , @xmath66 , and the disappearance of the chemical potential @xmath70 from the theory .
the solution of the linearized differential equations leads to the expressions @xmath242 ; @xmath243 . let us now plug this solution in the non - linear @xmath244 system and see if the solution survives or how could it get modified .
near @xmath245 the full temporal , radial and angular equations are @xmath246    @xmath247    @xmath248    respectively .
replacing the possible solution , @xmath242 this eliminates any derivative of @xmath117 and @xmath243 into these equations , we get in turn @xmath249 \longrightarrow\ , 3\,a - 2\,a + 5\,a^2 - 4\,a^2 = a+a^2 = ( 1+a)a = -x\,a \end{split}\ ] ]    @xmath250 \longrightarrow\ , -\frac{a+a^2}{(z-1 ) } = -\frac{(1+a)a}{(z-1 ) } = \frac{x\,a}{(z-1 ) } \end{split}\ ] ]    @xmath251 \longrightarrow\ , -\frac{2\,a}{(z-1 ) } + \frac{2\,a}{(z-1 ) } -\frac{4\,a^2}{(z-1 ) } + \frac{4\,a^2}{(z-1 ) } = 0 \,.\end{split}\ ] ]    therefore , quite surprisingly , the linear solution is still an exact solution of the non - linear quadratic differential equations and we conclude that @xmath252    are solutions .
the only difference is that the dimensionless chemical potential is at this order @xmath253 rather than @xmath128 as suggested by the linearized approximation .
in this appendix , a derivation of the differential equation for the chemical potential @xmath70 is presented .
we are going to derive this equation for a theory where only the background metric is considered for the integration measure and to raise and lower indices .
covariant conservation of the einstein tensor implies automatically the same for the rhs of our equations @xmath254
\end{split}\ ] ]    in the last equality , the covariant derivative is expressed in terms of the christoffel symbols of the theory . performing the sum in @xmath255 because @xmath256 appears
always in the perturbation and only exists in the temporal and radial components , and doing the same procedure with @xmath257 the following relation is obtained @xmath258 \,\end{split}\ ] ]    where we have taken into account that there is only radial derivatives .
the latter equation is a set of 4 equations , @xmath259 .
the only non trivial equation is the radial one , i.e. @xmath260 . dropping all the null christoffel symbols and the off - diagonal elements of the @xmath27 matrix , makes the equation reads as @xmath261   \end{split}\ ] ]