Switchbacks and the Bridge to Nowhere

Leonard Susskind, Ying Zhao

Introduction

In reference evidence was given for the duality between two time-dependent black hole phenomena:

the growth of the volume of an Einstein-Rosen bridge (ERB) for a two-sided ADS black hole.

the growth of computational complexity of the two-sided state of the dual gauge theories.

Although it has not been widely discussed, the concept of an ERB applies equally well to one-sided black holes as to the two-sided case. Furthermore the duality between ERB-volume and computational complexity may also be formulated for the one-sided case.

The plan of this paper is to first generalize the concept of an ERB to the case of a one-sided black hole in a pure state. At first sight such a generalization seems contradictory; ERBs are ordinarily considered to be the geometric connection formed between entangled systems according to the ER=EPR conjecture. Because a black hole in a pure state is not entangled with anything else it would seem that there is no such thing as a one-sided ERB. Nevertheless we will see that there is a very definite sense in which the one-sided black hole is attached to a “bridge-to-nowhere.” Moreover, as in the two-sided case, the bridge-to-nowhere grows with time, tracking the growth of complexity.

The second direction that we explore has to do with the definition of complexity. Modeling black hole evolution by random quantum circuits has been a useful guide for some purposes. But black holes are not random circuits; they evolve by Hamiltonian evolution. The definition of complexity in terms of quantum gates may be ideal for computer science, but not for black holes.

In an attempt to be quantitative about gate complexity, Nielsen and collaborators constructed a continuum approximation to gate complexity which involves a new kind of “complexity geometry.” Thus far the Nielsen geometrization of complexity has not proved useful for quantum-computer science, but it is just what we need for black holes. We will illustrate the value of the geometric approach by applying it to an interesting aspect of black hole evolution: the switchback effect associated with precursors and shockwave geometries . The switchback effect is also closely related to fast-scrambling .

The methods we describe apply to both one and two sided cases, but for definiteness we assume a single copy of a gauge theory and a one-sided black hole.

The Bridge to Nowhere

According to the ER=EPR hypothesis of an ERB is the wormhole that connects a black hole to whatever it is entangled with. The most studied case is the two-sided ADS black hole (really a pair of black holes), dual to a pair of entangled conformal field theories. When viewed from the outside the black holes on either side are stationary, but the ERB is strongly time-dependent. Its volume increases with time, and appears to track the growing complexity of the conformal field theory state.

As we will explain, a one-sided black hole although unentangled, is attached to a growing ‘bridge-to-nowhere. To minimize new terminology we’ll use ERB1 and ERB2 to indicate one and two-sided ERBs.

To define an ERB1 we begin with a conformal field theory on a spatial sphere cross time—S(D−1)×R.S(D-1)\times R. Initially the system is in the ground state dual to the ADS vacuum. At t=0t=0 a source is activated which injects enough energy to create a stable black hole. The most tractable situation is when the matter comes in as a spherically symmetric light-like shell. The Penrose diagram in figure 1 shows the process.

The ERB1 is defined by analogy with the two sided construction of . Pick a time tt on the boundary. This defines a boundary (D−2)(D-2)-sphere. Then consider the family of all (D−1)(D-1)-dimensional volumes (space-like hypersurfaces which end on that sphere). The surface of maximal volume defines the ERB1 at time t.t. In figure 2 the evolution of the ERB1 is shown as a series of blue spacelike curves which intersect the boundary at various times.

The equations governing the maximal surface are the same as for the two-sided black hole and the results are similar. The details are in the Appendix.

The spatial volume of these maximal surfaces is infinite for the obvious reason that the metric diverges at the ADS boundary. That divergence can easily be regulated, either by cutting off the integrals at some radial regulator coordinate or by subtracting the vacuum ADS value. Another way to accomplish the same thing is to only count the volume behind the horizon.

As in the ERB2 case, as t→∞t\to\infty the ERB1 asymptotically approaches a fixed hypersurface of infinite regulated volume, show as the green curve in figure 2. For finite tt it grows by a process of stretching which takes place just inside the horizon. Figure 3 shows embedding diagrams for both the two and one-sided ERB’s at some particular time.

A worthwhile point to note is that if there are two unentangled black holes then the embedding diagram is just two copies of the left side of figure 3 with no internal connection. The right side of the figure is an example of the slogan “Entanglement provides the hooks that hold space together.” In this case entanglement stitches the two one-sided ERBs together to form a two-sided ERB.

For t<0t<0 the ERB1 is merely a constant time slice through ADS. In general the portion of the ERB1 inside the shell is part of such a constant time slice, but after it passes through the shell it becomes modified. Once the black hole has formed the volume of the ERB1 grows linearly with time,

with AerbA_{erb} being the cross-sectional area. Within a numerical factor discussed in AA is the area of the horizon.

The computational complexity of the state also grows linearly with a rate equal to the entropy SS times the temperature TT of the black hole. Calling the complexity C,{\cal{C}},

The example that was discussed at length in was the minimal stable black hole of Schwarzschild radius equal to the ADS radius of curvature l.l.

