Hawking Radiation of Extended Objects
George Johnson John March-Russell



Multi-Antenna Channel Interpolation via Tucker Decomposed Extreme Learning Machine
    Han Zhang, Student Member, IEEE, Bo Ai, Senior Member, IEEE, Wenjun Xu, Senior Member, IEEE, Li Xu, Member, IEEE, and Shuguang Cui, Fellow, IEEE
	Manuscript received February 24, 2019; revised April 05, 2019; accepted April 23, 2019; date of current version April 25, 2019. The work was supported in part by NSF with grants CNS-1824553, DMS-1622433, AST-1547436, and ECCS-1659025.
	Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org. 
	H. Zhang and S. Cui are with the Department of Electrical and Computer Engineering, University of California, Davis, CA 95616, USA (e-mails: {hanzh, sgcui}@ucdavis.edu), Bo Ai is with the State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing 100044, China (e-mail: aibo@ieee.org), Wenjun Xu is with the Key Lab of Universal Wireless Communications, Ministry of Education, Beijing University of Posts and Telecommunications(e-mail: wjxu@bupt.edu.cn), and Li Xu is with the Institute of Computing Technology, Chinese Academy of Science, Beijing 100190, China(e-mail: lixu@ict.ac.cn).
    December 2017 (Updated Feb 2018 and Dec 2018)
========================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================





§ INTRODUCTION


Since Hawking's original discovery <cit.> in 1974 that black holes emit radiation with a temperature inversely proportional to their mass, there has been a wealth of research into the precise spectrum of this radiation. In particular, there have been detailed studies of the rate of emission of particles of different masses and spins, for black holes charged and uncharged, spinning and static, and in different dimensions.

Several different derivations of the Hawking temperature exist in the literature; whilst the original analysis considers the scattering of wavepackets off the horizon towards future null infinity, Hawking and Gibbons, for instance, approached the topic from the perspective of the Euclidean path integral for the gravitational field <cit.>. In 1999, Wilczek and Parikh showed that Hawking radiation can be viewed as a tunnelling process <cit.>, in which particles pass through the contracting horizon of the black hole, lending formal justification to the most intuitive picture we have of black hole evaporation.

All such studies, however, assume the emitted particles are pointlike. In nature, of course, many particles are extended objects. Bound states of QCD, for instance, have finite size. Solitons, such as monopoles in gauge theories or Q-balls in scalar field theories, are not pointlike either. Even non-composite particles, such as fundamental strings or even black holes themselves, can have finite extent.

In this paper we study the effects of the extended nature of such objects on the temperature of the black hole and the precise spectrum of emission. The finite size of these objects modifies their interaction with the gravitational field; in particular, they are subject to tidal forces. We expect that these forces may alter the process of tunnelling through the horizon, as well as the subsequent evolution of the particles in the black hole background. For small black holes, such as the primordial black holes speculated to have existed in the early universe <cit.>, these tidal forces should be large enough as to modify the spectrum of emission. In this paper we restrict attention to uncharged and nonrotating black holes, and consider only spinless emitted particles.

Deviation from the pure black-body spectrum predicted by Hawking's original calculation can be described by frequency-dependent greybody factors. They arise due to the interaction of the emitted particles with the background gravitational field, after having been produced outside the horizon. These greybody factors are one of several effects that modify the pattern of radiation from that of a purely Planckian spectrum, whose supposed absence of correlations is the source of the black hole information puzzle <cit.>.

This paper is divided into five sections. In Section 2, we discuss the general theory that determines the motion of extended objects in curved spacetimes. In particular, we establish an expansion of the action in powers of the curvature of the spacetime, that we can truncate to lowest order. In Section 3, we compute the temperature of a black-body that radiates extended objects, and show that it coincides with the usual Hawking temperature. In Section 4 we compute the greybody factors associated with emission of such particles, at least in the long wavelength regime. In this regime, only the s-wave solution is relevant. Finally, in Section 5, we discuss the results, and mention the implications for various finite size particles we observe, or hypothesise to exist, in nature.

We take c=ħ=1 throughout. We write M for the length GM = M/M_P^2.



§ THEORY OF EXTENDED OBJECTS


