Hot Spaghetti: Viscous Gravitational Collapse

Berndt M¨uller Department of Physics, Duke University, Durham, NC 27708-0305, USA

Andreas Sch¨afer Institut f¨urTheoretische Physik, Universit¨atRegensburg, D-93040 Regensburg, Germany

We explore the fate of matter falling into a macroscopic Schwarzschild black hole for the simplified case of a radially collapsing thin spherical shell for which the back reaction of the geometry can be neglected. We treat the internal dynamics of the in-falling matter in the framework of viscous relativistic hydrodynamics and calculate how the internal temperature of the collapsing matter evolves as it falls toward the Schwarzschild singularity. We find that viscous hydrodynamics fails when either, the dissipative radial pressure exceeds the thermal pressure and the total radial pressure becomes negative, or the time scale of variation of the tidal forces acting on the collapsing matter becomes shorter than the characteristic hydrodynamic response time.

I. INTRODUCTION Here we focus on the dynamics of matter falling into a black hole before it reaches the region of Planck scale space-time curvature. The study of the viscous heat- The fate of matter collapsing into a black hole [1] re- ing of matter collapsing into a black hole presented here mains one of the great riddles of theoretical physics. In remains fully within the framework of well established General Relativity the matter encounters a space-time physics. It makes use of the insight that quantum physics singularity within finite proper time after passing an imposes a lower bound on the ratio of the shear viscosity event horizon, and no further prediction of its evolution η to the entropy density s of the matter. The Kovtun- is possible. When quantum physics in the semiclassical Son-Starinets (KSS) conjecture [11] asserts that the low- approximation is added to this picture, the black hole est value of this ratio is 4πη/s ≥ 1 for matter governed eventually radiates away the mass of the collapsed mat- by any consistent relativistic quantum theory. The KSS ter in the form of Hawking radiation [2], but this effect bound is believed to be a strict bound for quantum field does not address the problem of what happens to the theories that admit an Einsteinian gravity dual, but is matter that has fallen into the black hole, creating the known to be modestly violated for holographic bound- information loss “paradox” [3]. ary theories of non-Einstein gravity models [12, 13]. We It is generally agreed that the resolution of both prob- do not consider the effects of bulk viscosity on the dy- lems, the singularity problem and the information loss namics of collapsing matter, because the bulk viscosity problem, requires a full quantum theory of gravity be- vanishes in the conformal limit, which is usually attained yond the semiclassical approximation. In one widely in elementary quantum field theories at high density or studied version of quantum gravity, loop quantum grav- temperature. ity, the interior space-time of the black hole “bounces” 4 While the effect of the shear viscosity is negligible in when the energy density reaches MP where MP is the Planck mass [4, 5]. In this way the singularity is avoided, the early stages of the collapse, it becomes relevant in the collapse is reversed, and the black hole radiates its the approach toward the space-time singularity as the mass away while the incipient event horizon dissolves gravitational shear forces acting on the collapsing mat- smoothly. A conceptual description of this process, in- ter grow rapidly. As we will show, the Navier-Stokes cluding a semiclassical model of the singularity-free evap- approximation to viscous hydrodynamics fails below a certain distance from the singularity, requiring the use of arXiv:1710.00136v2 [gr-qc] 11 Oct 2017 oration of the black hole, has been presented in ref. [6]. A simple extension of General Relativity that incorpo- the second-order, causal theory of relativistic hydrody- rates a limiting curvature has recently been proposed by namics. We obtain analytical solutions to the dynamical Chamseddine and Mukhanov [7], and its relationship to equations for shear stress and entropy density in a simple, loop quantum gravity was explored in ref. [8]. radially symmetric collapse model neglecting back reac- In string theory, much work has been devoted to an un- tion of the in-falling matter onto the gravitational field, derstanding of the microscopic origin of the Bekenstein which allow us to explore the effect of viscous heating entropy of black holes and how information can be en- all the way down to the region where effects of quantum coded at the event horizon. The dynamics of the collapse gravity become relevant. and ensuing evaporation process is not well understood. Our paper is organized as follows. First, we briefly Competing pictures of the black hole as a “fuzzball” of remind the reader of the effect of tidal forces on matter highly excited strings [9] or the event horizon as a “fire- falling into a black hole. We next summarize the equa- wall” [10] that would immediately melt additional incom- tions governing the gravitational collapse of matter in ing matter have been proposed, but their validity remains the framework of second-order relativistic fluid dynam- controversial and unresolved. ics. We then study the collapse of a thin shell into a 2