From here on the term black hole will refer to one satisfying 2.7 and 2.8Larger black hole in ADS will be discussed in a upcoming paper with Dan Roberts and Douglas Stanford .

The Geometry of Complexity

Classically the growth of the ERB (either one or two-sided) takes place forever. Quantum mechanically the global classical geometry must break down by the quantum recurrence time, and probably by the classical recurrence timeThe quantum recurrence time is doubly exponential in the entropy while the classical recurrence time is given by a single exponential.. These are extremely long times; much longer than any ordinary thermalization time or evaporation time. The question is what quantity in the dual gauge theory continues to grow for such long periods? The answer is quantum computational complexity. In this section we will review some geometric ideas about complexity which are due to Michael Nielsen and collaborators , and then apply them to fast scrambling, and to the switchback effect of .

A quantum circuit (QC) is a device composed of qubits and gates (called gg) whose purpose is to implement specific unitary transformations on an initial state of the qubits. Let the number of qubits be K.K. The gates are one or two-qubit special-unitary transformationsBy special-unitary we mean with determinant 1.1.. The gates can either occur in series—in each time step a single gate acts—or in parallel as explained by Hayden and Preskill . In the Hayden-Preskill parallel case the qubits are paired at each step, and K/2K/2 gates act simultaneously on the pairs. For our purposes in this section the distinction is not important and it is easier to describe the series case.

Given an initial state ∣Ψ(0)⟩|\Psi(0)\rangle after nn time-steps the state has been transformed to ∣Ψ(n)⟩,|\Psi(n)\rangle,

Here u(n)u(n) is an element of SU(2K).SU(2^{K}). Let us think of 3.3 as defining a path in the space SU(2K).SU(2^{K}). This is schematically shown in the left side of figure 4.

The path begins at u(0)=Iu(0)=I and ends at u(n).u(n). The rule for such paths is that every link corresponds to a gate and therefore displaces the endpoint by a one or two qubit operator.

With these concepts in hand we can define the complexity of a unitary operator uu as the smallest number of gates of any circuit that can yield uu as an outcome. That is to say, it is the number of links of the shortest allowable path connecting II and u.u.

In the past, random quantum circuits have been used to model black hole evolution but our real interest is in continuous Hamiltonian evolution. Part of the reason for this paper is to draw attention to an innovation of Nielsen and collaborators who introduced a continuum description of complexity. Their purpose was to construct an approximation to a quantum circuit that used Hamiltonian evolution and Riemannian geometry. However, the methods of seems well suited to the study of Hamiltonian systems of the kind that may represent black hole evolution.

What the authors of do is to introduce a Riemannian geometry on the space SU(2K).SU(2^{K}). The discrete time steps are replaced by infinitesimal transformations and the discrete path through SU(2K)SU(2^{K}) by a continuous path. This is shown on the right side of figure 4.

The path may be parameterized by a generator hh that acts as a possibly path-dependent HamiltonianThe notation hh is used to represent any generator composed out of the ZAZ_{A} basis of operators. The notation HH is reserved for the actual Hamiltonian governing the black hole evolution. Similarly the notation u(s)u(s) is used for a generic trajectory on SU(2K)SU(2^{K}) and UU is used for the actual evolution induced by H.H..

where PP indicates path-ordered product. Here ss is a path parameter that we’ll think of as time, and h(s)h(s) is a time dependent Hamiltonian. The construction yields a path through SU(2K)SU(2^{K}) where at every point the tangent to the path is determined by h(s).h(s). It is natural in this setting to define the complexity of uu to be the length of the shortest geodesic connecting II and u.u.

The only question is how to choose the metric on SU(2K).SU(2^{K}). There is a conventional metric on unitary groups that is invariant under conjugation u→ vuv†,u\to\ vuv^{{\dagger}}, but it does not have the right properties for defining complexity. What one wants is a right-invariant metric that incorporates a rule that all gates should be one and two qubit operators. The idea is to choose the geometry so that a high price (in length) is paid for all infinitesimal motions that are not generated by one and two qubit Hamiltonians. More precisely a large price should be paid for Hamiltonians which are not sums of one and two qubit terms.

Let us begin with the algebra of all Hermitian operators constructed from KK qubits. Each qubit has four operators consisting of three traceless Pauli operators and the unit operator. We label the four operators with an index ii running over four values. Let Xi(n)X^{i}(n) be the four operators for the nthn^{th} qubit (i=1,..,4i=1,..,4). The algebra is spanned by a basis consisting of products of KK factors, one XX for each qubit. There are 4K4^{K} basis operators of this type. If we want to consider only traceless operators then there the basis contains 4K−14^{K}-1 operators that we will indicate by ZA.Z_{A}. The index AA runs over all 4K−14^{K}-1 orthogonal directions in the algebra.

Let us consider a point u0u_{0} in SU(2K).SU(2^{K}). Consider the displacements

defined by left multiplication by Z.Z. The conventional SU(2K)SU(2^{K}) metric is equivalent to the following inner product structure.

