RobinsonTrautman spacetimes and gauge/gravity duality
Abstract
We study farfromequilibrium field theory dynamics using gauge/gravity duality applied to the RobinsonTrautman (RT) class of spacetimes and we present a number of new results. First, we assess the applicability of the hydrodynamic approximation to inhomogeneous plasma dynamics dual to RT spacetimes. We prove that to any order in a late time expansion it is possible to identify variables corresponding to the local energy density and fluid velocity. However, we show using numerical examples that this does not hold at the nonperturbative level; for sufficiently inhomogeneous initial data a local rest frame does not exist. Second, we preset a new class of holographic inhomogeneous plasma flows on the plane. The corresponding spacetimes are not of the RT type but they can be obtained from RT spacetimes with spatially compact boundaries by coordinate transformations which generate Poincaré patchlike coordinates with planar boundaries. We demonstrate the application of this procedure using numerical examples.
RobinsonTrautman spacetimes and gauge/gravity duality
\FullConferenceCorfu Summer Institute 2016 ”School and Workshops on Elementary Particle Physics and Gravity”
31 August  23 September, 2016
Corfu, Greece
1 Introduction
The description of quantum systems far from equilibrium remains a challenging exercise. The machinery of gauge/gravity duality has become an attractive tool in this arena, which is arguably at its most mature in addressing CFTs with gauge group where is large. In this context nonequilibrium dynamics in the CFT are encoded by the a theory of classical gravity in one higher dimension with a negative cosmological constant, . Often one also includes a nonextremal black hole which provides a thermal contribution.
There are by now many examples of the application of holography to nonequilibrium physics, with much attention being payed to sudden quenches of couplings or field theory sources, and the subsequent relaxation back to equilibrium (see for example [1, 2, 3, 4]). When there is no null infinity to which radiation can escape, and so – in the absence of external dissipative couplings – the spacetime settles down entirely due to the dissipation provided by the black hole horizon. These late time dynamics are in part described by quasinormal modes of the black hole. At sufficiently late times an effective theory of the slowest () modes of the black hole persists, which takes the simple form of relativistic hydrodynamics, and equilibrium is ultimately reached. Other examples involve the construction of states which never return to equilibrium, such as the recent interest in the Riemann problem for CFTs considered both hydrodynamics and holography [5, 6, 7, 8, 9, 10], and the construction of stationary inhomogeneous plasma flows such as [11]. In such cases the chief interest is the existence of steadystate regions which are also out of equilibrium.
Whilst the use of gauge/gravity duality provides a considerable computational advantage over a direct field theory computation, one is still faced with solving the Einstein equations. For generic initial data one may utilise the tools of numerical relativity. In these proceedings we consider instead a specialised class of metrics: the RobinsonTrautman (RT) spacetimes [12]. These are solutions to the Einstein equations for any , but we focus exclusively on the case . They are inhomogeneous and time dependent, and are governed by a volumepreserving geometric flow equation for the spatial metric at the boundary of spacetime, the Calabi flow. This geometric flow is a fourthorder nonlinear diffusion process in dimensions, and its relative simplicity allows for valuable analytic control in the absence of any bulk isometries.
Summary of previous results. The study of the holographic aspects of RT was initiated in [13, 14]. The RT metrics can be viewed as the nonlinear extension of algebraically special perturbations of Schwarzschild black holes. The algebraically special modes are purely dissipative, and when describe radiation escaping to null infinity. When the modes remain purely outgoing, and describe radiation which deforms the conformal boundary metric. As such, the Calabi flow equation determines the evolution of both the boundary metric and the interior of the spacetime, including expectation values of the stress tensor – explicit expressions for these quantities in terms of the Calabiflow variables are given in [13, 14].
The dissipation in RT is due to the outgoing radiation rather than absorption by the black hole. An interpretation of this is the dissipation that arises due to external couplings. As such the physics is somewhat different to the usual studies of out of equilibrium physics in holography. It is therefore important to analyse how RT behaves at late times, since its late time behaviour will no longer be governed by quasinormal modes of a black hole. In cases where the boundary metric is topologically , an analysis was made in [14], finding exponentially damped modes with the Schwarzschild solution reached at late times. This has its origins in the late time properties of the Calabiflow on [15, 16, 17, 18, 19, 20, 21, 22].
For noncompact boundaries, a study was initiated in [23] for the case . In general, at late times Schwarzschild is not reached, and instead the late time behaviour is governed by families of selfsimilar attractor solutions. For instance, depending on the boundary conditions chosen at infinity, a class of solutions may be approached that depends only on . These solutions relate to the Riemann problem in the context of the field theory duals to RT spacetimes. The late time attractor solutions are governed by an underlying Lifshitz symmetry (which is subsequently broken at finite ), and the approach of a general dynamical evolution towards them is powerlaw in time, with the powers determined by the spectrum of operators about this Lifshitz invariant fixed point.
Summary of new results in these proceedings. Sufficiently close to equilibrium, the theory of relativistic conformal hydrodynamics should provide a good approximation of the dynamics of the field theory. In particular, the expectation is that hydrodynamics should apply when all gradients are small compared with the temperature scale. At least, this regime will occur at sufficiently late times in unforced situations described through holography via the Einstein equations, since it is expected that a final stationary black hole will be reached, corresponding to equilibration. As discussed, this is the case for the RT spacetimes in the case, which approach Schwarzschild at late times.
However, for a fully nonlinear evolution of the Einstein equations, it is useful to assess the validity of the hydrodynamic approximation. Indeed there have been some surprises, with hydrodynamics applicable in regimes where it would naively be expected not to apply (for a recent example see [24]). On the other hand, nonlinear dynamics may also produce regions with no local rest frame and local energy density [25]. The RT spacetimes provide a convenient way to test these issues, given the relative ease of obtaining a genuinely nonlinear inhomogeneous time dependent bulk solution. We will discuss this in section 3.
One goal of out of equilibrium holography is to understand the dynamics of a plasma on the plane. In order to study plasma dynamics on the plane, one needs an evolving bulk spacetime where the conformal boundary is planar. One efficient way to construct such solutions is to take known solutions in globalAdS coordinates where the boundary is the Einstein static universe, and perform a coordinate transformation such that the boundary is Minkowski space. Such a chart covers only part of the original spacetime, and corresponds to a conformal transformation from the boundary point of view. For instance the case of SchwarzschildAdS perturbed by quasinormal modes was considered in [26].
We would like nontrivial plasma dynamics encoded by the RT solutions with boundary . One way to get this is to consider Calabiflow on as considered in [23]. Another approach is as described above, to consider an RT solution with conformal boundary of topology supplemented by a coordinate transformation which gives a bulk solution boundary with topology . The new coordinates can be thought of as Poincaré patchlike coordinates, but note that the new boundary metric of topology is deformed by the dynamics associated to the outgoing radiation, just as in the case. Solutions so obtained describe nonlinear plasma dynamics on the plane which do not themselves correspond to a planar RT evolution. The details of the transformations and some numerical examples are given in section 4.
We begin with a brief review of the RT solution in section 2.
Note. The work described here is part of a bigger program that was initiated by Ioannis Bakas, aiming to explore the connection between the RobinsonTrautmann spacetimes, Calabi flow and holography. Unfortunately, Ioannis passed away before we could complete this program. Ioannis was an outstanding physicist and a very good friend and he will be missed by the community. We would like to acknowledge his many contributions and we hope to be able to finish the work we started together.
2 Brief review of the RT metrics and Calabiflow
The RT metric can be written,
(1)  
(2)  
(3) 
is the Ricci scalar of the metric , is a function of and the two , denotes the line element for an Euclidean 2metric with constant scalar curvature so that is the plane or torus, is the unit round sphere and hyperbolic space. Additionally the solution has a constant mass parameter . The RT metric (1)(3) solves the Einstein equations with cosmological constant provided obeys the Calabiflow equation,
(4) 
which, after inserting (2) can be seen to be a nonlinear fourthorder diffusion equation for . Here and are the Laplacian and Ricci scalar for .
In the case the conformal boundary metric – reached as , where the outgoing coordinate becomes boundary time coordinate – is given by
(5) 
where . Explicit expressions for the expectation value of the dual field theory stress tensor can be found in [14] depending on up to six spatial derivatives of .
3 Validity of the hydrodynamic approximation
In this section we assess the validity of the hydrodynamic approximation applied to evolutions described by RT spacetimes. In particular, in order for a hydrodynamic description to hold, we must first identify the local fluid variables: an energy density and a 4velocity field . To do so we look for solutions of the eigenvalue problem,
(6) 
where is the expectation value of the field theory stress tensor, which is provided here by holography from the near boundary behaviour of RT solutions. In domains where (6) yields physical solutions for , we may compare with the expected result from conformal relativistic hydrodynamics, whose constitutive relations are
(7) 
where represents higher order derivative corrections to the ideal hydrodynamic behaviour, and where we work in the Landau frame . Note also that , which is due to the conformal invariance of the dual QFT.
In section 3.1 we prove that it is possible to obtain local hydrodynamic variables to all orders in a perturbative late time expansion. However, interestingly this does not continue to hold at the nonperturbative level, which we demonstrate with explicit numerical solutions in section 3.2. Similar results were reported in [25]. We also use the numerical solutions to demonstrate the expected convergence to hydrodynamics at late times as governed by the late time expansion.
For simplicity we focus the Calabi flow on with axial symmetry. Thus we are constructing RT spacetimes with a conformal boundary of the form . The solutions will be functions of boundary time and the polar angle .
3.1 Proof of solution to eigenvalue problem (6) to any order in a late time expansion
Details of the late time expansion for the axially symmetric case were given in [14]. Here we utilise these results to prove that the eigenvalue problem (6) can be solved in a late time expansion to any order. As in [14] we adopt a normalisation where in this section.
First we note that for the late time expansion of ,
(8) 
we can obtain the boundary metric and holographic stress tensor in that expansion. The index in parentheses indicates the order in the late time expansion, i.e. , and explicit expressions for the first few coefficients can be found in [14].
Our goal here is to show that we can always find hydrodynamic variables to solve
(9) 
order by order, thereby solving (6). We begin by expanding the hydrodynamic variables in a late time expansion, first the velocity,
(10) 
and similarly for the energy density, . Similarly, we expand the dissipative contribution to the stress tensor in this way,
(11) 
We begin by considering explicitly the solution to a few orders in the expansion. By explicit computation we have [14]
(12) 
then .
We now solve all first order quantities. Knowing (12), allows us to iteratively construct given the Landau frame condition, . At first order we have
(13) 
Using this result and the fact that vanishes [14] we then find that
(14) 
Then from the and components of (9) we obtain
(15) 
We now use to obtain
(16) 
which upon use of (14) implies . The component of (9) yields
(17) 
and the component of (9) yields
(18) 
So we managed to obtain all quantities, , , , to this order.
The general strategy adopted to solve the eigenvalue problem to first order can be extended to give a proof to all orders, by induction. We have all quantities of interest at order 0, i.e. , and . Next suppose that we have determined everything at order and below, i.e. , , and lower order terms. The following steps can be used to advance to order :