Schwarzschild black hole in a convenient comoving co- Whatever the mechanism driving the transformation ordinate system and solve the resulting equation for the is, during the final stages of its collapse the in-falling mat- shear stress in the limit of a constant ratio η/s = (4π)−1, ter forms a quark-gluon plasma, a fluid with the lowest where η denotes the shear viscosity and s is the entropy specific shear viscosity known and the greatest ability to density of the fluid and conclude with some speculations flow smoothly under conditions of large shear. Yet, as for about the implications of our result. any fluid, a hydrodynamic description of the quark-gluon plasma fails when it is forces to respond to changes on time or distance scales that are shorter than the relevant II. SPAGHETTIFICATION hydrodynamic scales. In this study we will attempt to an- swer two questions: One is whether the rapidly growing In the final phase collapsing matter is transformed into gravitational shear strain on the collapsing quark-gluon a quark-gluon plasma. The transformation can either plasma leads to a failure of the hydrodynamical descrip- be driven by compressional or viscous heating or, in the tion, and if so when; the other is, how hot the plasma absence of a sufficiently high entropy at the start of the becomes before it either enters a non-hydrodynamic do- collapse, will ultimately be driven by the action of tidal main or reaches the domain where quantum gravity ef- forces. The destructive shear effect of the tidal forces on fects become relevant. matter falling into a black hole is colloquially known as “spaghettification” [14]. The tidal force compresses in- falling matter in the transverse directions and stretches III. VISCOUS COLLAPSE it in the radial direction. The radial stretching first leads to the rupture of macroscopic solid objects, but can even In order to avoid the analytical complications of solv- pull apart atoms and eventually hadrons. ing the full nonlinear equations governing gravitational The radial tidal force over a distance ∆r on an object collapse, we consider the collapse of a thin spherical shell with mass m in the gravitational field of a spherically of matter into an existing black hole of mass M. We fol- symmetric mass M is: low the article by Adler et al. [16] on analytical models of gravitational collapse, which analyzes the collapse of mc2R fluid spheres. For spheres or shells of dust (pressureless ∆F = ∆r, (1) r3 matter) the particles fall freely in the Schwarzschild ge- ometry; when an equation of state p = αε with α > 0 where R = 2GM/c2 is the Schwarzschild radius, and G is used, the spherical shell needs to have a nonvanishing is Newton’s constant. Hadrons get disrupted when the surface tension for consistency. We will not be concerned tidal stretching force over the size of a hadron exceeds the here with this aspect but simply consider the fluid dy- QCD string tension κ ≈ 1 GeV/fm. For a very rough, namics in an interior region of the collapsing shell. Also order-of-magnitude estimate we set m ∼ Λ/c, ∆r ∼ ~/Λ, from now on, we adopt natural units ~ = c = 1 for con- κ ∼ Λ2c/~ where Λ ≈ 200 MeV/c is the non-perturbative venience. QCD scale parameter. This yields the condition For the effect of viscosity in gravitational collapse we follow Herrera et al. [17]. Adopting, in part, the notation ΛcR Λ2c ∆F ∼ ~ > , (2) of [17] we use the following form of the metric (note the r3 Λ ~ opposite overall sign): which is satisfied for ds2 = A2dt2 − B2dr2 − (Cr)2(dθ2 + sin2 θdφ2). (5)  1/3 ~ ΛR We also need the shear stress tensor π , which we write r < rh ∼ . (3) µν Λ ~ in the form For a 10 solar mass black hole, which has a Schwarzschild  1  6 πµν = Ω χµχν − hµν , (6) radius R = 30 km, we obtain rh ∼ 3 × 10 ~/Λ ≈ 3 nm. 3 Matter can still be described in this domain by stan- µ dard techniques of thermal quantum field theory in where χµ is a unit vector, χµχ = −1, along the radial curved space-time, because the hadronic disruption dis- direction and hµν is the projector on the hypersurface tance is much larger than the radial distance at which orthogonal to the four-velocity. The nonvanishing com- the space-time curvature reaches the Planck scale. This ponents of the shear strain are (see eq. (9) in ref. [17]): 3 −2 condition is satisfied as long as R/r  rP , where −35 2 1 rP ≈ 10 m is the Planck length. Inserting the ex- 1 2 3 σ1 = σ, σ2 = σ3 = − σ, (7) pression for the Schwarzschild radius R we obtain 3 3