It is invariant with respect to right multiplication by unitary matrices and defines a homogeneous geometry. It also has a discrete form of isotropy. Any two basis generators ZAZ_{A} and ZBZ_{B} can be related by conjugation by an SU(2K)SU(2^{K}) operation. However, in general normalized superpositions of the basis generators are not equivalent. An easy way to see that is to note that any ZAZ_{A} satisfies

but superpositions of the generators do not satisfy the second equation for any normalization.

However this metric is not quite what we want. We would like a more anisotropic geometry in which displacements are very costly for directions ZAZ_{A} corresponding to three or more qubit operators. To that end introduces projection operators P{\cal{P}} and Q{\cal{Q}} and a cost-parameter q>>1.q>>1. The projector P{\cal{P}} projects onto the subspace of one and two qubit ZZ-operators. Those directions are called “easy” directions and no large cost is attached to them. The projector Q{\cal{Q}} projects onto the orthogonal “hard” subspace of three and more qubit operators. The C{\cal{C}}-metric is then given by,

If AA and BB are both easy then the inner product is δAB.\delta_{AB}. If they are both hard the inner product is qδAB.q\delta_{AB}. Otherwise the inner product is zero. For q≫1q\gg 1 the metric defined by 3.7 and 3.14 was called the “standard” metric in . We will call it the C{\cal{C}}-metric to avoid confusing it with other more conventional metrics. When qq is large the C{\cal{C}}-metric introduces a large price in metrical distance for curves whose tangents project onto hard directions. The geometry of SU(2K)SU(2^{K}) defined by the C{\cal{C}}-metric will be referred to as C{\cal{C}}-geometry.

C{\cal{C}}-geometry is homogeneous but when qq is large it is extremely anisotropic. Geodesics strongly avoid hard directions. Reference shows that if q>4K,q>4^{K}, then the complexity is insensitive to its precise valueThe limit q→∞q\to\infty was discussed in ..

We may now define the complexity of an element uu of SU(2K).SU(2^{K}). Using the C{\cal{C}}-metric with large qq the complexity is defined by the distance between II and uu along the shortest geodesic connecting them.

There is a very important case which the authors of study in detail. If HH is a fixed easy HamiltonianHere we have in mind the actual Hamiltonian governing the evolution of the the black hole. , i.e., a Hamiltonian composed as a sum of one and two qubit terms, independent of s,s, then the curve

is a geodesic. In other words time-independent easy Hamiltonians generate geodesics. Moreover for ss smaller than some critical value scs_{c} the curve in 3.15 is the shortest geodesic connecting II and u(s).u(s). The length of such a geodesic is proportional to ss which shows that for some length of time the complexity grows linearly with time. This was assumed in but in it is proved.

Beyond the critical time scs_{c} equation 3.15 continues to define a geodesic but it is not the shortest one. The transition to a shorter geodesic occurs at or before the first conjugate point at a so-called cut pointThere may be more than one geodesic connecting two points in a Riemannian space. As one proceeds along the geodesic generated by 3.15, at some point a second geodesic may replace it as the shortest one. Such a point is called a cut point . At a conjugate point a continuous family of geodesics connect two points. If a conjugate point exists then a cut point must occur at or before the conjugate point. . One can see that a cut point must happen for sufficiently large s.s. The quantum recurrence theorem implies that that e−iHs≈Ie^{-iHs}\approx I for s∼e2K.s\sim e^{2^{K}}. Thus for very long times the geodesic 3.15 cannot be the shortest.

This double exponential time is probably overkill. One knows that there is a maximum complexity for circuits. It is of order eKe^{K} beyond which the complexity cannot increase. Therefore after a singly exponential time (classical recurrence time), shorter geodesics must replace 3.15. But there is no general argument that requires that they occur before the classical recurrence time. In previous work we assumed that complexity increases linearly until it reaches is maximal value eK.e^{K}. It is important to know whether this is true for the kinds of Hamiltonians which can describe black holes.

An important fact noted in is that most sectional curvatures of the C{\cal{C}}-geometry are negative. This means that neighboring geodesics tend to deviate exponentially. The qualitative reason is not too hard to understand; most classical Hamiltonians are chaotic which means that trajectories have a tendency to exponentially separate. This is also a property of geodesics in a negatively curved space. We will return to this in the next section on switchbacks.

2 State Complexity

Thus far we have discussed the complexity of unitary operators. One can also define complexity for states. To do so we need to first define simplicity. For a system of qubits that’s fairly easy. Intuitively a state is simple if it not entangled. In other words product states are simple. The simplest of all product states is one in which all the qubits are all in the same state. For example the state

An arbitrary normalized state ∣Ψ⟩|\Psi\rangle may be identified with a point on the homogeneous manifold CP(2K−1).CP(2^{K}-1).

The complexity of a general state ∣Ψ⟩|\Psi\rangle can be defined to be the complexity of the least complex unitary that can act on ∣0⟩|0\rangle to give ∣Ψ⟩.|\Psi\rangle.

