Towards the web of quantum chaos diagnostics

We study the connections between three quantities that can be used as diagnostics for quantum chaos, i.e., the out-of-time-order correlator (OTOC), Loschmidt echo (LE), and complexity. We generalize the connection between OTOC and LE for infinite dimensions and extend it for higher-order OTOCs and multi-fold LEs. Novel applications of this intrinsic relation are proposed. We also investigated the relationship between a specific circuit complexity and LE by using the inverted oscillator model and made a conjecture about their relationship. These relationships signal a deeper connection between these three probes of quantum chaos.


Introduction
Characterizing the nature of quantum chaos [1] in quantum many body systems can be challenging. This area of research is versatile and appears in many branches of theoretical and experimental physics. This has spurred a renewed interest in the quest for a quantum version of a classical chaos in the last few years. For recent developments, interested readers are referred to [2] and the references therein. Quantum chaos has found applications and received considerable attention across physical disciplines such as condensed matter physics, quantum information theory and high energy physics, in particular, in the context of black hole and holography [3]. Several diagnostic tools have been proposed to quantify it's diverse applications.
Over time, the endeavors to improve the current diagnostic gadgets and develop new ones have gone a long way. The out-of-time-order correlator (OTOC) [4,5] has been intensively used to examine the early time chaotic behaviour providing deeper understanding for long-standing problems. Loschmidt echo (LE), introduced as another powerful toolkit [6,7], has also played a pivotal role in demystifying the structure of (quantum) chaos. Very recently, a quantum information theoretic tool called quantum circuit complexity has joined the club of quantum chaos diagnostics [8][9][10][11]. In [12], the authors have shown that quantum complexity for a specific type of quantum circuit, namely, circuit complexity [13][14][15][16][17], can capture the chaotic features. Another particularly useful diagnostic is the spectral form factor (SFF) that interpolates between the quantum mechanical OTOC and the random matrix theory (RMT) [18]. In general, SFF is very hard to compute, apart from some special cases, due to various subtleties related to the spectral analysis of the chaotic Hamiltonian. Other measures characterizing the randomness of the quantum dynamics have been used in the studies of quantum chaos as well, such as unitary designs [19,20] and their corresponding frame potentials. Information theoretic tools such as tripartite mutual information [21,22] and entropy generation [23] are used to analyze more directly information propagation in chaotic systems.
There are serious indications that some of these proposed chaos quantifiers are related to each other. For instance, there had been a strong belief that the OTOC and LE are connected to each other due to the intrinsic nature of the echo of the OTOC [24][25][26]. Indeed, in [27] a major step has been taken to establish a direct link. It is worth mentioning that, previously there had been several attempts [26,28] to achieve the same goal, but all of them are resorted to some variants of the OTOC or specific choice of operators. In [12] a close connection between certain OTOC and complexity has been proposed. These results indicate a deeper connection between these diagnostics and provide motivations to explore it further.
In this paper we would like to establish a relationship between OTOC, LE and complexity to initiate a program that ultimately intends to investigate towards a complete web of quantum chaos diagnostics. To be specific, the aim of this paper is two-fold; (i) first, generalize the OTOC-LE connection of [27] for infinite dimensional system and extend it to k multi-fold and provide examples, (ii) second, explore the relationship between LE and complexity. To investigate the LE-complexity relation we use the inverted harmonic oscillator model for computing complexity and establish its connection to a particular type of LE. Finally we will comment and speculate on the possible ways to go beyond this example and highlight some future directions. Note that, in our discussion for the possible connection between the above mentioned diagnostics, we will not include another more recent diagnostic constructed from the retarded correlator for a Holographic system called "pole skipping" [29][30][31]. Pole skipping is related to the maximal Lyapunov exponent for maximally chaotic systems by analytic continuation. We would like to explore this connection and extend the web of the chaos diagnostics in a future work.

