Black Hole Evaporation, Quantum Hair and Supertranslations
Cesar Gomez Sebastian Zell
hep-th gr-qc hep-ph quant-ph




Polarization Transfer Observables in Elastic Electron-Proton Scattering at 𝐐^2= 2.5, 5.2, 6.8 and 8.5 GeV^2
    L. Zhu
    December 30, 2023
===========================================================================================================






 In a black hole, hair and quantum information retrieval are interrelated phenomena. The existence of any new form of hair necessarily implies the existence of features in the  quantum-mechanically evaporated  radiation.   Therefore, classical supertranslation hair can be only distinguished from global diffeomorphisms if we have access to the interior of the black hole. Indirect information on the interior can only be obtained from the features of the quantum evaporation. 
 We demonstrate that supertranslations (T^-,T^+) ∈ BMS_-⊗ BMS_+ can be used as bookkeepers of the probability distributions of the emitted quanta  where the first element describes the classical injection of energy and the second one is associated to quantum-mechanical emission. However,  the connection between T^- and T^+ is determined by the interior quantum dynamics of the black hole.  We argue that restricting to the diagonal subgroup is only possible for decoupled modes, which do not bring any non-trivial information about the black hole interior and therefore do not constitute physical hair.  It is shown that this is also true for gravitational systems without horizon, for which both injection and emission can be described classically. Moreover, we discuss and clarify the role of infrared physics in purification.








§ INTRODUCTION AND SUMMARY




  §.§.§ The Puzzle of Black Holes

A black hole is an extraordinary physical system. While in a classical theory, it is extremely simple for an outside observer, as a consequence of the no-hair-theorem (see <cit.>), its internal quantum complexity measured by the Bekenstein-Hawking-entropy <cit.> N=M^2/M_p^2 is enormous. Both properties are obviously interrelated. The black hole entropy appears because many different matter configurations can collapse into the same black hole geometry. The no-hair-theorem prevents an outside observer from resolving these differences which remain hidden behind the horizon. Quantum-mechanically, the black hole evaporates <cit.> and unitarity requires that along the evaporation process the black hole should deliver back the information which was classically hidden in its interior <cit.>. This means that although the classical metric has no hair, the evaporation products should have features which compensate for this lack of information. In other words, the quantum radiation emitted during the black hole evaporation should carry the same information which the classical no-hair-theorem prevents us to extract from the geometry. 

In the last twenty years it has become popular to use the AdS/CFT correspondence as strong indication of the unitarity of black hole evaporation. However, this hope will not be fulfilled until counting with the CFT dual of a small evaporating black hole has been achieved. More generally, we shall argue in this note that the solution to the evaporation problem requires to have a microscopic model of the black hole as a quantum system –  whether obtained from AdS/CFT or differently. In <cit.>, we have developed such a model that, among other things, indicates the existence of forms of quantum hair effects of order 1/N. Moreover, and in a model independent way, it is easy to see that taking into account the change of the black hole mass due to Hawking evaporation leads to deviations from featureless emission on precisely this order of magnitude <cit.>.



  §.§.§ Classical BMS-Hair

Recently, a new way to attack the problem has been suggested <cit.> based on asymptotic BMS-symmetries <cit.>. This approach has received widespread attention (see <cit.>). In particular, a potential new form of classical hair for a black hole has been proposed <cit.>.
The idea is simply the following. One starts with a black hole of mass M and injects an energy μ in the form of incoming radiation with some angular features.[All quantities will be properly defined at the beginning of section <ref>.]

This incoming radiation can be associated with a supertranslation in BMS_- which we denote by T^-. Classically, the resulting system is a black hole with total mass M+μ but supertranslated by T^-. One can do the same construction with identical μ but with different angular features, different supertranslations T_i^-, to obtain a family of different metrics all of them with the same ADM-charges. Thus, it seems that one can indeed define classical hair if all these metrics sharing the same ADM-charges are physically inequivalent. At the classical level, this means that those metrics are not just the same metric written in different coordinate systems, that they are not related by a globally defined diffeomorphism. 
 As we shall elucidate, the problem with this form of classical hair is that for an observer outside, there is no way to decide if all these metrics are different or simply the same metric in different coordinates. In order to decide that, the observer needs to have information about the interior of the black hole. In summary, defining hair by means of the classical gravitational memory associated to some incoming radiation is only operative if somehow we can have extra information about the memory effects in the interior of the black hole which is, in a different guise, the essence of the no-hair-theorem.


Fortunately, there is an indirect way to decide from the outside whether two black hole metrics defined by injecting the same amount of energy but with different angular features are physically different or not. We can just wait until the black holes emit some radiation and compare the radiation produced by the two black holes. For simplicity, we restrict ourselves in our discussion to a pure theory of gravity in which only gravitational radiation can be emitted. The corresponding process is depicted in figure <ref>, where we distinguish the classical, semi-classical and purely quantum contributions. The first thing to be noticed is that this test is purely quantum in the sense that only quantum-mechanically, the black hole can emit radiation. The second thing is that the information we can get on the emitted quantum radiation by actual measurements is necessarily encoded in the form of probability distributions. Thus, if those black holes defined by different T^i_- are indeed different, we should expect that the corresponding quantum probability distributions are also different. 





  §.§.§ Insufficiency of Calculation in Classical Background Metric

  It is natural to expect that this difference has a non-trivial projection on deviations from isotropy, that the emitted quanta carry angular features. Then we obtain quantum probability distributions P_i(θ,ϕ) from the measurement of the radiated quanta. We can use those to define a classical supertranslation T_i^+. By that we simply mean a classical supertranslation with the flow of emitted radiation determined by the measured quantum probability functions P_i(θ,ϕ). In this sense, the former experiment produces a set of couples (T_i^-, T_i^+) where the first supertranslation in BMS_- is classical and the second one in BMS_+ is determined by the quantum probability distribution. From this point of view, if the classical T_i^- really implants hair, then the quantum T_i^+ should be non-trivial, contain spherical harmonics with l≥ 2.  The crucial point is that this behavior cannot be achieved by the standard Hawking computation performed in a supertranslated Schwarzschild metric, as pair creation in the background vacua defined by the near horizon geometry. The reason is that the supertranslation acts as a diffeomorphism near the horizon and does not change the local geometry. Therefore, it does not suffice if the P_i only depend on the injected radiation and the geometry of the black hole. Instead, they must also depend on its internal dynamics. 

We can make the argument a bit more quantitative and assume that from the whole energy μ injected a fraction μ̃ is associated to angular features. This means that the part of the incoming classical flow ℱ_in with angular labels l≥ 2 contributes to ∫[2]Ωℱ_in with a value equal to μ̃.   Clearly, μ̃=0 would correspond to the injection of featureless radiation.[In this case, the associated supertranslation will not support angular features and will only project on the l=0,1 spherical harmonics.]

 In order to parametrize how the P_i depend on μ̃ and the internal structure of the black hole, we shall use the typical number of quantum constituents of the black hole. In this sense, we expect P_i(θ,ϕ; μ̃, N) where the label i refers to the dependence on the incoming T_i^- and where we identify the number of quantum constituents of the black hole with the entropy N. 

Then the natural dimensionless parameter measuring the dependence of P_i on the internal structure is μ̃/√(N) where √(N) is the black hole mass in Planck units. In this setup, angular features in the evaporation, finite N effects in P_i(θ,ϕ; μ̃, N), depend on the black hole microscopic model. Those will define a couple (T_i^-, T_i^+)  generically not in the diagonal subgroup BMS_0 of BMS_- ⊗ BMS_+.  Thus, the first important message of our note will be that the information about T_i^- cannot determine the quantum probability distribution T_i^+, we cannot predict the quantum probability distribution solely from the incoming radiation implanting the hair. 



  §.§.§ Subleading Soft Modes

 We can investigate how this situation changes in the semi-classical limit M→∞, N→∞, in which the Hawking computation becomes exact. In this case, the energy associated to features becomes zero so that angular features can only be encoded in zero-energy modes. The effective decoupling of these modes will lead to a P_i identical to the incoming radiation. This produces couples (T_i^-, T_i^+) in the diagonal subgroup BMS_0. In more concrete terms, the lim_N→∞ P_i(θ,ϕ; μ̃, N) will only capture local horizon physics or zero-energy modes.[It is important to stress that the pseudo Goldstone-Bogoliubov modes identified in <cit.> are not equivalent to near horizon diffeomorphisms and consequently are good microscopic candidates to describe the low energy effective changes of the microstate of the black hole during the process of absorption and evaporation.]



This brings us to our second point, namely how the actual features of the quantum probability distribution P_i depend on infrared physics.[See <cit.> for a recent suggestion for purification by infrared modes.] 

 We know that in gapless theories such as gravity, evaporation interpreted as a S-matrix process has a zero probability amplitude without any accompanying soft gravitons. In order to obtain a finite answer, one has to include the emission of a certain class of soft radiation, namely IR-modes. However, this fact by no means implies that this companion radiation should carry the angular features that we need to purify the evaporation.  On the contrary, we know from infrared physics that IR-radiation is only sensitive to the initial and final scattering states. It is independent of the details of the process or in our case of the microscopic details of the black hole, cannot resolve the microstate.

 Independently of the question to what extent IR-radiation and hard quanta are correlated, we can quantitatively estimate the amount of information we could lose when we integrate over unresolved IR-modes. From well-known results of infrared physics it follows that their number only grows logarithmically with the resolution scale ϵ, n_soft∼ -lnϵ. However, what we have discussed implies that the natural resolution scale of features should be ϵ∼ 1/N. Thus, the second important message of our note is that unresolved IR-modes cannot account for the  bulk of information in black hole evaporation, but could only contribute as a subleading logarithmic correction. The part which carries features is the part of the radiation that can be resolved and that depends not on the infrared divergences but on the inner structure of the black hole, or in scattering language, on the details of the scattering process. 
 A possible candidate is soft non-IR radiation, which is independent of infrared divergences. As it should be, non-IR radiation depends on the details of the scattering process so that it cannot be predicted without a microscopic theory of the black hole. 



  §.§.§ Summary and Outline

In summary, non-trivial hair can be only defined by couples (T_i^-, T_i^+) ∈ BMS_-⊗ BMS_+ where the element in BMS_- is classical  and carries some finite energy and the element in BMS_+ is defined as a bookkeeper of the quantum probability distribution of the radiated quanta. What concrete element T_i^+ is associated with a given T_i^- cannot be derived solely from the classical geometry, but depends on the internal quantum structure of the black hole.  This non-trivial mapping is precisely what makes the quantum hair informative. For a system with non-trivial dynamics, it is therefore impossible to restrict to a subgroup of BMS_- ⊗ BMS_+.
Predictivity on this quantum output can be only achieved in the zero-energy (or equivalently N=∞ limit) where we only get elements in the diagonal subgroup BMS_0.[In <cit.> and <cit.> it is suggested to constraint the potential values of T_i^+ using an infinite set of conserved charges. Imposing these conservation laws makes the corresponding S-matrix completely insensitive to the internal structure of the black hole and consequently, in the language we are using here, can only capture unobservable zero-modes.]