We will study the dynamics of extended objects in a background gravitational field  by approximating them as point-particles with some effective action that is a modification of the typical relativistic point-particle action. All the effects of tidal forces on the body can be captured in a systematic fashion through terms coupling to the curvature of the metric. If we assume the extended object is spherical, the problem of writing down terms compatible with rotational symmetry, worldline reparametrization invariance and general coordinate invariance is straightforward <cit.>. The first few terms in such an effective action, linear in the curvature, are 

    S = -m ∫dτ - c_1 ∫dτ  R - c_2 ∫dτ  R_μνẋ^μẋ^ν + ⋯

where c_1, c_2 are constants, R_μν is the Ricci tensor, and R the Ricci scalar. For Schwarzschild spacetime, which will be the case we consider in this paper, R and R_μν are zero. The lowest-order non-vanishing terms in the action are constructed from the Riemann tensor, and involve two powers of the curvature. We will employ this effective action:

    S = -∫dτ(m + c  E_μνE^μν + c' B_μνB^μν)  ,

where

    E_μν   = R_μανβẋ^αẋ^β , 
    
    B_μν   = ϵ_μαβγ R^αβ_δνẋ^γẋ^δ ,

and c,c' are assumed to depend on the size of the object. In general the action depends on four powers of the particle velocity. Assuming purely radial motion, we find for Schwarzschild spacetime that B_μν=0[We emphasise that for more general motion, B_μν is not zero.], and the effective action can be written

    S = -∫dτ(m + c 6M^2/r^6((1-2M/r)ṫ^2 - ṙ^2 - 2 ṫṙ√(2M/r))^2 )  .

We note that such objects no longer move on geodesics; the deviation from geodesic motion is due to the aforementioned tidal forces.



§ HAWKING TEMPERATURE FOR EXTENDED OBJECTS

In this section we study the semi-classical tunnelling of extended objects through the black hole horizon, in line with the analysis of <cit.>. See also <cit.> for further discussion. The rate of emission Γ will have exponential part given by

    Γ∼exp(-2  S)  ,

where S is the tunnelling action. According to the Planck law, particles with frequency ω should be radiated at a rate with exponential dependence Γ∼exp(-ω/T). We can hence read off the temperature at which the black hole radiates according to

    = ω/2   S .


We will study the process in Painlevé-Gullstrand coordinates. The time coordinate t, just as for Schwarzschild time, corresponds to the time measured by a stationary observer at infinity. Ignoring angular directions, the metric is

    ds^2 = -(1-2M/r)dt^2 + 2 √(2M/r)dt dr + dr^2  ,

where M is the mass of the black hole. 

We will firstly present a calculation of the imaginary part of the action for a massless point-particle tunnelling through the horizon of Schwarzschild spacetime in the spirit of <cit.>. We will then repeat the calculation for a massive particle, and finally present the calculation for extended objects. In all three cases we find that the temperature of the black hole is given by the usual Hawking formula =1/8 π M. We can argue on simple physical grounds that this must be the case. In particular, if a given black hole is in thermal equilibrium with a bath of particles of two different species, then by the zeroth law of thermodynamics, those two baths of particles must be in equilibrium with each other, and hence have the same temperature. It is nevertheless useful to see how the Hawking temperature falls explicitly out of the calculation in each case.

We can write the action for a particle moving freely in a curved background as

    S = ∫ p_μ dx^μ    with     p_μ = g_μνx^νσ ,

where σ is an affine parameter along the worldline of the particle, chosen so that p^μ coincides with the physical 4-momentum of the particle. For a massive particle, this requires that dσ = dτ/m, with τ the proper time.



 §.§ Case m=0
   

The radial dynamics of massless particles in Schwarzschild spacetime are determined by the equations

    (1- 2M/r)ṫ^2 - 2√(2M/r)ṙṫ - ṙ^2 = 0   , 
    (1- 2M/r)ṫ - √(2M/r)ṙ = ω .

The second equation is the geodesic equation corresponding to the time-independence of the metric; in terms of the momentum defined in Eq. (<ref>), it can be written p_t = ω, and so ω has the interpretation of the energy of the particle as measured at infinity. The first equation can be factorised to yield, for an outgoing particle, the equation

    rt = 1 - √(2M/r) .

We then have

    S    =  ∫ω dt+ ∫ p_r  dr 
       = ∫( √(2M/r)ṫ + ṙ) dr 
       = ∫ṫ( √(2M/r)  + rt) dr 
       = ∫ṫ dr  .

Using Eqs. (<ref>) and (<ref>) we can write ṫ in terms of ω to yield:

    S = ∫ω/1-√(2M/r) dr  .

The integrand has a pole at r=2M, the horizon. Choosing the prescription to integrate clockwise around this pole (into the upper-half complex-r plane), we find

    S = 4 π M ω ,

giving =1/8 π M, as expected.



 §.§ Case m >0

For massive point-particles, the dynamics are governed by

    (1- 2M/r)ṫ^2 - 2√(2M/r)ṙṫ - ṙ^2 = 1  ,
    (1- 2M/r)ṫ - √(2M/r)ṙ = ω .

Using these equations to solve for p_r is somewhat more complicated than before, but essentially the same. Taking Eq. (<ref>) and substituting for ṙ into Eq. (<ref>) yields a quadratic equation in ṫ which we can solve:

    ṫ^2 (r/2M-1) - ṫ(ω r/M) + (rω^2/2M +1) = 0  .

We find

    (1-2M/r) ṫ   = ω±√(2M/r)√(ω^2 - 1 + 2M/r) , 
    ṙ   = ±√(ω^2 - 1 +2M/r) .

From here we compute the imaginary part of the action thus:

    S    = ∫ p_r  dr
       = ∫( √(2M/r)ṫ + ṙ) dr 
       = ∫(ṙ + ω√(2M/r)/1-2M/r + ṙ) dr  .

Taking the positive sign for ṙ and substituting yields

    S = ∫(√(ω^2 - 1 + 2M/r) + ω√(2M/r)/1-2M/r) dr  .

Finally we integrate around the pole at r=2M as before. And as before, we find

    S = 4 π M ω .





 §.§ Case c>0

In order to calculate the tunnelling rate for particles of finite size, we first need to determine the modifications to the equations of motion that result from the additional term in the action, before substituting these back into the action. Matters are simplified greatly for radial motion; using the standard normalisation of the 4-velocity, the action is simply

    S = -∫dτ(m + c 6M^2/r^6) 
    ≡ -m∫dτ  f(r)  .
 
As per usual, we can use the 4-velocity normalisation in place of the equation of motion for r(τ), leaving us with only one equation to determine. The t-independence of the effective action makes this relatively simple (cf. (<ref>)):

    ω = ((1-2M/r)ṫ - √(2M/r)ṙ)f(r)  .

The manipulations of the previous section proceed here as before, but with ω replaced everywhere by ω/f(r). We thus compute 

    S = ∫ f(r)  p_t  dt + ∫ f(r)  p_r  dr  .

The first term has no imaginary part. The second term is given by

    ∫ f(r)  p_r  dr =
    ∫ f(r)(√((ω/f)^2 - 1 + 2M/r) + (ω/f)√(2M/r)/1-2M/r) dr  .

At the pole, the factors of f(2M) cancel, yielding

    S = 4 π M ω ,

as before.



§ GREYBODY FACTORS


We have found that a Schwarzschild black hole with mass M radiates extended objects at the Hawking temperature . We next wish to find the exact rate at which such objects are radiated. For a perfect black-body emitting massless particles, the power received per solid angle subtended per unit frequency is given by the Planck law:

    P(ω) = Aω^3/8 π^31/exp(ω/k_B T)-1 ,

where A is the projected area of the black-body. For a non-gravitating sphere of radius r, we would have A=π r^2. For a black-body which is not perfect, this relation is modified to

    P(ω) = σ(ω)ω^3/8 π^31/exp(ω/k_B T)-1 ,

where σ(ω) is known as a greybody factor. One finds that this greybody factor is precisely the absorption cross-section for the appropriate particle incident on the black hole. In particular, if the black hole is to be in thermal equilibrium with a bath of particles surrounding it, the rate at which it emits these particles must equal the rate at which particles in the bath fall into the black hole. The relation between the absorption cross-section σ and the absorption probability 𝒜 is given by the usual relation

    σ = π/k^2∑_ℓ(2 ℓ + 1 ) 𝒜(ℓ)   ,

where ℓ is the partial wave number and k is the particle's momentum, given by k = √(ω^2 - m^2). We will calculate the quantity 𝒜 by solving the field equations for the relevant particle in the black hole background, imposing purely ingoing boundary conditions at the horizon, corresponding to total absorption.

We will be interested exclusively in the low-momentum regime,  k M ≪ 1. This corresponds to the wavelength of the particle being much larger than the Schwarzschild radius of the black hole.  Analytic results for both massless and massive particles are well-known in this regime. Furthermore, the greybody factor in the low-momentum regime will be determined by the ℓ = 0 partial wave, for which there is a great simplification of the tidal terms in the action, seen in Eq. (<ref>).


An analytic expression for the greybody factor in the regime where the de Broglie and Compton wavelengths of the particle are much larger than the Schwarzschild radius of the black hole (that is, ω M ≪ 1) was found in <cit.>:

    σ(k) = 32 π^2 M^3 (2k^2 + m^2)√(k^2 + m^2)/ k^2(1 - exp(-2 π M 2k^2 + m^2/k) )^-1 .

For low-momenta massive particles, this becomes

    σ(k) ≃32 π^2 M^3 m^3/k^2 .

In the low-energy massless case, the limiting behaviour is rather different. Taking the m=0, k → 0 limit of (<ref>), we have

    σ(k) ≃ 16 π M^2  .





 §.§ Effective Field Theory of Extended Objects


Recall that the action for a radially infalling extended object can be written 

    S = -∫dτ(m + c 6M^2/r^6) .

Though the action in principle contains four powers of the velocity, for radial motion these appear precisely in the  combination of the 4-velocity to the fourth power (for more general motion this is not the case). We note that this amounts to the action for a relativistic point-particle, but with a position-dependent effective mass

    μ(r) = m + c6M^2/r^6 .

The corresponding classical field theory that describes spinless coherent collections of these particles is thus the Klein-Gordon theory with such an effective mass:

    ℒ = -1/2(∂_μϕ)(∂^μϕ) - 1/2μ(r)^2ϕ^2  .

We can formalise this correspondence as follows. First note that we can introduce an einbein e(x) into the action to modify it thus[Indeed, this is the more appropriate action one should use for massless particles.]:

    S = 1/2∫dσ( e^-1 g_μνẋ^μẋ^ν- eμ(r)^2 )  .

Upon quantising this theory of a point-particle we find the effective action <cit.>

    Γ[g_μν] = 1/2Trln(∇^2 - μ(r)^2)  ,

which is precisely the effective action of the Klein-Gordon theory (<ref>). The classical equation of motion derived from this Lagrangian is simply

    g^μν∇_μ∇_νϕ - μ(r)^2 ϕ = 0  .



 §.§ Solving The Equations of Motion

It remains to solve equation (<ref>) subject to appropriate boundary conditions. We will restrict our analysis to the ℓ = 0 partial wave solution of this equation, since the equation is derived from the assumption of radial particle motion. In Schwarzschild coordinates (now using t to mean Schwarzschild time), we decompose our field ϕ as 

    ϕ(x) = e^- i ω t R(r)  .

Then equation (<ref>) reads

    1/r^2(1-2M/r)∂_r (r^2(1-2M/r)∂_r )R(r) - μ(r)^2(1-2M/r)R(r) + ω^2 R(r) = 0  .

We can simplify this equation using the tortoise radial coordinate r_*:

    r_* = r + 2M ln(r/2M-1)         r_*r = (1-2M/r)^-1 .

Then the equation of motion becomes, with S = rR <cit.>,

    -Sr_* + V(r)S = ω^2 S     with     V(r) = (μ(r)^2 + 2M/r^3)(1 - 2M/r)  ,

where r is implicitly a function of r_*. To calculate the absorption probability, we impose purely ingoing boundary conditions at the horizon. When r = 2M, this equation becomes that of a simple harmonic oscillator, with solution

    R(r_*) = A_+ exp(i ω r_*) + A_- exp(-i  ω r_*)  .

Ingoing boundary conditions correspond to A_+ = 0. At large distances, the effects of the tidal terms are negligible. In the massless case, the solutions are spherical Bessel functions, whilst the long-range 2M/r term in the massive case yields Coulomb wavefunction solutions, which in the far-field region have the asymptotic form <cit.>:

    R(r) = B_+/rexp(i k r - i ηln (2 k r) + i θ) + B_-/rexp(-i k r + i ηln (2 k r) - i θ)  ,

where η=- M(2k^2+m^2)/2k and θ is an unimportant constant phase.

Starting with boundary conditions at the horizon, we numerically integrate this equation out to large r to extract the reflection coefficient

    ℛ = |B_+/B_-|^2  .

The absorption probability we desire is then simply

    𝒜 = 1 - ℛ ,

and from here the greybody factor can be computed according to Eq. (<ref>).




 §.§ Results


We present in this section the results of the numerical simulations. We first consider the regime where the mass and energy of the particle are small, in the sense that ω M ≪ 1 and m M ≪ 1. Thereafter we consider the large mass (mM ≫ 1) case.



  §.§.§ Low Mass Emission



The effect of the finite size of the object is always to reduce the greybody factor σ. We plot in Figure <ref> the behaviour of the greybody factor as a function of momentum k, in both the massive and massless case. In particular, we find that, irrespective of the wavelength of the particle, the effect of the tidal terms is to reduce the absorption probability by a constant factor which we call X(c):

    X(c) = σ(c)/σ(0) .




    
We find that for small c, X(c) is equal to one, as expected, whilst for large c, X(c) tends to zero, suppressing the greybody factor:

    X(0) = 1  ,
    lim_c→∞X(c) = 0  .

We wish to know how this quantity behaves as a function of the parameters m and M. We find that as a function of the dimensionless ratio

    s = c/M^3 ,

the quantity X(s) is independent of the masses m and M. The form of X(s) as a function of s is plotted in Figure <ref>. 



















The point at which the greybody factor is reduced to half, s_1/2, is given by 

    s_1/2≃ 5.9  .

We find that beyond this point, the tail of the function X(s) is well-modelled by an exponential. Numerically we determine that

    X(s) ∼exp(-s/ s_0)     with     s_0 ≃ 3.3  .


We can understand this scaling relation by rescaling Eq. (<ref>) in terms of a dimensionless radius ρ = r/M (and its corresponding tortoise coordinate). We have

    -Sρ_* +  ((mM + 6c/M^3ρ^6)^2 + 2/ρ^3)(1 - 2/ρ)S = ω^2 M^2 S   .

The dimensionless parameters in this equation are mM, ω M and c/M^3. The effective mass is a rapidly decaying function of ρ, and so we expect the short distance behaviour to depend only on c/M^3, and the large distance behaviour on m and ω. In principle, these parameters could be mixed together in a complicated way when we match in the intermediate regime. However, for small masses and energies, the absorption probability is expected to be very small (cf. Eqs. (<ref>), (<ref>) and (<ref>)). Equivalently, the coefficients B_+ and B_- in (<ref>) are almost equal and opposite. The matching of the two solutions in the intermediate regime thus amounts to a single continuity condition, relating the size of the small r solution (which depends only on c/M^3) and the size of the large r solution (which depends only on mM and ω M). This explains the factorisation of the cross-section into the product of the point-particle cross-section and the tidal effects.

We can also understand this relation another way, examining the effective potential that the particle moves in. Eq. (<ref>) resembles the 1D Schrödinger equation for a particle moving in a potential well. For the particle at infinity to reach the horizon, it must tunnel through a potential barrier, which exists for all values of c for sufficiently light particles. We plot this effective potential for certain values of s in Figure <ref>. We hence expect the absorption probability to have an exponential suppression of the form

    𝒜∼exp(-2 ∫dr_* √(V(r;c)-ω^2)) ≡exp(-2 I(c))  .

For small m, ω, we find

    I(c) ≃∫_2^∞dρ√(((6c/M^3ρ^6)^2 + 2/ρ^3)(1 - 2/ρ)^-1) .

We see again that this integral depends on c only  through the dimensionless parameter c/M^3. Furthermore, numerical evaluation of this integral as a function of c shows that it remains roughly constant until around c/M^3 ∼ 5, and then rises approximately linearly, with gradient 1/3.3=1/s_0. This explains the exponential suppression of the absorption probability. We plot this integral against s on the left of Figure <ref>.










  §.§.§ High Mass Emission


In the high-mass regime, the restriction kM ≪ 1 corresponds to consideration of only non-relativistic particles. We can understand the behaviour of the absorption probability in this regime by considering once again the behaviour of the effective potential. It was shown in <cit.> that for m greater than a critical value m_crit≃ 0.385/M, the effective potential becomes a monotonically increasing function of r. Above this mass, therefore, there should be a large enhancement of the absorption cross-section due to the disappearance of a barrier to tunnel through. However, for a sufficiently large tidal term, the potential once more develops a barrier, and we expect exponential suppression of the cross-section. 

For large m we find a simple analytic relation for the value of s=c/M^3 at which this barrier develops:

    s_1/2≃1/2(3/2)^11 m M  .

This relation is remarkably well borne out by numerical simulations for mM ≳ 1. We can again examine the integral I(c) of our effective potential in the regime where s is greater than s_1/2, and a barrier exists. As before, the integral depends linearly on the value of s, and so we expect exponential suppression of the cross-section for s above s_1/2. The gradient is approximately 1/3.3 as before, independent of the masses m and M. Thus 

    X(s) ∼exp((s-s_1/2)/s_0)     with     s_0 ≃ 3.3  .

We plot the barrier integral against s on the right of Figure <ref>.





















§ DISCUSSION


We have found that Hawking radiation of spatially extended objects always occurs at a lower rate to that of the corresponding pointlike objects, at least for scalar particles in the low-momentum regime. It is interesting to ask whether this effect is significant for extended objects in nature, such as protons, strings, GUT monopoles, or other objects. To estimate this, we need to find some relation between the physical size of these objects and the corresponding parameter c in the action.

The parameter c has dimensions of volume. For a given object, there are three characteristic length scales: its Compton radius 1/m, its gravitational radius m/M_P^2, and its physical size d. For highly gravitationally-bound objects (we can consider the possibility of black holes radiating black holes themselves), the physical size is on the order of the gravitational radius, and a matching calculation indicates that c ∼ M_P^2 d^5 <cit.>. For non-gravitational objects, we expect d and 1/m to be the only relevant scales. However, we also expect the effects of the tidal forces to be well-behaved in the massless limit[For many extended objects, such as protons and GUT monopoles, the Compton and physical radii are comparable, so this argument is unnecessary.], implying no dependence on the Compton wavelength: c ∼ d^3.

We compute in the following table the value of the black hole mass M (in Planck units) that results in an order two-fold suppression of the emission rate for various extended objects. The size and mass of the Q-ball are typical values taken from <cit.>.




Object      Mass m     Size d     M/M_P 

Fundamental String     0     10^-35 m     1 

Proton     1 GeV     10^-15 m     10^20 

Q-ball     10^8 GeV      10^-17 m    10^16 

GUT Monopole     10^16 GeV     10^-31 m     10^4 

Micro BH     m     m/M_P^2      m/M_P 






Finally, we must discuss the limitations of this analysis. Firstly, by reducing the theory of extended objects to an effective theory involving couplings to the curvature, we implicitly assume the object is rigid. For non-rigid objects, the point-particle approximation is inappropriate. However, this analysis should be valid whenever there is a sufficiently large mass gap in the spectrum of the extended object. For sufficiently low energies, we can expect the particle to remain in its ground state, and hence behave as a rigid object. We should also expect that for very large curvatures, it is not sufficient to take only the first term in the expansion of the action, and that higher-order terms become relevant.

Secondly, we have restricted attention to the long wavelength regime. In this regime the effects on the spectrum are dominated by s-wave radiation, for which there is a large simplification of the effective action. Furthermore, we assumed in our analysis that the emitted particles are spinless, which led us to an effective theory of a scalar field. We refer the analysis of higher spins and higher partial waves to future work.

Thirdly, we have only considered radiation from Schwarzschild black holes. It would be interesting to study how finite size effects modify the spectrum of rotating or charged black holes, or of black holes in different spacetime dimensions. For charged black holes radiating charged objects, for instance, there will be forces on the object as it escapes to infinity due to inhomogeneities in the electric field as well as those in the gravitational field, which would likely give rise to a rich range of phenomena.



§ ACKNOWLEDGEMENTS

GJ is supported by the Science and Technology Facilities Council (STFC).


JHEP