Loschmidt echo and OTOC
We start with an introduction of Loschmidt echo and the regular 4-point OTOC. Then we will discuss the general properties of these two quantities, as well as the intrinsic connection between them. Our first result is to generalize the link between the regular OTOC and LE to higher-order OTOCs and a echo quantity with multiple loops. This leads to a wide range of novel applications. The LE is formally defined as [7] where |ψ 0 is the initial state of a quantum system, H 1 and H 2 are two slightly different Hamiltonian, e.g., H 1 = H 0 is the unperturbed Hamiltonian, and H 2 = H 0 + V with V a small perturbation. One can interpret the LE in two ways. First, it can be considered as an "echo" process. It quantifies how much of the complex system is recovered upon applying an imperfect time-reversal, as sketched in Fig. 1. The other way is to interpret it as the overlap (the "distance") between two wavefunctions from the same initial state ("trajectories") evolving under slightly different dynamics. This is analogous to the classical notion of chaos, though in the latter case perturbations are applied to the initial condition in the classical phase space, while in the quantum case the perturbations are applied to the Hamiltonian. (Due to the fundamental unitary dynamics in quantum systems, any small perturbations on the initial wavefunctions remain unchanged during time evolution.) In this sense, the LE is related to the butterfly effect, so one can consider it as a diagnostic for chaos.
The regular 4-point OTOC is formally defined as Here the average is taken over a thermal state at inverse temperature β. W and V are two local operators on distinct local subsystems. W (t) ≡ e −i Ht W e i Ht is the Heisenberg evolution of the operator W . The OTOC has been extensively studied in various context and different variants of it has been proposed [32,33]. We note the following universal features of the OTOC: • When W and V are both Hermitian and unitary, the OTOC is related to the squared commutator Two local operator W and V commute at t = 0. The Heisenberg evolution converts W (t) into a global operator; the commutator hence fails to vanish and induces decay of the OTOC. For chaotic dynamics the OTOCs exhibits fast decays. • The OTOC has several decay regimes. At early stage before the Ehrenfest time scale (also known as the scrambling regime [34]), the decay of OTOC is manifested as an exponential growth, 1 − δe λt , where δ 1. In the scrambling regime, the exponential growth rate is conjectured to be bounded by the temperature, i.e., λ ≤ 2π/β from holography [34]. This type of decay certainly does not converge, and will switch to a pure exponential decay e − t (intermediate regime) before saturation. The scrambling time is determined by the small parameter δ: For fast scrambling systems, the scrambling time scale is proportional to the logarithm of the number of degrees of freedom [34]. Such a small parameter δ does not always exist. In many situations early scrambling of the above form, even when exists, can only last for a transient moment (see Refs. [35,36] for detailed discussions). In these cases the majority of the decay of the OTOC is in the pure exponential regime. It is recognized [37] that systems that can exhibit clear and well-separated early scrambling are typically in some kinds of semi-classical limit, or a have a large number of degree of freedoms, e.g., large-N gauge theories. For lattice models, e.g., spin chains with finite number of local degrees of freedom, only the intermediate exponential decay is visible. In the asymptotic regime (late time) the OTOC typically shows model-dependent power law behaviors, and saturates to universal values [38,39].
• The OTOC recovers the essential pieces of the classical notion of chaos in phase space. A heuristic way to see this is to look at the semi-classical limit [34] for the choice of operators W (t) = q(t) and V = p, where q and p are conjugate pair of variables. In the semi-classical limit, the commutator reduces to the Poisson bracket. This gives [q(t), p] → ih{q(t), p} ∼ ih∂q(t)/∂q(0), which grows as e λ L t with λ L the Lyapunov exponent. For various systems with classical counterparts, e.g., the kicked rotor [40] or the cat map [41], the decay rates of the OTOC were shown to match the classical Lyapunov exponents. • Different choices of operators of W and V share common features of their OTOCs. For complex enough systems, the OTOC is expected to be not sensitive to the particular form of the operators, especially when we are interested in extracting the universal characteristics. This makes it possible to get rid of the operator dependence by averaging over all operators of given subsystems. In the following sections we explore the consequences of this averaging procedure. We will show that after averaging over local operators, the OTOC takes the form of an LE. This relation was demonstrated [27] to hold in the intermediate (exponential decay) regime, and also works accurately in the scrambling (exponential growth) regime for a model of coupled inverted harmonic oscillators. Whether this relation generally holds in the scrambling regime is not clear yet. However, it is well-known [7] that the LE in the semi-classical limit can decay in a double exponential form, i.e., exp (−a exp (bt)), which, upon a is small, expands at early times into the universal exponential growth form of the OTOC as discussed in the previous points. This provides evidence for the OTOC-LE connection in the scrambling regime, and shed new light on the understanding of the early scrambling behaviors of the OTOC.
As noted in the previous section, the insensitivity to the choice of operators allows one to extract the universal features of OTOC by taking the average over a given set of operators. The averaging procedure has been considered for different variants of the OTOCs. When restricted to the original form with local operators, there exists a strong relation between the Loschmidt echo and the 4-point OTOC: Without losing the local structure of a many-body system, the supports of two operators W and V are chosen as two distinct subsystems A and B, where A is a small subsystem, while B is the complement of A to the total system, as illustrated in Fig. 2. We then take the average of the two operators over the set of all unitaries on the two fixed subsystems with the "largest randomness", i.e., with respect to the Haar measure. It has been demonstrated in [27] that the OTOC and LE are ultimately related as Here the Hamiltonian of the larger subsystem B plays the role of the unperturbed Hamiltonian; and the perturbation P naturally emerges from the interaction between the two subsystems. It is the projection of the interaction to the Hilbert space of the subsystem B (see Ref. [27] or the following derivations for the construction of the effective perturbations). This relations was shown to be valid in both the scrambling and the intermediate decay regime. In the following section, we will further generalize this result for 2k-point OTOC to the 2(k − 1)-fold Loschmidt echo, and also to the case where both of the operators W and V are supported on small local subsystems. For simplicity, we will restrict ourselves to the case of infinite temperatures. With proper regularization of the thermal state, this relation generalizes to finite temperature as well, using techniques developed in Ref. [27].
2.1 2k-OTOC and 2(k − 1)-fold echo We will start with a formal definition of the 2k-point OTOC.
Then we will demonstrate that it is linked to a LE with 2(k−1) forward and backward loops. The regular 4-point OTOC, probes the spreading of the local operator W over the entire system. The work of Refs. [19,42] suggests the study of the generalized 2k-OTOC, defined as where W ≡ W k−1 (t k−1 ) . . . W 1 (t 1 ) indicates the ordering of the operators in the correlator. Note that there are other definitions of the 2k-OTOC, such as the ones used in [19] to probe k-designs and the one connected to the spectral form factors in [20]. The operators in the correlator could be interpreted as either global as in [19,20]) or local operators as in [42]) for different purposes. We first focus on the most restricted choice, i.e., the operator W k 's are all local operators applying on distinct local (and small) subsystems, such that the 2k-OTOC mentioned in (5) probes the scrambling of multiple local perturbations. We choose V as an operator on the complement of the k − 1 local subsystems. The structure of this 2k-OTOC is illustrated in Fig. 3. We denote where U k = e i Ht k is the evolution operator. Consider the averaged 2k-OTOC with respect to the Haar integral at infinite temperature, Here d is the dimension of the total Hilbert space. The particular ordering of operators in the integrand allows us to perform the integral one-by-one, e.g., the inner-most integral for the W k−1 operators can be computed first (see Appendix A for Haar average over subsystems): where T r k−1 represents the partial trace over the subsystem k − 1, and d k is the dimension of the local Hilbert space supporting operator W k . Performing the integral for all the W operators gives Let us define the following: Using the same trick provided for the 4-point OTOC in [27], the partial traces in piece-A can be evaluated one-by-one (see Appendix A). For instance, the inner-most partial trace is Here N k−1 is the number of different P k−1 operators, which serve as the perturbations. The summation range over all of them. These noisy operators emerge from the interaction between the k−1's subsystem with the rest of the total system (see Appendix A for details). Note that the LHS of the above equation, after tracing over the (k − 1)'th subsystem, is an operator that involves not only the subsystem-V, but also subsystems 1, 2, . . . , k − 2. However, we assume that it only evolves (under noises) in subsystem-V and it does not "leak" to other subsystems.
, the above procedure for partial tracing can be repeated to all partial traces, which give the following expression for A. Before proceeding further, let us define the following, Then we get, Finally, As has been discussed before, U 1 V U † 1 is a global operator, while DV † D † is an operator with support on system-V only. Thus the trace in the above equation can be evaluated with two partial traces Tr = Tr V TrV , namely, Tr which is the dimension of the Hilbert space of the subsystem complementary to subsystem-V . The right hand side of the Eq. (10) continues as (11) P 1 , . . . , P k−1 are perturbations emerge from the tracing our the subsystems −1, . . . , k − 1; and they have, respectively. P 0 emerges from tracing out the subsystemV . For complex systems, the structure of these perturbation operators are not essential. Hence, we can eliminate the average over all the perturbations and treat each P i as a constant perturbation instead of a variable. Finally we get, This is the expected Loschmidt echo with 2(k − 1) loops.