But since the soft modes are decoupled once the infrared divergences of the theory are properly taken into account <cit.>,  they cannot lead to observable features.[This decoupling of soft modes is a quantum effect that should not be confused with the existence, for instance in asymptotically Minkowski space time, of a non-trivial family of asymptotically flat connections defining a representation of the BMS-group (see <cit.> and references therein). This multiplicity of classical inequivalent vacua is quantum-mechanically reabsorbed in the cancellation of infrared divergences.] 


 The outline of the paper is as follows. In section <ref>, we first recap some properties of BMS-gauge. In particular, we show how angular features of radiation define a supertranslation, which can be measured as a memory effect. Moreover, we discuss the role of soft modes. Then we use a combination of injected and emitted radiation of the same total energy to define Goldstone supertranslations as element (T^-,   T^+) ∈ BMS_-⊗ BMS_+. 
In section <ref>, we first concentrate on a gravitational system without horizon, which we shall call planet for concreteness, and show how we can use Goldstone supertranslations to change its angular distribution of mass. In doing so, the key point is that it is impossible to infer T^+ from T^- unless one knows the internal dynamics of the planet. Moreover, we highlight the importance of angular features by showing that it is impossible to determine the angular mass distribution of the planet without access to its interior. Subsequently, we apply Goldstone supertranslations to a black hole. We demonstrate how supertranslations can be used as bookkeeping tool for the emitted quanta. However, without knowledge of the microscopic dynamics of the black hole, they have no predictive power.  We also point out how we can use Page's time to estimate the magnitude of deviations from featureless evaporation.
After concluding in section <ref>, we provide a more detailed discussion of IR-physics in appendix <ref>. In appendix <ref>, we discuss the matching of the supertranslation field in advanced and retarded coordinates and finally we explicitly calculate a Goldstone supertranslation of a planet in appendix <ref>.





§ QUANTUM HAIR




 §.§ Recap of BMS-Gauge and Memory Effect




  §.§.§ Retarded Coordinates

We first recap some properties of BMS-gauge, which is defined by the four gauge conditions <cit.>

    g_11 = g_1A = 0  ,     g_AB = r^2 sin^2θ ,
 
where A,B, … = 2,3. Typically, BMS-gauge is used to study a spacetime asymptotically, for r→∞, but it is possible to extend the metric to the bulk by imposing the conditions (<ref>) to all orders in 1/r. In a typical situation, however, a metric in BMS-gauge does not cover the whole spacetime.

 A metric in BMS-gauge exists both in retarded time u, which is suited to describe outgoing radiation, and in advanced time v, which is suited to describe incoming radiation. The matching between these two metrics will be crucial for our treatment. Explicitly, an asymptotically flat metric in retarded time takes the form <cit.>:

    s^2=   (-1 + m^+_B/r + O(r^-2)) u^2 - (2+O(r^-2))  u  r 
       + r^2 	(γ_AB + C_AB^+ r^-1 + O(r^-2))  x^A  x^B + O(r^-2)  x^A  u ,

where the metric on the sphere has to fulfill the requirement g_AB = r^2sin^2 θ. Here m_B^+ is the Bondi mass, γ_AB the standard metric on the sphere and 

    C^+_AB = (2 D_A D_B - γ_ABD^2 )C^+

is determined by the supertranslation field C^+, where D_A is the covariant derivative on the sphere. It is helpful to expand the supertranslation field in spherical harmonics. Then the mode l=0 represents a time shift and the mode l=1 corresponds to spatial translations. Therefore, all modes with l≥ 2 define proper supertranslations.  Metrics with different values of C^+ are connected via asymptotic diffeomorphisms, the choice of the supertranslation field constitutes a residual gauge freedom of BMS-gauge. These diffeomorphisms are the famous supertranslations. Therefore, we can define a supertranslation T^+ by the change it induces in the supertranslation field:

    T^+ := Δ C^+  .

	

In order to analyze the effect of supertranslations, we will need the constraint equation G_00=8π G T_00, whose leading order reads in BMS-gauge:

    ∂_u m^+_B =  1/4GD^2(D^2+2) ∂_u C^+ -ℱ_out ,

where 

    ℱ_out = 1/8(∂_u C^+_AB)(∂_u C^+AB) + 4πlim_r→∞ (r^2T_uu)

is the total incoming null energy, composed of gravitational waves (first summand) and other forms of gravitating energy (second summand). 



  §.§.§ Advanced Coordinates

The situation in advanced coordinates is very similar. The metric takes the form

    s^2=   (-1 + m^-_B/r + O(r^-2)) v^2 + (2+O(r^-2))  v  r 
       + r^2 	(γ_AB + C_AB^- r^-1 + O(r^-2))  x^A  x^B + O(r^-2)  x^A  v ,

where the supertranslation field and the supertranslations in advanced coordinates are defined as in (<ref>) and (<ref>). The constraint equation becomes

    ∂_u m^-_B =  1/4GD^2(D^2+2) ∂_u C^- +ℱ_in ,

where ℱ_in is the incoming energy, in analogy to (<ref>).



  §.§.§ Measurement of the Supertranslation Field: Memory Effect

As already discussed, one can change the value of the supertranslation field by a diffeomorphism. Therefore, it follows by general covariance that the value of the supertranslation field cannot have in general any experimental implication. However,  since it corresponds to physical outgoing or ingoing radiation, the difference of the supertranslation field at different times does have experimental implications: It describes the memory effect caused by the radiation, a permanent displacement of test masses after the radiation has passed <cit.>. 

We will restrict ourselves to a simple situation in which we start with some stationary metric g^1_μν and we finish in a different stationary metric g^2_μν. In between, there is a radiation epoch, ℱ_in/out only has support during this time span. Asymptotically on 𝒥^±, the process defines a non-stationary metric interpolating between g^1_μν and g^2_μν which should be a solution to the Einstein equations.

 Since Birkhoff's theorem implies that we can set ∂_A m_B^±  = 0 in a stationary metric, we can single out the zero-mode from (<ref>) by integrating over the sphere:

    μ^+ = -∫u∫[2]Ωℱ_out/4π ,

where we first consider retarded time and μ^+ = m^+_B, 2- m^+_B, 1 is the total change of Bondi mass due to the radiation epoch. This formula shows explicitly that the Bondi mass m^+_B is monotonically decreasing, it measures the energy which has not yet left the bulk.  Defining the emitted energy with non-trivial angular distribution as Δℱ̃_out := ∫uℱ_out - μ^+, the constraint (<ref>) becomes

    0 =  1/4GD^2(D^2+2) T^+ - Δℱ̃_out .

 Thus, angular features in the outgoing radiation induces a supertranslation T^+=Δ C^+. Note that it is independent of the total emitted energy μ^+.

 In advanced coordinates, we get from the constraint (<ref>):

    μ^- = ∫v∫[2]Ωℱ_in/4π .
 
The advanced Bondi mass m^-_B is monotonically increasing, it measures the energy which has already entered the bulk. Defining Δℱ̃_in := ∫vℱ_in - μ^-, the constraint (<ref>) becomes

    0 =  1/4GD^2(D^2+2) T^- + Δℱ̃_in .

 This formula implies that an advanced supertranslation T^- tracks angular features in the incoming radiation.




 §.§ Goldstone Supertranslations

As already pointed out, we shall define hair on the basis of scattering processes where some injected gravitational energy is radiated back by the system. The hair will be encoded in the angular features of the injected radiation and the outgoing radiation. In this sense, we define hair as a typical response function. Through these formal scattering processes we define a map relating gravitational systems, black holes or planets, in different states sharing the same values for all the ADM-conserved quantities. We denote this induced map a Goldstone supertranslation since it relates states which are degenerate in energy. Note that this scattering definition of hair is tied to the mechanism of radiation whatever it could be. 



  §.§.§ Relationship to Antipodal Matching

 As a first step, it is important to discuss whether there are general constraints on this scattering process. Namely, it has been suggested in <cit.> that any gravitational S-matrix in an asymptotically flat spacetime must satisfy the following relation for an arbitrary initial quantum state |α⟩:

    S T^- |α⟩ = P(T^+) S |α⟩ ,
 
where (T^-, T^+) ∈ BMS_- ⊗ BMS_+ and P is the antipodal map on the sphere:

    P(T^+)(θ, φ) = T^+(P(θ, φ))  .

 Imposing this invariance implies that if a matrix element ⟨β| S |α⟩ is non-vanishing, then

    ⟨α| T^- |α⟩ = ⟨β|P(T^+) |β⟩ .

This means that the memory effect of the outgoing wave, parameterized by T^+, must match the memory effect of the incoming wave, parameterized by T^-, antipodally at each angle. Because of the constraints (<ref>) and (<ref>), this is equivalent to the statement that the outgoing energy Δℱ̃_out matches the ingoing energy Δℱ̃_in antipodally at each angle, in particular that Δℱ̃_out is fully determined in terms of Δℱ̃_in.

This criterion has a very interesting connection to IR-physics. As discussed in appendix <ref>, we know that in a gapless theory such as gravity most process in which no soft modes are emitted have zero probability <cit.>. In order to obtain a finite answer, one has to include the emission of a certain class of soft radiation, namely IR-modes. 
 The sole exception are processes for which the kinematical factor B_α,  β defined in <cit.> is zero. The crucial point is that this happens if and only if the ingoing energy matches the outgoing energy antipodally at each angle, as discussed in detail in <cit.>. Thus we conclude that[We will elaborate on this point in <cit.>.]

    ⟨α| T^- |α⟩ = ⟨β|P(T^+) |β⟩  ⇔    B_α, β=0  .

This means that restricting to processes which fulfill the condition (<ref>) is equivalent to only considering processes that are IR-finite even without including IR-emission.

A priori, there is nothing wrong with solely considering such processes. However, they only form a set of measure zero of those processes that occur in reality. Namely any realistic scattering is accompanied by the emission of soft IR-modes. Once we include soft IR-emission, we know that all processes – with an arbitrary non-zero value of B_α,  β – are IR-finite. Thus in reality, any process can occur, also ones that do not fulfill the antipodal matching condition (<ref>). For this reason, we will not restrict ourselves to processes that obey (<ref>).




  §.§.§ Role of Soft IR-Gravitons

 Since we consider processes that include the emission of soft IR-modes, it is natural to ask if those modes could carry information about the black hole state and if they could even suffice to purify black hole evaporation. This is only possible if two conditions are fulfilled. First, IR-modes would have to be sensitive to the microstate of the black hole. We expect this not to be the case since they only depend on the initial and final scattering state, but not on the details of the process.
While we leave the above question for future work, we now focus on the second condition, namely that the number of resolvable IR-modes would have to be big enough to be able to carry the whole black hole entropy.