Let’s see if we can define a metric on CP(2K−1).CP(2^{K}-1). Consider two neighboring points, i.e., normalized vectors in the qubit Hilbert space modulo overall phase. Call them ∣f⟩|f\rangle and ∣g⟩.|g\rangle. It is assumed that the difference is infinitesimal. There exist unitary transformations which connect the two states. The unitary operators are not at all unique but among them there is a unitary of minimal complexity. This means that it is the shortest distance from the identity using the complexity metric. We take the complexity of that unitary to define the distance between the two points f, gf,\ g on CP(2K−1).CP(2^{K}-1).

It is not hard to prove that the defining properties of a metric are satisfied:

We can now define the complexity of any state to be the minimum geodesic distance from the simple state to the state in question.

The choice of simple state is not unique. One could choose any unentangled state instead of

However all such states are close because one can go from one to another with at most KK gates. In metric terms they are all within distance KK from one another. Distance KK sounds large but in the context of the parallel computing model of Hayden and Preskill it corresponds to a single time-step.

The complexity of quantum states can be vastly larger than the complexity of classical states. As a model for a classical system let’s consider KK c-bits. Each c-bit is like a coin, it can be in one of two states—heads or tails. We will add the rule that states are identified under the global Z2Z_{2} that flips all coins.

The state (hhhhh....hhh)=(ttttt...ttt)(hhhhh....hhh)=(ttttt...ttt) is a simple state; you can describe it in two words, all equal. Instead of quantum gates we can allow configurations to be updated by simple single-coin flips—heads goes to tails; tails goes to heads. The complexity of a configuration like (htthtthhhththhhtt...)(htthtthhhththhhtt...) may be defined as the minimum number of flips needed to get to it from the simple state. Here are two facts:

The maximum complexity of any classical state is K/2.K/2. This is obvious; any configuration can be achieved by just flipping the initial configuration to the desired final configuration—one coin at a time.

The second fact is almost as obvious. It says that when KK is large almost all states have complexity close to K/2.K/2.

To put it another way, if one imagines cycling though the states by some rule that flips one coin at a time (the rule for which coin is flipped may depend on the configuration of the other coins) in a reversible manner, then it takes a time of order K/2K/2 to reach a maximally complex state. The complexity will hover around K/2K/2 for a long time but eventually come back to simplicity in a time 2K.2^{K}.

Something like that is true for quantum complexity except the magnitudes are much larger. First of all, within a given accuracy any state can be achieved using ∼eK\sim e^{K} gates. Secondly, almost all states have complexity of order eK.e^{K}. Finally, it typically takes an exponential time to get to most states from the simple state and a doubly exponential time to get back to a simple state. The set of states that you can get to in polynomial time is very special and not at all generic.

Now that we have an understanding of complexity we can state the conjecture relating complexity to the geometry of the black hole interior.

For a black hole of ADS radius and entropy KK the volume of the ERB1 is proportional to the complexity of the quantum state.

A more general formula can be based on the following hypothesis. The rate of growth of complexity (long before the recurrence time) should reflect the number of qubits and the energy per qubit. This translates to the rate

where SS and TT are the entropy and temperature of the system. Writing S=A/4GS=A/4G one finds

To specialize this general formula to the case R=R= we recall that the temperature in that case in T=1/lT=1/l giving 3.24.

This relation has been confirmed for the TFD state and a wide variety of shockwave deformations of the TFD.

The Geometry of Precursors and Switchbacks

The section makes heavy use of ideas due to Dan Roberts and Douglas Stanford which will appear shortly .

Precusors are useful probes that have been used to test the connection between complexity and ERB . In this section we will use the geometric complexity ideas of to analyze them. Precursors, and switchbacks are related concepts that apply to products of operators taken out of time order. A precursor is an operator of the form

where U(t)U(t) is the time evolution operator for an arbitrary system. In practice this will mean black hole satisfying 2.7 and 2.8. The reason they are interesting is that their complexity evolves in a very characteristic manner. Moreover when a precursor acts on a black hole, the evolution of precursor complexity closely parallels the evolution of the ERB. Our goal in this section is to relate the properties of precursors to properties of the C{\cal{C}}-geometry. We will then have three distinct ways to evaluate the complexity of precursors: two old and one new. The old ways are: using Einstein field equations to calculate the response of the ERB geometry to the action of a precursor ; and using random quantum circuits to evolve the precursor. The new way is to study the geodesics of C{\cal{C}}-geometry. We already know that the old ways agree with each other . With some simple assumptions about the C{\cal{C}}-metric we will see that the new way agrees with the old.

Two quantities will occupy our attention; the first is the complexity of the precursor labeled CW(t).{\cal{C}}_{W}(t). The other is the size of the precursor . Roughly speaking it is the a measure of how strongly the other qubits are causally influenced by the qubit WW after time t.t.

The magnitude of the effect is given by the norm of the commutator,

One way to define the size of the precursor is by,

We can also define a symmetric version of the size by summing over all pairs of qubits,

2 Quantum Circuit Analysis

We are interested in understanding the complexity and size of precursors by relating them to properties of the C{\cal{C}}-geometry, but it is helpful to first get an idea of how they work in the Hayden-Preskill quantum circuit (QC) model. We begin with the size sws_{w}.