Local OTOCs
In the previous section, we have established the OTOC-LE correspondence for one type of generalized 2k point OTOC. This 2k-OTOC involved many local operators, and a giant operator over the compliment system of the other local operators (see Fig. 3) for illustration. In this section, we consider a even more general case, where the OTOC only involves operators on small local subsystems. For simplicity, we only consider the case of the 4-point OTOC and infinite temperature. Generalization to the multi-point case follows from the techniques developed in the previous section. The four point OTOC at infinite temperature is with the operators W and V constrained to local subsystems (see Fig. 4 for illustration.). Averaging W and V over all unitaries on the corresponding subsystems, we get In the above equation, the reduced evolution of operator W (partial trace Tr V over the subsystem of V ) can be approximated with where d V is the dimension of the support of operator V . HV is the Hamiltonian of the subsystem excluding the support of V . Hence, the evaluation of the four point OTOC continues as We now use again the approximation for the reduced dynamics, i.e., where R labels the part the total system that excludes the subsystems of W and V , Note that a new emerged perturbation P W . Hence, This is a special type of LE, with four loops and the perturbations are local in each loop, i.e., P W and P V are local perturbation emerged from the contact with the subsystems of W and V , respectively. Let's redefine the unperturbed Hamiltonian as H = H R + P V + P W in the above expression, and extract a simplified form of this LE, i.e., In Sect. 2.4, we will present an application of this particular LE for detecting the butterfly velocity.