In contrast to the proposal made in <cit.>, we argue that generic properties of IR-physics imply that this is not the case. As follows from equation (<ref>),  the number of unresolved soft modes scales logarithmically with the IR-resolution scale. Thus, when we lower the energy scale of resolution from ϵ_1 to ϵ_2, the number of additional IR-modes that we can resolve is:

    n^res_soft∼ B_α,  βlnϵ_1/ϵ_2 ,

where B_α,  β∼ G s is determined by the energy scale s of the process. 
We apply this formula to the single emission of a Hawking quantum of energy r_g^-1. It will be crucial in this argument that Hawking radiation gets softer for bigger black holes. The worst resolution scale compatible with observing this process is ϵ_1=r_g^-1. The key point is that the resolution scale in this process cannot be arbitrarily good. Namely, it is set by the time-scale of the process, ϵ_2∼ t_b-h^-1. Since the life-time of a black hole scales as t_b-h∼ N r_g, we get 

    n^res_soft≲1/Nln N  .

Thus, after the black hole has evaporated by emitting N Hawking quanta, the maximal entropy contained in the soft IR-modes is

    S_soft≲ln N  .

Independently of the question whether IR-modes are strongly correlated with the Hawking quanta, this shows that they cannot account for the whole entropy of the black hole, but could only give a logarithmic correction. Of course, this leaves open the possibility that non-IR soft modes could account for the bulk of black hole information. However, since they are independent of IR-divergences and accompanying dressing tools, the results of infrared physics do not constrain them.




  §.§.§ Role of Zero-Energy Gravitons


 Finally, we briefly discuss the role of zero-energy gravitons. To this end, we consider the process of a Goldstone supertranslation in the limit of zero energy injected and zero energy radiated. This is equivalent to the scattering with a graviton of zero energy.  Since those carry no energy, they cannot emit IR-modes and therefore obey the antipodal matching condition (<ref>). This fact simply reflects the well-known decoupling of soft modes <cit.>. The physical interpretation of this phenomenon is that any bulk configuration is transparent for decoupled soft modes so that the energy profile of the outgoing wave is antipodally related to that of incoming energy. 

But when the emitted/injected radiation does not carry energy, μ^±=0, then the constraint (<ref>) (or respectively (<ref>)) implies that ∫[2]Ωℱ_in/out =0. Since ℱ_in/out represents real gravitational radiation, it follows by the requirement of positive energy that ℱ_in/out=0. Thus, only supertranslations with D^2(D^2+2) T^± = 0 can occur in such a zero-energy process. This means that only the angular modes l=0,1 are left, time- and space-translation. Hence zero-energy radiation cannot lead to a physical memory effect that is observable in finite time. In other words, it is impossible to measure a zero-energy graviton in finite time.


The upshot is that predicting T^+ from the knowledge of T^- is only possible for zero-energy modes. Those are, however, unphysical since they cannot be measured in finite time. So we will only consider processes of non-zero energy in our paper. As explained, it is not possible for them to constrain or even predict T^+ from T^- without detailed knowledge of the dynamics in the bulk. The response function, which determines T^+ in terms of T^-, is trivial only for modes of zero energy. 

 



  §.§.§ Physical Hair With Non-Zero Energy

So from here on, we consider the case where after we inject radiation ℱ_in of non-zero total energy μ, the system radiates back the same total amount of energy, but with a possibly different distribution ℱ_out.[We recall that we restrict ourselves for now to a pure gravitational radiation, which propagates along null geodesics. Therefore, all emitted energy is bound to reach future null infinity 𝒥^+.]

While such systems are of course special, we will see that black holes can be one of them. This is a zero-energy process in the sense that the total energy of the system does not change. Thus, this process, which is depicted in figure <ref>, constitutes a transformation between degenerate systems and therefore defines hair.


As far as we reduce ourselves to gravitational radiation, we can generically describe this process in terms of two supertranslations: At 𝒥^-, T^- is determined by the angular distribution Δℱ̃_in of incoming energy according to the constraint (<ref>) and at 𝒥^+, T^+ follows from the angular distribution Δℱ̃_out of outgoing energy via the constraint (<ref>).  Thus, the whole process is associated to an element (T^-,  T^+) ∈ BMS_-⊗ BMS_+. It describes a zero-energy transition which interpolates between two spacetimes of the same total energy, but contrary to the case of a zero-energy mode, this transformation is non-trivial and it is not decoupled. 

 It is crucial to note that for an asymptotic observer, T^- and T^+ are independent. Whereas one is free to choose T^- by preparing an appropriate incoming radiation, T^+ is sensitive to the properties of the system in the bulk. In other words, T^+ is a response of the system which does not only depend on the ingoing radiation, parameterized by T^-, but also on the state of the system and its particular dynamics, which are not entirely visible asymptotically. In particular, there is no reason why (T^-,  T^+) should be in any subgroup of BMS_- ⊗ BMS_+.






  §.§.§ Coordinate Matching



 In order to compare ingoing and outgoing radiation, T^- and T^+,  we need to relate the supertranslation field C^- in advanced coordinates to the supertranslation field C^+ in retarded coordinates.  
Namely, we assume that we are given a classical spacetime whose asymptotic behavior is fully known to us. Then it is possible to describe this spacetime both in advanced and retarded BMS-gauge. Given an advanced coordinate system g_μν^v, we want to know if there is a unique retarded coordinate system g_μν^u we can associate to it. If we have such a mapping, it determines the relation of the advanced supertranslation field C^-, defined as the r^1 part of g^v_AB, and the retarded supertranslation field C^+, defined as the r^1 part of g^u_AB.

Given g_μν^v, we therefore have to find a diffeomorphism 𝒟 such that g_μν^u:=𝒟(g_μν^v) is in retarded BMS-gauge. Then we can read off from g_μν^u the C^+ associated to C^-. However, we could have instead considered the diffeomorphism 𝒟'= T^+∘𝒟, where T^+ is a supertranslation diffeomorphism in retarded coordinates. Also 𝒟' transforms the metric in advanced BMS-gauge to a metric in retarded BMS-coordinates. Clearly, if T^+ is a nontrivial supertranslation, the supertranslation field in the resulting metric differs from the one in g^u. From this consideration it is obvious that the matching between the advanced and the retarded supertranslation field is in general not unique. 

The only hope we could have is that there is a natural way to identify C^- and C^+.  In a static situation, a natural prescription is to require that the spatial part of the two metrics matches, 
    g^u_AB = g^v_AB .

As is shown explicitly in appendix <ref> for the example of the Schwarzschild metric, we can achieve this by identifying C^+(θ, φ) = - C^-(θ, φ), as also proposed in <cit.>. Up to a sign, we match the supertranslation field angle-wise.   Consequently, the same matching holds for the supertranslations:

    T^+(θ, φ) = - T^-(θ, φ)  .


There are several reason why the coordinate matching (<ref>) is natural. First of all, the prescription (<ref>) comes from a simple intuition. For an observer in a static spacetime who lives on a sphere of fixed radius, the description of the sphere should be the same independently of the choice of time coordinate. More generically, it is possible to require that the action of advanced and retarded supertranslations is the same in the bulk. This was done in <cit.> for the cases of Schwarzschild and Minkowski.


Moreover, we can consider a detector at big radius which is sensitive to gravitational memory. Then we investigate a process of back scattering, in which the angular distributions of incoming and outgoing energy are identical at each angle. This corresponds to a wall in the bulk which reflects the wave without further modifying it. In this case, the memory effect the ingoing wave causes, parameterized by T^-, is exactly canceled by the memory effect of the outgoing wave, parameterized by T^+, so that there is no overall memory effect after the process. In that case, if we match T^- and T^+ at each angle as in (<ref>), it is possible to simply describe the overall memory effect as T^-+T^+.

However, it is crucial to stress that the coordinate matching (<ref>) does not have any constraining power on the physical process. It does not predict outgoing from ingoing radiation, but only shows how one and the same setup can be described in different coordinates. This is also evident from figure <ref>. The matching condition at i^0 only relates the absolute values of the supertranslation fields. In contrast, processes of non-zero energy solely determine a change of the supertranslation field, as is clear from equations (<ref>) and (<ref>). Thus, radiation of non-zero energy is independent of the coordinate matching.




§ APPLICATION OF HAIR




 §.§ Planetary Hair


In order to make the ideas presented above concrete, we discuss an explicit example, namely the  application of a Goldstone supertranslation to a certain class of planets. We start from a spherically symmetric nongravitational source T_μν, which sources a spherically symmetric spacetime g_μν with ADM-mass M. In such a spacetime, we want to realize a Goldstone supertranslation, we send in a wave with total energy μ and angular distribution Δℱ̃_in in such a way that after some time, the planet emits a wave of the same energy μ but with a possibly different angular distribution Δℱ̃_out. Of course, only a special class of planets behaves in that way.

We explicitly construct such spacetimes in appendix <ref>, to which we refer the reader for details of the calculation. First, we consider the incoming wave. As discussed, the angular distribution Δℱ̃_in of injected energy determines an advanced supertranslation T^-. As derived in equation (<ref>), we can use it to describe the change of the metric due to the injected radiation:

    δ g^v_μν = τ_v_0, v_1(v) s^-(r) (ℒ_ξ_v(T^-) g^v_μν + 2 μ G/rδ_μ^0 δ_ν^0 )  ,

where ℒ_ξ_v(T^-) g^v_μν is an infinitesimal supertranslation which changes the supertranslation field by a small amount T^-. Whereas the asymptotic supertranslation T^- only depends on the leading part of the incoming energy, it is crucial to note that the transformation (<ref>) also depends on a careful choice of the subleading components of the incoming wave.[Subleading terms are the 1/r^3-term in T_00 and the whole T_0A in (<ref>). If one does not insist that the wave acts as a supertranslation also in the bulk, one is free to choose the coefficient of one of the two terms. The other one is determined by energy conservation: T_μν^μν;μ=0.]

Only with a particular choice, the wave acts as a diffeomorphism not only asymptotically but also in the bulk outside the planet.

We observe that the effect of the wave is twofold. First, it adds the total mass μ to the planet and secondly, it supertranslates the metric by T^-. However, these effects are localized both in space and time. The function τ_v_0, v_1(v) describes the smooth interpolation between g^v_μν and g^v_μν + δ g^v_μν, we have τ_v_0, v_1(v<v_0)=0 and τ_v_0, v_1(v>v_1)=1. The function s^-(r) describes the absorption of the wave, namely absorption takes place whenever s^- '(r)<0. There is no absorption outside the planet, s^-(r>R)=1, where R is the radius of the planet, and the wave is fully absorbed before it reaches the center, s^-(r=0)=0. It will be crucial to note that the transformation s^-(r)ℒ_ξ_v(T^-) g^v_μν only acts as a diffeomorphism when s^- '(r)=0. 