In the Hayden-Preskill model, in each time step K/2K/2 gates act. The qubits are randomly paired and a gate acts on each pair. In the next step the pairing is randomized again and another K/2K/2 gates act. The time steps are taken to be thermal times equal to the inverse temperature. For the black holes we are considering the time steps are the ADS time l.l. Thus after time t,t, t/lt/l steps have occurred and Kt/2lKt/2l gates have acted.

At time t=0t=0 all qubit operators for different qubits commute. After one time-step, the first qubit W=q1W=q_{1} has paired with one other qubit (call it q2q_{2}), and has interacted with it via a two-qubit gate. Therefore WW no longer commutes with q2,q_{2}, but it commutes with all others. Let’s say that WW has infected q2.q_{2}. Thus after one step there are two infected qubits.

In the next step the two infected qubits interact with new partners. If we ignore the small probability (of order 1/K1/K) that WW and q2q_{2} interact again with each other, then four qubits will become infected, meaning that they will not commute with W.W. Moreover the non-zero commutators will be of order unity.

Clearly the epidemic of qubits infected by WW grows exponentially for a time. The size of the precursor—the number of infected qubits—will be denoted by sw.s_{w}. Since an infected qubit has an order one commutator with WW it follows that the size is of order

We can compute how sws_{w} grows in the QC model as follows: Suppose at any given instant the number of infected qubits is ss with s<K.s<K. In the next time step the average number of newly infected qubits will be,

ss is given by the cumulative sum of the Δs\Delta s up to time t.t.

A single step of the quantum circuit corresponds to a time of order Δt∼l.\Delta t\sim l. Smoothing out the time steps gives a differential equation for s(t),s(t),

For s≪Ks\ll K and K≫1K\gg 1 we see that s(t)s(t) grows exponentially with t/l.t/l. But once s∼Ks\sim K the growth saturates. The solution to 4.10 is,

Notice that the saturation takes place at about the scrambling time. Beyond that ss tends quickly to its asymptotic value s→K.s\to K.

The size ss is one term in the double sum 4.7 that defines S(t).{\cal{S}}(t). Therefore we multiply ss by KK to get,

S{\cal{S}} asymptotically tends to the value K2.K^{2}.

Now let us consider the complexity. After time tt the number of gates that has acted to form U(t)U(t) is p=K2tl.p=\frac{K}{2}\frac{t}{l}. We may write,

All together 2p+12p+1 gates have acted (We may consider WW to be a one-qubit gate). But the total number of gates is an over-estimate of the complexity for the simple reason that fewer gates would have gotten us to the same result for W(t).W(t). It is easy to see that only the gates that were involved in the spread of the epidemic were required, all the others having canceled in 4.14. This implies that the complexity is given by the sumWe thank Dan Roberts and Douglas Stanford for explaining this point. of the sizes for all time-steps earlier than t.t. We may write this in differential form,

The complexity can be obtained by integration using 4.11 and 4.15. We find,

where t∗=llog⁡Kt_{\ast}=l\log{K} is the scrambling time.

3 Switchback Effect

at large time 4.19 behaves likeThere is a tricky factor of 22 here. This expression is actually twice the complexity of the operator U(t).U(t). The origin of the factor of 22 is the fact that K/2K/2 gates occur in a unit time step.

The subtraction of Kt∗lK\frac{t_{\ast}}{l} in 4.20 is the switchback effect that was explained in . Let us take a moment to recall its importance. One might have thought that the precursor’s effect on the volume of an ERB would be proportional to the total number of gates that act in the product U(t)WU(−t).U(t)WU(-t). But complexity is not defined by the number of gates in a construction like equation 4.14. It is defined by the minimum number of gates in any quantum circuit that can produce W(t).W(t). Any gate in U(t)U(t) that does not act on an infected qubit might as well not have occurred since it will cancel with a corresponding gate in U(−t).U(-t). Thus these gates should be deleted in the complexity book keeping.

What is remarkable is that exactly the same cancelation occurs when computing the effect of a Shenker-Stanford shockwave on the volume of the ERB. Although the shockwave was created by acting with a precursor, the effect is a General Relativity effect having nothing obvious to do with quantum complexity.

Later we will see that exactly the same effect occurs in the Nielsen geometric way of thinking about complexity.

Finally we remark that for small tt both ss and C{\cal{C}} grow exponentially,

The qualitative behavior of C{\cal{C}} and ss are sketched in figure 5.

4 Complexity-geometry of Precursors

Our goal is to understand how the properties of complexity and size relate to C{\cal{C}}-geometry.

The precursor 4.1 defines a path through the C{\cal{C}}-geometry that begins with the evolution operator U(s)=e−iHsU(s)=e^{-iHs} sweeping out a geodesic connecting the unit operator with U(t).U(t). We will refer to it as a “vertical” geodesicThe terminology vertical and horizontal are just names with no technical significance.. Next, WW acts to add a short segment of unit length in a new “horizontal” direction. Finally U†U^{{\dagger}} acts and continues the curve more or less backward vertically along the direction of U.U.

We will now define a two dimensional sub-manifold of C{\cal{C}}-geometry associated with the precursor: call it a C2{\cal{C}}_{2}-geometry. U(t)U(t) and U(−t)U(-t) define segments of neighboring geodesics connected by a third short segment representing the one-qubit operator W.W. This is depicted in figure 6.