1/3 where r  rc = rp(M/MP ) , (4) 1 d B −22 σ = ln (8) with rc ≈ 10 m for a 10 solar mass black hole. A dt C 3 is the shear scalar. It is convenient to consider ρ as the comoving radial coor- The shear stress tensor can be expressed in terms of a dinate. We now transform the metric from the Painlev´e- single function Ω(r, t) as: Gullstrand form with coordinates (r, t) to comoving co- ordinates (ρ, t). Using the differential 2 1 π1 = Ω, π2 = π3 = − Ω. (9) 1 3 2 3 3 rρ dr = dρ + ψdt, (16) r The relaxation equation of the shear stress (eq. (50) in ref. [17]) simplifies in the limit considered here. We go we easily obtain beyond the simplest form of the Israel-Stewart theory considered in [17] and include the complete set of second- ρ r2 ds2 = dt2 − dρ2 − ρ2(dθ2 + sin2 θdφ2). (17) order terms obtained in the so-called 14 moment approxi- r ρ2 mation to kinetic theory [18] (see also [19]). The equation then reads Comparing with (5) we read off the following metric func- tions: dΩ 4 10  τπ + ΩΘ + Ωσ = −2ηAσ − ΩA, (10) dt 3 21 A = 1,B = pρ/r, C = r/ρ; (18) where η is the shear viscosity, τπ is the thermal relaxation with r(ρ, t) given by (15). We can now calculate the shear time of the shear stress, and scalar (8): √ 1 d 2 d 3 R Θ = ln(BC ) (11) σ = − (ln r3/2) = . (19) A dt dt 2r3/2 is the volume expansion scalar. In the comoving coordi- It may be tempting to simply adopt the N-S approxi- nates we will be using, A = 1, and the equations sim- mation, where the shear stress is given by (13). However, plify further. In the Navier-Stokes (N-S) limit, τπ → 0, this limit must fail at some point during the collapse, be- eq. (10) implies cause the shear stress in the fluid cannot instantaneously adjust to the ever more rapidly growing gravitational Ω = −2ησ (12) shear strain. We thus need to solve the relaxation equa- and thus tion for the shear stress (10) with A = 1. Since the relaxation time τπ is a function of the tem- perature T , the relaxation equation for the shear stress πµν = −2ησµν . (13) cannot be solved in isolation, but needs to be solved to- 2 2 µν gether with the equation governing the growth of the The amount of shear 3 σ = σµν σ and volume con- traction Θ in the collapsing fluid depend greatly on the entropy density. The equation for entropy production in details of the collapse. In the case of the collapse of a viscous fluid can be found, e.g., in ref. [15], eq. (28): a self-gravitating homogeneous dust cloud, there is no shear (σ = 0) but rapid contraction governed by a con- µ 1 √ µ µν 2 T s;µ = T √ ∂µ( gs ) = π σµν = − σΩ, (20) tracting FWG geometry. Joshi et al. [20] give an explicit g 3 expression for the case of dust (p = 0), but the general solution is more complicated. For a thin shell collaps- where as usual the semicolon denotes the covariant ing into a Schwarzschild black hole, there is both strong derivative and shear and volume contraction (Θ = −σ). Other collapse g = − det g = A2B2(Cρ)4 sin2 θ = ρr3 sin2 θ (21) scenarios probably lie in between these two extremes. µν Here we adopt the following simplified approach: We is the metric determinant. Since σ and Ω have opposite consider a pre-existing black hole into which a thin signs this implies entropy growth as dictated by the sec- shell of matter falls, which perturbs the metric only ond law of thermodynamics. In the comoving frame the marginally, so that we can assume the shell to fall in spatial entropy current ~s vanishes, and we have s0 = s. a Schwarzschild background with Schwarzschild radius The left-hand side of the entropy growth equation thus R = 2GM. The Schwarzschild metric in Painlev´e- takes on the form Gullstrand form (see eq. (50) of ref. [16]) is given by: µ T d  3/2  2 2 2 2 2 2 2 T s = r s . (22) ds = dt − (ψdt − dr) − r (dθ + sin θdφ ) (14) ;µ r3/2 dt with ψ = −pR/r. We assume that the shell falls along For an ultrarelativistic gas in the high temperature limit, 3 the geodesic of a free particle: we have s ∝ T and after a few manipulations we find: √ √ s d √ 3/2 3 3 3/2 µ
r + R t = R t ≡ ρ . (15) T s;µ = 3√ rT . (23) 2 2 0 r dt 4