Moreover, the transformation (<ref>) shows that we focus on planets which have a second very special property aside from the fact that they emit as much energy as they receive: Namely there is no transport of energy between different angles. This means that the mass of the planet does not redistribute after absorption (the same will be true after emission).  The fact that this assumption is unnatural and not true for generic systems will contribute to our conclusions.

As a second step, we consider the emission of a wave by the planet.  Of course, the properties of the emitted wave depend on the internal dynamics of the source T_μν. It is crucial to note we cannot resolve them in our purely gravitational treatment, we cannot predict what wave will be emitted. From the point of view of gravity, any emission process is possible as long as it respects energy-momentum-conservation. However, we can study the effect of a given emitted wave. As derived in equation (<ref>), it can be described in terms of the supertranslation T^+ induced by the angular distribution Δℱ̃_out of outgoing energy:

    δ g^u_μν = τ_u_0, u_1(u)s^+(r) (ℒ_ξ_u(T^+) g^u_μν - 2μ G/rδ_μ^0 δ_ν^0 )  .

As for the case of absorption, the emission has two effects: It decreases the total mass by μ and it supertranslates the metric by T^+. Moreover, it is localized in space and time in an analogous manner. 

We want to compare the planet before and after the Goldstone supertranslation, we are interested in the combined effect of the transformations (<ref>) and (<ref>). To this end, we have to specify a mapping between the advanced and retarded supertranslations. As explained in section <ref>, we employ the angle-wise matching (<ref>). Thus, we obtain the static final state of the planet:

    δ g^tot_μν    =  θ(r-R) ℒ_ξ_u(T^+-T^-) g_μν
        + θ(R-r)(s^+(r)ℒ_ξ_u(T^+)g_μν - s^-(r) ℒ_ξ_u(T^-) g_μν) .

 We get a planet which has the same ADM-mass but a different angular distribution of mass. This is clear from the fact that the transformation (<ref>) acts as a diffeomorphism only outside the planet. 
 
 Since we used in our computation a planet with the special property that its angular distribution of energy is frozen, we can read off the distribution from difference of energy distributions of the injected and emitted wave. In this case, T^- - T^+ encodes all information about the angular energy distribution of the planet in the bulk.[For the planet with frozen energy distribution, there is also a very literal way in which one can interpret the quantity T^- - T^+: One can imagine a gedankenexperiment where a source of light is located in the interior of the planet after the Goldstone supertranslation and we collect the light rays on the sky. The light sent from this common center point determines in this way a section at infinity described by the supertranslation field T^- - T^+.  Thus, the different redshift effects due to the inhomogeneities of the planet matter distribution define a supertranslated section in the sky as the one for which light rays originate from a common spacetime point. This is reminiscent of Penrose's concept of "good sections" <cit.>.]

However, this is no longer true for generic systems which exhibit non-trivial dynamics after absorption and emission. In that case, T^- and T^+ merely encode the initial state. Only with full knowledge of the theory which governs the internal dynamics of the planet, we can infer the state of the planet at a later time from the asymptotic data T^- and T^+. 



  §.§.§ The Role of Supertranslations

 In summary, we obtain the following key properties of a Goldstone supertranslation in the case of a planet: Outside the planet, it acts as a diffeomorphism. In particular, it does not change its ADM-mass. In contrast, it does not act as a diffeomorphism inside the planet where absorption takes place. Therefore, it is not a trivial global diffeomorphism but changes the spacetime physically. Thus, the Goldstone supertranslation encodes differences in the angular distribution among matter configurations degenerate with respect to the ADM-conserved quantities. 

It is crucial to discuss the role of supertranslations in this process:

	
  * For an asymptotic observer, (T^-, T^+) can be used as label for the angular features of ingoing and outgoing radiation. 
	
  * An asymptotic observer, however, cannot infer T^+ from T^-. This is only possible with knowledge of internal dynamics of the planet.
	
  * Thus, (T^-, T^+) is a bookkeeping tool but without detailed information about the interior, it does not have predictive power.

As we shall discuss in a moment, the same conclusions hold in the black hole case. The only difference is that the internal dynamics leading to emission are fully quantum mechanical for a black hole. This will mean that in any classical description, supertranslation cannot constrain or even predict black hole evaporation.

 Using the example of the planet, it is easy to convince ourselves that antipodal matching cannot play a role in processes of non-zero energy. Namely if it did, this would mean that the only planets which could exist would have the extremely special property that they emit all energy they receive from one side exactly on the other side.




  §.§.§ Hidden Angular Features

Finally, we discuss the transformation (<ref>) when we do not have access to (T^-, T^+), when we do not record ingoing and outgoing radiation but only compare the initial and final state of the planet. In that case, the planet possesses an interesting property, namely a special kind of no-hair-theorem. Concretely, we take the perspective of an observer who has no access to the interior of the planet and discuss the difference between two planets which have the same mass but a different angular mass distribution. As we have observed, the transformation (<ref>) acts as a diffeomorphism outside the planet. Therefore, an outside observer cannot distinguish the two following cases when he is given a supertranslated outside metric. First, it could be the result of the transformation (<ref>), where the planet was physically changed due to a Goldstone supertranslation. Secondly, however, one can also obtain the supertranslated metric by acting on the initial planet with a global diffeomorphism. In this case, clearly, the planet does not change. Thus, also for a planet, an outside observer is not able to resolve angular features. In order to decide whether two asymptotic metrics differing by a supertranslation describe two different distributions of matter or the same distribution of matter in different coordinates, one needs access to the whole spacetime, the interior of the source.

We conclude that generic gravitational systems posses physical angular features which are inaccessible for an outside observer. This strengthens our believe that the microstates of a black hole have a non-trivial projection on angular features. The only difference is that while the restriction to outside measurements was artificial in the case of the planet, an outside observer has in principle no access to the interior of a black hole. As we will discuss in the next section, he can therefore never decide whether a supertranslated metric corresponds to a physical change of the matter inside the black hole or to a global and therefore meaningless diffeomorphism. This is the reason for the classical no-hair theorem of a black hole and why we assign an entropy to the black hole and not to the planet.





 §.§ Black Hole Quantum Hair




  §.§.§ Supertranslations as Bookkeeping Device

Now we are ready to discuss the system of our interest, namely black holes. Since absorption and emission are of different nature in that case, we will discuss them separately. For absorption, we can proceed in full analogy to the planet and inject a wave with total energy μ and arbitrary angular distribution Δℱ̃_in. By Birkhoff's theorem, the spacetime outside the black hole is the same as for the planet so that the wave behaves identically. As in the case of the planet, the wave cannot be absorbed outside the horizon  and acts as a diffeomorphism everywhere outside the black hole and also on the horizon.

 For the planet, we observed that the knowledge of injected energy alone does not suffice to predict what radiation the planet emits. Instead, this can only be done with knowledge of the interior dynamics of the planet. Those, however, can be described classically in the case of the planet. For the black hole, the situation is even worse. Not only do we not have access to any interior dynamics, but these dynamics are also fully quantum. It is impossible to describe them even with full classical knowledge of the interior of the black hole. 
	
 Before we elaborate on this point, we first show how it is possible to use supertranslations as bookkeeping device for black hole evaporation. Unlike for the case of the planet, this is a non-trivial question since the evaporation products are generic quantum states.  In order to define an associated supertranslation, we shall proceed as follows. We consider an ensemble of quantum-mechanically identical black holes of mass M.[Experimentally, we can realize this by preparing identical quantum states in such a way that they collapse and form black holes.]
 
For each black hole, we wait until it has emitted exactly one Hawking quantum.  We only record their angular features, the deviation from an isotropic emission. This means that we assume that the microstates of the black hole have a non-trivial projection on angular features of the evaporation products. As explained in the introduction and illustrated for the case of the planet, we believe this assumption to be natural. Thus, we record the Hawking quanta using a filter for angular features, where we use one for each spherical mode (l,m). This defines a probability distribution for the angular features of the ensemble:
 
    P(l,m)  .

 Obviously, the probability distribution (<ref>) only contains a part of the quantum-mechanically available information. However, we will only focus on it since it can be described in terms of classical supertranslations. At this point, it is crucial to point out that the probability distribution (<ref>) does not originate from a mixed state but as a result of an ordinary quantum measurement. Thus, unlike in a description in terms of a mixed state, it is not associated to any fundamental loss of information.
 
 Since we need to recover a featureless emission in the semi-classical limit, it follows that 
 
    P(0,0)=1-ϵ ,
 
 where ϵ→ 0 in the semi-classical limit. This means that only a fraction ϵ of the emitted quanta carries features. For l≥ 2, we consequently get
 
    P(l,m) = ϵ A_l,m ,

 where ∑_l=2^∞∑_m=-l^m=l A_l,m =1. The information contained in the P(l,m) is purely quantum mechanical. At the semi-classical level, we have that P(l,m)=δ_l0 and in the classical limit, we have no emission at all. 
 
 Using the quantum probability distribution (<ref>), we can associate to every Hawking quantum an average energy flux:
 
    ℱ_out = ħ r_g^-1∑_l=0^∞∑_m=-l^m=l  P(l,m) Y_l,m ,

 where Y_l,m are the standard spherical harmonics. Just like for the case of the planet, where we considered a classical process of emission, we can use the flux (<ref>) to define a classical supertranslation T^+. Of course, this is only possible as long as ħ≠ 0 since the energy flux is zero otherwise.
 When we record the quantum-mechanically emitted energy ℱ_out, we can proceed in analogy to the planet and use the supertranslation fields T^- and T^+ to track the evolution of the black hole. Concretely, in order to perform a Goldstone supertranslation, we first inject an energy μ and then we wait until n_H=μ/(ħ r_g^-1) quanta have evaporated, as is depicted in figure <ref>. Then we end up with a black hole of the same mass as before the process. Of course, the sensitivity of the final state on the initial state is suppressed by μ/M but unitarity dictates that the dependence is never trivial.




  §.§.§ Insufficiency of Supertranslation Hair

 However, it is impossible to predict T^+ solely from the knowledge of T^-. The reason is that the wave that we inject acts as a diffeomorphism outside the horizon and also on the horizon. Therefore, the geometry outside the black hole is unchanged after the wave has passed. Since the semi-classical Hawking calculation is only sensitive to the geometry on the horizon and outside the black hole, its result cannot change as a result of a supertranslation diffeomorphism. Therefore, additional knowledge about the interior is required to predict T^+.
 
 We can make this argument more concrete by taking the perspective of an observer who lives in a Schwarzschild metric supertranslated by T^-. The observer has no record of how the black hole was formed and is only allowed to make experiment outside the horizon. Her goal is to determine the microstate of the black hole. More specifically, she wants to know if the black hole is in the bald microstate, whose evaporation products are featureless and in particular perfectly isotropic, or in a non-trivial microstate, whose evaporation products carry some angular features. By our definition of microstate, one way to do so is to wait till the black hole has evaporated and to determine the properties of the evaporation products. 
 
 The question we are asking is if there is another way to determine the microstate of a black hole. The answer is negative, for the following reason: When an outside observer finds herself in a black hole metric with supertranslation field T^-, this can happen because of two very distinct reason. Firstly, it could be the result of injecting a wave with a non-trivial angular distribution of energy into a black hole. In that case, the black hole is in a non-trivial microstate and T^- indeed characterizes the microstate. 
 
 However, there is a second way in which we can obtain a supertranslated Schwarzschild metric. Namely, we can consider a featureless microstate, whose evaporation products are isotropic, and apply a supertranslation diffeomorphism to this setup. In this way, we do not change the physical state of the black hole but only describe it in a different metric. Thus, T^- can also correspond to a featureless microstate described in different coordinates.
 
 Without access to the evaporation products, the only way to distinguish those two cases – injection of wave with angular features versus global diffeomorphism – is to enter the black hole. There, the wave acts non-trivially, not as a diffeomorphism, whereas the global diffeomorphism still does.
 Since the same exterior metric can correspond to both a trivial and a non-trivial microstate, the metric alone cannot suffice to predict the evaporation products.  From the outside, it is therefore impossible to distinguish classical supertranslation hair and global diffeomorphisms.

 In summary, as in the case of a planet, we can use (T^-, T^+) as a natural bookkeeping device for the black hole to track the angular features of ingoing and outgoing radiation. However, knowing T^- does not suffice to predict T^+, an observer outside the black hole cannot infer T^+ from T^-. This is only possible with a microscopic model of the interior dynamics of the black hole.



  §.§.§ Generalization to Evaporation

Having discussed how we can implant hair on a black hole with a Goldstone supertranslation, it is trivially to consider the case of pure evaporation. We obtain it if we just leave out the first part of the Goldstone supertranslation, namely the injection of a wave. Therefore, it suffices to consider 𝒥^+ as screen, where the constraint (<ref>) determines the retarded supertranslation field T^+ in terms of the angular distribution Δℱ̃_out. In that case, the metric outside the black hole changes according to (<ref>):

    δ g^u_μν = τ_u_0, u_1(u) (ℒ_ξ_u(T^+) g^u_μν - 2μ G/rδ_μ^0 δ_ν^0 )  .

This equation shows that the back reaction splits in two parts. First, energy conservation dictates that the mass of the black hole is reduced by the total emitted energy μ = ∫u∫[2]Ωℱ_out. This part of the back reaction is undebatable but does not suffice to ensure unitarity of the process. Fortunately, ℱ_out contains more information than just the emitted energy, namely the supertranslation T^+. Consequently, we obtain the back reacted black hole not only by reducing its mass, but by supertranslating it by T^+.  This approach is only valid if the supertranslation acts non-trivially in the interior of the black hole, because it is induced by a physical wave. But in that case, the ability to associate hair to a black hole is equivalent to the ability to purify its evaporation.






 §.§ A Comment on Page's Time


So far, we have not specified the magnitude of deviations from a thermal evaporation. We can estimate them by requiring that we reproduce Page's time in our approach. In its most basic formulation, Page's time is a direct consequence of describing the black hole evaporation in a Hilbert space of fixed dimension. In brief, if we keep the dimension of the full Hilbert space, which describes at any time both the black hole and the emitted radiation, fixed and equal to 2^N, then at t=t_P, which corresponds to the half life-time, the evaporation of ∼ N/2 quanta, there is no place to continue increasing the entanglement between the radiation and the black hole internal degrees of freedom. At this time, entanglement starts to decrease and information starts to be delivered. This makes clear why purification of black hole evaporation relies on fixing N and keeping it finite.

 Page's time can be defined as the time-scale for the emission of the order of N quanta. Therefore, we first consider an ensemble of N identical quantum mechanical black holes and for each of them, we record the first emitted quantum.
For a measurement on a single black hole, the standard deviation is

    σ_1 ∼ O(1)

since the quanta are distributed isotropically to leading order. However, when we average over N measurements, the standard deviation decreases as

    σ_N∼1/√(N) .

Features become visible as soon as their strength becomes bigger or equal than the uncertainty of the measurement. After Page's time we can therefore resolve features with the relative amplitude
	
    ϵ∼1/√(N) .

	In the formulation of the probability distribution (<ref>), this means that after O(N) measurement, those features becomes visible which are only carried by a fraction 1/√(N) of the quanta. 

 So far, we have only considered one emission for N identical black holes. If we consider instead O(N) emissions of a single black hole, the difference is that the probability distribution for each emission step is generically different. This is true because of the back reaction of the previously emitted quanta. However, the argument in terms of the resolution stays the same, after Page's time, we can still resolve those features which are only carried by a fraction ϵ∼ 1/√(N) of quanta. This argument provides evidence for the black hole N-portrait <cit.> where features are 1/N-effects with a resolution scale O(1/√(N)).[An interesting question that we shall not discuss in this note but that can be worth to mention is the possibility that a quantum computer designed using a Grover like algorithm <cit.> can reduce t_P from O(N) to O(√(N)).]

 In particular, in an S-matrix analysis along the lines of <cit.>, where black hole formation was studied as 2→ N-scattering process, angular features should appear as 1/N-correction to the leading amplitude.



§ CONCLUSION




The main message of this note is easy to summarize. Any form of black hole hair should imply the existence of features in the black hole evaporation products, in the emitted radiation. This obvious requirement immediately entails, given the intrinsically quantum nature of black hole radiation, that black hole hair should be defined quantum-mechanically and that such a definition is inseparable from the mechanism through which the black hole delivers, in the radiation, information about its internal structure.

  In this note, we have suggested to define hair on the basis of elementary processes of classical absorption followed by quantum emission. Moreover, we specialized to angular features in the radiation. This simplification has been done in order to use the asymptotic symmetry group and the corresponding supertranslations to parametrize both the incoming and the emitted radiation. Since we have complete control over the angular features of the injected radiation, we can define hair on the basis of the angular features of the quantum-mechanically emitted radiation. These features encode information about the internal structure of the black hole which can be measured by an external observer. In this sense, they provide an operative and intrinsically quantum definition of hair. 
  
  In principle, we can imagine two different sources of those features of the emitted radiation. The first one is a classical modification of the near horizon geometry that will modify the corresponding semi-classical Bogoliubov transformations. The second one is a real quantum interaction of the injected radiation with the quantum constituents of the black hole. The first possibility requires to define local changes of the horizon geometry that preserve all the ADM-charges. Thus, locally, they can be always tuned to be equivalent to a diffeomorphism. Therefore, they  cannot have observable consequences, classical supertranslations do not suffice to define observable black hole hair. So the only real possibility of quantum emitted radiation with features is having a non-trivial scattering between the injected radiation and the microscopic constituents of the black hole. This means that the features that define hair in the way we are suggesting depend on the microscopic quantum structure, which we can parametrize as a dependence on the black hole entropy N. Thus, the hair that we are defining vanishes in the limit N=∞. 
  
  As it is clear from the discussion, this way of addressing the definition of hair is what we can call an S-matrix approach, where by that S-matrix we simply mean the dynamics involved in the complex process of actual absorption and quantum emission. If we focus on angular features, we can encode the properties of the hair in terms of the commutators, as operators, of this S-matrix and the generators of the asymptotic symmetry group. Associating with the injected energy a supertranslation T^- in BMS_-, a way to approach the existence of hair is by considering the commutator [S, T^-]. Generically, the non-trivial hair will be associated with the symmetry generators that are broken since those are the ones that will create net differences between the angular features of the injected and emitted radiation. Although the infrared dynamics of gravity selects the zero-energy modes as natural symmetries of S, they are not able to tell us anything about the internal structure of the black hole since they are decoupled. Zero-modes are unable to encode observable features.
  
  What we have presented in this note is just the general framework to address the problem of quantum hair. In order to go further, it is necessary to use a concrete model of the black hole interior. The model in <cit.> provides, in principle, the tools to  address this questions in a quantitative way, something to which we hope to come back in the future. 










§ RECAP OF IR-EFFECTS



In this appendix, we shall collect some well-known facts about infrared physics which could be useful to clarify some controversial aspects on the meaning of soft modes. Some of these issues have been revisited recently in a series of papers <cit.>. 
	
	
  * In QED, asymptotic physical states associated with freely moving charged particles should be dressed in order to satisfy the Gauss law constraints. This dressing simply adds to the freely moving charge its companion electrostatic field, the non radiative part of the retarded Lienard-Wiechert-field behaving at large distances as 1/r^2. In quantum field theory, this dressing can be defined using a coherent state of off-shell photons <cit.> with dispersion relation ω(k) = k⃗v⃗ for v⃗ the velocity of the charged particle. This coherent state dressing contains an infinite number of k=0 modes and it is identical to the dressing operator defined in <cit.>. In scattering theory, one can define physical asymptotic states and an IR-safe S matrix using this dressing operator.
	
  * Alternatively, one can use no dressing. Then, in perturbative QED as well as in perturbative gravity, we find IR-divergences due to virtual photon/graviton loops. These, after a careful analysis of overlapping divergences, can be resummed and exponentiated <cit.>.  When we consider the transition from an initial state |α⟩ to a final state |β⟩, we obtain

    S_α, β^loop = ^B_α, βlnλ/Λ/2 S_α, β^0  ,

where S_α, β^0 is the amplitude without taking into account soft loops whereas S_α, β^loop contains them. Here Λ is a UV-cutoff that defines what is soft, λ is a IR-cutoff and B_α, β is a non-negative number, which only depends on the initial state |α⟩ and the final state |β⟩. It is zero if and only if the ingoing current in |α⟩ matches the outgoing current in |β⟩ antipodally at each angle. In the case of gravity, it scales as B_α, β∼ G s, where s is the energy of the process.[That this scaling also holds for graviton scattering at an ultra-Planckian center of mass energy was shown in <cit.>.]

  For B_α,  β≠ 0, soft loops clearly lead to a vanishing amplitude in the limit λ→ 0.

  * In order to cancel the IR-divergences due to virtual photons/gravitons, the Bloch-Nordsieck-recipe <cit.> requires to add a certain class of soft emission processes. Again the effects of emitting theses soft IR-modes of energies below ϵ can be resummed and exponentiated, yielding the rate <cit.>
 
    |S_α, β^full|^2 := ∑_γ  |S_α, βγ^soft|^2 =  ^B_α, βlnϵ/λ f(B_α, β)  |S_α, β^soft|^2  ,

 where f(B_α, β) is due to energy conservation and f(B_α, β)≈ 1 for small B_α, β. Combing the contribution from (<ref>) and (<ref>), one obtains a rate which is independent of the IR-cutoff λ and in particular finite for λ→ 0. 
  This cancellation leads to the connection, highlighted in <cit.>, between the soft photon theorem and the electrostatic coherent state dressing. In QED, we do not have new symmetries besides the decoupling of zero-energy photons. The same is true in perturbative gravity.	

  * In the correction factor ^B_α, βlnϵ/λ in (<ref>), the n^th summand of the exponential series comes from the emission of n IR-modes. Therefore, we can estimate the number of soft modes from the term which gives the biggest contribution in the series. This gives

    n_soft^unres∼ B_α, βlnϵ/λ .

We conclude that the number of unresolved soft modes only scales logarithmically with the
infrared resolution scale ϵ.





§ MATCHING IN SCHWARZSCHILD COORDINATES


In this section, we demonstrate explicitly how we can transform a Schwarzschild metric with non-trivial supertranslation field from advanced to retarded coordinates. In this way, we show how we can naturally identify the advanced supertranslation field C^- with the retarded one C^+. We start from the Schwarzschild metric g_μν^v ,0 in advanced coordinates without supertranslation field:

    d s^2 = -(1-2GM/r) v^2 + 2  v r +r^2 Ω^2  .
 
The corresponding generators of supertranslations are
 
	
    ξ_v^v =    f^-  ,  
    ξ_v^r =    -1/2 D^2 f^-  , 
    ξ_v^A =   f^-, A/r ,


which are characterized by an arbitrary function f^- on the sphere. Thus, the supertranslated metric is

    g_μν^v(f^-) = g_μν^v ,0 + ℒ_ξ_v(f^-) g_μν^v ,0 .


In retarded coordinates, the Schwarzschild metric g_μν^u ,0 without supertranslation field is:

    d s^2 = -(1-2GM/r) u^2 - 2  u r +r^2 Ω^2  .
 
The corresponding generators of supertranslations are

	
    ξ_u^v =    f^+  ,  
    ξ_u^r =   1/2 D^2 f^+  , 
    ξ_u^A =    -f^+, A/r ,


where it is important to note that the signs of ξ_u^r and ξ_u^A have changed with respect to (<ref>). The supertranslated metric is:

    g_μν^u(f^+) = g_μν^u ,0 + ℒ_ξ_u(f^+) g_μν^u ,0 .



The task now is to transform g_μν^v to retarded coordinates. As explained in section <ref>, there can in general not be a unique way to match the advanced and retarded supertranslation fields. However, a natural choice in a static metric is to require that the spherical metrics match: g_AB^v(f^-) = g_AB^u(f^+). Therefore, we use the diffeomorphism 𝒟_m defined by

    v = u + 2 _r_0^r 1/1-2GM/rr' - D^2 f^-/1-2GM/r - 2 f^-  .

Then it turns out that

    𝒟_m(g_μν^v(f^-)) = g_μν^u ,0 - ℒ_ξ_u(f^-) g_μν^u ,0 = g_μν^u(-f^-)   .

Thus, we identify

    f^+=-f^-  .

Up to a sign, the supertranslation field in advanced coordinates matches the retarded one angle-wise. With this choice, not only the spherical metrics match, but also the g_00-components, the Newtonian potentials.



§ EXPLICIT SOLUTION FOR GOLDSTONE SUPERTRANSLATION OF A PLANET





  §.§.§ Step 1: Absorption

The Goldstone supertranslation consists of two steps: First, an initially spherically symmetric planet absorbs as wave. As is well-known (see (9.3) in <cit.>), the metric of a static spherically symmetric spacetime can be cast in the general form

    d s^2 = - A(r) t^2 + B(r)  r^2 +r^2 Ω^2  ,
 
where all physical information is contained in the tt- and rr-components. Since we want to describe a planet, there should neither be a surface of infinite redshift, A(r) > 0  ∀ r, nor an event horizon, B(r) < ∞ ∀ r. Furthermore, asymptotic flatness implies that A(r) r→∞⟶ 1 and B(r) r→∞⟶ 1 sufficiently fast. Using the transformation
 
    v = t + _r_0^r  r' √(B/A) ,

we obtain the metric g_μν^v in advanced BMS-gauge:

    s^2 = - A  v^2 + 2 √(AB) v  r + r^2 Ω^2  ,

which is suited to describe incoming radiation. Note that this metric describes the whole spacetime and not only its asymptotic region, r→∞.

We will restrict ourselves to infinitesimal supertranslations. In advanced time, these are generated by
 
	
    ξ_v^v =    f^-  ,  
    ξ_v^r =    -1/2 r D_B ξ_v^B  , 
    ξ_v^A =    f^-, A_r^∞ r' (√(AB) r'^-2)  ,