Next, “ horizontal” geodesics are constructed connecting corresponding points along U(±t).U(\pm t). As tt is varied from t=0t=0 to t=∞t=\infty the horizontal geodesics sweep out a surface with boundaries. This surface is the C2{\cal{C}}_{2}-geometry for the precursor 4.1. Its metric is inherited from the original C{\cal{C}}-metric defined in section 3. Note that the length of the horizontal geodesics represent the complexity of the precursor at time t.t.

In principle, given a Hamiltonian and a WW it is possible to construct the metric on the two-dimensional geometry, but this would be very hard. Our goal will more modest; namely we will do a bit of guessing about the metric and then attempt to verify the properties of complexity and scrambling that we indicated earlier.

First consider whether the two-dimensional C2{\cal{C}}_{2}-metric is positively or negatively curved. It it were positively curved the geodesics U(±t)U(\pm t) in figure 6 would tend to meet after a finite time. This is impossible because the geodesic generated by a given Hamiltonian, and passing through a given point, is unique. Moreover, if the Hamiltonian leads to chaotic behavior, as we expect for a black hole, then the trajectories should exponentially diverge. This suggests that the curvature of the C2{\cal{C}}_{2}-metric is negative. Our first guess is that the curvature is negative and constant over the the portion of the C2{\cal{C}}_{2}-geometry of interest. As for the magnitude of the curvature we could leave it as a parameter to be fixed by comparison with 4.3. This is an exercise that we will leave to the reader. The result is that the two-dimensional curvature scalar R2,R_{2}, must be chosen to satisfy,

With the assumption of uniform negative curvature, the C2{\cal{C}}_{2}-geometry may be drawn as a piece of the Poincare plane. This is shown in figure 7.

5 Complexity Calculation

Let us start by calculating the complexity of U(t).U(t). According to our assumptions the complexity should be given by Kt/2l,Kt/2l, and should be equal to the geodesic length of the vertical segment representing U(t).U(t).

The vertical geodesic distance between the top of the rectangle (location of WW) and the point y0y_{0} is given by

This allows us to identify the relation between tt and y0,y_{0},

The naive expectation for the complexity of W(t)W(t) is just the sum of the geodesic lengths of the three blue segments in figure 7. Adding the two vertical lengths and ignoring the small contribution of WW, the naive complexity would be,

But 4.26 is based on a wrong identification of the complexity with distance along the broken curve. According to the arguments of the complexity should be identified as the geodesic distance between the two points at y0.y_{0}. Thus we must find the horizontal geodesic connecting corresponding points on the vertical geodesics. This is easy to do since all geodesics on the Poincare plane are circles centered on the xx axis. The geodesic passing through the two points at y0y_{0} is shown in green. It is given by,

Its length is easily calculated from the metric 4.23,

From 4.32 it is easily seen that CW{\cal{C}}_{W} behaves like

This is very similar to the results from the quantum circuit analysis summarized in figure 5. It is significant evidence that the complexity properties of precursors are well described by a C2{\cal{C}}_{2}-geometry with the simple property of uniform negative curvature.

The shift from tt to (t−llog⁡K)(t-l\log K) in the second equation of 4.35 is the switchback effect of whose importance we discussed earlier. In C{\cal{C}}-geometric context it occurs because complexity is defined by the minimum geodesic distance and not the length of the original path consisting of U(t),U(t), W,W, and U(−t).U(-t).

6 Size of Precursor

Now let us consider the size function sw.s_{w}. We don’t have a first-principles way to connect the size with C{\cal{C}}-geometry, but inspired by equation 4.16 we may write

Equation 4.39 is a remarkably simple result which has the geometric interpretation of the distance along a line of constant yy tangent to the top of the circle x2+y2=L2x^{2}+y^{2}=L^{2} as shown in figure 8. We’ll denote it by the symbol d(t).d(t). It should be understood that the time tt does not refer to the points at which the red segment intersects the vertical geodesics; it refers to the points in the figured labeled y0.y_{0}.

The red segments in figure 8 are not geodesics, and its not clear what is special about them from the geometric viewpoint. We would like to suggest an alternate interpretation. We can get a hint from dimensional analysis, where the dimensions refer to C{\cal{C}}-geometry, not space-time. From the form of the metric in 4.23 we see that the constant KK has units of length. Both complexity and size scale with one power of KK and are therefore also lengths in “complexity units”.

On the other hand the double sum S(t){\cal{S}}(t) in 4.7 scales like K2K^{2} and is therefore an area in complexity units. It is of course just the product sK.sK.

It follows from the properties of the Poincare plane that the area of the column above the red segment, from the red segment to y=∞,y=\infty, is also d(t)K2.d(t)\frac{K}{2}. Let us call that area A′.{\cal{A}}^{\prime}. Equation 4.41 becomes