Equations (10) and (23) yield two equations governing The N-S approximation is applicable when the viscous the evolution of the temperature and the shear stress. It correction to the stress tensor is small compared with the is useful to rewrite them in terms of the entropy density- thermal stress. Using (9) and the ultrarelativistic form 3 normalized shear stress ω = Ω/s, which has dimension of of the thermal equilibrium entropy, sth = aT , the radial energy: pressure assumes the form:   r dω d ln s 4 10 η 1 2 1 η R τπ + ω + ωΘ + ωσ = −2 σ − ω (24) P = T s + ωs = T s − 2 s. (30) dt dt 3 21 s r 4 th 3 4 th s r3 d √
2√ rT = − rσω. (25) Like other thermodynamic variables, the entropy factor s dt 9 in the second term contains thermal equilibrium and vis- The advantage of this rewriting is that it removes a large cous contributions. Since the second term already rep- part of the temperature dependence from the shear vis- resents the viscous contribution to the radial pressure, cosity, in some limits, all of it. Equations (24, 25) can be which is assumed to be small in the N-S limit, it is con- solved when the time evolution of the shear strain σ and sistent to approximate the entropy factor s in the second the ratio η/s are known. For the collapsing shell in the term by the leading contribution, the thermal entropy Schwarzschild geometry, σ is given by (19), and we will sth. This yields assume that η/s is a constant. r ! There are then two different limits of (10) that need to 1 η R P = s T − 8 . (31) be considered. One is the domain, at large r, where the r 4 th s r3 shear strain σ is small and slowly varying, and therefore the left-hand side of (10) can be neglected. This is the The condition for the validity of the N-S approximation domain of validity of the Navier-Stokes approximation. thus becomes The other is the domain where σ is large and varying r rapidly (though still on time scales large compared to η R T  8 . (32) τπ), where the second term on the right-hand side of (10) s r3 can be neglected. We shall denote the radius at which the transition between these two domain occurs by rm. According to (29), and neglecting the viscous contribu- Finally, as we shall see below, there is a third domain at tion, this condition is satisfied as long as the smallest radii where the rate of change of tidal forces exceeds the inverse hydrodynamic response time and the η r R 1 collapsing matter dissolves into a collection of nearly free- r  rm ≈ 8 . (33) s r0 T0 streaming quanta. We shall denote the radius where this transition occurs by rg. It is easy to confirm that in this region the viscous cor- rection to the temperature is, indeed, small.