where an arbitrary function f^- on the sphere determines the change of the supertranslation field.  We denote it by f^- instead of T^- in this appendix to avoid confusion with the energy-momentum-tensor of the wave. The minus-superscript indicates that we deal with a supertranslation in advanced coordinates. Our goal is to realize the infinitesimal diffeomorphism defined by (<ref>) in a physical process, outside the planet, we want to have the stationary metric g^v_μν before some time v_0 and after some point of time v_1, we want to end up in the stationary metric g^v_μν + ℒ_ξ_v(f^-) g_μν^v. For v_0<v<v_1, physical radiation interpolates between the two metrics. Inside the planet, the wave should be absorbed so that the transformation fades out and the metric around the origin remains unchanged. Adding as final ingredient a change of the Bondi mass μ, which is necessary to ensure the positive energy condition, we obtain

    δ g^v_μν = τ_v_0, v_1(v) s^-(r) (ℒ_ξ_v(f^-) g^v_μν + 2μ G/rδ_μ^0 δ_ν^0 )  ,

where 0≤τ_v_0, v_1(v) ≤ 1 parameterizes the interpolation, τ_v_0, v_1(v<v_0)=0 and τ_v_0, v_1(v>v_1)=1. The function s^-(r) describes the absorption of the wave by the planet. It has the property that it is monotonically increasing with s^-(0)=0 and s^-(∞)=1, where s^-(0)=0 ensures that the wave is fully absorbed before the origin and no black hole forms. Moreover, s^- '(r)≠ 0 is only permissible whenever the local energy density of the planet is non-zero. The magnitude of s^- '(r) determines how much absorption happens at r. It is crucial to note that the transformation s^-(r) ℒ_ξ_v(f^-) g_μν is a diffeomorphism only where s^-(r) is constant. Thus, the transformation (<ref>) acts as a diffeomorphism only outside the planet, but not inside. This reflects the fact that we want to obtain a physically different planet. A transformation which acts as a diffeomorphism everywhere could not achieve this.