Still, the lack of a clear geometric meaning to the red line, and the area above it, remain unsatisfying. A more attractive quantity is the area above the green geodesic shown in figure 9. Call it A.{\cal{A}}. How much different are A′{\cal{A}}^{\prime} and A?{\cal{A}}? Not very. They differ by the area of the two small triangular regions between the two curves. One might worry that these areas grow proportionally larger than A′{\cal{A}}^{\prime} when y0y_{0} gets small, but this doesn’t happen. The maximum ratio of A{\cal{A}} to A′{\cal{A}}^{\prime} is,

For times less than the scrambling time the ratio is close to unity.

We find it tempting to speculate that in the realistic Hamiltonian evolution that an appropriately defined version of S{\cal{S}} is accurately represented by the area A.{\cal{A}}.

The calculation of A{\cal{A}} is straightforward. It is given by,

For small time, one finds the area having the form,

For large tt the area has the limiting value,

compared with a′a^{\prime} which has the limiting behavior

We have constructed three functions to represent the size sw(t),s_{w}(t), namely 4.11, 4.40, and 4.47. The first came from analyzing the behavior of random quantum circuits, and the other two from C2{\cal{C}}_{2}-geometry. It is interesting to compare them. We can do this best by first changing the tt variable to u,u,

For the three cases 4.11, 4.40, 4.47 the function FF takes the form,

for 4.47. The factor of 2/π2/\pi was included in F3F_{3} in order to normalize the asymptotic value to 1.1.

The three functions are plotted in figure 10

While the functions are generally similar, 4.55 and 4.53 are extremely close over most of the range, while 4.54 differs substantially. We can see the same thing by plotting the ratios in in figure 11 ( F2/F1F_{2}/F_{1} (red) and F3/F1F_{3}/F_{1} blue).

It is not clear to us why F3F_{3} is a better fit to F1F_{1} than F2.F_{2}.

Questions

There are many questions that this paper leaves unanswered. We will list a few.

When Hawking radiation is accounted for, the bridge-to-nowhere will become a bridge-to-somewhere, the somewhere being the Hawking radiation that purifies the state of the black hole . What are the properties of such a bridge? Can it be described by a quantum circuit and its close relative, a tensor network ?

The trajectory U(t)U(t) swept out by an easy Hamiltonian is a geodesic, but not necessarily the shortest geodesic. For some period of time it will be the shortest, but that cannot last forever. Since complexity is bounded by by exponential in K,K, shorter geodesics must come into play by the classical recurrence time ∼eK.\sim e^{K}. This happens when cut points occur. The occurrence of a cut point may signal the breakdown of the classical general relativistic description of the ERB.

How long does it take for the first cut point to occur? Our guess is that it takes the classical recurrence time but there is no strong evidence for this.

We conjectured a very simple C2{\cal{C}}_{2}-geometry to describe precursor complexity, and it appears to capture the expected properties of precursors rather well. But there was no first-principles derivation of it from the C{\cal{C}}-geometry of . Can we derive its main properties directly from definition of C{\cal{C}}-geometry and the properties of black holes such as fast scrambling? In particular can we demonstrate that the curvature of C2{\cal{C}}_{2} is of order −1/K2?-1/K^{2}?

We restricted our discussion to black holes of radius equal to the ADS radius l.l. How much of the analysis generalizes to larger black holes? This question will partly be answered in a forthcoming paper .

Perhaps the biggest issue is how to understand the relation between general relativity and complexity. GR gives a description of ERBs, and certain aspects of the geometry closely parallel the growth of quantum complexity including the details of precursors and switchbacks. How much more of GR can be found hidden in the study of quantum complexity? We expect much more.

Acknowledgements

We thank Patrick Hayden and Nathaniel Thomas for explaining the work of Dowling and Nielsen. We also thank Michael Nielsen for correcting some misconceptions.

Special thanks goes to Dan Roberts and Douglas Stanford for sharing their insights about the behavior of ss and C{\cal{C}} from the study of quantum circuit models.

Support for this research came through NSF grant Phy-1316699 and the Stanford Institute for Theoretical Physics.

Appendix A Geometry of ERB1

In this appendix we will calculate the functional form of the ERB1 for the one-sided version of a BTZ black hole.

Consider a one-sided AdS black hole formed by gravitational collapse of a spherical shell of null matter. The metric is

The rr-coordinates, which characterizes the proper area around the center, is continuous across the shock wave. We also choose tt-coordinates to be continuous along the boundary r=∞r=\infty as a convention.

This coordinate doesn’t cover the horizon in BTZ part. Passing to null coordinates, let

then the matching condition along the shock wave requires both rr and vv coordinates to be continuous. If the shock wave comes in at boundary time t=0t=0, then it lies at vw=0v_{w}=0.

A.2 Extremal Surface

We look for extremal spacelike surface ending at boundary time tbt_{b}. When tb<0t_{b}<0, the surface is simply a disk in pure AdS with t=t= constant. Its area is

where rbr_{b} is the boundary radius cutoff.

After the shockwave, the surface will contain two parts, one disk inside pure AdS, and another part inside BTZ. We parametrize the surface by a function r=r(v)r=r(v). Then the area element is