Infinite dimensional generalization
The previous discussions focus on finite dimensional Hilbert spaces. 1 In this section we argue that the OTOC-LE connection can be generalized to infinite dimension. The key ingredient is the Haar integral for unitary operators U on an O is a trace-class operator and μ is the Haar measure.
Here we consider the right Haar measure, which, by def- This means that the Haar-averaged operator is proportional to the identity operator I. In finite dimensions, its trace can be computed as Haar measure is unique up-to a constant multiplication factor; and the unitary groups on finite dimensional Hilbert spaces have finite measures. This allows us to normalize the Haar measure by choosing dμ = 1. Under this convention, the averaged operator has the representation where d is the dimension of the Hilbert space. For infinite dimensions, the Haar measure is not normalizable, and hence the averaged operator is not trace-class anymore. However, we are interested in the case where the averaged operator is still bounded (the OTOC takes finite values). In this case, the Haar averaged operator can be fixed as a constant multiplied by the identity, c(O)I. The functional c must be linear and invariant under unitary transformation, i.e., c(O) = c(U † OU ). By Riesz representation theorem, it is determined, up-to a multiplication factor, to be the trace, i.e., c(O) ∝ Tr(O). We have the freedom to remove the prefactor by absorbing it into the definition of the Haar measure. Under this convention, the desired integral for the Haar average matches precisely with the one in finite dimensions. Once this infinite dimensional Haar integral is evaluated, the OTOC-LE connection follows in the same manner as in the finite dimensional case.
If we average the OTOC over a given group of unitaries {U g }, rather than performing the average over all unitary operators with respect to Haar measure, we can firmly say that the OTOC-LE connection holds as well, as long as the group average, up-to a constant multiplication factor which can be removed by re-scaling the measure, gives the same result as the Haar average, namely, In other words, the group {U g } is an analog of the unitary 1-design in finite dimensions.
As an example, consider the Heisenberg group {U (q 1 , q 2 ) = e i(q 1x +q 2p ) }, wherex andp are the canonical position and momentum operator, q 1 and q 2 are real numbers.
To show that the Heisenberg group is a unitary 1-design, we will need to prove, for any trace-class operator O, This is equivalent to showing that the above operator D in the position representation has elements

Application I: butterfly velocity
The OTOC is designed as a diagnostic for chaos. For chaotic systems, it decays rapidly and converges to a persistent small value. While for integrable systems the OTOCs typically exhibit oscillatory behaviors, with finite recurrent times. Another intriguing feature of the OTOC is that it can detect information propagation in systems with higher spatial degree of freedoms. For instance, for a 1-D chaotic spin chain with local interactions, if the operators are chosen as Pauli operators on distinct sites, e.g., W =σ z i and V =σ z j , the OTOC does not decay immediately. Rather, it stays constant for a finite amount of time. This is the time for the operator W (t), which is initially local on the i'th site, to propagate to the j'th site. The propagation is ballistic, with a velocity known as the butterfly velocity [44][45][46][47][48]. This effect is absent for a regular LE, which decays immediately even for local perturbations.
Here, as a first novel application of the OTOC-LE connection, we propose to use a multi-loop LE to probe the butterfly velocity in chaotic systems. As derived in Sect. 2.2. the fourloop LE relating a regular 4 point OTOC with local operators takes the form, Here P 1 and P 2 are both local perturbations. Regular twoloop LEs of the form e i(H +P 1 )t e −i Ht and e i(H +P 1 )t e −i(H +P2)t both show instant decays with no dependence on locality.
We apply the above four-loop LE to study the butterfly velocity of a 1-D spin chain system. The Hamiltonian is given by whereh is set to zero. 1/J measures the unit of time. The parameters are fixed as h x /J = 1.05 and h z /J = 0.5, for which the model is known to be chaotic [49][50][51]. In our simulation of the LE, we chose perturbations as P 1 = gσ z 1 , the Fig. 5 Comparison between the regular two-loop LE and the four-loop LE. The former shows immediate decay after perturbation, while the latter exhibits an initial plateau regime Fig. 6 The decay of the four-loop LE at various n, where n labels the site on which the second perturbation P 2 is applied. The first perturbation P 1 is always applied to the first site. Inset: The time at which the LE starts to decay as a function of n, the distance between the two sites. Red crosses are numerical data. Solid line is the best liner fit. Note that the slop of the linear curve is the inverse of the butterfly velocity (in the unit of J ) Pauli operator on the first site, and P 2 = gσ z n on the n'th site. g = 0.2 is the strength of the perturbations. The total number of spins is N = 12. Figure 5 compares the regular LE, which decays exponentially immediately after perturbation, and the four-loop LE (with n = 12), which exhibits an initial plateau regime.
To extract the butterfly velocity, we simulate the four-loop LE at various n's -the position of the second perturbation, and read out the time at which the LE starts to decay. In Fig. 6, it can be seen that the larger the distance between the two perturbations, the longer the initial flat regime is; namely, it takes a longer time for the local perturbation to propagate. The width of the plateau regime is proportional to the distance between the two sites (site-1 and site-n), indicating that the propagation is indeed ballistic.