Since we work with infinitesimal supertranslations, it is important that we stay within the regime of validity of this first-order approximation, that terms linear in f^- dominate. As it will turn out in the calculation, this is the case if max_(θ, φ) |f^-| ≪ v_1-v_0. This means that the time-shift induced by the supertranslation must be much smaller than the time-scale of the process, the supertranslation must be performed slowly. We will choose v_1-v_0 such that this is the case and so that we can neglect all higher orders in f^- when we calculate the Einstein equations.

We have to show that the transformation (<ref>) leads to a valid solution of the Einstein equations. Thus, if we calculate the Einstein G_μν and consequently the new energy-momentum-tensor T_μν, we have to demonstrate that this is a valid source. To this end, we have to perform two checks. First of all, it must be conserved, T_μν^μν;μ = 0. This is trivially true in our construction because of the Bianchi identity, G_μν^μν;μ = 0. Secondly, we have to show that T_μν fulfills an appropriate energy condition. For that purpose, we first note that this perturbation only depends on the local geometry, except for ξ^μ_v and s^-(r), which also depend on spacetime points at bigger radii. Thus, outside the planet, we have the same solution as in <cit.>, except for the fact the we perform our supertranslation slowly:
 
	
    T_00 =   1/4π r^2[μ - 1/4 D^2(D^2 + 2) f^- + 3M/2r D^2 f^-] τ'_v_0, v_1(v)  ,
    
    	T_0A =    3M/8π r^2 D_A f^- τ'_v_0, v_1(v)  ,


where we used that there is no absorption outside the planet: s^-(r>R)=1.
Obviously, the energy condition is fulfilled. At this point, we remark that leaving out all subleading parts, which are proportional to M, would also lead to a valid wave in the metric (<ref>), i.e., T_μν^μν;μ would also be true to all orders if one only considered the leading order of (<ref>). This means that we add the subleading parts to (<ref>) not because of energy conservation but since we want to realize the transition (<ref>) not only to leading order in 1/r, but to all orders.[Of course, energy conservation relates the two subleading parts of T_00 and T_0A. When we choose one, it determines the other.]

Fortunately, we do not either have to worry about the energy condition inside the planet. For a small enough perturbation, this is true since the energy condition inside a planet is not only marginally fulfilled. This means that s^-'(r) can be non-zero inside the planet: This corresponds to absorption of the wave by the planet.

Lastly, we have to show that the wave is still a valid solution after it has been partly absorbed. For the purpose of illustration, we model the planet as a sequence of massive shell with vacuum in between, T_μν=0. In that case, the only non-trivial question is whether (<ref>) fulfills the energy condition after it passed some or all of the shells. Therefore, we calculate the energy-momentum-tensor in this region. By Birkhoff's theorem, the local geometry corresponds to a Schwarzschild solution with diminished mass M̃ (where M̃ can be zero). It only depends on the matter which it has passed via ξ_v^μ. We parameterize the difference of ξ_v^μ and the vector field one would get in a pure Schwarzschild geometry of mass M̃ by 

    σ := _R_min^R_max r' ((√(AB)-1) r'^-2) ,

where  we have no matter for r>R_max and between r and R_min.[We use that AB=1 in a Schwarzschild geometry of arbitrary mass.]

Explicitly, this means that we can write

    ξ_v^A =    f^-, A(1/r + σ)  ,

where it is important that σ does not depend on r in our region of interest. Of course, σ=0 corresponds to the case when there is no matter outside.

With the help of Mathematica <cit.>, we compute:
 
	
    T_00 =   1/4π r^2[μ̃ - (1+σ r)(1/4 D^2(D^2 + 2) f̃^- - 3M̃/2r D^2 f̃^-)] τ'_v_0, v_1(v)  ,
    
    	T_0A =    [3M̃/8π r^2 D_A f̃^- +σ/16 π D_A (D^2 + 2)f̃^-] τ'_v_0, v_1(v)  ,
    
    	T_AB =    -σ r/8 π[(2D_A D_B - γ_ABD^2 )f̃^- ]τ'_v_0, v_1(v)  ,


where f̃^-=s^-(r)f^- is the supertranslation which is attenuated because of absorption in the outer shells. It is crucial to note that s^- '(r)=0 in this calculation since we are not inside one of the shells of the planet and likewise μ̃ =s^-(r)μ.
As we can estimate σ very crudely as σ < 1/R, we see that for sufficiently large μ, the energy condition is fulfilled. With a more accurate estimate, we expect that the freedom of choosing μ is not restricted when the wave passes a massive shell. In summary, we have shown that the metric (<ref>), which describes the dynamical transition from a spherically symmetric planet to a counterpart with nontrivial angular distribution of mass, is a valid solution.



  §.§.§ Step 2: Emission