where r′r^{\prime} is the derivative of rr with respect to vv. We can consider L=r2r′−fL=r\sqrt{2r^{\prime}-f} as a lagrangian, vv as time, and rr as canonical coordinate. Then the conjugate momentum is

In both the AdS and the BTZ region, the Lagrangian is independent of vv (time), so the corresponding energy is separately conserved

Look at the pure AdS part first. Smoothness at r=0r=0 requires 0=dtdr=dvdr−f0=\frac{\displaystyle dt}{\displaystyle dr}=\frac{\displaystyle dv}{\displaystyle dr}-f. So r′=fr^{\prime}=f at r=0r=0, and therefore the energy EAdS=0E_{AdS}=0. Then it follows that r′=fr^{\prime}=f on the entire surface inside AdS, which just says it’s a disk t=t= constant. Let’s label the surface by the rr-coordinates of its intersection with the shockwave rwr_{w}, then

At the intersection point, the momentum is

What’s the matching condition across the shock wave? In the corresponding mechanics problem, there is a change of Lagrangian at time v=0v=0. First we must require the coordinate rr to be continuous. Also, if we vary the end point along that v=0v=0 time slice, when the equation of motion is satisfied, the change of the action is prδrp_{r}\delta r. So the local extremal condition requires that momentum prp_{r} is also continuous across the shockwave.

We get that just across the shock wave in the BTZ part,

To simplify notation, let’s parametrize the crossing point by tan⁡θ=rwl\tan\theta=\frac{\displaystyle r_{w}}{\displaystyle l}, then

Solving this, and picking up the appropriate root, we get

This gives the shape of the extremal surface in BTZ part parametrized by θ\theta.

We still need to relate the parameter θ\theta to the boundary time tbt_{b}. At boundary, vb=tbv_{b}=t_{b}, so

To consider small time behavior, we can expand everything to linear order in π2−θ\frac{\displaystyle\pi}{\displaystyle 2}-\theta, then we find

At large time, look at the equation (A.76). θ\theta decreases as we increase tbt_{b}. The expression inside the square root stays positive and the integral is finite until we reach sin⁡θ=12\sin\theta=\frac{1}{2}. At sin⁡θ=12\sin\theta=\frac{1}{2}, i.e., rw=13lr_{w}=\sqrt{\frac{1}{3}}l, the integral is infinite as rr goes from 13l\sqrt{\frac{1}{3}}l to 12l\sqrt{\frac{1}{2}}l, which means that the surface ends at intersection of the boundary and singularity. This is the limiting surface at large time. In this case the expression inside the square root is a complete square, and we can do the integral of (A.74) analytically.

This limiting surface lies completely inside the horizon. In tt, rr coordinates, it’s

In the black hole region r−f(r)r\sqrt{-f(r)} is maximized on the surface r=rm=12l.r=r_{m}=\sqrt{\frac{1}{2}}l. The limiting surface we found in (A.79) lies entirely inside r=rmr=r_{m}, and approaches it exponentially in tt coordinate as we go to the boundary.

At large tbt_{b}, the growth of the surface is controlled by the r=rmr=r_{m} surface.

In more details, consider time and area integrals (A.76) and (A.75). rr goes from rwr_{w} to the boundary infinity. At large tbt_{b}, sin⁡θ\sin\theta becomes very close to 12\frac{1}{2}, and the denominators of both integrals become very small near r=rm=12lr=r_{m}=\sqrt{\frac{1}{2}}l. We can divide the integrals into three regions, (I) rr goes from rw=ltan⁡θr_{w}=l\tan\theta to rm−Δr_{m}-\Delta, (II) from rm−Δr_{m}-\Delta to rm+Δr_{m}+\Delta, and (III) from rm+Δr_{m}+\Delta to rb⟶∞r_{b}\longrightarrow\infty, where Δr\frac{\Delta}{r} is small. At large boundary time tbt_{b}, we expect the surface to get very close to the limiting surface. How close? The crossing point on the shockwave approaches its limit exponentially in tbt_{b},

Also, in the time integral (A.76), contributions from region (I) and (III) stay almost constant (Δ\Delta-dependent), while region (II) gives most contribution, increases like tbt_{b}, i.e, the surface spends most of its time around the r=rmr=r_{m} surface, roughly half inside and half outside. The behavior of the area integral (A.75) is similar. The contributions from region (I) and (III) are almost constant, while the contribution of region (II) goes like tb2l\frac{\displaystyle t_{b}}{\displaystyle 2l}. This is where the large time linear behavior comes out.

Here is a picture of the large time surface. As tbt_{b} increases, region (I) stays almost exactly at the limiting surface and gives a constant area. Region (II) stays around the Hartman-Maldacena surface and gets longer, giving an area linear in tbt_{b}. The position of region (III) changes since it’s pushed out by region (II). It roughly looks like

But its contribution to area remains again almost constant. We can see this by evaluating the time integral (A.76) in region (III), or, we can see this directly from the above expression (A.83). It’s essentially the same surface got pushed out in vv-direction as tbt_{b} increases. Its contribution to area will be divergent, goes like rbl\frac{\displaystyle r_{b}}{\displaystyle l}, but is tbt_{b}-independent anyway.

References