Multiplicities from black-hole formation in heavy-ion collisions

The formation of trapped surfaces in the head-on collision of shock waves in conformal and non-conformal backgrounds is investigated. The backgrounds include all interesting confining and non-confining backgrounds that may be relevant for QCD. Several transverse profiles of the shocks are investigated including distributions that fall-off as powers or exponentials. Different ways of cutting-off the UV contributions (that are expected to be perturbative in QCD) are explored. Under some plausible simplifying assumptions our estimates are converted into predictions for multiplicities for heavy-ion collisions at RHIC and LHC.

Data from heavy-ion collisions from RHIC and LHC have revolutionized our perception of strong coupling physics in QCD, and revealed the characteristics of the deconfined phase. They have also become the testing ground of novel techniques emerging in string theory that attempt to control strong coupling phenomena using a gravitational description. In this description, a heavy-ion collision is described as the process of black hole formation and decay, albeit in a five-dimensional theory of gravity including also other fields, notable a scalar, the dilaton, [1]. In this direction, several attempts have been made to analyse the scattering of highenergy sources, using shock waves in AdS, [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20]. 1 A direct outcome of this approach, is the estimate of lower bounds of the final multiplicity by using Penrose's idea of trapped surfaces.
In this paper, we analyse the formation of trapped surfaces in head-on collisions of shock waves in gravitational theories with more complicated bulk dynamics, and different types of vacuum solutions. An example are Einstein-dilaton theories with a scalar potential. Such theories have been argued, [26][27][28], to describe holographic physics that is closer to QCD than the AdS theory. A phenomenological theory, Improved Holographic QCD (IHQCD) has been constructed, [29,30], that agrees well with both zero temperature and finite temperature YM data.
The idea is to explore how different aspects of the dynamics affect the trapped surface that forms during the collision of shock waves. The shock waves used are gravity (spin-2) shock waves. The different factors are the following: • Different bulk geometries. There are several different possible geometries that that have been classified in [26,27]. They are characterized by their IR and UV (nearboundary) behavior.
In the UV the typical behavior is asymptotically AdS. There is however interest in different non-AdS asymptotics, as they can capture the physics of the collision when the trapped surface forms mostly in the IR part of the geometry.
The IR behavior can be split into three large classes: (a) Confining geometries. (b) Non-confining geometries. (c) Unacceptable geometries (that violate the Gubser bound 2 ). AdS 5 belongs to class (b). There is another special geometry that lies at the boundary of (a) and (b), the linear dilaton geometry.
• Different profiles of the transverse distributions of energy. So far two types have been analyzed, uniform distributions or power-like (aka GPY, [32]) profiles. We will add one more class namely exponential profiles, that are well localized in the transverse JHEP04(2012)065 plane. Such energy profiles are closer to what the targets are in heavy-ion collisions. Except the uniform distributions, the other profiles are characterized by a length scale that controls the size of the energy distribution. Typically GPY profiles are fuzzier as they stretch to larger distances, while exponential profiles are sharply localized.
• The option of cutting off the UV part of the bulk geometry, [32]. This is motivated by the fact that we are striving to emulate the QCD behavior that is perturbative in the UV. In a perturbative regime, by definition, multiplicities are small. It is in this regime that the geometry is expected to be unreliable in a holographic description. The simplest way to implement this is to stop the geometry at a point deemed to be the transition to the perturbative regime, and ignore the contributions above that point. This is not a controlled approximation, but it is expected to give useful hints. 3 A shock wave is an idealization that is a good approximation for bodies moving with a velocity near that of light. As one is boosting the gravitational field of a point source to the speed of light, the field is shrinking into a transverse cone, and in the limit it is localized at x − = 0. The gravitational scattering of two bodies, in that limit, concerns the determination of the space-time metric that results from two such initial shock waves, initially far apart from each other, and moving towards each other.
In a generic case of the gravitational scattering of two bodies, the space-time metric is the sum of the two metrics when they are far apart, but is deformed already as they are approaching. The important simplifying property of the shock-wave limit is that the (gravitational) interaction that deforms the metric of two shock waves from the simple sum of the two metrics, does so after they cross. Therefore there is a well defined "crossing point", adjusted usually at x + = x − = 0. In the first quadrant of the (x + ,x − ) plane, (x + < 0,x − < 0) the metric is the sum of the two shock wave metrics. In the last quadrant, (x + > 0,x − > 0), the metric is deformed because of the collision.
The power of Penrose's trapped surface argument is that sometimes one can find a trapped surface that does not depend on the details of the interacting metric after the collision, where the "after" was defined above as the crossing point of the shock waves. Because of this, the calculation of the trapped surface depends alone on the original form of the shockwave metrics. Therefore, if the shock-wave metrics can be found, then the problem is solvable. The Penrose ansatz that we employ in (3.2) captures the form of such a surface and has been used in all such studies in many contexts.
However, it is not guaranteed that the trapped surface must always be of that form, and in numerical examples it can form in the part of the geometry that is the fourth quadrant, and therefore depends non-trivially on the metric "after" the interaction. Obviously, without solving for that metric we do not have access to that region and the formation of trapped surface therein. Therefore, if a trapped surface of the type (3.2) exists, then it can provide a lower bound on the entropy. If it does not, then no statement can be made in this formalism, and an analysis of the full collision is necessary.