Solve the Landau frame condition for .

From the component of (9) solve for in terms of the unknown .

From the component of (9) solve for in terms of the unknown .

From the component of (9) solve for in terms of the unknown .

From the condition solve for .

From the component of (9) solve for .

From the component of (9) solve for .
With the above results we now have all quantities at order , i.e. , and . This completes the proof by induction.
Since this is a perturbative construction, will stay positive within perturbation theory, and we would need to resum the series in order to see any possibility of regions where .
3.2 Numerical examples
In this section we use numerical solutions of the Calabiflow to assess the hydrodynamic approximation, both for the existence of solutions to the (6) away from the late time regime (considered separately in section 3.1), and also to illustrate the approach to the late time regime itself. Note that existence of solutions to the (6) can be immediately assessed without performing any evolution, simply by looking at initial data at , which is sufficient to obtain , where we write for convenience. We pick a class of initial data which characterise deviations from equilibrium using axiallysymmetric () spherical harmonics,
(19) 
where is a constant. We have omitted the redundant azimuthal angle argument. We also restrict attention to even in order to have an antipodal symmetry. For we find that there is a critical value , such that for all spatial points admit real solutions to (6) whilst for there are points which do not. Similarly for the case we find a critical value .
Let us now turn to the time evolution for the case and . The angular profile of , where defined, is displayed in figure 1 for a handful of different time slices. Initially there are two ‘belts’ on the sphere where there is no real solution to (6). The fluid flows towards the poles in each hemisphere, with zero velocity on the equator and at the pole itself.
We would like to quantify the approach to equilibrium at late times in terms of the rest frame variables we have identified. To this end, once we have , we can construct two 3vectors where such that,
(20) 
and use them to define two pressures,
(21) 
Additionally we may construct a strain but this vanishes for the axially symmetric case of interest here. Here we take . We can identify from the numerical solutions, and we can obtain the same quantities analytically in a late time expansion by following the methodology outlined in section 3.1. This allows us to compare the numerical evolution of rest frame variables with the late time expansion. In the case of antipodal symmetry, the first 14 orders in the late time expansion are determined by a single mode, the quadrupole (), which at linear order takes the form,
(22) 
where is the constant late time value of . New data, , can enter at order corresponding to the contribution of a linear hexadecapole () mode,
(23) 
which goes in hand with a term due to the overlap of the two mode expansions [14].
In figure 2 we show the numerical evolution of and a comparison to the linear quadrupole mode, as well as its nonlinear extension up to order 15 in the late time expansion. For this particular numerical solution we find and the late time values of the fluid quantities are as expected, i.e. , , , . Late time deviations are governed by the amplitude of the linear quadrupole mode, which we find to be . Beyond order in the late time expansion the agreement with the fluid variables is not visibly improved by adding higher order terms in this case, and we have not been able to reliably extract the hexadecapole amplitude at order 15 from our numerics.
We now make some observations on the timescales indicated by the rest frame variables displayed in figure 2. The first observation is that the first order late time expansion (the dashed lines), which describes a hydrodynamic regime [14], provides a good description of the evolution earlier than the time at which the pressures become equal. Such behaviour was observed in other conformal systems [27, 28]. To better characterise this behaviour we construct a ‘hydrodynamisation’ time, , here defined to be the time after which both pressures agree with the order1 late time expression to better than 10% for all , i.e., we seek the time after which
(24) 
satisfied by both where and the superscript in parentheses denotes the firstorder expression in the late time expansion. For the numerical solutions shown in figures 1 and 2 we find , indicated by a vertical grey line in figure 2. The second observation is that the energy density reaches equilibrium earlier than the time at which the pressures become equal. To quantify this statement further, we introduce two timescales: an energy equilibration time, , and an isotropisation time . We define to be the time after which differs from its eventual equilibrium value (here ) by less than 10% for all , i.e.
(25) 
Note that for conformal systems, such as the one we are studying, is the same as the timescale characterised by when is within 10% of its equilibrium value, (here ). Next we define to be the time after which the pressures differ from one another by less than 10% for all , i.e.
(26) 
For the evolution shown in figures 1 and 2 we find and , indicated by vertical grey lines in figure 2. We note that .
4 Transformations to the Poincaré patch
In this section we seek to generate nontrivial plasma dynamics on , which requires a nontrivial evolving planar bulk spacetime. One way to achieve this is to consider directly RT metrics corresponding to Calabi flows on as considered in [23]. Another approach which we adopt here, is to perform a coordinate transformation from RT solutions which have a boundary of topology . The new coordinates correspond to a Poincaré patchlike slicing of the bulk spacetime. Transformations of this type were previously considered in their action on Schwarzschild perturbed by quasinormal modes to generate inhomogeneous plasma flows on the plane [26].
4.1 Generalities
We begin by considering the near boundary metric in FeffermanGraham form,
(27) 
where . The stress tensor is given by [29],
(28) 
Expressions for the corresponding to RT solutions were computed in [14].
The Einstein universe and are conformally related and the conformal transformation can be implemented via a bulk diffeomorphism. This diffeomorphism has been determined in the near boundary region in [30]. Here we present a derivation following the approach set out in [26]. Consider AdS as an embedding in defined by the hyperboloid,
(29) 
The metric for global AdS in FeffermanGraham coordinates can then be reached by
(30)  
(31)  
(32)  
(33)  
(34) 
giving
(35) 
where and where , so that, for the terms presented in the FeffermanGraham expansion (with radial variable ),
(36)  
(37)  
(38) 
Note that the here differs from the one used in earlier sections by an overall normalisation by . Similarly, the metric for planar AdS with boundary polar coordinates can be reached by
(39)  
(40)  
(41)  
(42)  
(43) 
giving
(44) 
and again, for a FeffermanGraham expansion (with radial variable ), we have
(45)  
(46)  
(47) 
where the tilde notation has been introduced to indicate the planar case.
In order to perform the transformation of interest, we just need to find the conversion from to by equating the ’s:
(48)  
(49)  
(50)  
(51) 
This maps the north pole of the to the origin of , and at the south pole is mapped to a point at infinity. At different , the amount of the covered by the planar chart decreases, until at only the north pole is covered.
4.2 Applying to RT
For axiallysymmetric RT solutions the global metric in FeffermanGraham gauge takes the form,
(53)  
(54)  
(55) 
Applying the coordinate conversions (52) directly to the metric specified by these coefficients however leads to terms in the complete line element of the form,
(56) 
i.e. outside the FeffermanGraham gauge. The development of these terms can be compensated for orderbyorder in . Specifically, in addition to we apply the supplementary coordinate transformation,
(57)  
(58)  
(59) 
where,
(60)  
(61)  
(62)  
(63) 
then the cross terms vanish and we are in FeffermanGraham gauge to high enough order in the expansion to extract the stress tensor. The boundary metric is given by,
(64)  
(65)  
(66)  
(67)  
(68)  
(69) 
When this reduces to Minkowski in polar coordinates. Expanding near the origin of boundary polar coordinates () we see that the coordinate system is still regular.
The stress tensor is given in terms of the original metric components, . Again here we denote corresponding functions of the new leading order coordinates by .
The original stress tensor in the planar coordinate system is covariantly conserved and traceless, i.e.
(70) 
Here we verify that given the transformations above these conditions imply that the new planarcase stress tensor given by is covariantly conserved with metric and traceless. To use the relations derived from (70), note that
(71)  
(72) 
for some function . With these relations it is a straightforward to check that,
(73) 
As a further technical note, in order to numerically compute the stress tensor in the planar case at some arbitrary value of we need to know the numerical solution, at,
(74) 
Fortunately both are bounded, and we simply have to ensure that we have evaluated the solution in the domain and . In practise we compute the solution numerically on a rectangular grid of points covering this domain, then interpolate in order to extract the solution at the exact value of required.
4.3 Examples
4.3.1 Schwarzschild black hole
Here we consider the simple case of the Schwarzschild solution corresponding to the trivial Calabi flow solution, , where we may work analytically. In global coordinates we have the boundary metric
(75) 
and stress tensor
(76)  
(77)  
(78) 
with other entries vanishing. And so, in the planar coordinates we have the metric
(79) 
and stress tensor
(80)  
(81)  
(82)  
(83) 
This situation is the 3d analogue of a flow considered in [26]. The evolution of the stress tensor is shown in figure 3, illustrating how a shell of plasma collapses towards the origin when , and then recedes for – interpreted as the black hole passing through the Poincaré patch in the bulk.
4.3.2 RT with an initial inhomogeneity
To illustrate the application to RT spacetimes we here consider the evolution of initial data of the form (19) with . In linear theory the addition of an perturbation to corresponds to the addition of a quadrupole radiation component. We have freedom to choose the value of at which we begin the RT evolution, which through the coordinate transformation determines how much of the plasma flow on the plane we are describing. In the cases considered here we consider initial data at and . A component of the holographic stress tensor on the plane, , is shown in figure 4 for .
5 Discussion
The RT spacetimes provide a convenient window into nonlinear, isometryfree solutions to the Einstein equations, including for . It is natural to try and utilise these solutions in the context of gauge/gravity duality, where they capture the dynamics of farfromequilibrium plasmas in certain field theories.
We first focussed on the applicability of the hydrodynamic approximation to the field theory duals to RT spacetimes, in the case where the conformal boundary is topologically . The boundary metric and stress tensor can be extracted from the near boundary behaviour of the solutions. We proved that local fluid variables – energy density and 4velocity fields – can be constructed to any order in a perturbative late time expansion. However, by analysing the nonperturbative regime using numerical solutions, we found that such variables cannot be identified everywhere for initial data which is sufficiently deformed from equilibrium. It would be interesting to try and resum the late time expansion to see these nonperturbative effects. At late times the solutions approach equilibrium at a rate governed by the longest lived mode identified in [14] corresponding to quadrupole radiation. We considered the nonlinear extension of this quadrupole mode to order 15 in a late time expansion, governed by a single amplitude , finding excellent agreement with the numerical data.
The regime described by linear quadrupole radiation is also a hydrodynamic regime [14]. The local fluid rest frame variables are well described by the linear quadrupole after a time, , which we quantified. We found where is the time at which the transverse and longitudinal pressures become equal, consistent with other observations of hydrodynamisation in conformal theories [27, 28]. We also constructed the quantity describing deviation of the energy density from equilibrium (or here equivalently, the deviation of average pressure from equilibrium), and found that .
We also generated a new class of inhomogeneous plasma flows on the plane by performing coordinate transformations which correspond to moving to region analogous to the Poincaré patch in AdS. These coordinate transformations correspond to covering part of the boundary of topology with a chart of topology . We also showed how this transformation could be applied efficiently to solutions obtained numerically in global coordinates. We gave examples with initial inhomogeneities and the corresponding evolution of the expectation value of the dual stress tensor. Such solutions provide a way of introducing nonlinear deviations from the inhomogeneous flows provided by conformal transformations of the Schwarzschild solution, which do not themselves correspond to a RT solution on the plane.
Acknowledgements
BW is supported by the NCCR under grant number 51NF40141869 “The Mathematics of Physics” (SwissMAP). KS is supported in part by the Science and Technology Facilities Council (Consolidated Grant “Exploring the Limits of the Standard Model and Beyond”). This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie SkodowskaCurie grant agreement No 690575.
References
 [1] P. M. Chesler and L. G. Yaffe, Horizon formation and farfromequilibrium isotropization in supersymmetric YangMills plasma, Phys. Rev. Lett. 102 (2009) 211601, [0812.2053].
 [2] K. Murata, S. Kinoshita and N. Tanahashi, Nonequilibrium Condensation Process in a Holographic Superconductor, JHEP 07 (2010) 050, [1005.0633].
 [3] V. Balasubramanian, A. Bernamonti, J. de Boer, N. Copland, B. Craps, E. KeskiVakkuri et al., Holographic Thermalization, Phys. Rev. D84 (2011) 026010, [1103.2683].
 [4] A. Buchel, R. C. Myers and A. van Niekerk, Universality of Abrupt Holographic Quenches, Phys. Rev. Lett. 111 (2013) 201602, [1307.4740].
 [5] D. Bernard and B. Doyon, Energy flow in nonequilibrium conformal field theory, J. Phys. A45 (2012) 362001, [1202.0239].
 [6] D. Bernard and B. Doyon, Timereversal symmetry and fluctuation relations in nonequilibrium quantum steady states, J. Phys. A46 (2013) K2001, [1306.3900].
 [7] H.C. Chang, A. Karch and A. Yarom, An ansatz for one dimensional steady state configurations, J. Stat. Mech. 1406 (2014) P06018, [1311.2590].
 [8] M. J. Bhaseen, B. Doyon, A. Lucas and K. Schalm, Far from equilibrium energy flow in quantum critical systems, Nature Physics 11 (2015) 509, [1311.3655].
 [9] I. Amado and A. Yarom, Black brane steady states, JHEP 10 (2015) 015, [1501.01627].
 [10] D. Bernard and B. Doyon, A hydrodynamic approach to nonequilibrium conformal field theories, 1507.07474.
 [11] P. Figueras and T. Wiseman, Stationary holographic plasma quenches and numerical methods for nonKilling horizons, Phys. Rev. Lett. 110 (2013) 171602, [1212.4498].
 [12] I. Robinson and A. Trautman, Some spherical gravitational waves in general relativity, Proc. Roy. Soc. Lond. A265 (1962) 463–473.
 [13] G. B. de Freitas and H. S. Reall, Algebraically special solutions in AdS/CFT, JHEP 06 (2014) 148, [1403.3537].
 [14] I. Bakas and K. Skenderis, Nonequilibrium dynamics and RobinsonTrautman, JHEP 08 (2014) 056, [1404.4824].
 [15] P. T. Chrusciel and D. B. Singleton, Nonsmoothness of event horizons of RobinsonTrautman black holes, Commun. Math. Phys. 147 (1992) 137–162.
 [16] B. Schmidt, Existence of solutions of the robinsontrautman equation and spatial infinity, Gen. Rel. Grav. 20 (1988) 65–70.
 [17] A. D. Rendall, Existence and asymptotic properties of global solutions of the robinsontrautman equation, Class. Quant. Grav. 5 (1988) 1339.
 [18] D. Singleton, On global existence and convergence of vacuum robinsontrautman solutions, Class. Quant. Grav. 7 (1990) 1333.
 [19] P. T. Chrusciel, Semiglobal existence and convergence of solutions of the RobinsonTrautman (twodimensional Calabi) equation, Commun. Math. Phys. 137 (1991) 289–313.
 [20] P. T. Chrusciel, On the global structure of RobinsonTrautman spacetimes, Proc. Roy. Soc. Lond. A436 (1992) 299–316.
 [21] X. X. Chen, Calabi flow in Riemann surfaces revisited: A new point of view, Intern. Math. Res. Notices No. 6 (2001) 275–297, [math/0009246].
 [22] M. Struwe, Curvature flows on surfaces, Ann. Scuola Norm. Sup. di Pisa CI. Sci. (5) I (2002) 247–74.
 [23] I. Bakas, K. Skenderis and B. Withers, Selfsimilar equilibration of strongly interacting systems from holography, Phys. Rev. D93 (2016) 101902, [1512.09151].
 [24] M. Attems, Y. Bea, J. CasalderreySolana, D. Mateos, M. Triana and M. Zilhao, Phase Transitions, Inhomogeneous Horizons and SecondOrder Hydrodynamics, 1703.02948.
 [25] P. Arnold, P. Romatschke and W. van der Schee, Absence of a local rest frame in far from equilibrium quantum matter, JHEP 10 (2014) 110, [1408.2518].
 [26] J. J. Friess, S. S. Gubser, G. Michalogiorgakis and S. S. Pufu, Expanding plasmas and quasinormal modes of antide Sitter black holes, JHEP 04 (2007) 080, [hepth/0611005].
 [27] P. M. Chesler and L. G. Yaffe, Holography and colliding gravitational shock waves in asymptotically AdS spacetime, Phys. Rev. Lett. 106 (2011) 021601, [1011.3562].
 [28] M. P. Heller, R. A. Janik and P. Witaszczyk, The characteristics of thermalization of boostinvariant plasma from holography, Phys. Rev. Lett. 108 (2012) 201602, [1103.3452].
 [29] S. de Haro, S. N. Solodukhin and K. Skenderis, Holographic reconstruction of spacetime and renormalization in the AdS / CFT correspondence, Commun. Math. Phys. 217 (2001) 595–622, [hepth/0002230].
 [30] K. Skenderis, Asymptotically Antide Sitter spacetimes and their stress energy tensor, Int. J. Mod. Phys. A16 (2001) 740–749, [hepth/0010138].