IV. NAVIER-STOKES LIMIT V. SECOND-ORDER VISCOUS We first consider the Navier-Stokes (N-S) limit (r > HYDRODYNAMICS rm): Since we are interested in the behavior approaching the Ω η singularity, we now consider the opposite case where the ω = = −2 σ, (26) s s growth rate of the shear strain σ is too fast for the shear stress to adjust instantaneously. In this domain (r < rm) and (25) takes the form it is possible to neglect the second term on the right-hand side of (24). Again using dt = −drpr/R at fixed comov- d √ 4 η √ η R rT
= rσ2 = . (27) ing position ρ, introducing the dimensionless shear stress dt 9 s s r5/2 ω¯ = ω/T , and inserting the Schwarzschild expressions for σ and Θ, the shear stress relaxation equation becomes: Using dt = −drpr/R at fixed comoving position ρ, we obtain dω¯ d ln(sT ) 9 ω¯  3 η τ +ω ¯ + = . (34) √ π dr dr 7 r r sT d √ η R rT
= − , (28) dr s r2 In order to make progress, we need to fix the temperature dependence of the relaxation time τ of the shear stress. which has the solution π Here we adopt the parametrization r rr η R 1 1  T (r) = T 0 + − . (29) η 0 τπ = β . (35) r s r r r0 sT 5

The proportionality constant has the value β = 5 in the VI. SOLUTION MATCHING conformal Boltzmann gas limit [18], and β = 2(2 − ln 2) in the strongly coupled super-Yang-Mills theory [21, 22]. As discussed above, the N-S approximation applies for 3 Specializing to the conformal limit s ∝ T , the equations large values of r where the gravitational shear strain is forω ¯ and T become: small. The limiting solution of the second-order theory that we derived in the previous section applies for small dω¯ 5 ω¯ 4 ω¯2 3 − + = , (36) values of r where the shear strain is large and growing dr 7 r 3 r βr rapidly. In order to match the two solutions we choose d ln T 1 ω¯ the radius r where the pressure anisotropy ξ coincides + = , (37) m dr 2r 3r for the two solutions. Using (9), ξ is in both cases given by the expression where we used the second equation to eliminate 8ω(r) d(ln T )/dr from the first equation. Remarkably, the two 1 1 2 T sT + Ω 1 + 3T (r) equations no longer depend on the value of η/s nor on ξ(r) ≡ 1 = 4 3 = . (45) T 2 1 sT − 1 Ω 4ω(r) the value R of the Schwarzschild radius. 2 4 3 1 − 3T (r) Equation (36) is solved by a constantω ¯ that satisfies the condition Using the expression (26) together with (19) and the ideal fluid limit of the temperature (29) we obtain for the N-S 4 5 3 solution: ω¯2 − ω¯ − = 0. (38) 3 7 β q 1 − 8 η R 1 s r0 T0r The negative solution of this equation for β = 5 is ξNS(r) ≈ q . (46) 1 + 4 η R 1 s r0 T0r s 15 152 9 ω¯ = − + ≈ −0.454; (39) For the causal viscous hydrodynamics solution with 56 56 20 ω(r) =ω ¯ T (r) the anisotropy is for the strongly coupled SUSY theory we find 1 + 8 ω¯ ξ (r) = 3 , (47) IS 1 − 4 ω¯ s 3 15 152 9 ω¯ = − + ≈ −0.698. (40) which is time independent. The two expressions coincide 56 56 8(2 − ln 2) when r We can now use (37) to solve for T (r): 3 η R 1 r = rm ≈ − . (48) ω¯ s r0 T0 r b T (r) = T 1 (41) 1 r It is reassuring that the coefficient −3/ω¯ in (48) is nu- merically of similar size as the coefficient in (33). with The parameters of the matched second-order solution are then related to those of the N-S solution as 1 ω¯ r b = − . (42) r0 2 3 T1 ≈ T0 . (49) rm The expression for ω(r) =ω ¯ T (r) grows more slowly for r → 0 than σ ∼ r−3/2, confirming our assumption that We still need to check whether the assumption underlying the second term on the right-hand side of (24) can be viscous hydrodynamics, that the viscous corrections to neglected in this domain. the stress tensor do not exceed the pressure, are satisfied. The solution (41) of second-order viscous hydrodynam- Using the results just obtained, we find for the shear ics assumes that the time scale of change of the gravita- stress scalar: tional shear strain σ is longer than the hydrodynamic 2 2 π1 = Ω = ωsT,¯ (50) response time τπ. The rate of change of σ is given by 1 3 3