JHEP04(2012)065
Our approach is to solve first the equations for the shock waves, in the presence of different types of metrics, with sources in the bulk. Such solutions determine also the transverse distributions. We then find the associated trapped surfaces, and the estimated lower bound on the generated entropy. 4 We will use the word entropy liberally in this paper to mean the area of the trapped surface. Strictly speaking this is a lower bound on the entropy, but it will always be obvious which entropy we are referring to. Moreover, entropy can be converted to multiplicity, and it is in this sense that we will use this two terms interchangeably.
Two characteristic examples are analyzed quantitatively and eventually compared to data.
1. The first is AdS 5 , with a transverse distribution having a localized exponential profile and a cutoff in the UV at r = 1/Q s with Q s the saturation scale (section 8.2.2). We will name this setup "AdS-Q s " for future reference.
2. The second is a simplified metric emulating the IHQCD solution of [26,27], implementing confinement and asymptotic linear glueball trajectories, while it is asymptotically AdS 5 in the UV. In this case we have again a localized exponential profile in the transverse plane and no UV cutoff, (section 8.2.3). We will name this setup, with a slight abuse of language, "IHQCD".
What we find is as follows: • Both, our analytical calculations and our numerical analysis show that most of the entropy comes from the UV part of the geometry provided that the geometry reduces to the AdS 5 space near the boundary. This is not in contradiction with the naive expectation that in a strongly coupled theory the multiplicity comes from the highenergy part of the phase space available during the collision.
• For uniform transverse distributions, the AdS 5 geometry, produces the least entropy among the geometries of section 5 . For non-confining geometries the entropy scales as s a for large s, with 1 3 ≤ a < 1 2 with a = 1 3 for AdS 5 . For scaling confining geometries 1 2 ≤ a < 1. The case a = 1 2 corrected by logs corresponds to the geometries related to IHQCD. We have assumed above that the other geometries persist up to the boundary, or equivalently the trapped surface forms in the regime in which they are valid. These results are summarized in table 2. However, as most realistic geometries are asymptotically AdS, the thrust of this result may be important only at intermediate energies, where IR geometries may take over.
• For geometries that have a mass gap, discrete spectrum and confinement, the allowed transverse distributions has a spectrum of scales that is in one to one correspondence with the discrete spectrum of 2 ++ glueballs. The associated entropy production, independently of transverse distribution, is less than AdS 5 . It is assumed that all such confining geometries are asymptotically AdS 5 . Table 1. The several cases analysed and compared. Two geometries have been compared: A non confining geometry (the AdS 5 ) and a confining (IHQCD) one. None has a cutoff. The transverse profiles correspond to the cases of being constant, GPY (falling-off as a power law, [32]) and fallingoff exponentially.
• We also find that the scattering of a distribution with transverse size associated with the lowest lying 2 ++ glueball does not seem to lead to a (Penrose type) trapped surface that forms before the collision of the shock waves. In such a case, numerical techniques are necessary to investigate the full metric after the shock-wave collision.
• There are geometries with a UV energy-independent cut-off that lead to an asymptotic ∼ log 2 (s) behavior for the entropy (see table 3).
• A general trend is that at equal total energy, the collision of distributions that have a larger transverse size leads to a larger entropy production. This implies that more dilute energy distributions produce more entropy at fixed total energy. In particular a uniform transverse energy distribution produces (at equal total energy) more entropy that one with power-law or exponential transverse distribution. Similarly a powerlaw transverse distribution at the sane total energy produces more entropy compared to an exponential transverse distribution with the same length scale.
• Comparing an AdS geometry, with an asymptotically AdS geometry that is confining in the IR we find that for all different types of transverse energy distributions that we have examined (uniform, power-law, or exponential), the AdS geometry generates substantially more entropy than the confining geometry, at the same (and large ) total energy, and transverse scale.
More to the point, table 1 defines the different contexts studied in this paper. We find the following inequalities between the various trapped surface areas (keeping total energy and transverse size fixed) • In confining backgrounds, the entropy production increases as the confinement scale Λ QCD decreases, provided that the total energy and the transverse size are kept fixed.

JHEP04(2012)065
• In the AdS-Q s setup the multiplicities grow with the atomic number almost linearly; in particular as ∼ A 17/18 [64].
• The lower bound on multiplicities for AdS-Q s is given in equation (8.4). For the IHQCD setup it is given in for the (8.6). Both formulae once fit to RHIC data, make the same (correct) prediction 2.76 TeV PbPb LHC data, [64]. A single parameter is used in these fits, namely the overall constant coefficient of the leading (large) s-dependence.
The implications of the conclusions above do not a priori apply verbatim to QCD because: 1. We are discussing glueballs only, and not nuclei.
2. The holographic backgrounds used are approximations to various gauge dynamics.
Only AdS is an exact background for sYM, but all backgrounds, including AdS used for QCD, are approximations.
3. The final entropy produced may be and is usually larger than the trapped surface area.
The first issue may not be of prime importance at high energy. The reason is that the energy released at mid-rapidity in heavy-ion collisions is expected to be mostly gluons. Quarks will provide corrections to this but they are not expected to change this picture drastically.
The second issue is of importance. Although by now bottom-up models of YM can provide reliable calculations matching lattice calculations, their structure in the UV is more shaky. It is not a priori clear where the geometric description breaks down, and such transitions may be at different places for different observables.
The third issue is also important for quantitative predictions. There are very few cases where the final entropy has been calculated numerically in the collision of shock waves, [11,16]. In [16] in particular it was shown that the released entropy at the end of the collision process is 60% larger that the bound found from the trapped surface calculation. Moreover conformal invariance in AdS makes this percentage to be independent of the collision energy.
This suggests that for collisions and bulk geometries where the majority of the trapped surface area come for the UV, AdS part of the geometry, we should still expect that the relative factor relating the final entropy release to the area of the trapped surface to be almost energy-independent.
The analysis in this paper is providing important and potentially general clues on multiplicity generation in high energy collisions using holography. It must be however be backed-up by a more reliable calculation of the gravitational evolution, but we leave this for future work. Another interesting issue is whether there are differences between the high-energy scattering of 2 ++ glueballs studied in this paper and 0 ++ glueballs. Although we do not have reasons to expect major differences, an analysis in this case should be done.

JHEP04(2012)065
We organize this write-up as follows: We begin in section 2 by showing the way one may build more realistic geometries that come closer to QCD. This is achieved by considering scalar gravity. In section 3, we state the equations that compute the entropy production of the collided shock-waves and hence of the gauge matter indicating that the calculation reduces to a boundary valued problem. In section 4 we present the subclass of geometries we investigate. Sections 5 and 6 refer to various geometries with uniform or non-uniform transverse dependence respectively. The entropy of the trapped surface is computed for each case. Section 7 deals with the simplest way of removing the weak coupling entropy production: At higher energy scales, where coupling is weak, the contribution to the entropy production should be less important. In the geometrical language of AdS/CFT this would imply that one should modify the way the fifth coordinate of AdS/CFT at the UV is treated. In section 8, we present our results and make our predictions. Finally, in section 9, we comment on various aspects of our investigations and conclusions. We omit intermediate details of our calculations for the three appendices at the end. In particular, appendix A serves as a practical introduction to the theory of trapped surfaces. Appendix B shows how one may localize the (five-dimensional) bulk sources and appendix C proves a useful equation that is needed for subsection 6.2.

Einstein-dilaton gravity
We start from the action: where Φ s is the scalar field dual to the YM coupling constant. We first find shock wave solutions in this theory of the form 2) with the asymptotically AdS boundary at r = 0. Compatibility of these equations implies that ∂ + Φ = 0. Eliminating Φ s using the equations of motion, the equation for φ is where we have introduced a stress-tensor J ++ . It should be stressed here that the solutions in 2.3 are exact solutions to the equations stemming from (2.1) modulo the presence of point-like of string like-sources in bulk along the lines of [38]. Potential backreaction corrections may appear after the collision, but the trapped surfaces we investigate can be traced to before the collision of the shock waves. It should be stressed that the source stress tensor J ++ is strictly conserved, ∂ − J ++ = 0, and remains so even if we make arbitrary distributions of shock waves (no x − dependence on the metric).

