We review recent developments encompassing the description of quantum chaos in holography. We discuss the characterization of quantum chaos based on the late time vanishing of out-of-time-order correlators and explain how this is realized in the dual gravitational description. We also review the connections of chaos with the spreading of quantum entanglement and diffusion phenomena.
Ministry of Science, ICT and Future PlanningNRF2017R1A2B4004810Gwangju Institute of Science and Technology1. Introduction
The characterization of quantum chaos is fairly complicated. Possible approaches range from semiclassical methods to random matrix theory: in the first case one studies the semiclassical limit of a system whose classical dynamics is chaotic; in the later approach the characterization of quantum chaos is made by comparing the spectrum of energies of the system in question to the spectrum of random matrices [1]. Despite the insights provided by the above-mentioned approaches, a complete and more satisfactory understanding of quantum chaos remains elusive.
Surprisingly, new insights into quantum chaos have come from black holes physics! In the context of so-called gauge-gravity duality [2–4], black holes in asymptotically AdS spaces are dual to strongly coupled many-body quantum systems. It was recently shown that the chaotic nature of many-body quantum systems can be diagnosed with certain out-of-time-order correlation (OTOC) functions which, in the gravitational description, are related to the collision of shock waves close to the black hole horizon [5–9]. In addition to being useful for diagnosing chaos in holographic systems and providing a deeper understanding for the inner-working mechanisms of gauge-gravity duality, OTOCs have also proved useful in characterizing chaos in more general nonholographic systems, including some simple models like the kicked-rotor [10], the stadium billiard [11], and the Dicke model [12].
In this paper we review the recent developments in the holographic description of quantum chaos. We discuss the characterization of quantum chaos based on the late time vanishing of OTOCs and explain how this is realized in the dual gravitational description. We also review the connections of chaos with spreading of quantum entanglement and diffusion phenomena. We focus on the case of d-dimensional gravitational systems with d⩾3, which excludes the case of gravity in AdS2 and SYK-like models [13–16]. (Another interesting perspective on the characterization of chaos in the context of (regularized) AdS2/CFT1 is provided by [17–19].) Also, due the lack of the author’s expertise, we did not cover the recent developments in the direct field theory calculations of OTOCs. This includes calculations for CFTs [20], weakly coupled systems [21, 22], random unitary models [23–25], and spin chains [26–30].
2. A Bird Eye’s View on Classical Chaos
In this section we briefly review some basic aspects of classical chaos. For definiteness we consider the case of a classical thermal system with phase space denoted as X=(q,p), where q and p are multidimensional vectors denoting the coordinates and momenta of the phase space. We can quantify whether the system is chaotic or not by measuring the stability of a trajectory in phase space under small changes of the initial condition. Let us consider a reference trajectory in phase space, X(t), with some initial condition X(0)=X0. A small change in the initial condition X0→X0+δX0 leads to a new trajectory X(t)→X(t)+δX(t). This is illustrated in Figure 1. For a chaotic system, the distance between the new trajectory and the reference one increases exponentially with time(1)δXt~δX0eλtor∂Xt∂X0~eλt,where λ is the so-called Lyapunov exponent. This should be contrasted with the behavior of nonchaotic systems, in which δX(t) remains bounded or increases algebraically [31].
Variation of a trajectory in the phase space under small modifications of the initial condition. For a chaotic system the distance between two initially nearby trajectories increases exponentially with time, i.e., δqt=δq0eλt.
The exponential increase depends on the orientation of δX0 and this leads to a spectrum of Lyapunov exponents, {λ1,λ2,…,λK}, where K is the dimensionality of the phase space. A useful parameter characterizing the trajectory instability is (2)λmax=limt→∞limδX0→01tlogδXtδX0,which is called the maximum Lyapunov exponent. When the above limits exist and λmax>0, the trajectory shows sensitivity to initial conditions and the system is said to be chaotic [31].
The chaotic behavior can be a consequence of either a complicated Hamiltonian or simply the contact with a thermal heat bath. This is because chaos is a common property of thermal systems. For the latter to make contact with black holes physics, we consider the case of a classical thermal system with inverse temperature β. If F(X) is some function of the phase space coordinates, we define its classical expectation value as(3)Fβ=∫dXe-βHXFX∫dXe-βHXwhere H(X) is the system’s Hamiltonian.
Classical thermal systems have two exponential behaviors that have analogues in terms of black holes physics: the Lyapunov behavior, characterizing the sensitive dependence on initial conditions, and the Ruelle behavior, characterizing the approach to thermal equilibrium [32, 33].
To quantify the sensitivity to initial conditions in a thermal system we need to consider thermal expectation values. Note that (1) can have either signs. To avoid cancellations in a thermal expectation values, we consider the square of this derivative (4)Ft=∂Xt∂X02β.The expected behavior of this quantity is the following [34] (5)Ft~∑kcke2λkt,where ck are constants and λk are the Lyapunov exponents. At later times the behavior is controlled by the maximum Lyapunov exponent F~e2λmaxt.
The approach to thermal equilibrium or, in other words, how fast the system forgets its initial condition can be quantified by two-point functions of the form (6)Gt=XtX0β-Xβ2,whose expected behavior is [34] (7)Gt~∑jbje-μjt,where bj are constants and μj are complex parameters called Ruelle resonances. The late time behavior is controlled by the smallest Ruelle resonance G~e-μmint.
3. Some Aspects of Quantum Chaos
In this section we review some aspects of quantum chaos. For a long time, the characterization of quantum chaos was made by comparing the spectrum of energies of the system in question to the spectrum of random matrices or using semiclassical methods [1]. Here we follow a different approach, which was first proposed by Larkin and Ovchinnikov [35] in the context of semiclassical systems, and it was recently developed by Shenker and Stanford [6–8] and by Kitaev [9].
For simplicity, let us consider the case of a one-dimensional system, with phase space variables (q,p). Classically, we know that ∂q(t)/∂q(0) grows exponentially with time for a chaotic system. The quantum version of this quantity can be obtained by noting that (8)∂qt∂q0=qt,p0P.B.,where {q(t),p(0)}P.B. denotes the Poisson bracket between the coordinate q(t) and the momentum p(0). The quantum version of ∂q(t)/∂q(0) can then be obtained by promoting the Poisson bracket to a commutator (9)qt,p0P.B.→1iħq^t,p^0where now q^(t) and p^(0) are Heisenberg operators.
We will be interested in thermal systems, so we would like to calculate the expectation value of [q^(t),p^(0)] in a thermal state. However, this commutator might have either signs in a thermal expectation value and this might lead to cancellations. To overcome this problem, we consider the expectation value of the square of this commutator (10)Ct=-q^t,p^02β,where β is the system’s inverse temperature and the overall sign is introduced to make C(t) positive. More generally, one might replace q^(t) and p^(0) by two generic Hermitian operators V and W and quantify chaos with the double commutator(11)Ct=-Wt,V02β.This quantity measures how much an early perturbation V affects the later measurement of W. As chaos means sensitive dependence on initial conditions, we expect C(t) to be ‘small’ in nonchaotic system and ‘large’ if the dynamics is chaotic. In the following we give a precise meaning for the adjectives ‘small’ and ‘large’.
For some class of systems the quantum behavior of C(t) has a lot of similarities with the classical behavior of 〈(∂q(t)/∂q(0))〉β. However, the analogy between the classical and quantum quantities is not perfect because there is not always a good notion of a small perturbation in the quantum case (remember that classical chaos is characterized by the fact that a small perturbation in the past has important consequences in the future). If we start with some reference state and then perturb it, we easily produce a state that is orthogonal to the original state, even when we change just a few quantum numbers. Because of that it seems unnatural to quantify the perturbation as small. Fortunately, there are some quantum systems in which the notion of a small perturbation makes perfect sense. An example is provided by systems with a large number of degrees of freedom. In this case a perturbation involving just a few degrees of freedom is naturally a small perturbation.
For some class of chaotic systems, which include holographic systems, C(t) is expected to behave as (see [30, 36] for a discussion of different possible OTOC growth forms)(12)Ct~Ndof-1fort<td,Ndof-1expλLtfortd≪t≪t∗O1fort>t∗,where Ndof is the number of degrees of freedom of the system. Here, we have assumed V and W to be unitary and Hermitian operators, so that VV=WW=1. The exponential growth of C(t) is characterized by the Lyapunov exponent (this is actually the quantum analogue of the classical Lyapunov exponent; the two quantities are not necessarily the same in the classical limit [21]; here we stick to the physicists long standing tradition of using misnomers and just refer to λL as the Lyapunov exponent) λL and takes place at intermediate time scales bounded by the dissipation time td and the scrambling time t∗. The dissipation time is related to the classical Ruelle resonances (td~μ-1) and it characterizes the exponential decay of two-point correlators, e.g., 〈V(0)V(t)〉~e-t/td. The dissipation time also controls the late time behavior of C(t). The scrambling time t∗~λL-1logNdof is defined as the time at which C(t) becomes of order O(1). See Figure 2. The scrambling time controls how fast the chaotic system scrambles information. If we perturb the system with an operator that involves only a few degrees of freedom, the information about this operator will spread among the other degrees of freedom of the system. After a scrambling time, the information will be scrambled among all the degrees of freedom and the operator will have a large commutator with almost any other operator.
Schematic form of C(t). We indicated the regions of Lyapunov and Ruelle behavior. C(t)~O(1) at t~t∗.
To understand how the above behavior relates to chaos, we write the double commutator as (13)Ct=-Wt,V02β(14)=2-2WtV0WtV0β,where we made the assumption that W and V are Hermitian and unitary operators. Note that all the relevant information about C(t) is contained in the OTOC: (15)OTOt=WtV0WtV0.The fact that C(t) approaches 2 at later times implies that the OTO(t) should vanish in that limit. To understand why this is related to chaos we think of OTO(t) as an inner-product of two states (16)OTOt=ψ2∣ψ1,where (17)ψ1=W-tV0β,ψ2=V0W-tβwhere β is some thermal state and we replace t→-t to make easier the comparison with black holes physics.
If [V(0),W(t)]≈0 for any value of t, the two states are approximately the same, and ψ1∣ψ2≈1, implying C(t)≈0. That means the system displays no chaos—the early measurement of V has no effect on the later measurement of W. If, on the other hand, [V(0),W(t)]≠0, the states ψ1 and ψ2 will have a small superposition ψ1∣ψ2≈0, implying C(t)≈2. That means that V has a large effect on the later measurement of W.
In Figure 3 we construct the states ψ1 and ψ2 and explain why ψ1∣ψ2≈0 for large t means chaos. Let us start by constructing the state ψ1=W-tV0β=e-iHtW0eiHtV0β. The unperturbed thermal state is represented by a horizontal line. We initially consider the state V0β, which is the thermal state perturbed by V. If we evolve the system backwards in time (applying the operator eiHt) for some time which is larger than the dissipation time, the system will thermalize and it will no longer display the perturbation V. After that, we apply the operator W, which should be thought of as a small perturbation, and then we evolve the system forwards in time (applying the operator e-iHt). The final results of this set of operations depend on the nature of the system. If the system is chaotic, the perturbation W will have a large effect after a scrambling time, and the perturbation that was present at t=0 will no longer rematerialize. This is illustrated in Figure 3. In contrast, for a nonchaotic system, the perturbation W will have little effect on the system at later times, and the perturbation V will (at least partially) rematerialize at t=0.
Construction of the state W-tV0β. For a chaotic system the perturbation V fails to rematerialize at t=0. In a nonchaotic system we expect the perturbation V to rematerialize at t=0.
We now construct the state ψ2=V0W-tβ=V0e-iHtW0eiHtβ. This is illustrated in Figure 4. We start with the thermal state β and then we evolve this state backwards in time eiHtβ. After that, we apply the operator W and then we evolve the system forwards in time, obtaining the state e-iHtW0eiHtβ. Finally, we apply the operator V, obtaining the state V0W-tβ. Note that, by construction, this state displays the perturbation V at t=0, while the state W-tV0β does not. As a consequence, the two states are expected to have a small superposition W-tV0∣W-tV0β≈0. This should be contrasted to the case where the system is not chaotic. In this case the perturbation V rematerializes at t=0, and the states ψ1 and ψ2 have a large superposition, i.e., W-tV0∣W-tV0β≈1.
Construction of the state V0W-tβ. By construction, this state displays the perturbation V at t=0.
In this construction we assumed the operators V(0) and W(-t) to be separated by a scrambling time, i.e., t>t∗. This is important because, at earlier times, the two operators, which in general involve different degrees of freedom of the system, generically commute. The operators manage to have a nonzero commutator at later times because of the phenomenon of operation growth that we will describe in the next section.
3.1. Operator Growth and Scrambling
The operators V and W act generically at different parts of the physical system, yet they can have a nonzero commutator at later times. This is possible because in chaotic systems the time evolution of an operator makes it more and more complicated, involving and increasing number of degrees of freedom. As a result, an operator that initially involves just a few degrees of freedom becomes delocalized over a region that grows with time. The growth of the operator W(t) is maybe more evident from the point of view of the Baker-Campbell-Hausdorff (BCH) formula, in terms of which we can write (18)Wt=eiHtW0e-iHt=∑k=0∞-itkk!HH,…H,W0….From the above formula it is clear that, at each order in t, there is a more complicated contribution to W(t). In chaotic systems the operator becomes more and more delocalized as the time evolves, and it eventually becomes delocalized over the entire system. The time scale at which this occurs is the so-called scrambling time t∗. After the scrambling time the operator W(t) manages to have a nonzero and large commutator with almost any other operator, even operators involving only a few degrees of freedom.
This can be clearly illustrated in the case of a spin chain. Let us follow [7] and consider an Ising-like model with Hamiltonian (19)H=-∑iZiZi+1+gXi+hZi,where Xi,Yi, and Zi denote Pauli matrices acting on the ith site of the spin chain. The above system is integrable if we take g=1 and h=0, but it is strongly chaotic if we choose g=-1.05 and h=0.5.
To illustrate the concept of scrambling, we consider the time evolution of the operator Z1. Using the BCH formula we can write the following. (20)Z1t=Z1-itH,Z1-t22!H,H,Z1+it33!H,H,H,Z1+…Ignoring multiplicative constants and signs we can write the above terms (schematically) as follows.(21)H,Z1~Y1H,H,Z1~Y1+X1Z2H,H,H,Z1~Y1+Y2X1+Y1Z2H,H,H,H,Z1~X1+Y1+Z1+X1X2+Y1Y2+Z1Z2+X1Z2+Z3Y1+Y1Z2Y2+Z1X2X1+X2Z3X1As the time evolves, higher order terms become important in series (20), and the operator Z1(t) becomes more and more complicated, involving terms in an increasing number of sites. For large enough t the operator will involve all the sites of the spin chain and it will manage to have a nonzero commutator with a Pauli operator in any other site of the system. In this situation the information about Z1 is essentially scramble among all the degrees of freedom of the system. As discussed before, this occurs after a scrambling time. Above this time the double commutator C(t) saturates to a constant value. This should be contrasted to what happens for an integrable system. In this case the operator grows, but it also decreases at later times. In the chaotic case, the operator remains large at later times [7].
3.2. Probing Chaos with Local Operators
In quantum field theories we can upgrade (11) to the case where the operators are separated in space(22)Ct,x=-V0,0,Wt,x2β.Strictly speaking, the above expression is generically divergent, but it can be regularized by adding imaginary times to the time arguments of the operators V and W. For a large class of spin chains, higher-dimensional SYK-models, and CFTs, the above commutator is roughly given by (23)Ct,x~expλLt-t∗-xvB,where vB is the so-called butterfly velocity. (Actually, vB represents the “velocity of the butterfly effect”. Here we continue to follow the tradition of using misnomers.) This velocity describes the growth of the operator W in physical space and it acts as a low-energy Lieb-Robinson velocity [37], which sets a bound for the rate of transfer of quantum information. From the above formula, we can see that there is an additional delay in scrambling due to the physical separation between the operators. The butterfly velocity defines an effective light-cone for commutator (22). Inside the cone, for t-t∗⩾x/vB, we have C(t,x)~O(1), whereas for outside the cone, for t-t∗<x/vB, the commutator is small, C(t,x)~1/Ndof≪1. Outside the light-cone the Lorentz invariance implies a zero commutator. The light-cone and the butterfly effect cone are illustrated in Figure 5.
Light-cone (gray region) and butterfly effect cone (dark gray region). Inside the butterfly effect cone, for t-t∗⩾x/vB, we have C(t,x)~O(1), whereas for outside the cone, for t-t∗<x/vB, the commutator is small, C(t,x)~1/Ndof≪1. Outside the light-cone the Lorentz invariance implies a zero commutator.
4. Chaos and Holography
In this section we review how the chaotic properties of holographic theories can be described in terms of black holes physics. Black holes behave as thermal systems and thermal systems generically display chaos. This implies that black holes are somehow chaotic. This statement has a precise realization in the context of the gauge/gravity duality. According to this duality, some strongly coupled nongravitational systems are dual to higher-dimensional gravitational systems. In the most known and studied example of this duality the N=4 super Yang-Mills (SYM) theory living in R3,1 is dual to type IIB supergravity in AdS5×S5. More generically, a d-dimensional nongravitational theory living in Rd-1,1 is dual to a gravity theory living in a higher-dimensional space of the form AdSd+1×M, where M is generically a compact manifold. The nongravitational theory can be thought of as living in the boundary of AdSd+1 and because of that is usually called the boundary theory. The gravitational theory is also called the bulk theory.
There is a dictionary relating physical quantities in the boundary and bulk description [3, 4]. An example is provided by the operators of the boundary theory, which are related to bulk fields. The boundary theory at finite temperature can be described by introducing a black hole in the bulk. The thermalization properties of the boundary theory have a nice visualization in terms of black holes physics. By applying a local operator in the boundary theory we produce some perturbation that describes a small deviation from the thermal equilibrium. The information about V(x) is initially contained around the point x, but it gets delocalized over a region that increases with time, until it completely melts into the thermal bath. In the bulk theory, the application of the operator V(x) produces a particle (field excitation) close to the boundary of the space, which then falls into the black hole. The return to the thermal equilibrium in the boundary theory corresponds to the absorption of the bulk particle by the black hole. Figure 6 illustrates the bulk description of thermalization.
Bulk picture of thermalization. The figure represents an asymptotically AdS black hole geometry. The boundary is at the top edge, while the black hole horizon is at the bottom edge. The black hole’s interior is shown in gray. The boundary operator V is dual to the bulk field ϕ. From the point of view of the boundary theory the perturbation produced by V is initially localized around the point x, but it gets delocalized over a region that increases with time. In the bulk description this is described by a particle (field excitation) that is initially close to the boundary and then falls into the black hole.
The approach to thermal equilibrium is controlled by the black hole’s quasinormal modes (QNMs). In holographic theories, the quasinormal modes control the decay of two-point functions of the boundary theory(24)VtV0β~e-t/tdwhere the dissipation time td is related to the lowest quasinormal mode (Im(ω)~td-1). From the point of view of the bulk theory the QNMs describes how fast a perturbed black hole returns to equilibrium. Clearly, the black hole’s quasinormal modes correspond to the classical Ruelle resonances. In holographic theories the dissipation time is roughly given by td~β.
Another important exponential behavior of black holes is provided by the blue-shift suffered by the in-falling quanta or, equivalently, the red shift suffered by the quanta escaping from the black hole. The blue-shift suffered by the in-falling quanta is determined by the black hole’s temperature. If the quanta asymptotic energy is E0, this energy increases exponentially with time (25)E=E0e2π/βt,where β is the Hawking’s inverse temperature. Later we will see that this exponential increase in the energy of the in-falling quanta gives rise to the Lyapunov behavior of C(t,x) of holographic theories.
4.1. Holographic Setup
The TFD State & Two-Sided Black Holes. In the study of chaos it is convenient to consider a thermofield double state made out of two identical copies of the boundary theory (26)TFD=1Z1/2∑ne-βEn/2nLnR,where L and R label the states of the two copies, which we call QFTL and QFTR, respectively. The two boundary theories do not interact and only know about each other through their entanglement. This state is dual to an eternal (two-sided) black hole, with two asymptotic boundaries, where the boundary theories live [38]. This is a wormhole geometry, with an Einstein-Rosen bridge connecting the two sides of the geometry. The wormhole is not traversable, which is consistent with the fact that the two boundary theories do not interact.
For definiteness we assume a metric of the form (27)ds2=-Gttrdt2+Grrrdr2+Gijr,xkdxidxj,where the boundary is located at r=∞, where the above metric is assumed to asymptote AdSd+1. We take the horizon as located at r=rH, where Gtt vanishes and Grr has a first order pole. For future purposes, let β be the Hawking’s inverse temperature, and SBH be the Bekenstein-Hawking entropy.
In the study of shock waves it is more convenient to work with Kruskal-Szekeres coordinates, since these coordinates cover smoothly the globally extended spacetime. We first define the tortoise coordinate (28)dr∗=GrrGttdr,and then we introduce the Kruskal-Szekeres coordinates U,V as follows. (29)U=+e2π/βr∗-t,V=-e2π/βr∗+tleft exterior regionU=-e2π/βr∗-t,V=+e2π/βr∗+tright exterior regionU=+e2π/βr∗-t,V=+e2π/βr∗+tfuture interior regionU=-e2π/βr∗-t,V=-e2π/βr∗+tpast interior regionIn terms of these coordinates the metric reads (30)ds2=2AUVdUdV+GijUVdxidxj,where (31)AUV=β28π2GttUVUV.In these coordinates the horizon is located at U=0 or at V=0. The left and right boundaries are located at UV=-1 and the past and future singularities at UV=1. The Penrose diagram for this metric is shown in Figure 7.
Penrose diagram for the two-sided black holes with two boundaries that asymptote AdS. This geometry is dual to a thermofield double state constructed out of two copies of the boundary theory.
The global extended spacetime can also be described in terms of complexified coordinates [39]. In this case one defines the complexified Schwarzschild time (32)t=tL+itE,where tL and tE are the Lorentzian and Euclidean times, and then one describes the time in each of the four patches (left and right exterior regions, and the future and past interior regions) as having a constant imaginary part. (33)tE=0right exterior regiontE=-β4future interior regiontE=-β2left exterior regiontE=+β4past interior regionThe Euclidean time has a period of β. The Lorentzian time increases upward (downward) in the right (left) exterior region, and to the right (left) in the future (past) interior.
Note that, with the complexified time, one can obtain an operator acting on the left boundary theory by adding (or subtracting) iβ/2 to the time of an operator acting on the right boundary theory.
Perturbations of the TFD State & Shock Wave Geometries. We now turn to the description of states of the form (34)WtTFDwhere W is a thermal scale operator that acts on the right boundary theory. This state can be describe by a ‘particle’ (field excitation) in the bulk that comes out of the past horizon, reaches the right boundary at time t, and then falls into the future horizon, as illustrated in Figure 8.
Bulk description of the state WtTFD. In blue is shown the trajectory of a ‘particle’ that comes out of the past horizon, reaches the boundary at time t producing the perturbation W, and then falls into the future horizon. From now on, we will refer to this bulk excitation as the W-particle.
If t is not too large, the state WtTFD will represent just a small perturbation of the TFD state and the corresponding description in the bulk will be just an eternal two-sided black hole geometry slightly perturbed by the presence of a probe particle. This is no longer the case if t is large. In this case there is a nontrivial modification of the geometry. A very early perturbation, for example, is described in the bulk in terms of a particle that falls towards the future horizon for a very long time and gets highly blue-shifted in the process. If the particle’s energy is E0 in the asymptotic past, this energy will be exponentially larger from the point of view of the t=0 slice of the geometry, i.e., E=E0e2π/βt. Therefore, for large enough t, the particle’s energy will be very large and one needs to include the corresponding back-reaction.
The back-reaction of a very early (or very late) perturbation is actually very simple—it corresponds to a shock wave geometry [40, 41]. To understand that, we first need to notice that, under boundary time evolution, the stress energy of a generic perturbation W gets compressed in the V-direction and stretched in the U-direction. For large enough t we can approximate the stress tensor of the W-particle as (35)TVV~PUδVax→,where PU~β-1e2π/βt is the momentum of the W-particle in the U-direction and a(x→) is some generic function that specifies the location of the perturbation in the spatial directions of the right boundary. Note that TVV is completely localized at V=0 and homogeneous along the U-direction. Besides, even if the W-particle is massive, the exponential blue-shift will make it follow an almost null trajectory, as shown in Figure 9.
Bulk description of the state W-tTFD. An early enough perturbation produces a shock wave geometry. The effect of the shock wave (shown in blue) is to produce a shift U→U+h(t,x→) in the trajectory of a probe particle (shown in red) crossing it.
The shock wave geometry produced by the W-particle is described by the metric (36)ds2=2AUVdUdV+GijUVdxidxj-2AUVht,x→δVdV2,which is completely specified by the shock wave transverse profile h(t,x→). This geometry can be seen as two pieces of an eternal black hole glued together along V=0 with a shift of magnitude h(t,x→) in the U-direction. We find it useful to represent this geometry with the same Penrose diagram of the unperturbed geometry, but with the prescription that any trajectory crossing the shock wave gets shifted in the U-direction as U→U+h(t,x→). See Figure 9.
The precise form of h(t,x→) can be determined by solving the VV-component of Einstein’s equation. For a local perturbation, i.e., a(x→)=δd-1(x→), the solution reads(37)ht,x→~GNe2π/βt-μx→,withμ=2πβd-1Gii′rHGtt′rH,where, for simplicity, Gij has been assumed to be diagonal and isotropic.
Interestingly, the shock wave profile contains information about the parameters characterizing the chaotic behavior of the boundary theory. Indeed, the double commutator has a region of exponential growth at which C(t,x→)~h(t,x→). From this identification, we can write (38)ht,x→~e2π/βt-t∗-x→/vBwhere (the leading order contribution to) the scrambling time scales logarithmically with the Bekenstein-Hawking entropy (39)t∗~β2πlog1GN~β2πlogSBH,while the Lyapunov exponent is proportional to the Hawking’s temperature. (40)λL=2πβ.The butterfly velocity is determined from the near-horizon geometry. (Here we are assuming isotropy. In the case of anisotropic metrics the formula for vB is a little bit more complicated. See, for instance, Appendix A of [42] or Appendix B of [43].)(41)vB2=Gtt′rHd-1Gii′rH.
4.2. Bulk Picture for the Behavior of OTOCs
In this section we present the bulk perspective for the vanishing of OTOCs at later times. In order to do that, we write the OTOC as a superposition of two states (42)OTOt=TFDW-tV0W-tV0TFD=ψout∣ψin,where the ‘in’ and ‘out’ states are given by the following. (43)ψin=W-tV0TFD,ψout=V†0W†-tTFDThe interpretation of a vanishing OTOC in terms of the bulk theory is actually very simple. Let us go step by step and construct first the state V0β. This state is described by a particle that comes out of the past horizon, reaches the boundary at t=0, and then falls back into the future horizon. See the left panel of Figure 10.
Left panel: bulk description of the state V0TFD. The V-particle comes out of the past horizon, reaches the boundary at time t=0 producing the perturbation V, and then falls into the future horizon. Right panel: bulk description of the ‘in’ state ψin=W-tV0TFD. The W-particle (shown in blue) produces a shock wave along V=0. The trajectory of the V-particle (shown in red) suffers a shift and the perturbation V is produced at the boundary with some time delay.
Now the ‘in’ state can be obtained as (44)ψin=W-tV0TFD=e-iHtW0eiHtV0TFD.This amounts to the following: evolving the state V0TFD backwards in time, applying the operator W, and then evolving the system forwards in time. The corresponding description in the bulk is shown in the right panel of Figure 10. From this picture we can see that the perturbation W produces a shock wave that causes a shift in the trajectory of the V-particle, which no longer reaches the boundary at time t=0, but rather with some time delay. The physical interpretation is that a small perturbation in the asymptotic past (represented by W) is amplified over time and destroys the initial configuration (represented by the state V0TFD).
The bulk description of the ‘out’ state can be obtained in the same way. As this state displays the perturbation V at t=0, the V-particle should be produced in the asymptotic past in such a way that, after its trajectory gets shifted as U→U+h, it reaches the boundary at the time t=0 producing the perturbation V.
Comparing the bulk description of the state |ψin〉 (shown in the right panel of Figure 10) with the description of the state |ψout〉 (shown in Figure 11) we can see that these states are indistinguishable when h(t,x→) is zero, but they become more and more different for large values of h(t,x→). As a consequence, the overlap C(t)=ψout∣ψin is equal to one when h=0, but it decreases to zero as we increase the value of h.
Bulk description of the ‘out’ state ψout=V0W-tTFD. The W-particle produces the shock wave geometry. The trajectory of the V-particle is such that, after suffering the shift U→U+h(t,x→), it reaches the boundary at time t=0, producing the perturbation V.
The exponential behavior of h(t,x→) implies that an early enough perturbation can produce a very large shift in the V-particle’s trajectory, causing it to be captured by the black hole and preventing the materialization of the V perturbation at the boundary. See Figure 12. This should be compared with the physical picture given in Figure 3.
Bulk description of the state W-tV0TFD for the case where t≳t∗. The V-particle’s trajectory undergoes a shift, and it is captured by the black hole. The perturbation V never forms, and the corresponding state has no superposition with the ‘out’ state V0W-tTFD, resulting in a vanishing OTOC.
The physical picture of the process described in Figure 12 is quite simple. The state V0TFD can be represented by a black hole geometry in which a particle (the V-particle) escapes from the black holes and reaches the boundary at time t=0. The state W-tV0TFD is obtained by perturbing the state V0TFD in the asymptotic past. This corresponds to the addition of a W-particle to the system in the asymptotic past. This particle gets highly blue-shifted as it falls towards the black hole. The black hole captures the W-particle and becomes bigger. The V-particle fails to escape from the bigger black hole and never reaches the boundary to produce the V perturbation. This physical picture is illustrated in Figure 13.
Left panel: bulk description of the state V0TFD. The V-particle (whose trajectory is shown in red) escapes from the black hole (shown in gray), reaches the boundary at time t=0, and then falls back towards the horizon. Right panel: bulk description of the state W-tV0TFD for the case where t≳t∗. The W-particle (whose trajectory is shown in blue) gets blue-shifted and increases the black hole size as it falls into it. The V-particle fails to escape from the larger black hole, and the perturbation V never forms at the boundary.
The precise form of the above OTOC can be obtained by calculating the overlap ψout∣ψin using the Eikonal approximation [8], in which the Eikonal phase δ is proportional to the shock wave profile δ~h(t,x→). The OTOC can be written as an integral of the phase eiδ weighted by kinematical factors which are basically Fourier transforms of bulk-to-boundary propagators for the V and W operators.
The result for Rindler AdS3 reads (the below result assumes ΔW≫ΔV) (45)Viϵ1Wt+iϵ2Viϵ3Wt+iϵ4Viϵ1Viϵ3Wiϵ2Wiϵ4=11-8πiGNΔW/ϵ13ϵ24∗e2π/βt-x→/vBΔVwhere ΔV and ΔW are the scaling dimensions of the operators V and W, respectively, and ϵij=i(eiϵi-eiϵj). For this system β=2π and vB=1. This formula matches the direct CFT calculation (the CFT perspective for the onset of chaos has been widely discussed in [44]; other references in this direction include, for instance, [45–48].) obtained in [20]. It can also be derived using the geodesic approximation for two-sided correlators in a shock wave background [5, 20].
Expanding the above result for small values of GNe2π/βt-x→/vB, we obtain(46)OTOt=1-8πiGNΔVΔWϵ13ϵ24∗e2π/βt-x→/vB,and, since h(t,x→)~GNe2π/βt-|x→|/vB, the above result implies (47)Ct,x→~ht,x→.The above result is valid for small (in AdS/CFT the Newton constant is related to the rank of the gauge group of dual CFT as GN~1/Na, where a is a positive number that depends on the dimensionality of the bulk space time (cf. section 7.2 of [49]); our classical gravity calculations are only valid in the large-N limit (that suppresses quantum corrections) so it is natural to consider GN as a small parameter) values of GN, or for any value of GN, but for times in the range td≪t≪t∗, where t∗=β/2πlog1/GN.
Despite being true in the Rindler AdS3 case, the proportionality between the double commutator and the shock wave profile has not been demonstrated in more general cases. However, the authors of [8] argued that, in regions of moderate scattering between the V- and W-particle, the identification C(t,x→)~h(t,x→) is approximately valid.
At very late times, the behavior of the OTO(t) is expected to be controlled by the black hole quasinormal modes. Indeed, in the case of a compact space it is possible to show that (48)Ct~e-2iωt-t∗-R/vB,withImω<0,where R is the diameter of the compact space and ω is the system lowest quasinormal frequency [8].
4.2.1. Stringy Corrections
In this section we briefly discuss the effects of stringy corrections to the Einstein gravity results for OTOCs. We start by reviewing the Einstein gravity results from the perspective of scattering amplitudes. In the framework of the Eikonal approximation, the phase shift suffered by the V-particle is given by (49)δ=-PVht,x→~GNs,where we used the fact that h(t,x→)~GNPU and introduced a Mandelstam-like variable s=2A(0)PUPV. In a small-GN expansion the double commutator C(t) and the phase shift δ scale with s in the same way, namely, (50)Ct~GNs,where s~β-2e2π/βt.
The string corrections can be incorporated using the standard Veneziano formula for the relativistic scattering amplitude A~sδ. The phase shift can then be schematically written as an infinite sum (51)δ~∑JGNsJ-1,where each term corresponds to the contribution due to the exchange of a spin-J field. In Einstein gravity the dominant contribution comes from the exchange of a spin-2 field, the graviton. In string theory, we have to include an infinite tower of higher spin fields. Naively, it looks like these higher spin contributions will increase the development of chaos. However, the resummation of the above sum actually leads to a decrease in the development of chaos. The string-corrected phase shift has a milder dependence with s, namely, (52)δ~GNsJeff-1,with the effective spin given by [8] (53)Jeff=2-dd-1ls24lAdS2where ls is the string length, lAdS is the AdS length scale, and d is the number of dimensions of the boundary theory. As a result, the string-corrected double commutator grows in time with an effective smaller Lyapunov exponent (54)Cstringt~e2π/β1-dd-1ls2/4lAdS2t,and this leads to a larger scrambling time. (At small scales, the string-corrected shock wave has a Gaussian profile, and the concept of butterfly velocity is not meaningful. It was recently shown, however, that at larger scales is possible to define a string-corrected butterfly velocity. The result for N=4 SYM theory reads [50] vB=2/31+23ζ3/161/λ3/2, where λ is the ’t Hooft coupling, which can be written in terms of string length scale as λ=(lAdS/ls)4.)(55)t∗string=t∗1+dd-1ls24lAdS2.
The above discussion implies that for a theory with a finite number of high-spin fields (J>2) chaos would develop faster than in Einstein gravity. These theories, however, are known to violate causality [51]. It is then natural to speculate that the Lyapunov exponent obtained in Einstein gravity has the maximal possible value allowed by causality. This is indeed true and this is the topic of the next section.
4.2.2. Bounds on Chaos
One of the remarkable insights that came from the holographic description of quantum chaos is the fact that there is a bound on chaos—the quantum Lyapunov exponent is bounded from above, while the scrambling time is bounded from below. A distinct feature of holographic systems is that they saturate these two bounds.
Let us follow the historical order and start by discussing the lower bound on the scrambling time. In black holes physics the scrambling time defines how fast the information that has fallen into a black holes can be recovered from the emitted Hawking radiation. (This assumes that half of the black hole’s initial entropy has been radiated [52].) In the context of the Hayden-Preskill thought experiment, the scrambling time is barely compatible with black hole complementarity [52], since a smaller scrambling time would lead to a violation of the no-cloning principle. This led Susskind and Sekino to conjecture that black holes are the fastest scramblers in nature; i.e., they have the smallest possible scrambling time [53]. The lower bound on the scrambling time of a generic many-body quantum system can be written as (56)t∗⩾CβlogNdofwhere C(β) is some function of the inverse temperature. In the case of black holes this function is simply given by C(β)=β/2π.
The scrambling time defines a stronger notion of thermalization and should not be confused with the dissipation time. In fact, for black holes, one expects the dissipation time to be given by the black hole quasinormal modes (this is true in the case of low dimension operators) td~β, while the scrambling time is parametrically larger t∗~βlogNdof. This brings us to the second bound on chaos: for systems with such a large hierarchy between the scrambling and the dissipation, time is possible to derive an upper bound for the Lyapunov exponent [54]: (57)λL⩽2πβ.One should emphasize that this bound does not depend on the existence of a holographic dual. It can be derived for generic many-body quantum systems under some very reasonable assumptions.
The fact that black hole always has a maximum Lyapunov exponent led to the speculation that the saturation of the chaos bound might be a sufficient condition for a system to have an Einstein gravity dual [9, 54]. In fact, there have been many attempts to use the saturation of the chaos bound as a criterion to discriminate holographic CFTs from the nonholographic ones [20, 44–48, 55, 56]. It was recently shown, however, that this criterion, though necessary, is insufficient to guarantee a dual description purely in terms of Einstein gravity [57, 58].
Since vB defines the speed at which information propagates, it is natural to question whether this quantity is also bounded. From the perspective of the boundary theory, causality implies (58)vB⩽1,meaning that information should not propagate faster than the speed of light. Indeed, the above bound can be derived in the context of Einstein gravity by using Null Energy Condition (NEC) and assuming an asymptotically AdS geometry (this derivation uses an alternative definition for vB, which is based on entanglement wedge subregion duality [59]) [60]. This is consistent with the expectation that gravity theories in asymptotically AdS geometries are dual to relativistic theories. In contrast, for geometries which are not asymptotically AdS, the butterfly velocity can surpass the speed of light [42, 60], which is consistent with the non-Lorentz invariance of the corresponding boundary theories.
If we further assume isotropy, it is possible to derive a stronger bound for vB [61] (59)vB⩽vBSch=d2d-1,where vBSch is the value of the butterfly velocity for an AdS-Schwarzschild black brane in d+1 dimensions. This is also the butterfly velocity for a d-dimensional thermal CFT.
The above formula shows that, for thermal CFTs, vB does not depend on the temperature. However, if we deform the CFT, vB acquires a temperature dependence as we move along the corresponding renormalization group (RG) flow. In fact, by considering deformations that break the rotational symmetry, it was noticed that the butterfly velocity violates the above bound, but remains bounded from above by its value at the infrared (IR) fixed point, never surpassing the speed of light [62–64]. The above bound can also be violated by higher curvature corrections, but vB remains bounded by the speed of light as long as causality is respected. (For instance, in 4-dimensional Gauss-Bonnet (GB) gravity, the butterfly velocity surpasses the speed of light for λGB<-3/4, but causality requires λGB>-0.19[65, 66]. Moreover, it was recently shown that, unless one adds an infinite tower of extra higher spin fields, GB gravity might be inconsistent with causality for any value of the GB coupling [51].) The violation of the bound given in (59) by anisotropy or higher curvature corrections is reminiscent of the well-known violation of the shear viscosity to entropy density ratio bound [67–72].
4.3. Chaos and Entanglement Spreading
The thermofield double state displays a very atypical left-right pattern of entanglement that results from nonzero correlations between subsystems of QFTL and QFTR at t=0. The chaotic nature of the boundary theories is manifested by the fact that small perturbations added to the system in the asymptotic past destroy this delicate correlations [5].
The special pattern of entanglement can be efficiently diagnosed by considering the mutual information I(A,B) between spatial subsystems A⊂QFTL and B⊂QFTR, defined as (60)IA,B=SA+SB-SA∪B,where SA is the entanglement entropy of the subsystem A and so on. The mutual information is always positive and provides an upper bound for correlations between operators OL and OR defined on A and B, respectively [73], (61)IA,B⩾OLOR-OLOR22OL2OR2.
The thermofield double state has nonzero mutual information between large (for small subsystems, the mutual information is zero) subsystems of the left and right boundary, signaling the existence of left-right correlations. These correlations can be destroyed by small perturbations in the asymptotic past, meaning that initially positive mutual information drops to zero when we add a very early perturbation to the system.
Interestingly, the vanishing of the mutual information can be connected to the vanishing of the OTOCs discussed earlier. If, for simplicity, we assume that OL and OR have zero thermal one point function, then the disruption of the mutual information implies the vanishing of the following four-point function (62)OLORW=TFDWR†OLORWRTFD=0,which is related by analytic continuation to the one-sided out-of-time-order correlator introduced earlier. (To obtain an OTOC with operators acting only on the right boundary theory, one just needs to add iβ/2 to time argument of the operator OL in the above formula.)
The disruption of the mutual information has very simple geometrical realization in the bulk. The entanglement entropies that appear in the definition of I(A,B) can be holographically calculated using the HRRT prescription [74, 75] (63)SA=AreaγA4GN,where γA is an extremal surface whose boundary coincides with the boundary of the region A. There is an analogous formula for SB. Both γA and γB are U-shaped surfaces lying outside of the event horizon, in the left and right side of the geometry, respectively. There are two candidates for the extremal surface that computes SA∪B: the surface γA∪γB or the surface γwormhole that connects the two asymptotic boundaries of the geometry. See Figure 14. According to the RT prescription, we should pick the surface with less area. If γA∪γB has less area than γwormhole, then I(A,B)=0, because Area(γA∪γB)=Area(γA)+Area(γB). On the other hand, if γwormhole has less area than γA∪γB, i.e., Area(γwormhole) < Area(γA)+Area(γB), then we have a positive mutual information (64)IA,B=14GNAreaγA+AreaγB-Areaγwormhole>0.Now, an early perturbation of the thermofield double state gives rise to a shock wave geometry in which the wormhole becomes longer. As a consequence, the area of the surface γwormhole increases, resulting in a smaller mutual information. It is then clear that the mutual information will drop to zero if the wormhole is longer enough. The length of the wormhole depends on the strength of the shock wave, which, by its turn, depends on how early the perturbation is producing it. Therefore, an early enough perturbation will produce a very long wormhole in which the mutual information will be zero. The fact that the shock wave geometry produces a longer wormhole (along the t=0 slice of the geometry) is clearly seen if we represent the shock wave geometry with a tilted Penrose diagram. See, for instance, Figure 3 of [76].
Illustration of the entangling surfaces in the t=0 slice of a two-sided black brane geometry. The U-shaped surfaces (γA and γB) are represented by blue curves. The surface stretching through the wormhole is given by the union of the two red surfaces γwormhole=γ1∪γ2. In the left panel we represent the unperturbed geometry, in which the two horizons coincide. In the right panel we represent the geometry in the presence of a shock wave added at some time t0 in the past. In this case the size of the wormhole is effectively larger, and the two horizons no longer coincide.
The mutual information I(A,B) decreases as a function of the time t0 at which we perturbed the system. For t0≳t∗, the mutual information decreases linearly with behavior controlled by the so-called entanglement velocity vE [63](65)dIA,Bdt0=-dSA∪Bdt0=-vEsthAreaA∪B,where sth is the thermal entropy density and Area(A∪B) is the area of A∪B (or the volume of the boundary of this region). The two-sided black hole geometry with a shock wave can be thought of as an additional example of a holographic quench protocol [63], and the time-dependence of entanglement entropy can be understood in terms of the so-called ‘entanglement tsunami’ picture. See [77] for field theory calculations and [78–82] for holographic calculations. However, it was recently shown that the entanglement tsunami picture is not very sharp. See [59] for further details. In [81, 82], the entanglement velocity was conjectured to be bounded as (66)vE⩽vESch=dd-11/2-1/d2d-11-1/d,where vESch is the entanglement velocity for a (d+1)-dimensional Schwarzschild black brane or, equivalently, the value of vE for a d-dimensional thermal CFT. This bound can be derived in the context of Einstein gravity assuming an asymptotically AdS geometry, isotropy, and NEC [61]. Just like in the case of vB, the entanglement velocity in thermal CFTs does not depend on the temperature. But vE acquires a temperature dependence if we deform the CFT and move along the corresponding RG flow [63, 64]. In these cases, vE violates the above bound, but it remains bounded by its corresponding value at the IR fixed point, never surpassing the speed of light.
One can also prove that the entanglement velocity is also bounded by the speed of light. (See [83, 84] for a discussion about small subsystems.) This can be done by using positivity of the mutual information [85] or using inequalities involving the relative entropy [86]. More generally, the authors of [59] conjecture that vE⩽vB, which implies the bound vE<1 in the cases where vB is bounded. However, both [85, 86] assumed that the theory is Lorentz invariant. In the case of non-Lorentz invariant theories (e.g., noncommutative gauge theories) the entanglement velocity can surpass the speed of light. This has been verified both in holography calculations [42] and in field theory calculations [87].
Finally, we mention that other concepts from information theory can also be used to diagnose chaos in holography. It has been shown, for instance, that the relative entropy is also a useful tool to diagnose chaotic behavior [88]. For a connection between chaos and computational complexity, see, for instance [89, 90].
4.4. Chaos and Hydrodynamics
Recently, there has been a growing interest in the connection between chaos and hydrodynamics [91–99]. Here we briefly review some interesting connection between chaos and diffusion phenomena.
A longstanding goal of quantum condensed matter physics is to have a deeper understanding of the so-called ‘strange metals’. These are strongly correlated materials that do not have a description in terms of quasiparticles excitations and whose transport properties display a remarkable degree of universality. In [100, 101] Sachdev and Damle proposed that such a universal behavior could be explained by a fundamental dissipative timescale (67)τP~ħ2πkBT,which would govern the transport in such systems.
Interestingly, the Lyapunov exponent defines a time scale τL=1/λL, and the upper bound on λL translates into a lower bound for τL that precisely coincides with τP(68)τL⩾ħβ2πkB,where we reintroduced ħ and the Boltzmann constant in the expression for the bound on the Lyapunov exponent. (In systems of units where ħ and kB are not equal to one, the bound on the Lyapunov exponent reads λL⩽2πkB/ħβ.) Holographic systems saturate the above bound, and this explains the universality observed in the transport properties of these systems.
A prototypical example of universality is the linear resistivity of strange metals. In [102], Hartnoll proposed that the linear resistivity could be explained by the existence of a universal lower bound on the diffusion constants related to the collective diffusion of charge and energy (69)D≳ħv2kBT,where v is some characteristic velocity of the system. As D is inversely proportional to the resistivity, systems saturating the above bound would display linear resistivity behavior. (See [103] for a recent successful holographic description of linear resistivity at high temperature.)
One should think of (69) as a reformulation of the Kovtun-Son-Starinets (KSS) bound [104] (70)ηs⩾14πħkB,which also relies on the idea of a fundamental dissipative timescale τL~ħ/(kBT) controlling transport in strongly interacting systems. Naively, the observed violations of the KSS bound would seem to indicate the existence of systems in which the bound (68) is violated. The bound (69) saves the idea of a fundamental dissipative timescale by introducing an additional parameter in the game, namely, the characteristic velocity v. The fact that η/s can be made arbitrarily small in some systems corresponds to the fact that the characteristic velocity is highly suppressed in those cases. See [91] for further details.
In [91, 92] Blake proposed that, at least for holographic systems with particle-hole symmetry, the characteristic velocity v should be replaced by the butterfly velocity. More precisely (71)Dc⩾CcvB2τL,where Dc is the electric diffusivity and Cc is a constant that depends on the universality class of theory. This proposal was motivated by the fact that both Dc and vB are determined by the dynamics close to the black hole horizon in the aforementioned systems. Despite working well for systems where energy and charge diffuse independently, this proposal was shown to fail in more general cases [93, 105–108]. This is related to the fact that, in more general cases, the diffusion of energy and charge is coupled, and the corresponding transport coefficients are not given only in terms of the geometry close to the black hole horizon. Hence, there is no reason for these coefficients to be related to the butterfly velocity, which is always determined solely by the near-horizon geometry.
There is, however, a universal piece of the diffusivity matrix that can be related to the chaos parameters at infrared fixed points. This is the thermal diffusion constant [94](72)DT⩾CTvB2τL,where CT is a universality constant (different from Cc). This proposal was shown to be valid even for systems with spatial anisotropy [109]. The above relation is not well defined when the system’s dynamical critical exponent z is equal to one, but it can be extended in this case (we thank Hyun-Sik Jeong for calling our attention to this) [110].
Finally, we mention that there is an interesting relation between chaos and hydrodynamics that manifests itself in the so-called ‘pole-skipping’ phenomenon. See [95–97] for further details.
5. Closing Remarks
The holographic description of quantum chaos not only has provided new insights into the inner-workings of gauge-gravity duality, but also has given insights outside the scope of holography: some examples include the characterization of chaos with OTOCs, the definition of a quantum Lyapunov exponent, and the existence of a bound for chaos.
The success of this new approach to quantum chaos explains the growing experimental interest that OTOCs have been received. Indeed, several protocols for measuring OTOCs have been proposed, and there are already a few experimental results. See [111] and references therein.
Finally, one of the remarkable features of quantum chaos is level statistics described by random matrices. The fact that this is present in the infrared limit of the SYK model [112–114] suggests that it should also be present in quantum black holes (we thank A. M. García-García for calling our attention to this), although this has not yet been verified [115].
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
It is a pleasure to thank A. M. García-García, S. Nicolis, S. A. H. Mansoori, and N. Garcia-Mata for useful correspondence. I also would like to thank Hyun-Sik Jeong for useful discussions. This work was supported in part by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF2017R1A2B4004810) and GIST Research Institute (GRI) grant funded by the GIST in 2018.
UllmoD.TomsovicS.2014http://www.lptms.u-psud.fr/membres/ullmo/Articles/eolss-ullmo-tomsovic.pdfMaldacenaJ.The large N limit of superconformal field theories and supergravity199822231252[International Journal of Theoretical Physics, vol. 38, article 1113, 1999]10.4310/ATMP.1998.v2.n2.a1MR1633016GubserS. S.KlebanovI. R.PolyakovA. M.Gauge theory correlators from non-critical string theory19984281-210511410.1016/S0370-2693(98)00377-3MR16307662-s2.0-0347811692WittenE.Anti de Sitter space and holography199822253291MR1633012Zbl0914.5304810.4310/ATMP.1998.v2.n2.a2ShenkerS. H.StanfordD.Black holes and the butterfly effect2014140306710.1007/JHEP03(2014)067MR3190979ShenkerS. H.StanfordD.Multiple shocks2014141204610.1007/JHEP12(2014)046Zbl06566194RobertsD. A.StanfordD.SusskindL.Localized shocks201515035110.1007/JHEP03(2015)051MR3340732Zbl1388.83694ShenkerS. H.StanfordD.Stringy effects in scrambling2015150513210.1007/JHEP05(2015)132MR3358871KitaevA.Hidden correlations in the hawking radiation and thermal noise, talk given at Fundamental Physics Prize SymposiumProceedings of the Stanford SITP seminars2014RozenbaumE. B.GaneshanS.GalitskiV.Lyapunov exponent and out-of-time-ordered correlator’s growth rate in a chaotic system201718808680110.1103/PhysRevLett.118.086801RozenbaumE. B.GaneshanS.GalitskiV.Universal level statistics of the out-of-time-ordered operatorhttps://arxiv.org/abs/1801.10591Chávez-CarlosJ.López-Del-CarpioB.Bastarrachea-MagnaniM. A.StránskýP.Lerma-HernándezS.SantosL. F.HirschJ. G.Quantum and classical lyapunov exponents in atom-field interaction systems2019122024101PolchinskiJ.RosenhausV.The spectrum in the Sachdev-Ye-Kitaev model20161604110.1007/JHEP04(2016)001MR3512315MaldacenaJ.StanfordD.Remarks on the Sachdev-Ye-Kitaev model2016941010600210.1103/physrevd.94.106002MR3751002MaldacenaJ.StanfordD.YangZ.Conformal symmetry and its breaking in two-dimensional nearly anti-de Sitter space201612C1041210.1093/ptep/ptw124MR3628687KitaevA.SuhS. J.The soft mode in the Sachdev-Ye-Kitaev model and its gravity dual2018180518310.1007/jhep05(2018)183MR3815004AxenidesM.FloratosE.NicolisS.Modular discretization of the AdS2/CFT1 holography2014140210910.1007/JHEP02(2014)109AxenidesM.FloratosE.NicolisS.Chaotic information processing by extremal black holes2015249154201210.1142/S0218271815420122MR3375705AxenidesM.FloratosE.NicolisS.The quantum cat map on the modular discretization of extremal black hole horizons201878541210.1140/epjc/s10052-018-5850-9RobertsD. A.StanfordD.Two-dimensional conformal field theory and the butterfly effect201511513131603StanfordD.Many-body chaos at weak coupling2016161000910.1007/JHEP10(2016)009MR3577478PlamadealaE.FradkinE.Scrambling in the quantum Lifshitz model20181806606310210.1088/1742-5468/aac136NahumA.VijayS.HaahJ.Operator spreading in random unitary circuits20188202101410.1103/PhysRevX.8.021014KhemaniV.VishwanathA.HuseD. A.Operator spreading and the emergence of dissipative hydrodynamics under unitary evolution with conservation laws20188303105710.1103/PhysRevX.8.031057RakovszkyT.PollmannF.von KeyserlingkC.Diffusive hydrodynamics of out-of-time-ordered correlators with charge conservation20188303105810.1103/PhysRevX.8.031058LuitzD. J.Bar LevY.Information propagation in isolated quantum systems201796202040610.1103/PhysRevB.96.020406BohrdtA.MendlC. B.EndresM.KnapM.Scrambling and thermalization in a diffusive quantum many-body system201719606300110.1088/1367-2630/aa719bHeylM.PollmannF.DóraB.Detecting equilibrium and dynamical quantum phase transitions in ising chains via out-of-time-ordered correlators2018121101680110.1103/PhysRevLett.121.016801LinC.MotrunichO. I.Out-of-time-ordered correlators in a quantum Ising chain20189714430410.1103/PhysRevB.97.144304XuS.SwingleB.Accessing scrambling using matrix product operatorshttps://arxiv.org/abs/1802.00801CenciniM.CecconiF.VulpianiA.2009SingaporeWorld ScientificMR2567372PolchinskiJ.Chaos in the black hole S-matrixhttps://arxiv.org/abs/1505.08108García-MataI.SaracenoM.JalabertR. A.RoncagliaA. J.WisniackiD. A.Chaos signatures in the short and long time behavior of the out-of-time ordered correlator20181212121060110.1103/PhysRevLett.121.210601StanfordD.Many-body quantum chaosProceedings of the Seminar at the school, IAS PiTP 2018: From Qubtis to Spacetime2018https://video.ias.edu/PiTP/2018/0723-DouglasStanfordLarkinA. I.OvchinnikovY. N.Quasiclassical method in the theory of superconductivity196928612001205KhemaniV.HuseD. A.NahumA.Velocity-dependent Lyapunov exponents in many-body quantum, semiclassical, and classical chaos2018981414430410.1103/PhysRevB.98.144304RobertsD. A.SwingleB.Lieb-robinson bound and the butterfly effect in quantum field theories2016117909160210.1103/PhysRevLett.117.091602MaldacenaJ.Eternal black holes in anti-de Sitter2003030402110.1088/1126-6708/2003/04/021MR1989847FidkowskiL.HubenyV.KlebanM.ShenkerS.The black hole singularity in AdS/CFT2004040201410.1088/1126-6708/2004/02/014MR2046795DrayT.'t HooftG.The gravitational shock wave of a massless particle1985253117318810.1016/0550-3213(85)90525-5MR789737SfetsosK.On gravitational shock waves in curved spacetimes1995436372174510.1016/0550-3213(94)00573-WMR1316141FischlerW.JahnkeV.PedrazaJ. F.Chaos and entanglement spreading in a non-commutative gauge theory2018181107210.1007/JHEP11(2018)072Zbl06995335BaggioliM.PadhiB.PhillipsP. W.SettyC.Conjecture on the butterfly velocity across a quantum phase transition20181807049MR3862004PerlmutterE.Bounding the space of holographic CFTs with chaos2016161006910.1007/JHEP10(2016)069MR3577538FitzpatrickA. L.KaplanJ.A quantum correction to chaos2016160507010.1007/JHEP05(2016)070MR3521789CaputaP.NumasawaT.Veliz-OsorioA.Out-of-time-ordered correlators and purity in rational conformal field theories2016201611113B0610.1093/ptep/ptw157MR3628681TuriaciG. J.VerlindeH.On CFT and quantum chaos2016161211010.1007/JHEP12(2016)110MR3602310CaputaP.KusukiY.TakayanagiT.WatanabeK.Out-of-time-ordered correlators in a T2n/Zn CFT2017964046020MR3852975TaylorM.Generalized conformal structure, dilation gravity and SYK2018180101010.1007/jhep01(2018)010MR3764265GrozdanovS.On the connection between hydrodynamics and quantum chaos in holographic theories with stringy corrections20192019110.1007/JHEP01(2019)048CamanhoX. O.EdelsteinJ. D.MaldacenaJ.ZhiboedovA.Causality constraints on corrections to the graviton three-point coupling2016160202010.1007/JHEP02(2016)020Zbl1388.83093HaydenP.PreskillJ.Black holes as mirrors: quantum information in random subsystems20070709912010.1088/1126-6708/2007/09/120MR2342311SekinoY.SusskindL.Fast scramblers200808101006506510.1088/1126-6708/2008/10/065MaldacenaJ.ShenkerS. H.StanfordD.A bound on chaos2016160810610.1007/JHEP08(2016)106MR3564142MichelB.PolchinskiJ.RosenhausV.SuhS. J.Four-point function in the IOP matrix model2016160504810.1007/JHEP05(2016)048MR3521767PadmanabhanP.ReyS.-J.TeixeiraD.TrancanelliD.Supersymmetric many-body systems from partial symmetries---integrability, localization and scrambling2017170513610.1007/JHEP05(2017)136MR3664112de BoerJ.LlabrésE.PedrazaJ. F.VeghD.Chaotic strings in AdS/CFT20181202020160410.1103/PhysRevLett.120.201604BanerjeeA.KunduA.PoojaryR. R.Strings, branes, schwarzian action and maximal chaoshttps://arxiv.org/abs/1809.02090MezeiM.StanfordD.On entanglement spreading in chaotic systems2017170506510.1007/JHEP05(2017)065MR3662866QiX. L.YangZ.Butterfly velocity and bulk causal structurehttps://arxiv.org/abs/1705.01728MezeiM.On entanglement spreading from holography2017170504610.1007/JHEP05(2017)064MR3662867GiataganasD.GürsoyU.PedrazaJ. F.Strongly coupled anisotropic gauge theories and holography20181211210.1103/PhysRevLett.121.121601JahnkeV.Delocalizing entanglement of anisotropic black branes20181801102MR3762544ÁvilaD.JahnkeV.PatiñoL.Chaos, diffusivity, and spreading of entanglement in magnetic branes, and the strengthening of the internal interaction2018180913110.1007/JHEP09(2018)131Zbl1398.81104CamanhoX. O.EdelsteinJ. D.Causality constraints in AdS/CFT from conformal collider physics and Gauss-Bonnet gravity2010100400710.1007/JHEP06(2010)099BuchelA.EscobedoJ.MyersR. C.PaulosM. F.SinhaA.SmolkinM.Holographic GB gravity in arbitrary dimensions20101003111BriganteM.LiuH.MyersR. C.ShenkerS.YaidaS.Viscosity bound violation in higher derivative gravity2008771210.1103/PhysRevD.77.1260061260062-s2.0-45249087420BriganteM.LiuH.MyersR. C.ShenkerS.YaidaS.Viscosity bound and causality violation200810019160110.1103/PhysRevLett.100.191601CamanhoX. O.EdelsteinJ. D.PaulosM. F.Lovelock theories, holography and the fate of the viscosity bound2011110512710.1007/JHEP05(2011)127MR2834530ErdmengerJ.KernerP.ZellerH.Non-universal shear viscosity from Einstein gravity2011699430130410.1016/j.physletb.2011.04.009RebhanA.SteinederD.Violation of the holographic viscosity bound in a strongly coupled anisotropic plasma201210802160110.1103/PhysRevLett.108.021601JahnkeV.MisobuchiA. S.TrancanelliD.Holographic renormalization and anisotropic black branes in higher curvature gravity2015150112210.1007/JHEP01(2015)122WolfM. M.VerstraeteF.HastingsM.CiracJ. I.Area laws in quantum systems: mutual information and correlations2008100707050210.1103/PhysRevLett.100.070502MR2385818RyuS.TakayanagiT.Holographic derivation of entanglement entropy from AdS/CFT20069618160210.1088/1126-6708/2006/08/045MR2249925HubenyV. E.RangamaniM.TakayanagiT.A covariant holographic entanglement entropy proposal2007070706210.1088/1126-6708/2007/07/062LeichenauerS.Disrupting entanglement of black holes201490404600910.1103/PhysRevD.90.046009CalabreseP.CardyJ.Evolution of entanglement entropy in one-dimensional systems20050504P0401010.1088/1742-5468/2005/04/p04010MR2141033Abajo-ArrastiaJ.AparícioJ.LópezE.Holographic evolution of entanglement entropy2010101114910.1007/JHEP11(2010)149AlbashT.JohnsonC. V.Evolution of holographic entanglement entropy after thermal and electromagnetic quenches201113404501710.1088/1367-2630/13/4/045017HartmanT.MaldacenaJ.Time evolution of entanglement entropy from black hole interiors2013130501410.1007/JHEP05(2013)014MR3080927LiuH.SuhS. J.Entanglement tsunami: Universal scaling in holographic thermalization201411201160110.1103/PhysRevLett.112.0116012-s2.0-84892378290LiuH.SuhS. J.Entanglement growth during thermalization in holographic systems20148906601210.1103/PhysRevD.89.066012KunduS.PedrazaJ. F.Spread of entanglement for small subsystems in holographic CFTs201795808600810.1103/physrevd.95.086008MR3813492LokhandeS. F.OlingG. W. J.PedrazaJ. F.Linear response of entanglement entropy from holography2017171010410.1007/jhep10(2017)104MR3731109CasiniH.LiuH.MezeiM.Spread of entanglement and causality2016160707710.1007/JHEP07(2016)077MR3547353HartmanT.Afkhami-JeddiN.Speed limits for entanglementhttps://arxiv.org/abs/1512.02695Sabella-GarnierP.Time dependence of entanglement entropy on the fuzzy sphere2017170812110.1007/jhep08(2017)121MR3697361NakagawaY. O.SárosiG.UgajinT.Chaos and relative entropy20181807002MR3862051MagánJ. M.Black holes, complexity and quantum chaos20181809043MR3871256MansooriS. A. H.QaemmaqamiM. M.Complexity growth, butterfly velocity and black hole thermodynamicshttps://arxiv.org/abs/1711.09749BlakeM.Universal charge diffusion and the butterfly effect in holographic theories2016117909160110.1103/PhysRevLett.117.091601BlakeM.Universal diffusion in incoherent black holes201694808601410.1103/physrevd.94.086014MR3737175DavisonR. A.FuW.GeorgesA.GuY.JensenK.SachdevS.Thermoelectric transport in disordered metals without quasiparticles: The Sachdev-Ye-Kitaev models and holography20179515155131BlakeM.DavisonR. A.SachdevS.Thermal diffusivity and chaos in metals without quasiparticles20179610106008MR3867816GrozdanovS.SchalmK.ScopellitiV.Black hole scrambling from hydrodynamics20181202323160110.1103/PhysRevLett.120.231601BlakeM.LeeH.LiuH.A quantum hydrodynamical description for scrambling and many-body chaos201820181010.1007/JHEP10(2018)127Zbl1402.81163GrozdanovS.SchalmK.ScopellitiV.Kinetic theory for classical and quantum many-body chaos201999110.1103/PhysRevE.99.012206BlakeM.DavisonR. A.GrozdanovS.LiuH.Many-body chaos and energy dynamics in holography2018181003510.1007/JHEP10(2018)035Zbl1402.83053HaehlF. M.RozaliM.Effective field theory for chaotic CFTs2018181011810.1007/JHEP10(2018)118Zbl1402.81223DamleK.SachdevS.Nonzero-temperature transport near quantum critical points199756148714873310.1103/PhysRevB.56.8714SachdevS.1999Cambridge, UKCambridge University PressMR1785382HartnollS. A.Theory of universal incoherent metallic transport2015111546110.1038/nphys3174JeongH.NiuC.KimK.Linear-T resistivity at high temperature2018181019110.1007/JHEP10(2018)191KovtunP. K.SonD. T.StarinetsA. O.Viscosity in strongly interacting quantum field theories from black hole physics2005941110.1103/PhysRevLett.94.1116011116012-s2.0-18044397667LucasA.SteinbergJ.Charge diffusion and the butterfly effect in striped holographic matter2016161014310.1007/JHEP10(2016)143MR3578568Zbl1390.83164BaggioliM.GoutérauxB.KiritsisE.LiW.-J.Higher derivative corrections to incoherent metallic transport in holography2017170317010.1007/JHEP03(2017)170MR3650647NiuC.KimK.-Y.Diffusion and butterfly velocity at finite density2017170603010.1007/JHEP06(2017)030MR3678391MokhtariA.MansooriS. A. H.FadafanK. B.Diffusivities bounds in the presence of Weyl corrections201878559160410.1016/j.physletb.2018.09.020JeongH.-S.AhnY.AhnD.NiuC.LiW.-J.KimK.-Y.Thermal diffusivity and butterfly velocity in anisotropic Q-lattice models2018135MR3762506DavisonR. A.GentleS. A.GoutérauxB.mpact of irrelevant deformations on thermodynamics and transport in holographic quantum critical stateshttps://arxiv.org/abs/1812.11060SwingleB.Unscrambling the physics of out-of-time-order correlators201814109889902-s2.0-8505431870210.1038/s41567-018-0295-5CotlerJ. S.Gur-AriG.HanadaM.PolchinskiJ.SaadP.ShenkerS. H.StanfordD.StreicherA.TezukaM.Black holes and random matrices20171705118Erratum: [Journal of High Energy Physics, vol. 1809, article 002, 2018]García-GarcíaA. M.VerbaarschotJ. J. M.Spectral and thermodynamic properties of the Sachdev-Ye-Kitaev model2016941212601010.1007/jhep04(2018)146MR3804351García-GarcíaA. M.VerbaarschotJ. J. M.Analytical spectral density of the Sachdev-Ye-Kitaev Model at finite N2017966066012SaadP.ShenkerS. H.StanfordD.A semiclassical ramp in SYK and in gravityhttps://arxiv.org/abs/1806.06840