The second step is to describe the emission of the wave by the planet. Thus, our initial metric is the one after absorption, as determined by equation (<ref>):

    δ g^v_μν = s^-(r) (ℒ_ξ_v(f^-) g^v_μν + 2μ G/rδ_μ^0 δ_ν^0 )  .

As we want to consider emission, our first step is to transform it to retarded coordinates. Intuitively, it is clear that it should be possible to describe a slightly asymmetrical planet also in retarded coordinates. While it is generically hard to write down the corresponding diffeomorphism which connects the two metrics, we can use that the metric of a planet does not differ from Schwarzschild in the exterior region. Therefore, we can use the diffeomorphism (<ref>) to obtain

    g_μν^u = 	g_μν^u, 0 + s^-(r)(ℒ_ξ_u(-f^-) g_μν^u, 0 + θ(R-r)dev)  ,

where g_μν^u, 0 is the metric of the initial, spherically symmetric planet in retarded coordinates. This means that we apply a supertranslation in retarded coordinates which is defined by the function f^- used to defined the advanced supertranslation. The function dev accounts for the fact that we do not know the continuation of the matching diffeomorphism (<ref>) to the interior of the planet. Therefore, g_μν^u might deviate slightly from BMS-gauge but only in the interior. We expect, however, that the matching diffeomorphism can be continued such that dev=0. Finally, we want to point out that g_μν^u(r=0) = g_μν^u, 0(r=0) since s^-(r=0)=0, the wave does not reach the center and the mass distribution of the planet is still spherically symmetric around r=0.

The case of the planet provides us with another justification why the matching (<ref>) is natural. With this identification, both the metric (<ref>) in advanced coordinates and the metric (<ref>) in retarded coordinates cover the whole manifold. Extrapolating the results of <cit.>, where finite supertranslations of Schwarzschild and Minkowski are discussed, we expect that for any other matching, for any other value of the supertranslation field, this is no longer the case. If this is true, the requirement that the BMS-coordinate system covers the whole manifold singles out a unique value of the advanced supertranslation field as well as a unique value of the retarded supertranslation field, and therefore a coordinate matching.

Next, we want to describe how the metric (<ref>) emits a wave. This wave should realize a supertranslation described by f^+, which is generically different from f^-:

    δ g^u_μν = τ_u_0, u_1(u)s^+(r) (ℒ_ξ_u(f^+) g^u, 0_μν - 2μ G/rδ_μ^0 δ_ν^0 )  ,

where we used that ℒ_ξ_u(f^+) g^u_μν = ℒ_ξ_u(f^+) g^u, 0_μν to first order in f^+ and f^-. Thus, working only to first order simplifies our calculations significantly since we can simply use the calculations for the absorption. The wave (<ref>) becomes:
 
	
    T_00 =   1/4π r^2[μ̃ - (1+σ r)(1/4 D^2(D^2 + 2) f̃^+ - 3M̃/2r D^2 f̃^+)] τ'_u_0, u_1(u)  ,
    
    	T_0A =     -[3M̃/8π r^2 D_A f̃^+ +σ/16 π D_A (D^2 + 2)f̃^+] τ'_u_0, u_1(u)  ,
    
    	T_AB =    -σ r/8 π[(2D_A D_B - γ_ABD^2 )f̃^+ ]τ'_u_0, u_1(u)  .


As for the absorption, we have shown that we can realize the transformation (<ref>) with a physical wave.

Finally, we analyze the joint effect of absorption and emission. Combining the transformations (<ref>) and (<ref>), we get total total change of the metric:

    δ g_μν^tot    = θ(r-R) ℒ_ξ_u(f^+-f^-) g_μν^u, 0
        + θ(R-r)(s^+(r)ℒ_ξ_u(f^+) g_μν^u, 0 - s^-(r) ℒ_ξ_u(f^-) g_μν^u, 0 + dev)  ,

	where we used retarded coordinates. As desired, the mass of the planet stays invariant. Moreover, δ g_μν^tot acts as a diffeomorphism outside the planet, namely it is the difference of the advanced supertranslation, described by f^-, and the retarded supertranslation, described by f^+. If we furthermore assume that the term dev, which reflects our incomplete knowledge of the matching between advanced and retarded coordinates in the planet, is zero, we see that the metric does not change for f^-=f^+. We obtain a trivial transformation if the angular energy distribution of ingoing and outgoing radiation is angle-wise the same.





 

§ ACKNOWLEDGEMENTS


We are happy to thank Gia Dvali for many stimulating discussions and comments. We also thank Artem Averin, Henk Bart and Raoul Letschka for discussions. The work of C.G. was supported in part by Humboldt Foundation and by Grants: FPA 2009-07908 and ERC Advanced Grant 339169 "Selfcompletion". 


tocsectionReferences




1
	
	noHair
	P. T. Chrusciel, J. Lopes Costa, and M. Heusler, Stationary Black Holes: Uniqueness and Beyond,
	http://dx.doi.org/10.12942/lrr-2012-7Living Rev. Rel. 15 (2012) 7,
	http://arxiv.org/abs/1205.6112arXiv:1205.6112 [gr-qc].
	
	
	bekenstein
	J. D. Bekenstein, Black holes and entropy,
	http://dx.doi.org/10.1103/PhysRevD.7.2333Phys. Rev. D7 (1973)
		2333–2346.
	
	hawkingRadiation
	S. W. Hawking, Particle creation by black holes,
	http://dx.doi.org/10.1007/BF02345020Commun. Math. Phys. 43
		(1975) no. 3, 199–220.
	
	page
	D. N. Page, Information in black hole radiation,
	http://dx.doi.org/10.1103/PhysRevLett.71.3743Phys. Rev. Lett. 71 (1993)  3743–3746,
	http://arxiv.org/abs/hep-th/9306083arXiv:hep-th/9306083
		.
	
	

NPortrait
G. Dvali and C. Gomez, Black Hole's Quantum N-Portrait,
http://dx.doi.org/10.1002/prop.201300001Fortschr. Phys. 61
	(2013)  742–767,
http://arxiv.org/abs/1112.3359arXiv:1112.3359 [hep-th].


1/NHair
G. Dvali and C. Gomez, Black Hole's 1/N Hair,
http://dx.doi.org/10.1016/j.physletb.2013.01.020Phys. Lett. B 719 (2013)  419–423,
http://arxiv.org/abs/1203.6575arXiv:1203.6575 [hep-th].


giaThermal
G. Dvali, Non-Thermal Corrections to Hawking Radiation Versus the
		Information Paradox,
http://dx.doi.org/10.1002/prop.201500096Fortsch. Phys. 64
	(2016)  106–108,
http://arxiv.org/abs/1509.04645arXiv:1509.04645 [hep-th].


StromingerMemory
A. Strominger and A. Zhiboedov, Gravitational Memory, BMS
		Supertranslations and Soft Theorems,
http://dx.doi.org/10.1007/JHEP01(2016)086J. High Energy Phys. 2016 no. 1, 86,
http://arxiv.org/abs/1411.5745arXiv:1411.5745 [hep-th].


HawkingFirst
S. W. Hawking, The Information Paradox for Black Holes,
http://arxiv.org/abs/1509.01147arXiv:1509.01147 [hep-th].



BMS
H. Bondi, M. Van der Burg, and A. Metzner, Gravitational waves in general
	relativity. VII. Waves from axi-symmetric isolated systems,
http://dx.doi.org/10.1098/rspa.1962.0161Proc. R. Soc. London, Ser. A
	269 (1962) no. 1336, 21–52.

R. K. Sachs, Gravitational waves in general relativity. VIII. Waves in
	asymptotically flat space-time,
http://dx.doi.org/10.1098/rspa.1962.0206Proc. R. Soc. London, Ser. A
	270 (1962) no. 1340, 103–126.

R. Sachs, Asymptotic symmetries in gravitational theory,
http://dx.doi.org/10.1103/PhysRev.128.2851Phys. Rev. 128 (1962)
	no. 6, 2851–2864.



BMSotherPenna

R. F. Penna, BMS invariance and the membrane paradigm,
http://dx.doi.org/10.1007/JHEP03(2016)023J. High Energy Phys. 2016 no. 3, 023,
http://arxiv.org/abs/1508.06577arXiv:1508.06577 [hep-th].
BMSotherDieter
G. Dvali, C. Gomez, and D. Lüst, Classical Limit of Black Hole Quantum
		N-Portrait and BMS Symmetry,
http://arxiv.org/abs/1509.02114arXiv:1509.02114 [hep-th].
BMSotherFlanagan
E. E. Flanagan and D. A. Nichols, Conserved charges of the extended
		Bondi-Metzner-Sachs algebra,
http://dx.doi.org/10.1103/PhysRevD.95.044002Phys. Rev. D 95
	(2017) no. 4, 044002,
http://arxiv.org/abs/1510.03386arXiv:1510.03386 [hep-th].
BMSotherDonnay
L. Donnay, G. Giribet, H. A. Gonzalez, and M. Pino, Supertranslations and
		Superrotations at the Black Hole Horizon,
http://dx.doi.org/10.1103/PhysRevLett.116.091101Phys. Rev. Lett. 116 (2016) no. 9, 091101,
http://arxiv.org/abs/1511.08687arXiv:1511.08687 [hep-th].
BMSotherBlau
M. Blau and M. O'Loughlin, Horizon Shells and BMS-like Soldering
		Transformations,  http://dx.doi.org/10.1007/JHEP03(2016)029J.
	High Energy Phys. 2016 no. 3, 29,
http://arxiv.org/abs/1512.02858arXiv:1512.02858 [hep-th].
BMSotherBarnich
G. Barnich, C. Troessaert, D. Tempo, and R. Troncoso, Asymptotically
		locally flat spacetimes and dynamical nonspherically-symmetric black holes in
		three dimensions,
http://dx.doi.org/10.1103/PhysRevD.93.084001Phys. Rev. D 93
	(2016) no. 8, 084001,
http://arxiv.org/abs/1512.05410arXiv:1512.05410 [hep-th].
BMSotherArtem
A. Averin, G. Dvali, C. Gomez, and D. Lust, Gravitational Black Hole Hair
		from Event Horizon Supertranslations,
http://dx.doi.org/10.1007/JHEP06(2016)088J. High Energy Phys. 2016 no. 6, 88,
http://arxiv.org/abs/1601.03725arXiv:1601.03725 [hep-th].
compereMatching
G. Compère and J. Long, Vacua of the gravitational field,
http://dx.doi.org/10.1007/JHEP07(2016)137J. High Energy Phys. 2016 no. 7, 137,
http://arxiv.org/abs/1601.04958arXiv:1601.04958 [hep-th].

compereFinite
G. Compère and J. Long, Classical static final state of collapse with
		supertranslation memory,
http://dx.doi.org/10.1088/0264-9381/33/19/195001Class. Quant. Grav.
	33 (2016) no. 19, 195001,