JHEP04(2012)065
3 Trapped surfaces Trapped surfaces are created when two shocks like the one of (2.2) 5 which moves along x − collide. 6 In terms of metrics before the collision, one then has Associated with the shock-wave φ 1 in (2.2), we parametrize (half of the) trapped surface S 1 by where ψ 1 remains to be determined. 7 It is also useful to rescale the functions φ 1 and ψ 1 by defining Ψ 1 satisfies the following differential equation where AdS 3 is defined with respect to the metric We point out that once (2.3) is solved, then ( AdS 3 −A(b(r)))Φ 1 provides a source term for ( AdS 3 −A(b(r)))Ψ 1 . The missing ingredient is the boundary conditions that are given by for some curve C which defines the boundary of the trapped surface and where both, S 1 and (the associated surface to φ 2 ,) S 2 end. The entropy is then bounded below by the area of the surface obtained by adjoining the two pieces of the trapped surface associated with each of the shocks as 5 We call it φ1 to distinguish it from the second that moves along x + that we call φ2. In our case the shocks will be taken identical and hence the subscripts will be soon dropped. 6 A more complete set of notes on the theory of trapped surfaces may be found in appendix A. 7 We are looking for Penrose-type trapped surfaces. Each (one of the two) piece(s) of Penrose surface moves along the light cone and surrounds the point-like source that creates the shock; this is precisely what the left equation of (3.2) describes. The trapped surface closes at the intersection of the two pieces. This happens at x + = x − = 0. See also [32], section 8 for a relevant discussion.

JHEP04(2012)065
where the (generalized) curve C defines the boundary of the trapped surface S 1 and S 2 which are identical; thus the overall factor of 2. The integral with respect to the transverse coordinates gives x 2 ⊥ (r) when considering a head-on collision. Typically, r C 1 , r C 2 and x ⊥ carry the information of the shock φ. 8 We conclude by pointing out that the un-scaled version of the trapped surface (see (3.3)) for a head-on collision is given by while the entropy is still given by (3.7). In the absence of transverse dependence one ignores x 2 ⊥ from (3.7) which measures entropy/transverse area in this case.

Shock geometries
We have analyzed a number of bulk geometries by examining a class of different scale factors, that classify the non-conformal behavior of Eistein-dilaton gravity models. In terms of their behavior in the far IR (r → ∞, or r → r 0 ), we have the following cases [26,27,29]: The AdS case corresponds to a = 1. This corresponds to quasiconformal geometries, with no confinement, continuous spectrum and a mass gap, with potential asymptotics as Φ s → ∞, V ∼ e QΦs , Q < 4 3 .
2. Confining backgrounds that are scale invariant in the IR, [57], with b(r) ∼ (r 0 − r) a , a > 1 3 . In this case r 0 is finite and signals the position of an IR singularity that satisfies the Gubser bound for a > 1 3 . They have a discrete spectrum of glueballs and a mass gap. The potential asymptotics as Φ s → ∞ are V ∼ e QΦs , Q > 4 3 .
3. Confining backgrounds with b(r) ∼ e −(Λr) a , a > 0. They have a discrete spectrum and a mass gap. The potential asymptotics as Φ s → ∞ are V ∼ e

Shocks with uniform transverse space dependence
In this case the shock φ(x + , r) can be determined and the trapped surface ends at r = r H with r H determined from 9 , and s the center of mass energy of the collision. The area of the trapped surface is 8 We drop the subscripts 1, 2 from φ's and ψ's from now on. 9 The delta function in ( 5.1) may be in principle replaced by any function of x + .

JHEP04(2012)065
We may therefore estimate the energy dependence of the trapped area for different bulk geometries.
• For non-confining scaling theories b ∼ r −a , with a ≥ 1 we obtain A trap ∼ s 3a−1 6a . The AdS case corresponds to a = 1. This agrees with previous estimates for the AdS case, [32].
at high energy. The exponent varies between 1 2 and 1. • at high energy.
Taking into account that S trap ∼ A trap we conclude that in all cases the entropy production is larger than AdS at high enough energies.

Non-uniform transverse dependence
We assume the shock(s) have non-trivial transverse dependence and we solve the homogeneous equation (2.3) by separating variables. We obtain the following set of differential equations (assuming rotational symmetry on the transverse plane) The first equation yields while the solution to g(r) depends on the scale factor b(r). Strictly speaking, the solutions corresponding to K 0 do not exactly solve the homogeneous differential equation for f k . They induce a delta function source and they satisfy It is evident that J ++ from (2.3) consists from an appropriate linear combination of K 0 (kx ⊥ )g k (r)'s. In the remaining of this section we analyze several b(r)'s.

Power-like transverse distributions
We will study solutions to (6.1) where the transverse distributions fall-off as a power in the transverse plane. Such solutions were considered for AdS in [32]. In the remaining of this section we analyze several different types of scale factors b(r) and we derive the boundary that specifies the trapped surface as a function of the energy (see equations (6.7) and (6.13)). Then using (3.6) one may compute the entropy bound S; the final answers for S are summarized in section 8 (see tables 2, 3, 4).

Power-law b(r)
We begin by studying the power-like scale factor When a > 1 3 and r 0 finite, this corresponds to a confining geometry, with IR singularity at r = r 0 , [26,27]. This can be resolved by uplifting to higher dimensions, [57]. It can also include non-confining geometries, if r 0 = 0 and a < −1. a = −1 corresponds to AdS. The boundary for a > 1/3 is at r → ∞ while for a < 0 is at r → 0.
Assuming a point-like where P µ ν is the associated Legendre polynomial and where the tilde onΦ is to remind us that the rescaling of φ 10 is not by a factor of b(r) to the first power (see (3.3)) and will be omitted from now on. The presence of the combination q, natural in AdS space, [32] for more general power-like metrics, can be explained, as such metrics are conformally AdS, and can be lifted to AdS metrics in a higher dimensional space-time, [57].
The equation for the trapped surface satisfying the first boundary condition is given by where Φ − is the solution to the homogeneous differential equation satisfied by Φ which is finite at q = 0. Q µ ν is the associated Legendre polynomial of the second kind while q C denotes the value of q at the boundary of the trapped surface. The second boundary condition in (3.8) yields where the quantity on the right of the arrow denotes the high energy limit. We note the agreement with [32] when a = −1, r 0 = 0 and r ′ = L. The S trap is given by (6.8b) 10 And likewise, the rescaling of ψ.

JHEP04(2012)065
The limits of integration are found from (6.7) when x ⊥ = 0. Performing the integrations and keeping the leading contribution in E taking into account the range of a for each case, we finally arrive at As a cross check we note that (6.9b) for a = −1 reproduces the result of [32]. We remark that the trapped surface for a > 1/3 exists because r 0 is finite.