Application II: shockwave and Loschmidt echo
We further illustrate the relation OTOC-LE with a brief application from AdS/CFT correspondence. We shall examine this   Fig. 7. This subsection is primarily based on the following work [52][53][54][55][56][57][58][59] and references therein. We shall also shed some light into the possible link that may exist between OTOCs, LE and quantum complexity which we will elaborate on in the next section.

Echo evolution, precursors and black holes
Let's consider two entangled black holes connected by an Einstein-Rosen bridge, aka wormhole Fig. 8. The holographic description of the wormhole volume is quantified by the complexity of the quantum state of the dual pair of CFTs at time t. For a given thermofield double state, we can evolve, for instance, the left side back in time for a time t L = −t w and then apply a simple localized precursor perturbation W L that adds a thermal quantum; a localized packet of energy in the left side. Then, we evolve this state forward in time, t L = t w (see Fig. 9). Due to the fact that the quantum state loses its memory, the left target state has to differ from the left initial state.
In AdS spacetime, this energy source starts to warp the spacetime near the horizon by creating a gravitational shockwave which expands away from the source, and remains Fig. 9 The operator W L creates an infalling quantum at |t w | >> t * , where t * is the scrambling time. It undergoes a huge blue shift as it moves at the speed of light toward the horizon highly energetic for most of its worldline. This kick changes the geometry and leads to a larger wormhole compared to the initial one. More precisely, the wormhole owes its growth to the gravitational back-reaction on the shape of the geometry (or alternatively to the hop and displacement of the trajectories crossing the shockwave Fig. 10). This corresponds, on the CFT side, to a greater decay for the correlation between the two different sides of the thermofield-double state. Consequently, one infers that there exists a correspondence between the amount of energy the shockwaves produce and the rate of decay of the correlations. In what follows, we shall provide a rough picture as to how this might play out.