σ 2 r 2  r 3/2 while the thermal equilibrium pressure is given by Peq = τg = = − = R . (43) 1 dσ/dt 3 dr/dt 3 R 4 sT . For our two limiting cases of the relaxation time coefficient, we have: The constraint τπ < τg is satisfied when r > rg with 2 ω¯ ≈ −0.303 [β ≈ 5], (51) √ 3 3 3β η R 2 −b 2 rg = b . (44) ω¯ = −0.465 [β = 2(2 − ln 2)], (52) 2 s T1r1 3 6 which means that the positivity of the radial pressure is Inserting this result into the previous equation, we obtain violated in both cases. This result shows that the validity 3   2 −b of the formalism of second-order viscous hydrodynamics rg 3β is, at best, marginal, and that higher-order corrections ≈ , (55) rm 16 should be taken into account. However, it is worth noting that the second-order theory represents a vast improve- which is slightly less than unity for the two values of β ment over the N-S approximation, which fails completely considered here (β = 5 and β = 2(2 − ln 2) ≈ 2.61). upon approach to the singularity as we saw in the previ- It is instructive to consider some semi-realistic exam- ous section. ples. We choose the following parameters: R = 30 km, In Figure 1 we show the complete solution for T (r) for r0/R = 1, β = 2(2 − ln 2), η/s = 1/4π, and consider the parameters r0/R = 1, β = 2(2−ln 2), and T0 = 20/R. the following two initial temperatures: T0 ≈ 0.25 meV The figure shows the temperature in viscous hydrody- (temperature of the cosmic background radiation) and namics (solid red line) in comparison with that found in T0 = 100 MeV (10 percent of the proton rest mass), ideal hydrodynamics (dashed blue line). The arrows in- assuming substantial pre-heating of the matter during dicates the matching radius rm and the radius rg where its accretion into the vicinity of the event horizon. The the hydrodynamic description fails. Figure 1 shows that matching radius rm and temperature Tm = T (rm) are the effect of viscosity on the temperature in the domain then found as shown in Table I. Both the radial dis- of validity of the N-S approximation is small. T0 rm Tm 0.25 meV 0.5 mm 1.1 eV 100 MeV 1.3 fm 2.8 × 108 GeV

TABLE I: Radial distance rm and temperature Tm at the distance where the Navier-Stokes approximation fails, because the radial pressure of the collapsing matter becomes negative. The hydrodynamic model fails shortly afterwards.

tances and the temperature fall well into the domain where quantum gravity effects can be ignored. The question remains what happens at distances r < rg . rm where the hydrodynamical description fails. One possibility is that the hydrodynamical model needs to be replaced by a kinetic description of interacting quanta in phase space for r < rg that can account for the fact that the relevant response times of the medium are shorter FIG. 1: The temperature function T (r) (in units of R−1) for than the hydrodynamical response times. Another possi- the selected parameters. The arrows shows the matching ra- bility is that the continuous fluid model fails as the radial dius rm and the radius rg where viscous hydrodynamics fails. pressure becomes negative for r < rm and medium me- The solid (red) line shows the solution for viscous hydrody- chanically ruptures (cavitates) in the radial direction into namics; the dashed (blue) line shows the solution for ideal increasingly smaller domains. Yet another possibility hydrodynamics. is that the increasingly anisotropic quark-gluon plasma spontaneously develops growing domains of partially co- herent color fields that modify the transport equations Quite generally, one finds that rg ≈ rm implying that second-order hydrodynamics only slightly extends the ex- [23, 24]. Possibly more than one of these scenarios is tent of validity of hydrodynamics beyond the domain of relevant. validity of the Navier-Stokes approximation. This can be seen as follows: Setting r1 = rm in (44) we can write the equation as VII. CONCLUSION