Exponential b(r)
We will take the scale factor b(r) in this case to be This corresponds to a marginal confining geometry, that of a linear dilaton, [26,27]. Rescaling φ = e 3r 2R Φ(u(r, x ⊥ )) 11 and assuming a point-like J ++ as in the power-low case, we obtain where r ′ is the r coordinate of the source. The trapped surface equation is while the rescaled boundary condition 4 R 2 e 3r/R u(Ψ ′ (u)) 2 C = 8 finally yields where the quantity on the right of the arrow denotes the high-energy limit. The trapped surface is found by working as in subsection 6.1.1. In this case also, the trapped surface exists either when r is allowed to extend in the interval [−∞, ∞] or when a UV cut-off is placed (see discussion below (6.8)). Allowing for r to take negative values, so that the geometry is extended towards the boundary, the entropy is computed from (6.13) numerically: One, solves this equation with respect to x 2 ⊥ = x 2 ⊥ (r, E) and integrates as in (3.7). The two limits of integration are found numerically from (6.13) again by setting x ⊥ = 0. The analysis shows that S trap ∼ s 1.66 log 1.17 (s). (6.14) The result indicates a larger growth of S trap than the uniform profile case (see third bullet of section 5). This, does not contradict our general conclusions as here the exponential metric is assumed to be valid up to the UV boundary. There is a way to derive a simple formula 11 And likewise, rescaling ψ.

JHEP04(2012)065
for S trap when a UV cut-off (see section 7) is placed. If for instance, we cut the surface at some r > 0 then equation (6.13) implies that as E increases, u C increases logarithmically: This result looks similar with Froissart bound estimates for cross sections [2,3]. We do not know if there is a connection.

Localized transverse distributions
We built a single-mode shockwave φ k out of the solutions of (6.1) demanding to be square integrable for all x ⊥ . The result is where E is the energy carried by the shock. g 1,2 (r), the radial solutions, are dimensionless, because φ k has dimensions of length. The overall constant ensures that the gauge stress tensor, which according to the AdS dictionary is given by when is integrated in space, it yields the total energy E. The trapped surface is defined by and yields for some coefficients C k to be determined. The K 0 's do not participate as they induce source terms. 12 If we assume that at small r, g 1 → r 4 Π(r) with Π(r) a regular function with at least one real root, and that g 1 has no multiple roots then, as it is proved in appendix C, only the coefficient C 1 k is non-trivial. This assumption turns out to be true for the cases we are studying in this paper. The boundary condition for the trapped surface then satisfies (6.20) The high-energy limit of the trapped surface is defined by JHEP04(2012)065