http://arxiv.org/abs/1602.05197arXiv:1602.05197 [gr-qc].
BMSotherEling
C. Eling and Y. Oz, On the Membrane Paradigm and Spontaneous Breaking of
		Horizon BMS Symmetries,
http://dx.doi.org/10.1007/JHEP07(2016)065J. High Energy Phys. 2016 no. 7, 65,
http://arxiv.org/abs/1605.00183arXiv:1605.00183 [hep-th].
BMSotherDonnay2
L. Donnay, G. Giribet, H. A. Gonzalez, and M. Pino, Extended Symmetries
		at the Black Hole Horizon,
http://dx.doi.org/10.1007/JHEP09(2016)100J. High Energy Phys. 2016 no. 9, 100,
http://arxiv.org/abs/1607.05703arXiv:1607.05703 [hep-th].
BMSotherCai
R.-G. Cai, S.-M. Ruan, and Y.-L. Zhang, Horizon supertranslation and
		degenerate black hole solutions,
http://dx.doi.org/10.1007/JHEP09(2016)163J. High Energy
	Phys. 2016 no. 9, 163,
http://arxiv.org/abs/1609.01056arXiv:1609.01056 [gr-qc].








HPS1
S. W. Hawking, M. J. Perry, and A. Strominger, Soft Hair on Black
		Holes,  http://dx.doi.org/10.1103/PhysRevLett.116.231301Phys.
	Rev. Lett. 116 (2016) no. 23, 231301,
http://arxiv.org/abs/1601.00921arXiv:1601.00921 [hep-th].


HPS2
S. W. Hawking, M. J. Perry, and A. Strominger, Superrotation Charge and
		Supertranslation Hair on Black Holes,
http://dx.doi.org/10.1007/JHEP05(2017)161J. High Energy Phys. 2017 no. 5, 1–33,
http://arxiv.org/abs/1611.09175arXiv:1611.09175 [hep-th].




criticalGas
G. Dvali and C. Gomez, Black Holes as Critical Point of Quantum Phase
		Transition,  http://dx.doi.org/10.1140/epjc/s10052-014-2752-3Eur.
	Phys. J. C 74 (2014)  2752,
http://arxiv.org/abs/1207.4059arXiv:1207.4059 [hep-th].





stromingerNew
A. Strominger, Black Hole Information Revisited,
http://arxiv.org/abs/1706.07143arXiv:1706.07143 [hep-th].





porrati1
M. Mirbabayi and M. Porrati, Shaving off Black Hole Soft Hair,
http://dx.doi.org/10.1103/PhysRevLett.117.211301Phys. Rev. Lett. 117 (2016) no. 21, , http://arxiv.org/abs/1607.03120arXiv:1607.03120 [hep-th].

sever
B. Gabai and A. Sever, Large Gauge Symmetries and Asymptotic States in
	QED,  http://dx.doi.org/10.1007/JHEP12(2016)095J. High Energy
	Phys. 2016 no. 12, 95, http://arxiv.org/abs/1607.08599arXiv:1607.08599 [hep-th].

raoul
C. Gomez and R. Letschka, Memory and the Infrared,
http://arxiv.org/abs/1704.03395arXiv:1704.03395 [hep-th].

mischa
M. Panchenko, The infrared triangle in the context of IR safe S
	matrices,  http://arxiv.org/abs/1704.03739arXiv:1704.03739
		[hep-th].

porrati2
R. Bousso and M. Porrati, Soft Hair as a Soft Wig,
http://arxiv.org/abs/1706.00436arXiv:1706.00436 [hep-th].

R. Bousso and M. Porrati, Observable Supertranslations,
http://arxiv.org/abs/1706.092801706.09280 [hep-th].




ashtekar
A. Ashtekar, Geometry and Physics of Null Infinity,
http://arxiv.org/abs/1409.1800arXiv:1409.1800 [gr-qc].






memoryFirst
Y. B. Zel'dovich and A. G. Polnarev, Radiation of gravitational waves by a
	cluster of superdense stars,  Sov. Astron. 18 (1974)  17.

wald
S. Hollands, A. Ishibashi, and R. M. Wald, BMS Supertranslations and
	Memory in Four and Higher Dimensions,
http://arxiv.org/abs/1612.03290arXiv:1612.03290 [gr-qc].






StromingerBMSInvariance
A. Strominger, On BMS Invariance of Gravitational Scattering,
http://dx.doi.org/10.1007/JHEP07(2014)152J. High Energy Phys. 2014 no. 7, 152,
http://arxiv.org/abs/1312.2229arXiv:1312.2229 [hep-th].


stromingerLecture
A. Strominger, Lectures on the Infrared Structure of Gravity and Gauge
	Theory,  http://arxiv.org/abs/1703.05448arXiv:1703.05448
		[hep-th].



YFS
D. R. Yennie, S. C. Frautschi, and H. Suura, The infrared divergence
		phenomena and high-energy processes,
http://dx.doi.org/10.1016/0003-4916(61)90151-8Annals Phys. 13
	(1961)  379–452.


weinberg
S. Weinberg, Infrared photons and gravitons,
http://dx.doi.org/10.1103/PhysRev.140.B516Phys. Rev. 140 (1965)
	B516–B524.





semenoff
D. Carney, L. Chaurette, D. Neuenfeld, and G. W. Semenoff, Infrared
	quantum information,  http://arxiv.org/abs/1706.03782arXiv:1706.03782 [hep-th].


	usNew
	C. Gomez, R. Letschka, and S. Zell, Infrared Divergences and Quantum Coherence,  http://arxiv.org/abs/1712.02355arXiv:1712.02355 [hep-th].







penrose
R. Penrose, Relativistic Symmetry Groups,  in NATO Advanced
		Study Institute – Group theory in non-linear problems Istanbul, Turkey,
		August 7 - 18, 1972, pp. 1–58.
1972.
In R. Pensore, Collected Works Volume II, Oxford University Press.



grover
L. K. Grover, Quantum Mechanics Helps in Searching for a Needle in a
	Haystack,  http://dx.doi.org/10.1103/PhysRevLett.79.325Phys. Rev.
	Lett. 79 (1997) no. 2, 325–328,
http://arxiv.org/abs/quant-ph/9706033arXiv:quant-ph/9706033.


2toN
G. Dvali, C. Gomez, R. Isermann, D. Lüst, and S. Stieberger, Black
		hole formation and classicalization in ultra-Planckian 2 → N
		scattering,
http://dx.doi.org/10.1016/j.nuclphysb.2015.02.004Nucl. Phys. B 893 (2015)  187–235,
http://arxiv.org/abs/1409.7405arXiv:1409.7405 [hep-th].



stromingerMasslessQED
T. He, P. Mitra, A. P. Porfyriadis, and A. Strominger, New Symmetries of
		Massless QED,  http://dx.doi.org/10.1007/JHEP10(2014)112J. High
	Energy Phys. 2014 no. 10, 112,
http://arxiv.org/abs/1407.3789arXiv:1407.3789 [hep-th].


mohd
A. Mohd, A note on asymptotic symmetries and soft-photon theorem,
http://dx.doi.org/10.1007/JHEP02(2015)060J. High Energy Phys. 2015 no. 2, 60, http://arxiv.org/abs/1412.5365arXiv:1412.5365 [hep-th].


campiglia
M. Campiglia and A. Laddha, Asymptotic symmetries of QED and Weinberg's
	soft photon theorem,  http://dx.doi.org/10.1007/JHEP07(2015)115J.
	High Energy Phys. 2015 no. 7, 115,
http://arxiv.org/abs/1505.05346arXiv:1505.05346 [hep-th].

stromingerQED
D. Kapec, M. Pate, and A. Strominger, New Symmetries of QED,
http://arxiv.org/abs/1506.02906arXiv:1506.02906 [hep-th].




stromingerRevisited
D. Kapec, M. Perry, A.-M. Raclariu, and A. Strominger, Infrared
	Divergences in QED, Revisited,
http://arxiv.org/abs/1705.04311arXiv:1705.04311 [hep-th].





FK
V. Chung, Infrared Divergence in Quantum Electrodynamics,
http://dx.doi.org/10.1103/PhysRev.140.B1110Phys. Rev. 140
	(1965) no. 4, B1110–B1122.

T. W. B. Kibble, Coherent Soft-Photon States and Infrared Divergences.
	II. Mass-Shell Singularities of Green's Functions,
http://dx.doi.org/10.1103/PhysRev.173.1527Phys. Rev. 173 (1968) no. 5,
	1527–1535.

T. W. B. Kibble, Coherent Soft-Photon States and Infrared Divergences.
	III. Asymptotic States and Reduction Formulas,
http://dx.doi.org/10.1103/PhysRev.174.1882Phys. Rev. 174 (1968)
	no. 5, 1882–1901.

T. W. B. Kibble, Coherent Soft-Photon States and Infrared Divergences.
	IV. The Scattering Operator,
http://dx.doi.org/10.1103/PhysRev.175.1624Phys. Rev. 175 (1968)
	no. 5, 1624–1640.

P. P. Kulish and L. D. Faddeev, Asymptotic conditions and infrared
		divergences in quantum electrodynamics,
http://dx.doi.org/10.1007/BF01066485Theor. Math. Phys. 4 (1970)
	745–757.







	veneziano
	M. Ciafaloni, D. Colferai, and G. Veneziano, Emerging Hawking-like
		Radiation from Gravitational Bremsstrahlung Beyond the Planck Scale,
	http://dx.doi.org/10.1103/PhysRevLett.115.171301Phys. Rev. Lett. 115 (2015) no. 17, 171301, http://arxiv.org/abs/1505.06619arXiv:1505.06619 [hep-th].

	M. Ciafaloni, D. Colferai, F. Coradeschi, and G. Veneziano, Unified
		limiting form of graviton radiation at extreme energies,
	http://dx.doi.org/10.1103/PhysRevD.93.044052Phys. Rev. D 93
		(2016) no. 4, 044052, http://arxiv.org/abs/1512.00281arXiv:1512.00281 [hep-th].

		A. Addazi, M. Bianchi, and G. Veneziano, Glimpses of black hole
				formation/evaporation in highly inelastic, ultra-planckian string
				collisions,  http://dx.doi.org/10.1007/JHEP02(2017)111J. High
			Energy Phys. 2017 no. 2, 111,
		http://arxiv.org/abs/1611.03643arXiv:1611.03643 [hep-th].



BN
F. Bloch and A. Nordsieck, Note on the Radiation Field of the
		electron,
http://dx.doi.org/10.1103/PhysRev.52.54Phys. Rev. 52 (1937)
	54–59.





hobson
M. P. Hobson, G. P. Efstathiou, and A. N. Lasenby, General Relativity: An
	Introduction for Physicists.
Cambridge University Press, 2006.
	
	mathematica
	Wolfram Research, Inc., Mathematica 11.0,  2016.
	

	