3 −b √ We have explored the fate of matter falling into a  r  2 3β η R g = , (53) macroscopic Schwarzschild black hole for the simplified 3/2 rm 2 s T1rm case of a radially collapsing thin spherical shell for which the back reaction of the geometry can be neglected. We where T1 = T (rm). On the other hand, combining (33) have treated the internal dynamics of the infalling matter and (49) yields in the framework of viscous relativistic hydrodynamics and calculated how the internal temperature of the col- η √ lapsing matter evolves as it falls toward the Schwarzschild T r3/2 ≈ 8 R. (54) 1 m s singularity. We found that viscous hydrodynamics fails 7 when either, the dissipative radial pressure exceeds the In memoriam thermal pressure and the total radial pressure becomes negative, or the time scale of variation of the tidal forces This article is dedicated to the memory of Walter acting on the collapsing matter becomes shorter than the Greiner, our teacher, advisor and mentor during early characteristic hydrodynamic response time. We denoted parts of our careers. Walter was an inspiring scientist and the first radius by rm, and the second radius by rg, and teacher, who taught all that came into close contact with found that rg . rm. him that the curiosity of the scientific mind should know We found that the temperature T of a shell of matter no bounds. In this, he led by example. It was Walter’s collapsing into a black hole grows as a power of the radius mantra that scientific inquiry should not impose limits parameter on its approach toward the quantum grav- on itself by sticking to generally accepted approaches to ity domain that shields the classical singularity. Since a given problem or established ways of thinking about it. the temperature of the matter remains far below the For much of his life he was fascinated by the question Planck scale, it can be described by conventional meth- how the vacuum – empty space – responds to external ods of thermal quantum field theory. We found, however, forces, such as strong fields. In his later years as a scien- that the pressure anisotropy becomes large in the small- tist, Walter was particularly curious about black holes, r region where the Navier-Stokes approximation fails. i.e. strong gravitational fields. He and his collaborators As even the validity of causal second-order viscous hy- explored theories that would avoid the existence of black drodynamics in this domain is questionable, more work holes [29, 30]. While we do not follow this idea here, we will need to be done before definitive conclusions can be recognize that the question, what kind of objects black reached. An approach that treats the anisotropy of the holes exactly are, still remains an area of inquiry. stress tensor in all orders, such as anisotropic hydrody- In our article we address a very specific question: how namics [25, 26] or a kinetic approach [27], may be better matter falling into a Schwarzschild black hole is viscously suited to describe this region. heated by the tidal forces it encounters. In doing so, Another interesting approach, especially for the treat- we connect with two of Walter’s lifelong interests: black ment of a collapsing, strongly coupled quark-gluon holes and relativistic hydrodynamics of matter under ex- plasma, would be a holographic description of the col- treme conditions. We believe that in doing so the article lapse. In such an approach the space-time metric (17) would, at the very least, have captured his interest. of the black hole represents the boundary geometry of a Acknowledgments: This work was supported in part by five-dimensional space-time with a black brane, and the grants from the U.S. Department of Energy (DE-FG02- hydrodynamical evolution is recovered from the infrared 05ER41367). AS thanks the Nuclear and Particle Physics near-boundary fluctuations around the black-brane back- Directorate of Brookhaven National Laboratory for their ground [28]. Many other future research directions come hospitality. We thank Ulrich Heinz for helpful comments to mind: One would be the numerical solution of the full on the correct form of relativistic second-order viscous set of Einstein equations combined with the relaxation hydrodynamics and Bill Zajc for several useful comments equation for the dissipative stress tensor, including the on an earlier version of our manuscript. We also thank back reaction of the metric to the collapsing matter. By the organizers of the FIAS International Symposium on modifying the equation for the gravitational field in such Discoveries at the Frontiers of Science for the invitation a way that a singularity is avoided, it would be possible to present a lecture. to describe the entropy generation during a gravitational bounce.