Discrete spectra
We consider a confining scale factor b which asymptotes to AdS in the UV and behaves as in IHQCD in the IR, Although the precise solution, that fits YM data is slightly different, [40], we will use the one in (6.22) as it does not affect the high energy asymptotics. This solution has been discussed in appendix G of [26,27]. The length scale R is the analogue of the QCD scale. We may parametrize it as R ∼ Λ −1 QCD . Unlike non-confining cases, in confining bulk geometries the physical spectrum of glueballs, namely the fluctuations of the metric and the scalar dilaton, is discrete, [26,27]. The radial equations that define both the shockwave profile, in equation (2.3), as well as the equation for the trapped surface, (3.8) involve the same radial equation as the one that determines the spectrum of spin-2 fluctuations (2 ++ glueballs).
When it comes to determine the presence of a trapped surface, a similar condition appears: there is no trapped surface unless the solutions to (3.8) are normalizable (see (6.24)). This is directly obvious in shockwaves determined by a single graviton wavefunction and can be shown in general.
As normalizability of the wavefunction is important for the existence of the trapped surface, the spectrum of shockwave profiles is discrete, with the transverse momentum determined from the 2 ++ glueball masses m n by |k n | = m n . Therefore, the transverse profiles of the shockwave distributions, not surprisingly, are determined by the normalizable 2 ++ (graviton) wavefunction.
For the metric (6.22), when k takes the particular subset of values the radial solution of equation (6.1) reduces to a (finite) polynomial that behaves as r 4 at small r and is normalizable. Normalizability in the radial direction is defined when the corresponding eigenfunction g 1 (k n r) satisfies b(r) 3 |g 1 (k n r)| 2 dr < ∞.
The set of the normalizable eigenfunctions is given by where L (2) n are the (finite) associated Laguerre polynomials of degree n. This is the qualitative behavior in any background that is confining with a discrete spectrum of glueballs and a mass gap. The values in (6.23) coincide with the mass-spectrum of 2 ++ glueballs.
We now analyze the example (6.22) and study the formation of a trapped surface via collision of such quantized transverse distributions. The case n = 0 gives a shock JHEP04(2012)065 x ⊥ R r Figure 1. The (closed) boundary x ⊥C (r C , E) of the trapped surface for the first excited state corresponding to b(r) = L/re −r 2 /R 2 for fixed energy E. There is a family of such curves; one for each value of E which moves upwards for larger E. The section x ⊥ = 0 is a part of the boundary as should because the source has coordinate x ⊥ = 0. The surface does not extend beyond r > R which implies that there is no entropy production for energies greater than 1/R, a scale which is naturally identified with Λ QCD . that behaves as g 1 ∼ r 4 for all r and hence it has only a single root (see discussion after (6.21). We therefore conclude that this ground state mode produce no Penrosetype trapped surface before the collision of the shock waves. This does not preclude the formation of a trapped surface once the interacting metric after the collision is used but our tools do not allow us to investigate this possibility.
We next considering the n = 1 case which yields the two independent solutions The solution g 1 has the right asymptotics in order to generate a closed trapped surface for all energies. The shock given by g 1 has r 4 behavior at small r in accordance with the expectation value of the gauge theory stress-energy tensor. The boundary along the r direction is given by r C 1 = 0 and r C 2 = R. Indeed, when r takes these two values, the numerator of the left-hand-side of (6.20) becomes zero and compensates the zero of the denominator when x ⊥ = 0 as figure 1 depicts. The (two pieces of the transverse) exponential profile trapped surface (without any cut-offs) yield 2 (n, ER; y)dy, ERκ 2 5 g 1 (y 12(n + 2)) 8πL 3 24(n + 2) =x ⊥ I 0 (x ⊥ ) (6.27)

JHEP04(2012)065
where y 0 is the highest root of g 1 (y 12(n + 2)) = 0 and n denotes the n th -excitation. For the case at hand, n = 1 while g 1 is given by (6.26) with y = r/R. It is evident that the entropy, unlike [32], depends not only on the transverse size R/ 12(n + 2) but in addition on the the confinement scale 1/R. In fact, these two parameters may be varied independently.
We may compare with the analogous AdS calculation where the transverse profile is taken to be the same as here, with characteristic scale k. The confining theory has another scale 1/R that can be traded with varying the integer n. We find that for any n the area of the trapped surface in AdS is always larger than that in the confining background.

Accommodating asymptotic freedom
As was pointed out in [38], the UV (small r) should not contribute importantly to the entropy production S. The reason is that by definition for perturbation theory to be valid, the generated particle multiplicities must be small. Many examples are known when large multiplicities imply the breakdown of perturbation theory, with the sphaleron case the most prominent one, [60].
We expect that in the QCD UV, at some point perturbation theory takes over, and this is defined as the regime in which the generation of entropy is small compared to that generated from lower scales. This transition we will approximate as an abrupt transition: we will assume that this is a radial position r U V , below which we can use the gravitational description, while above it standard perturbation theory takes over. We will neglect the perturbative contribution to the entropy as we expect it to be small. Therefore in this approximation we will introduce a UV cutoff r = r U V in the trapped surface that will simulate the emergence of weak coupling in the UV. The position of r U V must be determined, and at this stage it appears as an additional phenomenological parameter.
A related question is to what extend geometries with varying coupling constant like IHQCD implement the fact that interactions are weak near the UV, as we would have expected from QCD. The answer is that in asymptotically AdS backgrounds, even as the string coupling e φ → 0 in the UV, the graviton interactions remains strong. It would have been probably different if the geometry becomes AdS in the string frame. However in this case, (a) the boundary geometry is singular (b) Such a case cannot be a solution to a gravitational action with two derivatives only.
In the sequel we will cutoff the geometry in the UV and explore the result.

Energy-independent cut-off
The work [38] suggests that in the high-energy limit, the entropy, for the geometries studied, should be given by

JHEP04(2012)065
where C 1 is an energy independent constant. 13 When the energy dependence of x ⊥ is of the form x 2 ⊥ = i c i (E)c ′ i (r), as it is usually the case, the last equation reduces to In particular, the discussion above and equation (6.7) imply that the geometries (6.1.1) with a UV (constant) cut-off yield We remark that for a > 1, the UV cut-off is placed on the upper bound of the integral in (7.1) as the boundary theory is at r = ∞ in this case. Evidently, this procedure, modifies the (center of mass) energy dependence (that we denote by s) of S. In the case of [38], it reduces to S ∼ s

Energy-dependent cut-off and the saturation scale
We will now consider the shocks with uniform transverse dependence for simplicity.
In all examples we have analyzed, in the high-energy limit the trapped surface produces the following entropy where b 0 (r) is the asymptotic form of b for small r. r C = r C (E) is the (lower) boundary of the trapped surface and it is determined by the boundary conditions (see e.g. (6.20)). It is generically energy-dependent. We believe, that this behavior is much more general. Equation (7.4) implies that most of the entropy originates in the UV part of the trapped surface. As the energy becomes larger this part enters into the region of the asymptotic freedom where the coupling is small and where we do not expect a large multiplicity to be produced. Therefore, as we have argued in the beginning of this section, we will impose asymptotic freedom by cutting-off the trapped surface at some r 0 as in [32]. It is natural to expect that r 0 may be energy-dependent. We propose as a natural cut-off the saturation scale Q s (see figure 2) by identifying r 0 ∼ 1/Q s i.e.
where √ s N N is measured in GeV and denotes the c.m. energy in nucleon-nucleon collisions.
In this equation A is the atomic number while typically λ lies in [0.1, 0.15] for energies at 13 In [38] the authors consider in addition rC 2 (E) → C2 in order to remove the IR contribution to the entropy production. In our case, this has already been taken into account by the confining IR geometry.

JHEP04(2012)065
Freedom IR Entropy Production Figure 2. The entropy production occurs between the scales Λ QCD ∼ 1/R and Q s which in the geometric language corresponds to 1/Q s < r < R. The Λ QCD has already been incorporated by the usage of a confining scale factor b(r/R) as in (6.22). The weak coupling regime on the other hand can not be implemented rigorously in the framework of the AdS/CFT duality because string theory is not weakly coupled at high energies where 't Hooft coupling is small. The phenomenological Q s cut-off is proposed as a first approximation. RHIC (and LHC) depending from the nature of the nuclear matter (pp, AA etc.) participating in the collision [43][44][45][46][47][48][49][50][51][52][53][54][55] and √ s N N is measured in GeV's. The interval for the values of λ is obtained from fitting data from independent processes, e.g. Deep Inelastic Scattering (DIS) while they are close to the results predicted by analytical calculations [47]. In the uniform transverse dependence case and assuming that at small r the geometry maps to AdS 5 then cutting the surface at 1/Q s yields It is interesting to note the Q 2 s dependence for any scale factor b(r) that has the asymptotic behavior b 0 = L/r! The same dependence of entropy with Q 2 s was found recently in [56] for uniform transverse nuclei in the context of a different approach. In addition, multiplicity densities dN/dη are proportional to ∼ Q 2 s as these are predicted on theoretical grounds from perturbative methods [47]. η denotes the pseudorapidity while dN/dη ∼ Q 2 s describes the data quite well [48,50,51].
However, when non-uniform transverse dependence is present, we can not reach such a general conclusion as the one provided by (7.7) (compare with (8.4) for instance).

Behavior at large s
Using the asymptotic form of the boundary conditions, equations (6.7), (6.13), and equations (3.7), (6.27), (7.1), (7.2) and (7.5) wherever appropriate we analytically or numerically compute the entropy production for all the cases we have considered. The results are summarized in tables 2, 3 and 4.

JHEP04(2012)065
We remark that: (a) It makes sense to cut the surface at r s ∼ 1/Q s if and only if r s > r C 1 placing restrictions on the allowed values of a in the power-low case (see caption of table 3 and footnote 16).
(b) The power-law shocks for a = −1 (see [32]) when the surface is cut at 1/Q s yields S trap ∼ r ′ Q s (r ′ A) 1/3 s 1/6 ∼ A 17/18 s 0.24 (for AA; see (8.4)) where r ′ is the radius of the nucleus A. This energy dependence describes data well (see plots of figure 3). In addition, an almost linear dependence with the number of participants (see figure 5), in PbPb collisions at 2.76 TeV has been observed [64]. The result for 2.76 TeV concerns the ALICE experiment at LHC and is, up to this time, a preliminary result (see figure 5).
(c) The case of (6.22) has been analyzed for the exponential (in the transverse direction) profile corresponding to k 1 . The results found fit satisfactorily the RHIC data 14 up to 200 GeV (see plots of figures 4).
(d) The geometry corresponding to (6.10) is also interesting. In this case a trapped surface may be formed when a cut-off is placed in the UV region. The resulting entropy then behaves as S trap ∼ log 2 (s) (see [2,3]) at large s even when the cut-off of the trapped surface is energy-dependent. We remind the reader that the geometry (6.10) is a "marginal" case corresponding to continuous spectra with a mass gap, [26,27].
(e) We have numerically derived a set of inequalities about the entropy production. In all the cases analyzed and compared (see below), the energies are taken identical and large: 15 (a) For an exponential transverse profile corresponding to the geometry of (6.22) we find S k 1 trap > S k 2 trap > S k 3 trap > . . .; ER = fixed. Taking into account that k n sets the transverse size, we conclude that more dilute transverse distributions at fixed energy, result in more entropy.
(b) For an exponential transverse profile corresponding to the confining geometry of (6.22) we find that as ER → ∞ with E/k fixed, the entropy increases and becomes that of the geometry b = L/r (AdS 5 ) (see (9.1)).
(c) For the geometry b = L/r we compare the (trapped) entropy of a shock with exponential (transverse) profile (exp shock) with a shock having a power-like transverse profile, [32] (GPY shock). We assume that both of the transverse profiles fall-off for x ⊥ > 1/k; the first falls-off exponentially (as K 0 (kx ⊥ )) while the second as a power, 1/(x 2 ⊥ + 1/k 2 ) 3 . We find S GY P trap > S exp trap and we conclude (again) that more dilute energy distributions produce more entropy at the same total energy.

JHEP04(2012)065 b(r)
Confining Sources S with S with S with no cut cut at cut at Table 2. Classification of trapped surfaces and entropy production S = S trap at high-energies of shocks with uniform transverse distribution. In all the cases, the last column assumes that b(r) reduces to L/r at small r. The first column of the table displays the large r asymptotics. The cut-off at the last two columns refers to a cut-off of the UV region of the trapped surface. b(r) Con-Sources S with S with S with fining no cut cut at a cut at  . We conclude that confined matter produces less entropy than conformal matter at infinite coupling. This is accord with basic intuition. We point out however, that the difference in the two entropies is subleading at high energy because most of the entropy is produced at the UV where the two geometries coincide (L/re −r 2 /R 2 ≈ L/r for small r.).

Relating S trap with multiplicities
The trapped surface analysis does not give the produced entropy but it provides a lower bound S trap ≤ S prod. .    [20,200] GeV. In the same energy interval, c 1 ranges in [300,775]. An accurate fitting is found (see left plot of figure (4)).
introduction and commented on how much each of them is expected to affect the connection with the data. In particular we have argued that for collisions and bulk geometries where the majority of the trapped surface area comes from the UV, AdS part of the geometry, we should still expect that the relative factor relating the final entropy release to the area of the trapped surface to be almost energy-independent.
We must also quantify the relation of the total multiplicity and the produced entropy. The total entropy is given by the the number of charged particles N ch times ∼ 3/2 to account for the neutral particles multiplied by ∼ 5, that is the entropy per particle [32,58,59]. Hence, We have analyzed the two cases: (a) The case of the AdS geometry with a cut-off at the UV (at ∼ 1/Q s where Q s the saturations sale). (b) The case of the confining IHQCD-like geometry (6.22) for n = 1 (first excitation) without any UV cut-off. It should be stressed here that a single parameter is used in these fits, namely the overall constant coefficient of the leading (large) s-dependence.

Multiplicities for the AdS-Q s setup
This setup has an AdS metric, a GPY-like transverse profile and a UV cutoff at r ∼ 1/Q s with Q s given by (7.6). The gravity parameters are chosen according to [32] as L 3 /G 5 ≈ 1.9. We need in addition the following relations where A is the atomic weight of the participating nuclei in the collision and √ s N N the center of mass energy/nucleon. 16 16 We drop the subscript NN from √ sNN from now on for simplicity. We will cut-off the trapped surface at r s = c/Q s , where c a positive constant. 17 One then computes (at large s) where we normalized the formula of S trap with A Au = 197 for AuAu collisions. In normalizing, we used the fact that r ′ is the transverse size of the colliding glueball 18 (beam) assuming that it satisfies the empirical law (r ′ /r Au ) = A 1/3 /A
The limits of the r integration are found by setting x ⊥ = 0 (see figure 1). The excellent fitting of plot 3 with RHIC data is achieved for A=A Au = 197 and (1/c)×(overall coefficient 19 )≈ 1.54 taking r Au = 8 fm. 20 It is pleasing that both, the overall coefficient and c can simultaneously be of order one. The extrapolation to higher energies is done in section 8.3 and the right plot of figure 3.

Multiplicities for the IHQCD setup
This setup involves a metric b = L/re −r 2 /R 2 without a UV cut-off and an exponential transverse profile. In this case, S trap is given by (6.27) for n = 1 and has to be solved 17 As mentioned earlier, it makes sense to cut at c/Qs iff c/Qs > rC 1 = r ′ /2 A √ sG5r ′ /L 3 −1/3 (r ′ is the transverse nuclear size). For c ∼ 1, A any (reasonable) value (see (7.6)), r ′ ≈ (A/AAu) 1/3 × 8 fm and √ s ≥ 20 measured in GeV, this condition is satisfied.  Figure 4. Both figures concern the numerical analysis of the surface (resulted from the first excitation) of (6.27) for lower and higher energies (dashed plots). There is not a UV cut-off in this case.
In the left figure, the green plot is given by (8.6) and exhibits the accurate approximation of the numerical plot found for A=A Au . The overall coefficient of the numerical plot has been chosen in order to fit the RHIC data. The right plot is the extension of the approximating plot at higher energies for A= A P b and hence, according to (8.6), it is given by N ch = 79.8 ( √ s) 0.451 log 0.718 (562 √ s).
The points for √ s = 2.76, 5.5 and 7 TeV are inserted.
numerically. The parameters of gravity were determined by matching lattice data in [40] and we will use these in the sequel. Choosing N c = 3 we find where √ s is measured in GeV. In order to go higher in the energies, a more refined numerical analysis is needed. Hence, for the present work, we use the fitted curve given by equation (8.6) in order to predict multiplicities for higher energies. We do this in section 8.3 for the energies to be reached by LHC taking into account that for Pb A= 207. The corresponding plot is the one on the right of figure 4. We remark that cutting the surface at some UV cut-off does not improve the fitting. In particular, for energies higher than 200 GeV we either find very low multiplicities (for a constant cut-off) or very large multiplicities (for an energy dependent cut-off).
We close this section by noticing that in order to obtain a (more realistic) value for in the case of the trapped surface given by (6.27), the black hole ansatz for the scale factor of (6.22) should be solved. We postpone a detailed analysis with the precise IHQCD profiles for later work and choose the particular values mentioned below equation (8.2) in order to fit the data.

Multiplicities at LHC energies
In the previous section we have analyzed (a) the geometry of AdS 5 with a UV cut-off and (b) the first excited state of (6.22) without any UV cut-off. We have seen that both fit the RHIC data up to 200 GeV satisfactorily (see plots of figures 3 and 4). For the aforementioned geometries (a) and (b), we may extrapolate them at higher energy in order to assess what they predict for the multiplicities at the energies reached by LHC. For high-multiplicity 22 proton-proton (pp) central collisions where A= 1 we find N p ch ≈ 70 for 0.9 TeV, N p ch ≈ 110 for 2.36 TeV, N p ch ≈ 190 for 7 TeV and N p ch ≈ 260 for 14 TeV. It is pointed out that extracting experimental results (for total N p ch ) from ATLAS [61] and CMS [62,63] (for these energies) is not a trivial task. It involves model dependent procedures and Monte-Carlo simulations.
• Geometry (b) The numerical analysis corresponding to (8.6) (see right plot of 4) for PbPb colisions predicts N ch ≈ 19100 for 2.76 TeV, N ch ≈ 27000 for 5.5 TeV and N ch ≈ 30500 for 7 TeV. 23 In the above results, we have assumed exactly a zero impact parameter and hence these are zero centrality processes. This implies that we might predict slightly larger multiplicities than the upcoming data. The difference between the two predictions (at a given value of √ s) for the two cases considered above is ∼ 6% for PbPb collisions. It is interesting that the first preliminary results from ALICE at LHC for √ s = 2.76 TeV and PbPb collisions give N ch ≈ 17000 ± 1000 for 0 − 5% centralities [64]. In order to make a better estimation, we consider the right plot of figure 7 in [64] which shows that the maximum number of participants for the particular measurements, is approximately 380. Hence, plugging A= 380/2 = 190 and √ s = 2.76 TeV in (8.4) results N ch = 17300 which is extremely close to the (preliminary) measurements (see figure 5).

Outlook and discussion
The initial motivation of this work was to estimate the produced entropy of colliding shocks in different circumstances and obtain a qualitative answer for the results, in view of applications to heavy-ion collisions. In particular, the goal was to compute S trap and hence 21 All the energies that are mentioned in this subsection refer either to √ sNN or to √ spp. 22 The notion of centrality in pp collisions is trickier to define. The best definition is to select high multiplicities in the final state. The definition of high multiplicity is ambiguous, but may still defined. 23 A prediction for pp collisions using (8.6) is less reliable for the moment because the errors induced from factors of Ap/AAu = 1/197 are larger than in PbPb collisions. To deal with this, a more refined numerical work is needed.

JHEP04(2012)065
(estimate) the particle multiplicities (see section 8.2.1). We have analyzed various dilatongravity geometries, with or without transverse dependence, with or without confinement and with or without UV cut-offs. The cut-offs were assumed energy dependent or energy independent. We conclude the following: • Both our analytical calculations and our numerical analysis have shown that most of the entropy comes from the UV part of the geometry provided that the geometry is asymptotically AdS.
• There are geometries with a UV energy-independent cut-off that lead to an asymptotic ∼ log 2 (s) behavior for the entropy (see table 3).
• For uniform transverse distributions, the AdS 5 geometry, produces the least entropy among the geometries of section 5. This result assumes that the non-conformal geometries in 5 survive until the boundary. They are therefore not asymptotically AdS. A similar result would be valid if the associated geometries are asymptotically AdS, but the energies such that the trapped surfaces do not penetrate the AdS region. These results are summarized in table 2.
• We have constructed exact shock solutions, with non-trivial transverse distributions and point-like bulk sources, and have computed the entropy of the trapped surface (see table 3). Inserting an energy-independent UV cut-off (see fifth column of table 3), and choosing a ≈ 1 in the first line, 24 results in S ∼ ( √ s) 1/2 (see also (7.3a)). The particular energy dependence seems to describe RHIC data as equations (8.4), (8.6) and figures 3 and 4 suggest. A similar power-law dependence is obtained assuming a Landau hydrodynamical behavior [65] after the collisions . A characteristic of the geometries corresponding to a = 1 and −1/3 is that they do not reduce to AdS 5 at the UV unlike the a = −1 geometry.
• For geometries that have a mass gap and confinement, the entropy production, independent of transverse distribution, is subleading to that in AdS at the same total energy and transverse scale. Here it is assumed that all geometries are asymptotically AdS.
• The entropy production decreases when higher glueballs collide (see 6.2.1). This implies that more dilute transverse energy distributions produce more entropy at fixed total energy (ER = fixed, where R can be identified with 1/Λ QCD ).
• We will denote the trapped surface areas as follows: S GYP;AdS,k trap for the GYP-like shockwave in (6.5) in AdS 5 where k is the scale of the transverse profile, S k;Ads trap for an exponential transverse profile in AdS 5 as in (6.16), and S kn trap for an exponential transverse profile with scale k n for the confining metric (6.22). In all of the above the total energy, and the transverse scale are kept the same. It shows the total (charged) multiplicity (N ch ) as a function of A (see (8.4)) at fixed energy s. The red points are data taken from reference [64] with the error bars included: the horizontal coordinate of these points is taken to be N part /2 where N part is the number of nucleons (of the two Pb nuclei) participating in the collision at the fixed value of √ s = 2.76 TeV. The agreement of our theoretical prediction with the data seems to improve as N part increases, that is as the collision becomes more central.
We find the following inequalities for large E: We conclude that the entropy production increases as the confinement scale (∼ 1/R see subsection 6.2.1) decreases provided that the transverse size (with respect to the energy) is kept fixed.
• Equation (8.4) suggests that multiplicities should grow with A as A 17/18 . In fact, an almost linear dependence of N ch with the number of participants was observed in the recent ALICE experiments performed at 2.76 TeV [64] for PbPb collisions. According to figure 5, our result shows a similar behavior. The agreement becomes better as the number of participants increases, that is as the collision becomes more central which is the case that we assumed in this paper. We remark that the ALICE results are still preliminary.
• We have derived two formulae for multiplicities, for the AdS-Q s setup in (8.4) and the IHQCD setup in (8.6). We normalize the multiplicities using RHIC data and we then compare with known LHC data. We find that in both cases they are very close to the 2.76 TeV heavy-ion LHC data, [64] (see subsection 8.3).
• An important puzzle of our analysis is that the scattering of the lowest lying 2 ++ glueballs does not seem to lead to a trapped surface. This does not exclude however that such a surface would form later.
For future work, we propose to use the shock of equation (B.4) and estimate the produced entropy. This shock has the advantage of localizing the bulk sources and in JHEP04(2012)065 addition, it retains the ingredients of a mass gap and confinement. Also, more accurate calculations that will specify the gravity dimensionless parameter L 3 /G 5 (see e.g. (8.5)) could be made. Finally, it is important to understand the kind of black holes which give rise to the shocks (after boosting) and which in turn, mimic the nuclear matter in the colliders. This matter is best described by finite size black holes [1,67]. In particular, in [1] a description of the collisions in term of what is known about such black holes was described. It would be interesting to make this picture precise, and in this numerical work will be required.
Note added. According to the resend work [68] when the energy is large enough, the compact dimensions do not contribute to the entropy dependence as a function of the energy. Now, from the one hand, we know that for curved backgrounds the entropy S behaves as S ∼ E D−3 D−2 where D is the number of dimensions of space-time. 25 This implies that a smaller D (but larger than 5) results to a smaller entropy. On the other hand, for flat backgrounds, the reality is opposite: S ∼ E D−2 D−3 (see [41]). At high energies and for the flat background case with extra dimmensions, the entropy increases faster (than in the absence of the extra dimensions) as only the extended dimensions participate (D is smaller). This fact has been postulated in [69] and was later shown in [68].
Returning back to the AdS backgrounds, It would be interesting to investigate the scenario of extra dimensions and study how the presence of the S 5 or any other compact manifold (that surrounds the AdS background) would change the entropy growth and hence multiplicities.

JHEP04(2012)065
where ψ 1 will have to be determined. Before finding the differential equation satisfied by ψ 1 it seems necessary to perform the following coordinate transformation which eliminates the δ(x + ) from (A.1). In particular we use the transformation which transforms (A.1) into It is (sometimes) useful to exchange the functions φ 1 and ψ 1 with Φ 1 and Ψ 1 defined via The trapped surface associated with the shock φ 1 can be found solving the equation where a few explanations are in order: (a) The vector l (+1) ν 26 can be generally expressed by which implies that We emphasize that (i) the solution, B = 0 does not satisfy the rest requirements of l (+1) ν and hence it is discarded. (ii) Both of the subscripts ( i ) in the ∇ i 's of (A.12) are lower indices. We point out that condition (A.12) does not specify A and B uniquely but it only specifies their ratio Hence a choice for A and B is which is proportional (by the factor 1 b(r) 2 → r 2 L 2 ) with (101) of [32]. 27 One may wonder that this arbitrariness in choosing A and B may invalidate the procedure that determines the trapped surface and in particular equation (A.6). We show later that this is not the case and that (A.6) is invariant under rescalings of the vector l , l (+1) Explanation (b). One has then to compute the covariant derivative of l (+1) µ with respect to the metric (A.4) forming the tensor ∇ ν l (+1) µ . We do not display explicitly this tensor.
Explanation (c). The final step is to contract ∇ ν l (+1) µ with h µν . This yields where AdS 3 is defined with respect to the metric ds 2 = b(r) 2 dx 2 ⊥ + dr 2 .
(A.17) 27 Up to a factor of 2 for B which is traced to the authors different convention for the light-cone coordinates.

JHEP04(2012)065
Associated with a second shock φ 2 moving along x − there exists the surface S 2 constructed from a second vector l (+2) µ , the analogue of l (+1) µ . From symmetry considerations, l (+2) µ can be immediately found from l (+1) µ (see (A.15)) and hence reads Boundary conditions on C = S 1 ∩ S 2 . Ψ 1 | C = Ψ 2 | C (continuity of S on C) and . Assuming for simplicity identical shocks, we have Φ 1 = Φ 2 ≡ Φ and by symmetry Ψ 1 = Ψ 2 ≡ Ψ = 0 on C. Hence the boundary conditions we get finally are The shocks we have considered in subsection 6.2 correspond to a non-localized bulk source in the r direction (see (6.3)). In order to localize this source we begin from (6.3) and consider appropriate linear combination of solutions. This is not hard as one merely has to use the completeness relation of g k which schematically has the form dk C(k)g k (r)g k (r ′ ) = δ(r − r ′ ). (B.1) This would imply that the shock φ k should be given by φ = dk C(k)K 0 (kx ⊥ )g k (r)g k (r ′ ). (B.2)

JHEP04(2012)065
In particular the completeness relation ∞ n=0 n! Γ(n + a + 1) x ′a e −x ′ L (a) n (x ′ )L (a n )(x) = δ(x − x ′ ) (B.3) for the Laguerre polynomials implies that the correctly normalized shock 28 φ should be given by 29 In arriving to (B.4) we have used K 0 x ⊥ R 12(n + 2) d 2 x ⊥ = 2πR 2 12(n + 2) , L (2) n (0) = (n + 2)! n!2 , The idea is to note that the boundary of the surface should contain the point (r C , x ⊥C = 0) because the source is located at x ⊥C = 0 while both the terms (∂ r ψ k ) and (∂ x ⊥ ψ k ) should be finite at x ⊥C . We also state the following relations Beginning from (∂ x ⊥ ψ k ) we have that where we have used (C.1) in order to simplify the last expression. It is evident that g 1 (kr C ) should have (at least) two distinct real roots in order to have a trapped surface (see figure 1). For reasonable shocks, the g 1 decays for small r as r 4 and hence the one root is at r C = 0. By assumption of the claim (see subsection 6.2) there exists one more real root call it r C 2 and as a result the following is true lim r C →r C 2 g 1 (kr C ) = 0 (C.3) 28 Such that the gauge tensor T++ integrates to E. 29 For the AdS5 geometry, the sources are localized choosing C(k) = kJ2(kr ′ ). The shock then becomes that of [32]. This can be shown using results of reference [39].

JHEP04(2012)065
We now compute ∂ r ψ k again for x ⊥C → 0 and hence at r c → r C 2 . We have lim x ⊥C →0,r C →r C 2 ∂ r (ψ k ) = Eκ 2 5 L 3 k 2 lim x ⊥C →0,r C →r C 2 where in the second equality we have used (C.3). The simplified expression for ∂ r (ψ k ) implies a logarithmic divergence when x ⊥ → 0 because lim r C →r C 2 g ′ 1 (kr c ) = 0 by the hypothesis of the claim that there are no multiple roots for g 1 . 30 The fact that the righthand side of (C.4) is independent on the arbitrary coefficients C 1 k , C 2 k would yield to the naive conclusion that the whole second term of (6.19) should be absent. However, looking more carefully inside the sum of the numerator of (C.4), one realizes that the terms for g 1 (kx ⊥C ) = g 1 (k ′ x ⊥C ) do not actually participate in the sum and hence they are allowed. Thus the trapped surface equation is given by yielding to the trapped surface boundary determined by the condition (6.20) completing the proof of the claim.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.