Single shock
Let us apply the echo evolution e −i H L t w W L e i H L t w to the thermofield-double state: where for any single sided operator (precursor) W , W L = W ⊗I and W R = I⊗W . The operator W L (t ω ) is a Schrodinger picture operator acting at time t = 0. The effect of W L (t ω ) amounts to adding at t ω a thermal quantum to the left side. Note that despite the fact that the thermal quantum, being localized low energy perturbation, created by W L (t ω ), feeds a tiny bit of energy to the black hole, for most of its worldline it is astronomically energetic shockwave. The two-sided correlator is found to be (see for instance Fig. 10 In the absence of the shockwave the size of the maximal spacelike slice formed behind the horizon is null (it goes through the intersection point of the bifurcate horizon). Adding the shockwave allows the maximal slice to acquire a considerable volume as indicated by the red surfaces The transpose "T" is in the energy basis. The state is subject to the so-called operator pushing property by which we mean with |max being a maximally entangled state. Using (29), one can push V T from the right to the left. One therefore can convert the two-sided correlator to a one-sided correlator, i.e., correlation between the 2 sides after perturbation The negative time is not profoundly significant. Evolving the system according to negative time is expected to have the same behaviour as evolving it with positive time. As a matter of fact, the behaviour of the OTOC can be generic for different local operators. In this case for any thermal state the aforementioned argument may be extracted from The generalization of the previous claim to general temperature, for which the obtained OTOCs are thermally regulated, is straightforward [59].

Multiple shocks
The lesson one can draw out from the previous single shock case is that the more shockwaves (energy) you feed the black hole with, the greater the decay of the OTOCs becomes. To create two shockwaves one needs to consecutively repeat the Fig. 11 Multifold echo with 2(k − 1) loops. Each green insertion represents a tiny perturbation. The arrows point toward the order in which the precursors W i 's apply process introduced above twice, that is, from which one can derive correlation between 2 sides after perturbing twice Notice that we can get rid of the negative time because of the same previously mentioned reasons. This makes our claim true for the case of two shockwaves. Upon the application of multiple (k − 1) operators W (t) on |T F D (pictorially presented in Fig. 11), One can write from which one ought to obtain 2k-OTOC where k > 1, i.e., correlation between 2 sides after perturbing (k-1) times

The triangle links
As claimed above the shockwave has a large effect on the geometry. Without the shockwaves the volume of the maxi-mal slice behind the horizon at t = 0 is null (it goes through the bifurcate horizon). Upon the creation of the shockwaves the spatial maximal slice, representing the wormhole connecting the two-sided entangled black holes, gains a significant volume. Roughly speaking, one can anticipate that the correlation exponentially decays with the size (length) L(t) of the wormhole. Thus, 2 Using (30), it yields From the OTOC-LE connection, we have This implies, We find that It was conjectured by Susskind et al. that quantum complexity C (precisely introduced in Sect. 3) is related to the size (length/volume) V of the wormhole connecting the two entangled black holes, i.e., C = V Gl Ad S . Combining all these together we end up with Then conjugating (40) both sides by T F D| we get, Finally we arrive at the following, This derivation involves only one single shockwave. However, one can incorporate multitude of shocks for which the one-fold LE is superseded by multi-fold LE and the complexity associated with one localized precursor is replaced by whereL(t) is the stretched length of the ERB (wormhole) behind the horizon.

Perspectives from infinite dimensional continuous variable systems
A relation between an operator's distribution in phase space and OTOCs in continuous variable (CV) system has been established in [60]. Consider an operator that spreads in phase space having width/volume V. The OTOC was found to be [60] To derive (45) we shall introduce a few definitions and quantities. We begin by defining the displacement operator, the analog of the Pauli operator in discrete variables, for a simple harmonic oscillator (single mode CV system) Such shifts operators, being e.g., elements of the Heisenberg group, form a complete basis and act on a coherent state in phase space. For N -mode CV system they read with x x x = (q 1 , q 2 , . . . , q N , p N ) being the vector of quadrature operators. These N -mode displacement operators satisfy Tr(D(ξ ξ ξ)D(ξ ξ ξ )) = π N δ(ξ ξ ξ + ξ ξ ξ ), The CV OTOC is defined to be where the so-called displacement operator takes the following form From the above decomposition, which is allowed by the completeness of displacement operators, one can infer that scrambling in the CV system is featured by the growth of the Wigner characteristic χ [ξ ξ ξ 2 ; D(ξ ξ ξ 1 ; t)] given by Now, given and making use of the formulae presented above leads to (45). This shows that a larger width of the phase space results in greater decay of the OTOC. One should be able to directly relate the increase of the operator volume in the phase space with the size of the wormhole in the two entangled black holes model studied above rendering the connection between OTOC-LE and complexity more rigorous. This correspondence may be achieved by matching the norm of the displacement vector with the AdS radius, i.e., |ξ | ∼ 1/l AdS (cf. Eq. (36)). It is worth emphasising that averaging the OTOCs over ensembles of displacement operators may enable us to measure a coarse-grained spread of a time-evolved operator in phase space. This may allow one to gain better understanding into the link between various diagnostics. For more about using the average OTOCs as probes for finer-grained aspects of operator distribution, we refer the reader to [60]. In Sect. 3 we shall pursue a slightly different path to establish such a connection between the three diagnostics.
Due to some very recent progress one is now able to test these predictions. Based on the formalism developed in [60] and a generalization of quantum teleportation mechanism, a detailed experimental blueprint has been put forward [61]. The proposed experimental protocols can be potentially generalized to include (multi)-shockwaves. Along similar lines, Brown et al. have recently set a long-term goal of studying models of quantum gravity in the lab [62,63], which could mimic in particular the two entangled black hole setup. More precisely, they put forward holographic teleportation protocols that can be readily executed in table-top experiments. These quantum-teleportation-inspired experimental protocols are malleable to be devised such that they could potentially include multiple shockwaves, from which one would be able to test our predictions in experiments.
In summary, chaos is a keynote ingredient intimately related to the onset of thermalization. It has been shown that a signal of chaos is encoded in the behaviour of OTOCs and exponential growth of the commutators. In the gravity side, this growth manifests itself as a near horizon higher energy scattering/collision; semi-classically controlled by a shockwave geometry. A particle propagating at the speed of light can be described as a null-like delta function of the stressenergy tied to the horizon of the black hole. This source triggers what is called a gravitational shockwave that moves far away from the source. Going through the shockwave results in a kick in a certain null-like direction. Such an extra kick modifies the geometry, leading to a kind of decorrelation between the two sides of the perturbed thermofield double when the kick becomes strong. The boost the particle has gone through after falling for a time t ω is proportional to the strength of the kick that scales like the proper energy ∼ e 2π t ω /β . This dependence on the exponential time is one quantifier/measure of quantum chaos in the given system, with λ L = 2π β -being the quantum Lyapunov exponent of the black hole. This picture can be generalized for multiple shockwaves, which implies a clear connection between quantum chaos/fast scrambling and the recently established transversability of the wormhole.

Loschmidt echo and complexity
We will start with a brief review of the circuit complexity by using the Nielsen's method [13,14]. Given a reference state |ψ s=0 , a target state |ψ s=1 , and a set of elementary gates {exp(i M I )} where the {M I } are group generators, the goal is to build the most efficient circuit U (s) that starts at the reference state and terminates at the target state: where U (s) is the path-ordered operator where H (s ) is the Hamiltonian and can be written as The coefficients Y I are the control functions that dictates which gate will act at a given value of the parameter. The control function is basically a tangent vector in the space of unitaries and satisfy the Schrodinger equation Then we define a cost functional F(U,U ) as follows: Minimizing this cost functional gives us the optimal circuit. There are different choices for the cost functional [14]. In this paper we will consider

Introducing complexity
Recently complexity has been demonstrated as an equally powerful and computationally simpler quantity in some cases than OTOC to diagnose the chaotic behaviour of a quantum system [33? ]. In [10], by using the SYK model, the author shows that the complexity grows linearly for exponential time for a chaotic system. They also argued that the complexity would eventually saturate. In these papers, authors have used Nielsen's circuit complexity [13] for the Hamiltonian evolution operator. Furthermore, yet another notion of complexity, namely "Krylov Complexity," has been used to diagnose quantum chaos [64]. Although it has been argued that the early time growth of this Krylov Complexity can be thought of as an indicator of quantum chaos [65,66]. However, lately, counterexamples have also been provided, demonstrating that this early time growth of Krylov Complexity may happen even for free theories [67]. But we like to emphasize that none of these studies captures indicators like Lyapunov exponents, scrambling time of quantum chaos using the complexity. The authors of this paper have used a simple quantum mechanical model to study Nielsen's circuit complexity in the space of wavefunctions. They have demonstrated that it has the potential to provide information about scrambling time, Lyapunov exponents for some simple quantum mechanical systems [12,[68][69][70][71][72]. In the rest of this paper, we will follow this approach. Since all three of these quantities -LE, OTOC and Complexity -are providing similar information about the chaotic system, it is natural to anticipate that these three quantities are related to each other. In the previous sections we have established that the sub-system LE and averaged OTOC are the same. Therefore, to establish the relationship between the three quantities, we only need to explore the connection between LE and complexity.
To make progress in this direction, we will use the complexity for a particular quantum circuit from the inverted oscillator model: Classically, the inverted harmonic oscillator has an unstable fixed point and is not a chaotic system in the strict sense. Nonetheless, it has been used as a powerful toy model for studying quantum chaos in various quantum field theories [73][74][75][76][77], mostly because it is an exactly solvable system. The oscillator can be tuned to the regular and chaotic regime by changing the value of λ, i.e., for λ < m 2 the oscillator is simple, while for λ > m 2 the oscillator is inverted and chaotic.
In [12], it was shown that the appropriate quantum circuit in this regard is the one where the target state |ψ 2 is obtained by evolving a reference state |ψ 0 forward in time by Hamiltonian H and then backward in time with slightly different Hamiltonian H + δ H as follows For the inverted harmonic oscillator model the authors in [12] showed that the complexity evaluated by using the covariance matrix method [17,78] for the above mentioned target state with respect to the reference state |ψ 0 is given by whereω(t) is the frequency of the doubly evolved Gaussian target state which has the following form and In the last expression, = √ m 2 − λ is the frequency associated with the perturbed/slightly different Hamiltonian H = 1 2 p 2 + 2 2 x 2 and λ = λ + δλ with δλ very small. We make this perturbation by hand.
Note that the quantum circuit involving two time-evolutions with slightly different Hamiltonians is crucial for extracting the chaotic nature of the quantum system. Complexity for any target state will not capture similar information as OTOC. For example, the complexity of a target state which is forward evolved only once will not capture the scrambling time for the chaotic system as illustrated in Fig. 12. However, there is an alternative quantum circuit that will have the same complexity when evaluated by the covariance matrix method. In that circuit both the reference and target states are basically evolved states but with slightly different Hamiltonians from some other state. Once again this particular circuit also involves two evolutions.

LE-complexity connection
It was shown in [12] that complexity of the above mentioned target state (62) can capture equivalent information such as scrambling time and Lyapunov exponent as the OTOC for an inverted oscillator. In this paper, we want to make this statement more precise by using the fact that averaged OTOC is the same as (very close to) the sub-system LE. In Sect. 2.4 of this paper, we have proved this for the Heisenberg group. In the current section, we will use an explicit example from the Heisenberg group, namely the inverted oscillator to demonstrate that LE for the full system and complexity are very close quantities.
It is noteworthy that the construction procedure of this quantum circuit is conceptually similar to the LE, where one basically computes the overlap between these above mentioned states. Complexity simply offers us a different measure for the distance which is a more powerful measure for understanding various properties of quantum systems [79][80][81][82].
As shown in Fig. 13, both of the time evolutions of the LE and complexity exhibit two regimes of growth, i.e., an intermediate regime where complexity grows linearly in time, while the LE decays as a pure exponential function (Fig. 13, left); an early regime (scrambling [34]) where complexity grows exponentially and the LE decays as a double exponential. The growth pattern of complexity and the LE are suggestive of the following relationship between LE and complexity, and their universal forms: In the early scrambling regime, complexity and the LE have the form Figure 13 confirms that during the scrambling stage, these two quantities are remarkably close. The same Lyapunov exponent λ, which is a system characteristic, can be extracted from both of these two quantities. The time scales of the early scrambling are also the same. Note also that the double exponential decay of the LE, when expanded to first order of , reassembles the conjectured universal form of scrambling of the out-of-time ordered correlators of the form 1 − e t . Hence, in the scrambling regime, we conjecture a universal relation between complexity and the LE: While in the intermediate regime, we have observed similar relations between the their growth rates, i.e., indicated by the evolution forms Note that though they exhibit the same growth rate , we expect this to be an artifact of the harmonic oscillator model. The exponential decay of the LE is standard feature and it is well-known [6,7] that its decay rate is a perturbationdependent quantity, rather than a universal characteristic of the system. We can easily generalize this particular construction of quantum circuit to relate it with 2(k − 1)-fold LEs. The trick is to insert a pair of evolutions (forward and followed by a backward) for each fold of the echo. For example, for the In Fig. 14 we show a few of the higher fold-LE (− log [LE], to be precise) and the corresponding generalization of complexities. For each pair we see a clear match between complexity and the − log [LE] during the linear portion. The right panel of Fig. 14 displays the early time behaviour of these two quantities, which is similar to findings for the single fold case. We would like to stress that we do not have a concrete algebraic proof to establish the relationship at this point, therefore, it is just a conjecture and a concrete prove; we leave it for a future work. Note that the sub-system LE that we have used in the previous sections can be quite close to full system LE, when the sub-system associated with the LE is much larger that the other one. We will conclude this section by making the assertion that these three diagnostics of chaos -averaged OTOC, LE and a particular type of complexity -are not only carrying similar information about the underlying quantum system, but also have some direct connection with each other.

Discussion
In this paper we have extended the proof that the averaged (Haar average over unitaries) OTOC is the same as the LE (for a sub-system) as in [27] to higher point averaged OTOC and LE for finite dimensional system. Moreover, we have also generalized the proof for Haar average to infinite dimensional case. We have shown that the OTOC-LE relation holds in other averaging scenarios as well, e.g., the Heisenberg group average, as long as the given group is a unitary 1-design. We argue that if the sub-system for this LE is much larger than the other sub-system, this LE would be essentially the same as the LE of the full system. Furthermore, for an explicit example in the Heisenberg group we showed graphically that LE for the full system and complexity for some special type of quantum circuit is the same. Finally, we have extended this result for multi-fold LE and corresponding extensions of the complexity. These different results suggest that these three diagnostics of a chaotic quantum system, namely averaged OTOC, LE and complexity are secretly the same. However, we do not have a concrete proof at this point. Tying complexity as an alternative probe to OTOC or LE also provides a geometric meaning to the chaotic behaviour of a quantum system.
To give a proof-of-principle argument for the similarity between complexity and LE, we have used the inverted oscillator as a toy model. This is, however, a rather special example and not a realistic chaotic system. Also, we used graphical techniques to establish our result. To claim that our particular complexity and LE (and hence averaged OTOC) are basically the same probe for understanding a quantum chaos will require a rigorous algebraic proof by using more 'realistic' systems like the maximally chaotic SYK model and its many variants (see, for example, [83][84][85] and references therein).
Another possible extension of our work is to explore subsystem complexity in a system with N-inverted oscillators. This would help us make the connection between these quantities more rigorously. N -particle system with all-to-all two-body interactions, when the subsystem S A refers to a single particle, the energy scales of S A , S B , and the coupling between them, are on the order of 1, N 2 and N , respectively. The interaction can be decomposed as Here we are free to chose the operators {V i A } Hermitian and orthonormal, with respect to the Hilbert-Schmidt inner product, i.e., where d A is the dimension of the Hilbert space of S A . The operators V i B on S B are also Hermitian, but their (Hilbert-Schmidt) norms are fixed as equal to the norms of H B . Thus, the parameter λ qualifies the relative strength of the coupling compared to H B . We are interested in the reduced dynamics of an operator B on the subsystem S B , after the trace-out procedure, namely, This can be thought of as a decoherence process, i.e., the total system is prepared in an initial product state I A ⊗ B, where the subsystem S B has a "density matrix" B, and the subsystem S A , up-to normalization, is in a thermal state with infinite temperature. The "quantum state" B will become "mixed" with time evolution due to the presence of the couplings to subsystem S A . When λ 1, the above evolution of B(t) can be expanded to the second order of λ. This corresponds to the Born-Markov approximation, which leads the effective master equation for B(t) to a Lindblad form. It is known that in this case the effective master equation can be simulated with the evolution of B under H B without coupling to other systems, but subjects to a stochastic field with the correlations given by The approximation in the last step is due to the large energy hierarchy: the time scale of the dynamics of the subsystem S A is much larger than that of B(t) under consideration. Alternatively, this can be thought of as taking the zeroth order the H A . As a consequence, the noise field l i (t) can be taken as random constant valued, ±1, at equal probability. The reduced dynamics of the B operator is then given by averaged over the stochastic field. Note that the pre-factor d A appears from the normalization of I A . As the noise field are random ±1, each realization of the stochastic field F in the above solution of B(t) always appears as random combination of V i B 's. Suppose that are totally N realizations, the noisy evolution of B(t) is then (A12)