[1] J. R. Oppenheimer and H. Snyder, Phys. Rev. 56, 455 [9] K. Skenderis and M. Taylor, Phys. Rept. 467, 117 (2008) (1939). [arXiv:0804.0552 [hep-th]]. [2] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [10] A. Almheiri, D. Marolf, J. Polchinski and J. Sully, JHEP Erratum: [Commun. Math. Phys. 46, 206 (1976)]. 1302, 062 (2013) [arXiv:1207.3123 [hep-th]]. [3] G. ’t Hooft, Nucl. Phys. Proc. Suppl. 43, 1 (1995). [11] P. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. [4] L. Modesto, Class. Quant. Grav. 23, 5587 (2006) [gr- Lett. 94, 111601 (2005) [hep-th/0405231]. qc/0509078]. [12] M. Brigante, H. Liu, R. C. Myers, S. Shenker [5] L. Modesto and I. Premont-Schwarz, Phys. Rev. D 80, and S. Yaida, Phys. Rev. D 77, 126006 (2008) 064041 (2009) [arXiv:0905.3170 [hep-th]]. [arXiv:0712.0805 [hep-th]]. [6] S. Hossenfelder, L. Modesto and I. Premont-Schwarz, [13] R. Brustein and A. J. M. Medved, Phys. Rev. D 79, Phys. Rev. D 81, 044036 (2010) [arXiv:0912.1823 [gr- 021901 (2009) [arXiv:0808.3498 [hep-th]]. qc]]. [14] S. W. Hawking, A Brief History Of Time, (Bantam [7] A. H. Chamseddine and V. Mukhanov, JCAP 1703, 009 Books, New York, 1988). (2017) [arXiv:1612.05860 [gr-qc]]. [15] U. W. Heinz, nucl-th/0512049. [8] N. Bodendorfer, A. Sch¨afer and J. Schliemann, [16] R. J. Adler, J. D. Bjorken, P. Chen and J. S. Liu, Am. arXiv:1703.10670 [gr-qc]. J. Phys. 73, 1148 (2005) [gr-qc/0502040]. 8

[17] L. Herrera, A. Di Prisco, E. Fuenmayor and O. Troconis, [arXiv:0712.2916 [hep-th]]. Int. J. Mod. Phys. D 18, 129 (2009) [arXiv:0804.3584 [23] M. Asakawa, S. A. Bass and B. M¨uller, Phys. Rev. Lett. [gr-qc]]. 96, 252301 (2006) [hep-ph/0603092]. [18] G. S. Denicol, H. Niemi, E. Molnar and [24] M. Asakawa, S. A. Bass and B. M¨uller,Prog. Theor. D. H. Rischke, Phys. Rev. D 85, 114047 (2012) Phys. 116, 725 (2007) doi:10.1143/PTP.116.725 [hep- Erratum: [Phys. Rev. D 91, no. 3, 039902 (2015)] ph/0608270]. doi:10.1103/PhysRevD.85.114047, 10.1103/Phys- [25] W. Florkowski and R. Ryblewski, Phys. Rev. C 83, RevD.91.039902 [arXiv:1202.4551 [nucl-th]]. 034907 (2011) [arXiv:1007.0130 [nucl-th]]. [19] G. S. Denicol, U. W. Heinz, M. Martinez, J. Noronha and [26] D. Bazow, U. W. Heinz and M. Strickland, Phys. Rev. C M. Strickland, Phys. Rev. D 90, no. 12, 125026 (2014) 90, 054910 (2014) [arXiv:1311.6720 [nucl-th]]. doi:10.1103/PhysRevD.90.125026 [arXiv:1408.7048 [hep- [27] W. Florkowski, M. P. Heller and M. Spalinski, ph]]. arXiv:1707.02282 [hep-ph]. [20] P. S. Joshi, N. Dadhich and R. Maartens, Phys. Rev. D [28] P. Kovtun, D. T. Son and A. O. Starinets, JHEP 0310, 65, 101501 (2002) [gr-qc/0109051]. 064 (2003) [hep-th/0309213]. [21] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets [29] P. O. Hess, L. Maghlaoui and W. Greiner, Int. J. Mod. and M. A. Stephanov, JHEP 0804, 100 (2008) Phys. D 19, 1217 (2010). [arXiv:0712.2451 [hep-th]]. [30] W. Greiner, PoS FFP 14, 036 (2016). [22] M. Natsuume and T. Okamura, Phys. Rev. D 77, 066014 (2008) Erratum: [Phys. Rev. D 78, 089902 (2008)]