Barrier from chaos: operator entanglement dynamics of the reduced density matrix

It is believed that thermalization drives the reduced density matrix of a subsystem to approach a short-range entangled operator. If the initial state is also short-range entangled, it is possible that the reduced density matrix remains low-entangled throughout thermalization; or there could exist a barrier with high operator entanglement between the initial and thermalized reduced density matrix. In this paper, we study such dynamics in three classes of models: the rational CFTs, the random unitary circuit, and the holographic CFTs, representing systems of increasing quantum chaoticity. We show that in all three classes of models, the operator entanglement (or variant of) exhibits three phases, a linear growth phase, a plateau phase, and a decay phase. The plateau phase characterized by volume-law operator entanglement corresponds to the barrier in operator entanglement. While it is present in all three models, its persistence and exit show interesting distinctions among them. The rational CFTs have the shortest plateau phase, followed by the slowest decay phase; the holographic CFTs mark the opposite end, i.e. having the longest plateau phase followed by a discontinuous drop; and the random unitary circuit shows the intermediate behavior. We discuss the mechanisms underlying these behaviors in operator entanglement barriers, whose persistence might serve as another measure for quantum chaoticity.


Introduction
Quantum entanglement has played a central role in elucidating many recent progresses of theoretical physics, ranging from being an order parameter for topological orders [1,2] to an explicit probe for bulk geometries in AdS/CFT correspondence [3,4].Away from the static settings, the dynamics of quantum entanglement is also a crucial and universal ingredient in understanding systems out of equilibrium [5,6].It encodes the information about thermalization and entropy generation in the process of regaining equilibrium, and probes some surprising aspects of the systems such as pre-thermalization [7][8][9][10], many body quantum chaos [11,12], quantum scars [13][14][15], and in the context of AdS/CFT even the black hole interior behind the horizon [16].
While the entanglement properties are usually defined and studied with respect to a state |ψ , one can also study the entanglement properties of an operator U : H → H that maps within a Hilbert space H. Operationally we can pick a set of orthonormal basis in H and its dual H : {|i } and { j|}, and write U = ij U ij |i j|.The operator state can then be obtained by 1

|U
A convenient example to familiarize the set up is to imagine the operator being the thermal density matrix ρ β , and the corresponding operator state of ρ β being the thermal field double (TFD) state.The operator entanglement then refers to the entanglement of the operator state |U , defined with respect to a subsystem A ⊂ H ⊗ H [17][18][19][20].
In this paper, we study a particular class of dynamical operator entanglement: those associated with the reduced density matrices obtained from quenched states 2 3 .Let us lay out the set up step by step.Given a quenched state |ψ(t) that starts with only short-range entanglement at t = 0, and is defined on the total system ABC, we first obtain the reduced density matrix on the subsystem AB by tracing out H C : For the interest of this paper, we take C to be semi-infinite that surrounds the finite subsystem AB like a heat-bath.The wave-functional of the corresponding (unnormalized) operator state |ρ AB (t) on the doubled subsystem Hilbert space H AB ⊗ H AB (see Fig. 1(a)) is then given by the density matrix element (ρ AB (t)) ij = i|ρ AB (t)|j as [24,25]: where {|i } is an orthonormal basis in H AB .
To normalize the state, we need the norm of the |ρ AB (t) , which is the purity of the state ρ AB (t)2 = ρ AB (t)|ρ AB (t) = tr(ρ 2 AB (t)) (1.4)It is in general less than 1 and t-dependent when region AB is entangled with C. We are interested in computing the operator Rényi entropy of the normalized operator state |ρ AB (t) / ρ AB (t) on the subsystem A ⊗ A: where we use S op n (A, ρ AB ) to denote the operator Rényi entropy of ρ AB , and S n (AB, ψ) as the state entanglement of the quenched state |ψ .As |ψ(t) evolves in time t, so does ρ AB (t) and consequently S op n (A, ρ AB (t)).Finding out the t-dependent dynamics of S op n (A, ρ AB (t)) is the main task of this paper.
An important motivation to study the dynamics of such operator entanglement comes from the following question that is both theoretically interesting and of potential practical significance: is it possible to simulate the density matrix efficiently in the matrix product operator (MPO) framework throughout the thermalization process?Since large operator entanglement indicates large bond dimension in a MPO (Fig. 1(a)), the question can be quantified as weather or not to have low operator entanglement all the way to thermalization.
Various proposals [20,[26][27][28] have been made for directly simulating the reduced density matrix to study the evolution to the thermal states.The operator entanglement, crudely speaking, measures the logarithm the bond dimension to represent the reduced density matrix by a matrix product operator (Fig. 1(a)).In the process of quenched thermalization, the initial state is short-range entangled, and as a result the operator state is also short-range entangled with a low operator entanglement.On the other hand, if the system thermalizes, the final reduced density matrix will be close to an identity operator at high temperature.Thus while the state entanglement increases and saturates to volume-law, the operator entanglement is again very small.Assuming this feature of low operator entanglement persists in between, numerical algorithms that use variational approaches [26][27][28] at the low bond dimension manifold should be enough to reproduce all the correlation functions faithfully.
The question is whether the feature of low operator entanglement persists or not.In general, the time-dependence of operator entanglement may result in intermediate "bumps".If the values at the bumps are comparable to system size, then we should think of them as barriers in operator entanglement, i.e. barriers for efficient simulation of the density matrix.
In this paper we address this question in three classes of models, the two dimensional rational CFTs, the random unitary circuit, and the holographic CFTs.They represent systems with increasing quantum chaoticity [29].All of the three models give a "growthplateau-drop" pattern as shown in Fig. 1(b).The plateau after the initial growth is of volume-law value, and thus is a barrier that prevents efficient simulation using the MPO representation.The results are briefly summarized in the following outline of the paper.
In section 2, we study the problem in 2-dimensional CFTs.The computation was performed in Ref. [30] under the implicit assumption that the theory is rational.We redo Total system has length .The Entanglement quantities are displayed in units of the corresponding state equilibrium entropy density s (n) eq in each system, see section 5 for more detailed explanations.
the computation and allow the operator to be the powers of the reduced density matrix, i.e. |ρ α AB , α ∈ N. We find that the story is similar to what is raised up in Ref. [31]: the rational CFT assumption entails additional OPE regimes among wist operators, and a quasi-particle interpretation which gives the result in Fig. 1(b).For irrational CFTs, the corresponding OPE regime disappears and explicit calculations become difficult.
In section 3, we perform the computation in a random unitary circuit using tools developed in Refs.[32][33][34].The random unitary circuit is a recently proposed quantum circuit model that is believed to capture the universal features of chaotic dynamics with local interactions.The circuit averaged entanglement in this model is mapped to a statistical mechanical problem of domain walls.At different stages of the evolution, three different types of the domain wall configurations dominate, which gives rise to the three continuous segments in Fig. 1(b).
The holographic CFTs occupy the maximally chaotic corner of irrational CFTs, in which explicit calculations become plausible again via AdS/CFT.In section 4, we study the dynamics of entanglement entropy in holographic CFTs for the operator state |ρ 1/2 AB .This coincides with the reflected entropy considered in Ref. [35], where it was found to closely relate to the entanglement wedge cross section [36] in the static case.Using the picture derived in [35], the dynamics of operator entanglement in this case is driven by the interplay between two HRT surfaces homologous to AB and A respectively in the bulk of the original state.The three different dominant configurations of the HRT surfaces then give rise to the three phases for the operator entanglement.In particular, the holographic plot in Fig. 1(b) is obtained from explicit calculations for AdS 3 /CFT 2 in appendix A, which extends the results of section 2 into regimes inaccessible by OPE analysis.
While in all three models there exists a plateau barrier phase characterized by volumelaw operator entanglement, thus denying the possibility of efficiently simulating these density matrices throughout thermalization, the length/persistence of these barriers seem to be correlated with how chaotic the systems are, suggesting it as another measure for quantum chaoticity.We elaborate this point further in the discussion section 5.

Two Dimensional CFTs
In this section, we investigate the operator Rényi entropy in two dimensional conformal field theories defined on an infinite line.The partition into subsystems ABC is illustrated in Fig. 2 .
Figure 2: The space partition in the 1+1d CFT.We take subsystem A, B to be contiguous with size A and B .In general we take A to be smaller: A ≤ B .The total system size is taken to be much larger than = A + B , so that thermalization will take place at late time.

Replica structure and twist operators
When computing state Rényi entropy, one can work in an orbifold theory by introducing pairs of twist operators that is equivalent to the replica structure [37,38].In two dimensional CFTs the computation becomes especially tractable since the twist operators become quasi-local with known conformal dimensions.This has resulted in many progresses for computing state entanglements in 2d CFTs both in the ground state and quench scenarios [6,38].
The standard procedure to compute the state entropy is to represent the state by a Euclidean path-integral with an open boundary.The reduced density ρ A is then obtained by sewing along the compliment A c .Quantities such as tr(ρ n A ) then connect neighboring copies of ρ A by gluing the upper and lower rims of the open slits along A, leading to a partition function on a branched manifold of the form a staircase geometry.This is equivalent to the insertion of twist fields at the entanglement cut.
The operator Rényi entropy works in an analogous picture albeit with more complicated permutation structures.We consider the more general state |ρ α AB / ρ α AB , α ∈ N, i.e. operator state associated with the α-th power of the reduced density matrix.
The reduced density matrix ρ AB is represented by a path-integral with a slit along A ∪ B. The power of the reduced density matrix ρ α AB is given by gluing the upper and lower rims of the slits in a staggered manner out of α copies of ρ AB , thus leaving two rims, one from the first copy and one from the α-th copy.Placing its complex conjugate produces the unnormalized operator density matrix, see Fig. 3(a).
To compute the corresponding n-th operator Rényi entropy, we generalize Eq. (1.5): . . .It proceeds in a few steps.The first term in Eq. (2.1) involves partial tracing in Fig. 3(b) over B, i.e. gluing along B within each operator density matrix.The resulting path-integral is then replicated n times, whose slits along A are then glued staggeringly, see Fig. 3(b).The second term in Eq. (2.1) comes from the norm of the operator state, see Fig. 3(c), where the trace is performed in A ∪ B.
According to the descriptions above, one can work out the equivalent twist fields to be inserted (more details can be referred to Ref. [35]): Fig. 3(d) and Fig. 3(e) demonstrate the twist field insertions that reproduce Fig. 3(b) and Fig. 3(c) respectively.Our set up corresponds to the limit A and B being adjacent to each other, thus the two twist operators in the middle are fused via OPE to give Define the scaling dimension of size n cyclic permutation to be h n , then the dimensions of the twist fields are:

Boundary state set up
We set up the quench dynamics by considering the following initial state: where |B is a conformal boundary state.The imaginary evolution makes the state normalizable with short-range entanglement at length of order β.We follow Cardy and Calabrese [38] to set up the calculation for operator entanglement.This has been done for the operator state |ρ AB / ρ AB in [30].We repeat the analysis of [30] for general α ∈ N, and compute S op n (A, ρ α AB (t)).We shall emphasize that the results are only valid for rational CFTs.A similar problem also occurs in the two interval entanglement [31], where the behaviors of the four point function distinguish the integrable and chaotic models.We will clarify this as we proceed.
For CFTs, it is easier to first wick rotate to Euclidean time t → −iτ : and analytically continue back later.Both |ψ(τ ) and ψ(τ )| can be constructed from Euclidean path integrals with boundaries.The procedure for computing the operator Rényi entropy discussed in subsection 2.1 simply amounts to inserting the appropriate twist operators when gluing the two path integrals for |ψ(τ ) and ψ(τ )|.The setup is shown in Fig. 4: we have a strip geometry separated by the thermal length and three twists fields inserted for the unnormalized state |ρ AB and two twist fields for the normalization ρ α AB (Fig. 4).
The operator entanglement we are interested is the ratio of correlation functions where We proceed by conformally mapping to the upper half plane using The conformal boundary condition of |B then means that the correlation functions satisfy the same ward-identities as holomorphic correlation functions on the full complex plane, but with doubled number of points coming from the images: where z 4 = z1 , z 5 = z1 , z 6 = z3 are the images z 1 , z 2 and z 3 .When continued to real time, their coordinates become: (2.10)

Evaluation via OPE analysis
In taking the high temperature limit β → 0, the holomorphic conformal invariants of the correlation functions are pushed towards the boundary of moduli space for generic t.As a result the evaluation of correlation functions are controlled by particular OPE channels and the corresponding singularities, provided that they exist.In this subsection we proceed by assuming that all singularities arising from "chiral" OPEs exist, a point we shall come back to comment in the next subsection.As t evolves, the conformal invariants bounce from one boundary of moduli space to another, and the dominant OPE channels switch accordingly.As we summarize in the following, there are 4 phases of Eq. (2.9) with different OPE structure (see Fig. 5), they correspond to the segments of the "growth-plateau-drop" pattern in Fig. 1(b), also see review in Ref. [31] 2.3.1 0 < t < A /2 During this regime, the OPE structure is as follows: The fusion proceeds in the following steps: for the numerator and for the denominator.Thus we can evaluate the numerator and denominator respectively as: under which Eq. 2.9 evaluates to It gives the linear growth with the rate same as that of the state entanglement.

2.3.2
A /2 < t < B /2 In this regime, the OPE structure switches to: Now the fusion proceeds in the following steps: for the numerator.The denominator stays the same.The numerator changes to: Eq. 2.9 now evaluates to This gives the operator entanglement: It contains a time-independent volume-law term, and corresponds to the plateau barrier.

2.3.3
B /2 < t < /2 In this regime, the OPE structure switches to: Now the fusion proceeds in the following steps: for the numerator, while the denominator remains unchanged.So the numerator further changes to: Eq. 2.9 now evaluates to This gives the operator entanglement: which corresponds to the drop phase.

t > /2
In this last regime, the OPE structure settles down to: The fusion process decompose into two clusters: In this phase, the denominator also undergoes a phase transition, whose fusion process now takes the form: The numerator and denominator is then evaluated to be: Eq. 2.9 now evaluates to (2.25)This gives the short-range operator entanglement: The results for all four regimes are thus summarized as below (see Fig. 6)

OPE singularities and quasi-particle picture
The 4 segments of the results are fixed by the 4 "holomorphic" fusion channels of the 6-point correlation function Eq. (2.9), see Fig. 5, and the corresponding singular behaviors therein.However, as was pointed out in [31], in generic CFTs the class of singular behaviors is only a subset of those produced by chiral OPEs as in the last subsection.The holomorphic correlation function with images inserted only captures the conformal transformation property of the boundary state; the OPE occurs in the UHP and is controlled by both the holomorphic and anti-holomorphic coordinates (z, z).In particular, there are only two types of OPEs in the UHP: The first line Eq.(2.28) corresponds to ordinary OPE (z 1 → z 2 , z1 → z2 ) in the interior of UHP; the second line Eq.(2.29) corresponds to the dominant channel when O 1 approaches the boundary z 1 ↔ z1 , and fuses into a boundary operator (identity in this case).
One can check that of the 4 channels shown in Fig. 5, three of them consists solely of legitimate OPEs in the UHP and therefore whose singularities exist in generic CFTs: On the other hand, the fusion channel in the third phase ( B /2 < t < /2) of Fig. 5 cannot be decomposed into a series of legitimate OPEs in the UHP.Therefore the corresponding singularities giving rise to the linear decline phase in Fig. 6 do not exist in generic CFTs.It is only in rational CFTs that the existence of these singularities is necessary in order to satisfy the crossing symmetry of correlation function, see [31].We thus conclude that the linear decline behavior for the phase ( B /2 < t < /2) is a consequence of working with rational CFTs.
The nature of these singularities characterizing rational CFTs becomes more transparent if instead of the boundary state |B , we consider the thermal field double (TFD) state defined in the doubled CFT 1 ×CFT 2 : where |n are eigenstates of H.The reduced density matrix is thus defined on (AB) 1 ×(AB) 2 , and the corresponding operator state is defined on a further double copy of (AB) 1 ×(AB) 2 .The stripe geometry in Fig. 4 becomes a cylinder of circumference β, with another set of twist operators inserted at the images.The operator Rényi entropy is then related to the full 6-point correlation function on the complex plane: As conjectured and argued in [31], the TFD state can be used to reproduce the singularities in the boundary state.In this case, the singularities in ( B /2 < t < /2) takes the form of two "tangled" light cone singularities.To see what this means, we compute the following two conformal invariants: The first line Eq.(2.32) indicates that (z 2 , z2 ) is approaching the light cone "tip" of 1 and 3, i.e. simultaneously approaching the light cones of both 1 and 3; the second line Eq.(2.33) indicates that (z 6 , z6 ) is approaching the light cone tip of 3 and 5.Each one corresponds to a light cone singularity discussed in [31], see Fig. 7.However, we cannot isolate them into two separate "clusters" as in [31]; they are tangled together, and this makes it impossible to factorize the 6-point function into a product of lower-point functions.In principle, such light cone limits of correlation functions in chaotic 2d CFTs, including finite c corrections, can be explicitly studied via the structure of fusion matrix [39][40][41].We leave this for future studies.The time dependence in Eq. (2.27) describes rational CFTs.As a consistency check, it can be exactly reproduced by the quasi-particle picture, which is believed to model the essential physics regarding entanglement evolution in integrable models.We present this in the thermal field double set up, see Fig. 8.The initial state is represented by the middle line.After the quench, there are coherent pairs of quasi-particles generated and they move to the opposite directions.As time goes on, the overlap regions on the line, which are the area that emitted quasi-particles at time t that now entangled between region A and B , is proportional to the operator entanglement.It gives rise to the 4 segments in Eq. (2.27) and the symmetric trapezoid in Fig. 6.

A B A B
t

Random Unitary Circuit
In this section, we introduce the random unitary circuit model and the corresponding domain wall picture for the entanglement.Using this, we compute and interpret the change of the quench operator entanglement as the succession of three dominant domain wall configurations.
Random unitary circuit is a tensor network that models chaotic evolution with local interactions.The structure of the network is shown in Fig. 9(a), where the horizontal direction represents the system: a one-dimensional chain with q bits of local degrees of freedom.The dynamics evolves along the vertical direction.Each block is an q × q independent Haar random unitary matrix.The network is constructed to involve only local interactions that entangle (almost maximally) nearby sites in one time step.And the locations switch by one lattice spacing in alternating steps as the system evolves.The circuit retains minimally the features of evolution by a local Hamiltonian -the local generation and propagation of large entanglement -yet dispenses particular forms of interactions.The average over quenched randomness is what provides analytic handles on the dynamics of chaos.Recently, there have been fruitful results based on this model and its variants, about the scaling behaviors of out-of-time-order correlator [42][43][44], entanglement growth [32-34, 45, 46], hydrodynamic long tails in U (1) conserved systems [44], etc.Similar non-unitary random circuit were also proposed to reconstruct the holographic space time as an arguably more explicit tensor network structure.The Ryu-Takayanagi (RT) surface of entanglement, transition of entanglement scalings, etc. were considered [47][48][49] in this setting.
Among these results, the random unitary circuit enables a "domain wall" picture to understand the universal behavior of entanglement.It is motivated by the minimal cut picture of a general tensor network.In this picture, the number of bonds the minimal cut traverses through the tensor network gives an upper bound of the bond dimension to represent the states, and therefore is also an upper bound for the entanglement [50].The domain wall picture is a refinement of the minimal cut by taking unitarity into account [33].They can be viewed as world-lines of objects carrying permutation data.Unitarity puts causality constraints on their trajectories and possible forms of interactions.More explicitly, computing the circuit averaged entanglement is mapped to a statistical mechanical problem of domain walls, where the minimal free energy gives the entanglement.
The objects carrying permutation data can be understood from the underlying replica structure defined by tr(ρ n A ), see Fig. 9(c) for n = 2.They perform permutation operations at the entanglement cut, which is analogous to the twist fields in the field theory definition of tr(ρ n A ) using path integrals.In random unitary circuit, the random average yields an effective statistical mechanics model, whose solutions extend the permutations from the entanglement cut into the interior of the tensor network.In this sense, the domain wall gives a geometric picture to the entanglement of the random unitary circuit, in much the same way as the RT surface [51] in holography.In fact, domain walls of permutation in non-unitary random circuit were formulated to be the tensor network analog of the RT surface in attempts [47,48] to understand the geometric structure of the space time.In a different language, the formulation has also been studied in the context of machine learning [52].

The effective model
We now review the effective statistical mechanical model in the random unitary circuit.Let us begin by looking at the interior of the tensor network of the quantities tr(ρ n A ) shown in Fig. 10(a), for illustration we pick n = 2.The overlapped blue and red squares are n = 2 copies of the Haar random matrix U q×q and its complex conjugate U * q×q .Each block with U (blue in Fig. 10(a)) and U * (red in Fig. 10(a)) is replicated n times in the network.Random average gives rise to "contractions" between U and U * among the copies.Each contraction can be represented as a permutation element in the n-th symmetric group, see Fig. 10(b).The permutation elements are what emerge as the dynamical degrees of freedom after random average, and we call them "spins".The effective model is represented by the statistical mechanics of these spins4 .Boundary conditions are imposed at both the top and bottom boundaries of the tensor network.They play different roles: the boundary condition on the top comes from the replica structure for computing the Rényi entropy; the boundary condition on the bottom encodes the initial state, see Fig. 9(c).(12) (12) For the purpose of studying entanglement, it is more convenient to describe the configurations in terms of domain walls rather than spins, because the former fits better with the boundary conditions to be imposed and unitarity constraints.Different spin domains are characterized by different channels of contractions between the n copies of U and U * , the domain walls can be understood as the twists of those contractions.For example, in order to evaluate the n-th Rényi entropy, the contractions operation on the top boundary of Fig. 9(c) are identities acting on the subsystem A, and permutations (12 • • • n) acting on the subsystem B respectively.The domain wall is a twist between them, inserted at the entanglement cut.If the spins in the left domain are all σ 1 and those in the right domain σ 2 , then the domain wall is a permutation τ = σ −1 1 σ 2 .In a general tensor network, these permutations may become other tensors when propagating from the top boundary into the interior.But in random unitary circuits, the domain walls (or permutations) are the only degrees of freedom that survive the random average.For instance the average purity tr ρ 2 A will only contain wandering domain walls as in Fig. 11.The Rényi entropy thus is given by the freedom energy of the domain walls subject to the boundary conditions at both the top and bottom boundaries.For the product initial state we are considering, we have free boundary condition at the bottom.
This picture can be extended to more generic entropic quantities such as the operator Rényi entropy we are interested, which corresponds to imposing more complicated boundary conditions, see Fig. 12.In this case, the top boundary can host more generic domain walls.They can decomposed into elementary domain walls, which are defined to consist only of single transpositions.The elementary domain walls are the building blocks for the underlying dynamics.For example, the domain wall (123) consists of two elementary domain walls.The locality of the interaction and unitarity of the evolution restricts the allowed behaviors of elementary domain walls (see Fig. 13).When propagating downward, an elementary domain wall must remain inside the light cone and can not go back.When two elementary domain walls meet, they combine according to the group multiplication rule; and the same applies when a composite domain wall split.
(e) As mentioned, the domain walls need to minimize the free energy while satisfying the boundary conditions imposed at the top and bottom of the circuit.An elementary domain wall will cost s eq = ln q amount of free energy for each discrete time step it takes.Therefore at the leading order in large q expansion, the minimizing configuration is simply given by free elementary domain walls taking the shortest paths 5 .At higher orders in large q, elementary domain walls can have two types of interactions.One is the statistical interaction between two non-commutative domain walls, such as those with transpositions ( 12) and (13).It corrects the entanglement velocity in order O(1/ log q).For the special case of Rényi index n = 2, this effect is absent.Another weaker interaction at O(1/q 4 ln q) order appears between commutative domain walls like (12) and (34).The suppression of domain wall interaction at q = ∞ effectively reduce the quench disorder to annealed disorder.We can thus equating e −(n−1)Sn and e −(n−1)Sn .

Domain walls in operator entanglement
In this section, we study time evolution for the operator entanglement of the random unitary circuit.We proceed by listing the domain wall configurations that minimize the free energy at different stages of the evolution.We take n = 2 to simplify the domain wall structure, and work at leading order in q → ∞ so that the elementary domain walls can be treated independently.This brings in the reduction where the two terms corresponds have the initial domain wall configurations shown in Fig. 14.
To warm up, let us first consider the evolution of the 2nd Rényi entropy of the state, which is subtracted in Eq. (3.1).When t is small, the two domain walls at the entanglement cut cost the same amount of energy when going down in any direction within the light cone.Due to entropy consideration [32,33], macroscopically they go down vertically.The free energy increases linearly with rate v E ln q ≈ ln q.At time t > 2 , the two domain walls can meet and annihilate at the intersection points of their light cones.The free energy of this configuration is ln q, which will be the minimal when t 2 .This switch from the configuration 1 to configuration 2 in Fig. 14(b) gives rise to the linear growth and saturation behaviors of entanglement in a quenched state.
On the other hand, the first term in Eq. (3.1) has six elementary domain walls to begin with.Three types of (macroscopic) configurations will occur in their life time.Configuration 1 is the beginning stage, where all domain walls going vertically down6 when t is small.The growth rate is 6 ln q.In configuration 2, the four domain walls from the boundaries of region A first meet and fuse into two domain wall.These two go down vertically together with the other two domain walls, giving rise to a growth rate 4 ln q.In configuration 3 the middle domain walls split into two domain walls on the left and right, which then annihilate the remaining domain walls, giving the saturation value 2 ln q = 2S 2 (AB).By equating their free energies, we can determine the respective transition time: linear growth saturation .Their free energy is twice of a single pair of domain walls, due to their independence in the large q limit.The two configurations for the two phases of state quench entanglement are shown on the right of the brace: In the linear growth phase, the two domain walls microscopically going down vertically; in the saturation phase, the two domain walls meet and annihilate.
t In Fig. 15(a) we draw the time evolution of both terms in Eq. (3.1).The 1st term grows faster than the 2nd term in configuration 1 and then at the same rate in configuration 2. It eventually come to a tie with the state entanglement when reaching configuration 3. The operator entanglement is the difference between these two curves, which we draw in Fig. 15(b).It contains three stages of development: linear growth, plateau and linear decrease.
We can generalize these results to integer n and α assuming the q → ∞ limit.The domain wall at the top boundary becomes identical to the permutation data in the twist fields in Eq. (2.2).There are again two terms The first term is a generalization of Fig. 14(a): on the top boundary, both the cuts between A,C and B,C host (2α − 1)n elementary domain walls, the cut between A, B host 2(n − 1) elementary domain walls.Its dynamics also contain the 3 configurations in Fig. 14(a).
The second term is a generalization of Fig. 14(b), the difference is only that the number of elementary domain walls is (2α − 1)n rather than 2 at the entanglement cut.From the geometry, one can see immediately that the transition time scales do not dependent on n and α, as long as the order of the configurations to appear is the same as n = 2, α = 1.This can be confirmed by an explicit computation (3.4) Further computation shows that the plateau value of S op n (A, ρ α AB (t)) is also n and α independent.
We then claim that for a given integer n ≥ 2 and α ≥ 1, the q → ∞ collapses all the operator entanglement S op n (A, ρ α AB (t)) to the curve in Fig. 15(b).The results above are derived for the large q random unitary circuit.In a realistic system or even random circuit with finite q, we argue that such a generalized domain wall picture still exist, and the we can incorporate the non-universal feature into a coarse grained line tension function -the free energy of the domain wall per unit length.
We can not justify this assumption to the full extent, even when restricted to chaotic systems.Instead we motivate this assumption from the following aspects.
Microscopically, the domain walls are the twist of the permutations.These permutations, as we draw in Fig. 10, represents different ways to contract among copies of U and U * .These pairwise contraction gives a real deterministic value given in this case by formula shown in Fig. 10(b).There will contributions other than these contractions if we take a unitary gate in Fig. 10(a).But these will usually be complex numbers, can consequently dephase in time in chaotic systems.This is also the reason why the domain wall picture works so well in the random average of the random unitary circuit: the random average makes the dephasing exact.
Macroscopically, we can generally assume that the entanglement growth of a given cut can be written as a integral of a local quantity in space time.In the examples where the minimal cut works, it would be the integral of the increment of entanglement along the cut.Because of the translational invariance in space and time (true in the random average sense in the random circuit), we imagine that such minimal cut should be a line (or membrane in higher dimension), and the growth rate should depend only on the slope of this line.This motivates the line tension function E n (v) for Rényi index n and slope v.
Ref. [33] motivates the existence of such line tension function by first assuming the entanglement growth rate ∂ t S to be dependent only the a function of the slope Γ(∂ x S) at the leading order in coarse graining.The assumption only uses the local information of the entanglement (the derivative) which is reasonable in a system with only local interactions.The authors then deduce that the line tension function is the Legendre transform of Γ(x).By constraints of entanglement (e.g.subadditivity imposed on Γ(x)), they further derive that E n (v) has to satisfy several properties where v B is the butterfly velocity that characterizes the speed of operator spreading and many-body chaos [33].When entanglement saturates in Fig. 15, the free energy cost is 2 × E n (v)/v, which is minimized at slope v = v B rather than v = 1, though the saturation value is still the maximal entanglement.
In the configurations we considered in Fig. 14, the domain walls either go vertically or at slope v B (which maximizes the line tension), we can then only use E and E(v B ) = v B to estimate the transition points.
The random unitary circuit at the n = 2 case is a known example that the line tension function works well in practice even for q = 2.This is because after random average, the finite q correction only introduces weak interactions between the domain walls in Fig. 14.In leading order, we have [32,42,43,45] where v E is a perturbative series containing higher O(1/q 8 ) term due to the interactions of the domain walls.
Using these data, we can determine the transition points.
We look at the symmetric case A = B .When q = 2, we have t 1 > t 2 , which means that configuration 2 will be skipped.There is only one transition time from configuration 1 to configuration 3 This is verified in the numerical results in Fig. 16, where = 28, A = B = 6 and q = 2.
As we can see the plateau phase found in large q is absent, instead only a peak occurs at about t 3 ≈ 12.5.It therefore suggests that the persistence of the plateau regime is a signature of the strongly chaotic system, and likely will only occur with large on-site degrees freedom.

Holography
In this section, we look at holographic CFTs which represent systems that are maximally chaotic in terms of saturation of the chaos bound [29].In AdS/CFT, the entanglement entropy is conveniently captured geometrically by the area of the RT/HRT surfaces [3,4].To study the time evolution for operator entanglement in holographic CFTs, our strategy is to start from a bulk description for operator states, whose asymptotic boundaries are twocopied subsystems, and then apply the HRT prescription.The results are obtained using related pictures derived in Ref. [35].We first review the bulk picture for the stationary operator states, and then generalize to the dynamical case for generic α.After that we study the operator entanglement entropy for α = 1/2, which is the case considered in [35], and turns out to be the natural case for constructing the bulk operator state.We will see that the operator entanglement of the reduced density matrix undergoes three phases in the bulk, which is partly driven by the background geometry representing the density matrix operator state.

Stationary states
We begin by introducing the prescription derived in Ref. [35] using the language of Euclidean path-integral, which works explicitly for stationary states, or more generally states that are time reflection symmetric under: t → −t.
The goal is to construct the bulk geometry dual to the operator state |ρ AB , associated with the reduced density matrix: We assume that the global state |ψ in the boundary CFT has a smooth bulk geometry.Furthermore, it is stationary with a reflection symmetry t → −t and thus admits a Euclidean path-integral definition.The reduced density matrix element between two field configurations φ 1,2 on AB is then given by: where Φ denotes collectively the relevant quantum fields, and f (Φ) denotes the statecreating sources for |ψ .The path integral is over configurations that has a branch-cut along AB, and the field configurations along the branch-cut are fixed to be φ 1 and φ 2 from above and below respectively.Alternatively, this defines the (un-normalized) wave-functional of the operator state |ρ AB on the doubled subsystem {AB, A B }: In order to find the bulk dual of |ρ AB , let us consider computing its norm: For stationary states this is given by a boundary Euclidean path integral over a boundary manifold ∂M doubly branched across AB (see Fig. 17).Through the dictionary of AdS/CFT, this fixes the boundary conditions for the corresponding Euclidean path-integral over bulk fields and geometries, which matches the CFT result.Let us further suppose that the bulk Euclidean path integral is dominated by a particular saddle point contribution, denoted by (M, g, ϕ) where g, ϕ are respectively the metric and the collection of bulk fields: For simplicity we omit writing ϕ in specifying bulk solutions from now on.On the other hand, path-integrals over such branched manifolds are related to the refined 2nd Renyientropy [53].For holographic theories and assume replica symmetry, the dominant saddle B can be constructed as in Ref. [53] and takes the form of two identical "wedges" (see Fig. 17): where each wedge B1,2 can be obtained by solving under the original single-copied boundary condition, the back-reacted geometry upon inserting a "cosmic brane" γ AB 2 that anchors at the boundary along the entangling surface ∂(AB), and has the tension This creates a conical deficit of angle: ∆φ = π, so each wedge extends an angular range       Σ(h) also contains the cosmic brane γ AB 2 .The reflection symmetry characterizing the stationary state acts by reflecting across Σ and interchanging the two wedges: B1 ↔ B2 .Semi-classically, we can interpret Σ(h) as the bulk spatial geometry at t = 0 dual to the operator state |ρ AB .
Let us now describe the prescription operationally (see Fig. 18).To find the bulk geometry dual to the stationary operator state |ρ AB , we start from the bulk dual of the original global state |ψ , insert a cosmic brane γ AB 2 of tension T = 1 4G N anchoring towards the boundary at ∂(AB), let it back-react and settle down.The spatial geometry dual to |ρ AB is then obtained by identifying two copies of the back-reacted spatial geometries across the settled-down position of the cosmic brane γ AB 2 .

Dynamical states
We can extend the prescription to the time-dependent case, i.e. |ψ is a dynamical state.In this case, the (un-normalized) density matrix elements, i.e. the wave-functional of the operator state, need to be constructed using the Schwinger-Keldysh formalism [54,55].At the instant t = t 0 , it involves a path-integral over two Lorentzian time-folds L + and L − , ranging from t = (−∞, t 0 ) and t = (t 0 , −∞) respectively and glued at t = t 0 across C (see Fig. 19) : where f denotes the state-creating sources.Now let us follow the previous route and compute the norm: where ∂M is a boundary (Lorentzian) manifold branched over two Schwinger-Keldysh contours [56], i.e. four time-folds L ± 1,2 , and glued at t = t 0 (see Fig. 19).The gluing creates a local singularity on the branched manifold, in the form of a non-standard Rindler temperature: T Rindler = 1/4π.
Similar to before, we can push this path-integral into the bulk via AdS/CFT, and suppose further that there exists a dominant contribution to this from the Lorentzian saddle-point geometry (B(t 0 ), g) [56,57]: Strictly speaking, the existence of such a dominant Lorentzian saddle stands on shakier ground compared to the Euclidean counterpart.We shall not dwell on this too much in this paper other than acknowledging the potential subtlety, and proceed assuming such a saddle exists.The cosmic-brane construction has a natural generalization to the dynamical case.In particular, let us still assume that the replica symmetry on the boundary, which in the dynamical case interchanges the two Schwinger-Keldysh two-folds, remains a symmetry of the bulk saddle.There exists a co-dimension 2 bulk hypersurface γ AB 2 of fixed points under this symmetry, which ends on the boundary at ∂(AB) at t = t 0 .The bulk saddle B(t 0 ) takes the form of two identical bulk Schwinger-Keldysh two-folds, glued at a bulk Cauchy surface Σ(t 0 ) that contains γ AB 2 (see Fig. 20):

m k D H w K R c = " > A A A C W 3 i c d Z F P S 8 M w G M b T z j + z T t 0 U T 1 6 C Q / A 0 W h X 0 5 t C L x w n u D 2 x l p G m 6 x S V N S V J h l H 0 H r / r N P P h d T L c e n K s v B B 5 + z / u S N 0 + C h F G l X f f L s i t b 2 z u 7 1 T 1 n v 3 Z w e F R v H P e U S C U m X S y Y k I M
A K c J o T L q a a k Y G i S S I B 4 z 0 g 9 l j 7 v f f i F R U x C 9 6 n h C f o 0 l M I 4 q R N q g 3 S q Z 0 7 I 3 r T b f l L g t u C q 8 Q T V B U Z 9

y w 7 k e h w C k n s c Y M K T X 0 3 E T 7 G Z K a Y k Y W z i h V J E F 4 h i Z k a G S M O F F + t l x 3 A S 8 M C W E k p D m x h k v 6 e y I L A m 6 6
O N J T t c Y R V 2 r O g / + 8 H J Z 6 o c o v K r W U 2 X 5 K w l I v J 1 o I p s r c Y a q j O z + j c Z J q E u P V q 6 K U Q S 1 g H j Q M q S R Y s 7 k R C E t q g o F 4 i i T C 2 n y H Y y L 3 / g a 8 K X p X L e + 6 d f V 8 0 2 w / F O F X w R k 4 B 5 f A A 7 e g D Z 5 A B 3 Q B B q / g H X y A T + v b r t i O X V u 1 2 l Y x c w L W y j 7 9 A X U Q t m I = < / l a t e x i t >   Furthermore, being a saddle point, the bulk geometry near γ AB 2 should be locally smooth.For Lorentzian manifolds this implies that the local Rindler temperature takes the standard value T Rindler = 1/2π.This refers to the branched four-folds, and as a result, each of the two-folds B1,2 (t 0 ) = B+ 1,2 (t 0 ) ∪ B− 1,2 (t 0 ) has a local Rindler temperature near γ AB 2 of T Rindler = 1/π.This is the back-reacted bulk solution obtained by inserting a cosmic brane of tension T 2 = 1 8G N ending on ∂(AB) at t = t 0 , with the boundary being the Lorentzian two-fold L + ∪ L − that defines the density matrix in Eq. (4.9).The settled-down location of the cosmic brane is γ AB 2 .The bulk solution B1 (t 0 ) obtained this way then defines the leading order description of the bulk state dual to the operator state |ρ AB (t ) .We can understand the instantaneous spatial geometry Σ(t 0 ) as being obtained by identifying two copies of the back-reacted spatial geometries ending on AB at t = t 0 across the cosmic brane γ AB 2 .It is worthwhile to clarify the notions of dynamics.The operator state |ρ AB (t 0 ) "evolves" as we change t 0 .However, this evolution is not unitary, as manifested by the t 0 -dependent norm N (t 0 ).As a result, such t 0 -dependence is not encoded by a single Lorentzian bulk solution; instead, t 0 serves as a "label" for the one-parameter family of bulk saddles B1 (t 0 ), t 0 ∈ R .On the other hand, for a fixed t 0 one can evolve the state |ρ AB (t 0 , s) = U (s)|ρ AB (t 0 ) using some unitary operator U (s) and study the dynamics in

< l a t e x i t s h a 1 _ b a s e 6 4 = " q F M W z E S u R + a R Y J r p X F U 7 J x L Y h U I = " > A A A C W 3 i c d Z F P S 8 M w G M b T z j + z T t 0 U T 1 6 C Q / A 0 i n o z a E X j x P c H 9 j K S N N 0 j S Z N S V J h j H 0 H r / r N P P h d T L c e n K s v B B 5 + z / u S
, where K ψ AB is the half-side operator modular Hamiltonian, and is related to the original state Hamiltonian K Ψ AB = − ln ρ AB by: Notice that U (s) is not the full modular-flow, which acts as a symmetry of the operator state.In the vicinity of γ AB 2 (t 0 ), the action of U (s) can be approximated geometrically by two identical rindler boosts about γ AB 2 (t 0 ), with s being V W r l 5 j 1 7 5 A e s 0 0 9 / / L + F < / l a t e x i t > B 1 (t 0 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " M O p Z y P T Z d W t h X f H y E u l q P q v x 8 A w = "     the rapidity [58,59].This effectively generates a time-like evolution in the near-γ AB 2 (t 0 ) region of the bulk operator state.

Generalization to |ρ α
AB (t) One can generalize the construction discussed before, and consider operator states of the form considered in section 2: where we take α as an integer to begin with.By making the same assumptions (e.g.replica symmetry, existence of a dominant saddle, etc) it is easy to see that the construction for the bulk dual proceeds the same, with the only change being the tension of cosmic-brane: Next, we can analytically continue α away from integral values.A special case occurs at the value α = 1 2 .In this case T α = 0, i.e. the cosmic brane is tension-less, and the location of the brane is given by the HRT surface in the original geometry [4].More importantly, since the cosmic brane does not generate time-dependent back-reaction, the "parent" bulk geometry for constructing the dual of the operator state |ρ 1 2 AB (t) remain the same as we evolve in t.The only t-dependent ingredient is the HRT surface across which we make the identification.AB (t) , where we take the global state |ψ to be pure B-state e −Hβ/2 |B .In a holographic CFT d , this is dual to half side of a two-sided external black hole with the metric: Furthermore it contains an end-of-the-world (EoW) brane cutting through the interior [60].
The EoW brane can be understood as arising from extending the boundary condition on the CFT into the bulk [61].A particular type of boundary condition in the bulk is that the HRT suface can end perpendicular onto the EoW brane.We are interested in the operator entanglement entropy: This coincides with the reflective entropy considered in Ref. [35], where for static cases it was shown to be equal to twice of the entanglement wedge cross section, the latter has been conjectured as a holographic realization of the entanglement of purification [36].Via the HRT prescription this is obtained by finding the extremal surface γ A (t) in bulk geometry dual to the operator state.For α = 1 2 , the bulk dual is given by the original dual (i.e.half-sided eternal black brane with the EoW brane) identified across the the HRT surface γ AB (t).This allows us to extract S op A (t) directly in the original bulk geometry by modifying the rule of finding the HRT surface γ A (t) as follows: • A HRT surface γ A (t) in the original geometry remains a candidate HRT surface • A candidate HRT surface can also end perpendicularly on γ AB (t).

• S op
A (t) is given by twice the smallest area among all candidate HRT surfaces.In other words, S op A (t) in this case equals twice of the dynamical version of the entanglement wedge cross section for A inside AB.Let us focus on the case where both AB and A are flat stripes, ranging between (− A , B ) and (− A , 0) respectively along the spatial y direction, where 0 The first step is to determine the HRT surface γ AB (t) in the parent one-sided black brane geometry, on which the HRT surfaces γ A (t) of the operator state can end.The time dependence of γ AB (t) has been studied in [16], and exhibits two phases in the high temperature limit β 1: • For v d t < /2, γ AB (t) takes the form of two planes at fixed values of y = − A and B respectively, extending into the bulk and ending perpendicularly on the EoW brane, with the total proper area growing with time Area γ AB (t) ∝ 8π β v d t; • For v d t > /2, γ AB (t) is a connected extremal surface that extends in the y direction, most of which lie very close to the bifurcating surface ρ = 0 and has a volume-law proper area Area γ AB (t) ∝ 4π β .
where the velocity [2(d−1)] 1−1/d depends on the boundary CFT dimension d [16].On top of this evolution, we expect three phases for the behavior of the HRT surface γ A (t), summarized as in Fig. 21: • For v d t < f A with a proportionality constant f to be determined, γ A (t) is similar to each component of γ AB (t) at the same time, and takes the form of a plane at fixed y = 0, extending into the bulk and ending perpendicularly on the EoW brane.During this phase, the entanglement entropy grows linearly in time: • For f A < v d t < /2, γ A (t) latches onto one of the components of γ AB (t).During this phase the entanglement entropy is approximately time-independent and exhibits volume-law scaling: • For v d t > /2, due to the transition of γ AB (t), γ A (t) undergoes a discontinuous transition, changing from volume-scaling to being short-ranged.A (t) grows linearly in t; Center: for f LA < v d t < /2, after a phase transition for γ A (t) itself, S op A (t) settles down to be volume-law; Right: for v d t > /2, due to the phase transition of the underlying bulk state, i.e. the identification surface γ AB (t), γ A (t) is forced to jump discontinuously to be short-ranged. t AB (t)) As we can see, the decline phase is absent for holographic theories, instead the plateau barrier persists until the very end of thermalization.In appendix (A), we calculate the operator entanglement S op A (t) explicitly in AdS 3 /CFT 2 , where v 2 = 1.After subtracting the usual UV divergent contribution, the results are summarized below, and plotted in Fig. 22: This can be compared with the rational CFT results by sending n → 1, α → 1/2 in Eq. (2.27).In particular, we can identify explicitly that f = 1/2, agreeing with the CFT results.We make a comment about this.The constant f is determined by the HRT surface γ A (t) in the second phase, which has a proper area that scales with the volume ∼ A , see Fig. 21.Naïvely an HRT surface that scales with volume-law would have its most part sticking to the bifurcating surface, like (half of) γ AB (t) in the third phase, see Fig. 21.Let us call surfaces like this the thermal type.It is easy to see that the thermal type would give the incorrect answer f = 1.
So γ A (t) has smaller proper area than the thermal type by a factor of f , which is 1/2 in AdS 3 /CFT 2 .Geometrically this can be understood as follows.The thermal type for γ A (t) would emerge if we extremize all surfaces that end perpendicularly on the full time-like hyperplane y = − A ; while the actual HRT surface comes from extremizing those that end perpendicularly on γ AB (t), which is a particular Cauchy slice of the y = − A hyperplane.The extremizing procedure yielding the thermal type can be represented by a maxmin procedure [62]: foliate the relevant spacetime region into Cauchy surfaces; minimize on each Cauchy surface; then maximize the minimal values over all Cauchy surfaces.In this picture, the actual HRT surface corresponds to a particular minimal configuration on the Cauchy surface whose boundary is γ AB (t).In the maxmin construction, it has smaller proper area than the thermal type.
As a result, γ A (t) does not stick to the bifurcating surface, and in fact passes through the horizon.One can then investigate what fraction of γ A (t) is in the interior of the black brane.In appendix A.4 we check this explicitly, and found that the entire volume-scaling part of γ A (t) is from behind the horizon, see Fig. 23.More explicitly, the portion of of γ A (t) that is behind the horizon covers the interval in y coordinate: Such an extensive amount of entanglement from the black hole interior is only accessible from the operator entanglement.

Discussion
In this paper, we study the operator entanglement of the reduced density matrix of a quenched state in three representative systems: random unitary circuit with local interactions, rational CFT and holography.We find that in all those cases, the operator entanglement grows linearly from short-range value to an extensive value, stays constant for B A t < l a t e x i t s h a 1 _ b a s e 6 4 = " f F 6 q + 6 O x o t X L Z T B k t l X X 6 A 4 2 / u T 0 = < / l a t e x i t > `A < l a t e x i t s h a 1 _ b a s e 6 4 = " z P / J a V h 8 8 h m b R f g L x s C w N Q U A g b E = " > A A A C X H i c d Z H N S s N A F I U n 8 a + 2 V l s F N 2 6 C R X B j S a q g O 6 t u X F a w P 9 C G M p n c t E M n m T A z E U r o Q 7 j V J 3 P j s z h p s 7 A 2 X h g 4 f O d e 5 s 4 Z L 2 Z U K t v + M s y t 7 Z 3 d v d J + u X J Q P T y q 1 Y 9 7 k i e C Q J d w x s X A w x I Y j a C r q G I w i A X g 0 G P Q 9 2 Z P m d 9 / A y E p j 1 7 V P A Y 3  some time thus forming a plateau, and then drops to short-range value.These three phases has been summarized in Fig. 1(b), which we reproduce below in Fig. 24 for the readers' convenience.To answer the question raised in the introduction as motivation, our results suggest that for generic non-integrable systems such as the random unitary circuits or the holographic CFTs, there exists an operator entanglement barrier of volume-law value during the thermalization of a quenched state.Hence algorithms that use a low entanglement representation can not faithfully reproduce the reduced density matrix throughout the process.This is a natural expectation for chaotic systems.On the other hand, the rational CFTs also share the operator entanglement barrier despite the evolution Hamiltonian being integrable.The reason lies in our choice of the initial global quenched state |ψ = e −βH/4 |B , which is a conformal invariant boundary state perturbed by the energy density operator T tt .Such perturbation always exists in CFTs, and it has been shown in [63,64] that the reduced density matrix ρ A (t) under the global quench converges to the thermal density matrix 1 Z e −βH as t → ∞. 7 The fate of thermalization for the initial states we considered puts the rational CFTs in equal footing as the other two chaotic models for comparison.However, it is also known that other irrelevant perturbations of the boundary state (for example by the descendants of the stress tensor) could lead to a Generalized Gibbs ensemble (GGE) [65] at late times [63].The breakdown of thermalization could evade the decay phase, unless the chemical potentials are all tuned to be small.It would be interesting to study how the operator entanglement behave by quenching from the general initial states.
Strictly speaking, the explicit computations we performed in different models are for different operator entanglement quantities: in rational CFT and random unitary circuit we obtained the results for general n-th operator Rényi entropy, and for operator state with general α; in holographic CFTs what we computed is the reflective entropy, i.e. operator entanglement entropy n → 1 with α = 1 2 .That these results plotted on the same figure are showing the same plateau values therefore requires some technical explanation.First of all, the time scales of the systems are calibrated so that the effective speeds of light controlling the causality of local interactions are all set to 1. Furthermore, we have rescaled the entanglement quantities by their corresponding n-dependent equilibrium Rényi entropy density.More specifically, what we plot is the ratio S op n (A, ρ α AB (t))/( s eq ), where s (n) eq is the equilibrium Rényi entropy density of the quenched state at t → ∞.More importantly, one might worry that transition time could depend on the Rényi index n, and therefore the comparison between different models in Fig. 24 may be meaningless.We acknowledge this potential subtlety.However we point out that for integer n and α, both the CFT and random unitary circuit results are explicitly given: they are independent of α and n only enters as a pre-factor s (n) eq (in fact, it is also n independent in random circuit when q → ∞. 8 ).Therefore one can analytically continue them to n = 1, α = half and compare with holographic CFTs, and this should agree with Fig. 24.However, an honest comparison for generic n would require considering the cosmic-brane back-reactions in holography, which is beyond our current computational control.It is possible that Fig. 24 fails to accurately capture the case for generic n (for example, the transition to the decay phase in holography maybe continuous due to back-reactions.).We leave such investigations for the future.Now let us discuss the physics behind the three segments of the "growth-plateau-decay" pattern.
During the linear growth phase, we observe that the growth rates (in units of the equilibrium entropy density) are the same among the three models.This is very similar the the linear growth behaviors of quenched state entanglement, where the rates are also equal among the models despite different entanglement dynamics.However, we should caution that the linear growth of the operator entanglement is not solely a consequence of the corresponding state entanglement increase for the subsystem A. During this phase, the operator state |ρ α AB (t) ∈ H ⊗ H is building up entanglement between the subsystems A ∪ B and A ∪ B , thus is drifting away from the form of a product state in H ⊗ H.As a result, the operator entanglement growth rate is not twice of, but is equal to, the state rate.The linear growth stops at t = A /2.This coincides with the thermalization time for A, after which the operator entanglement barrier emerges.This suggests that the transition between the two phases is possibly related to the thermalization of region A.
The plateau phase represents an intriguing "equilibrium" in the course towards thermalization.Despite the operator entanglement taking constant value, the operator state is certainly not static, for example the state Rényi entropy of AB is still linearly increasing.This excludes the possibility of describing the reduced density matrix by any kind of static (equilibrium) ensemble.We can ask a few questions regarding this phase.Firstly, is the calibrated plateau value universal?For example, could it represent the saturation to the maximal operator entanglement subject to the given constraints on the density matrix?Secondly, how universal is the plateau phase?For example, does it survive in chaotic systems away from the large q or large N limits we are considering?And finally, what are the implications for the observables?For example, does it implies any properties in the local correlation functions within region A?
The existence of the decay phase for operator entanglement is dictated by the quenched thermalization of these models.In the high temperature limit, the thermalized reduced density matrix is approximately proportional to identity.Viewed as an operator state, it is therefore short-range entangled with correlation length of order β.However, the specific behaviors are highly model dependent; in fact, the models we study have different mechanisms for triggering the decay phase.
The linear growth phase and the early plateau phase are the same among all three classes of models.It is the point of exit from the plateau and the subsequent decay phase that serve to differentiate among them.A similar story for the two-interval entanglement/mutual information in the quenched state is discussed in [31], where it is also found to encode the distinction between the chaotic and rational CFTs.However, if we consider the same quantity for the random unitary circuits, at large q in large separation it exhibits the same behavior as the holographic CFTs9 .In contrast, by probing the later time behavior of the operator entanglement, the three models are fully distinguished.
Phenomenologically we observe that how late the system exits the plateau phase is correlated with how chaotic it is: the rational CFTs exit at t = B 2 ; the random unitary circuits exit at t = 2 − A 4 = B 2 + A 4 ; the holographic CFTs exit at t = 2 = B 2 + A 2 .By the underlying quenched dynamics, the reduced density matrix is fully thermalized at t = 2 , after which the operator state becomes short-range entangled and the operator entanglement has to drop to short-range value.Therefore the exit has to occur before t = 2 .This is saturated by a class of maximally chaotic systems: the holographic CFTs.It is therefore tempting to suggest that the persistence of the operator entanglement barrier, calibrated using the standard quenched state e −βH/4 |B , can serve as an alternative measure for the quantum chaoticity of the system.
To better check and understand the possible relation between the duration of the op-erator entanglement plateau and quantum chaoticity, it could be helpful in the future to study more examples of intermediate models, such as chaotic CFTs that are not holographic.Here based on the two chaotic models we study, it is still interesting to compare their mechanisms of exiting the plateau phase.In both models they are driven by geometrical transitions of the domain walls and the HRT surfaces respectively.The difference is as follows.In the random unitary circuits, the exit corresponds to a first order transition of the effective statistical mechanics, the operator entanglement is the free energy itself and thus is continuous across the exit.In particular, all the domain walls have the same tension in contributing to the free energy.In the holographic CFTs, the exit is driven by a phase transition of the bulk dual of the operator state, i.e. γ AB (t).The operator entanglement (or more specifically, reflective entropy) is given by the area of the HRT surface in the bulk, which changes by responding to the bulk dual transition, and thus is not continuous across the exit.In particular, contrary to the case in the random unitary circuits, γ AB (t) and γ A (t) play highly distinct roles in determining the operator entanglement.This distinction between the two models is not present when computing two-interval mutual information.This explains why the mutual information cannot distinguish between these two models as pointed out before.In the future, it would be interesting to connect this distinction with the different levels of quantum chaoticity between the two models.

Figure 1 :
Figure 1: Operator entanglement of the reduced density matrix.(a) Viewing the reduced density matrix as a matrix product operator.The operator entanglement is roughly logarithm of the bond dimension.The physics indices σ and σ are the indices in the operator Hilbert space.(b) The time dependence operator entanglement or reduced density matrix in 3 representative systems.Subsystem A is contiguous and has size A. Total system has length .The Entanglement quantities are displayed in units of the corresponding state equilibrium entropy density s

Figure 3 :
Figure 3: The Euclidean path integrals for evaluating the operator Rényi entropy,the two intervals represent A and B respectively.(a): The operator state |ρ α (t) (left) and the operator density matrix |ρ α AB ρ α AB | (right).(b): The non-trivial part of the operator Rényi entropy, where the way to connects layers is different in A and B. (c): The normalization of the operator density matrix, which is mapped to the state Rényi entropy.(d)/(e): The corresponding twist field insertion for (b)/(c).

Figure 4 :
Figure 4: Locations of the twist field in an infinite stripe.(a) The correlation function of the numerator in Eq. (2.6)(b) The correlation function of the denominator in Eq. (2.6).ψ(τ )| and |ψ(τ ) are defined by the upper and lower portion of the Euclidean path-integrals.

Figure 6 :
Figure 6: Results of rational conformal field theory.
t e x i t s h a 1 _ b a s e 6 4 = " A C e k 7 G x V 9 w S r B / M I Z X O i D Q z 8 X V c = " > A A A C V n i c d Z H N S s N A F I U n 0 d p a / 1 p d u g k W w V V J q q A 7 i 2 5 c t m B / o A 1 l M r l p h 0 4 y Y W Y i 1 N A n c K v v p i 8 j T t o s r I 0 X B g 7 f u Z e 5 c 8 a L G Z X K t r 8 M c 2 e 3 t F e u 7 F c P D o + O T 2 r 1 0 7 7 0 j j 6 M T + P b L J n l d a t p 5 D N n a K P M 2 g + K 1 b V 8 < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " m 4 B w 6 5 D c c l j 8 i D 0 8 s 5 P 7 f g L d F t 9 n w r x v N p 5 t 6 q / 2 2 C r 8 I z s E F u A I + u A U t 8 A j a o A M I m I F 3 8 A E + n W + 3 4 J b d y q r V d d Y L O w U b 5 Z 7 9 A G d z t 9 A = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " A C e k 7 G x V 9 w S r B / M I Z X O i D Q z 8 X V c = " > A A A C V n i c d Z H N S s N A F I U n 0 d p a / 1 p d u g k W w V V J q q A 7 i 2 5 c t m B / o A 1 l M r l p h 0 4 y Y W Y i 1 N A n c K v v p i 8 j T t o s r I 0 X B g 7 f u Z e 5 c 8 a L G Z X K t r 8 M c 2 e 3 t F e u 7 F c P D o + O T 2 r 1 0 7 7 k i S D Q I 5 x x M f S w B E Y j 6 C m q G A x j A T j 0 G A y 8 + W P m D 1 5 A S M q j 0 j j 6 M T + P b L J n l d a t p 5 D N n a K P M 2 g + K 1 b V 8 < / l a t e x i t > z < l a t e x i t s h a 1 _ b a s e 6 4 = " m 4 B w 6 5 D c c l j 8 i D 0 8d Q N m E b 7 h Q / s = " > A A A C X H i c d Z H N S g M x F I U z Y 9 X a W m 0 V 3 L g J F s F V m a m C 7 i y 4 c V n B / k A 7 l C S T t q H J Z E g y Q h 0 K v o J b f T I 3 P o u Z t g t r x w u B w 3 d u y M 2 5 O O Z M G 8 / 7 c t y d w u 7 e f v G g V D 6 s H B 1 X a y d d L R N F a I d I L l U f I 0 0 5 i 2 j H M M N p P 1 Y U C c x p D 8 8 e M r / 3 Q p V m M n o 2 8 5 g G A k 0 i N m Y E G Y t 6 Q 4 x U + r o Y V e t e w 1 s W 3 B b + W t T B u t q j m n M / D C V J B I 0 M 4 U j r g e / F J k i R M o x w u i g N E 0 1 j R G Z o Q g d W R k h Q H a T L e R f w 0 p I Q j q W y J z J w S X / f S D E W t k s g M 9 Ub H A m t 5 w L / 5 2 U w 1 w t 1 9 l C u p e 3 0 U x r m e h k x U n K d 5 w 4 S M 7 4 L U h b F i a E R W f 1 q n H B o J M y S h i F T l B g + t w I R x W w w k E y R Q s T Y f Z R s 5 P 7 f g L d F t 9 n w r x v N p 5 t 6 q / 2 2 C r 8 I z s E F u A I + u A U t 8 A j a o A M I m I F 3 8 A E + n W + 3 4 J b d y q r V d d Y L O w U b 5 Z 7 9 A G d z t 9 A = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " A C e k 7 G x V 9 w S r B / M I Z X O i D Q z 8 X V c = " > A A A C V n i c d Z H N S s N A F I U n 0 d p a / 1 p d u g k W w V V J q q A 7 i 2 5 c t m B / o A 1 l M r l p h 0 4 y Y W Y i 1 N A n c K v v p i 8 j T t o s r I 0 X B g 7 f u Z e 5 c 8 a L G Z X K t r 8 M c 2 e 3 t F e u 7 F c P D o + O T 2 r 1 0 7 7 k i S D Q I 5 x x M f S w B E Y j 6 C m q G A x j A T j 0 G A y 8 + W P m D 1 5 A S M q j

Figure 8 :
Figure 8: Quasi-particle interpretation of results in Eq. (2.27).Left: Quasi-particle in thermal field double states.Initially each pair of degrees of freedom in the left and right form an EPR pair.The entanglement is proportional to the overlap of past the light cones.Right: the overlap of the past light cones of A (red) and B (green) in different stages.Figure shows those turning points t = 0, A 2 , 1 2 ( − A),

Figure shows those turning points t = 0
Figure 8: Quasi-particle interpretation of results in Eq. (2.27).Left: Quasi-particle in thermal field double states.Initially each pair of degrees of freedom in the left and right form an EPR pair.The entanglement is proportional to the overlap of past the light cones.Right: the overlap of the past light cones of A (red) and B (green) in different stages.Figure shows those turning points t = 0, A 2 , 1 2 ( − A), 1 2 in Eq. (2.27)

Figure 9 :
Figure 9: Random unitary circuit: (a) its tensor network structure, where each blue block represents a q × q Haar random unitary matrix (b) the time dependent density matrix (c) computing tr (ρ n A ) , n = 2 by contracting different copies.

Figure 10 :
Figure10: Permutations in the circuit average.The blue and red blocks represents the independent Haar random matrices and their complex conjugate.(a) The circuit structure for n = 2. Red circle encloses 2n = 4 layers of unitaries, which can be regarded as an 8n-leg tensor that has 4n legs at the bottom as input and 4n legs on the top as output.(b) Random average over the 2n layers in each block generates a sum over permutation elements with weights.

Figure 11 :
Figure 11: Statistical model in terms of domain walls.Left: Domain wall configurations for the 2nd Rényi entropy.Right: Coarse grained domain wall configuration.Red and blue represent two domains of permutation extended from the boundary twist imposed by the definition of entanglement.

Figure 12 :
Figure 12: Domain walls in state and operator Rényi entanglement.The spin boundary conditions at the top boundary of the tenor network are shown for (a): operator R'enyi entanglement and (b): state entanglement.(c)/(d): The corresponding domain walls at the top boundary for operator/state entanglement.Notice the resemblance of domain walls to the twist fields in Eq. (2.2).

Figure 13 :
Figure 13: Domain wall rules: (a) domain wall must not go outside the light cone.(b) "bubble" is forbidden.(c) legitimate configuration: an elementary domain wall travel inside the light cone.(d) legitimate configuration: two domain walls meet and annihilate.(3) domain wall split according to multiplication rule.

Figure 14 :
Figure 14: Space time domain wall structures for the two terms in Eq. (3.1).According to our convention, the domain wall starts from the top boundary and go downward.(a) Domain wall configurations for the operator Rényi operator entanglement of unnormalized state (1st term in Eq. (3.1)).There are six incoming (elementary) domain walls on the top boundary (see Fig. 12(c)).Three configurations that possibly occur at different stages of evolution are shown on the right of the brace.(b) Domain wall configurations for the state entropy S2(AB) of region AB.There are four incoming (elementary) domain walls on the top boundary (see Fig.12(d)).Their free energy is twice of a single pair of domain walls, due to their independence in the large q limit.The two configurations for the two phases of state quench entanglement are shown on the right of the brace: In the linear growth phase, the two domain walls microscopically going down vertically; in the saturation phase, the two domain walls meet and annihilate.

Figure 15 :
Figure 15: S op 2 (A, ρAB) when A = B = 2 .(a).Operator entanglement as the difference of two terms in Eq. (3.1).As analyzed in the text, the solid curve has transitions between 3 configurations while the dashed curve representing the state entropy has only linear growth and saturation phases.(b).S op 2 (A, ρAB) as a function of time.
AB i < l a t e x i t s h a 1 _ b a s e 6 4 = " h C O X T g j S 0 p 8 8 a J T Q d G 5 r 6 8 b L H Z E = " > A A A C a H i c d Z H L T g I x F I Y 7 4 w 3 x B r o w x s 0 E Y u K K z K C J 7 k T c u M R E L g l D S K d z g I Z O O 2 k 7 G j L y K G 7 1 m X w F n 8 I O s B D B k z T 5 8 / 3 n p K d / g 5 h R p V 3 3 y 7 I 3 N r e 2 d 3 K 7 + b 3 9 g r w r B S X c + N d S / F U l P C Y J r 3 E w U x J m M 8 h K 6 R H E e g e u l s 9 6 l z Y U j o D I Q 0 h 2 t n R n 9 P p E E Q m a 4 I 6 5 F a 4 j h S a h I F / 3 k Z X O u F K r t o r a X M 9 i M I 1 3 o Z 0 U I w t c 7 t J n p w 2 0 s p j x M N n M x f N U i Y o 4 W T p e 6 E V A L R b G I E J p K a Y B w y w h I T b f 4 m b y L 3 / g a 8 K l r V i n d V q T 5 d l 2 v 1 R f g 5 d I 5 K 6 B J 5 6 A b V 0 C N q o C Y i 6 B W 9 o w / 0 a X 3 b B f v U P p u 3 2 t Z i 5 g Q t l V 3 6 A b / 8 u 5 4 = < / l a t e x i t > B1 < l a t e x i t s h a 1 _ b a s e 6 4 = " M w o R b N 6 r F t 4 d S m R O i P V N / b y S e F 5 i p L K l e 0 m m g 9 T T w / v M y c a 3 n 6 + y i t Z Y 2 6 c f g r / U y B a U U e p 3 b j X F 4 3 0 9 4 G M U I I Z u / a h g L G 6 W d N W / 7 X A F D M T W E M s V N M T Y b U 0 U Z m v 8 p m M r d v w W v k l a t 6 t 5 U a y + 3 l X p j U X 6 e l M k F u S Y u u S N 1 8 k S e S Z M w k p B 3 8 k E + c 9 / W q V W 2 z u e j V m 6 x c 0 K W Y F 3 9 A J d j v W I = < / l a t e x i t > B2 < l a t e x i t s h a 1 _ b a s e 6 4 = " P x Q 9 O e T + a U z t Z n 9 p F D A G 3 6 8 M 6 Y Y 7 U v j S K b w + l M i d E e 6 b + 3 k g 8 L z B T W V K 9 p N N A 6 2 n g / e d l 4 l r P 1 9 l F a y 1 t 0 o / B X + t l C k o p 9 D q 3 G + P w v p / w M I o R Q j Z / 1 T A W N k o 7 a 9 7 2 u Q K G Y m o I Z Y q b Y m w 2 p o o y N P 9 T M J W 7 f w t e J a 1 a 1 b 2 p 1 l 5 u K / X G o v w 8 K Z M L c k 1 c c k f q 5 I k 8 k y Z h J C H v 5 I N 8 5 r 6 t U 6 t s n c 9 H r d x i 5 4 Q s w b r 6 A Z l V v W M = < / l a t e x i t > ⌃(h) < l a t e x i t s h a 1 _ b a s e 6 4 = " L f y l H d u N O P W U g P t 7 m 6 M G N r q N O 7 Y = " > A A A C X n i c d Z F P S 8 M w G M b T + m 9 O 5 z a 9 C F 6 K Q 5 i X 0 U 5 B b w 6 9 e J z o / s h W R p q m a 1 j S l C Q V R t m n 8 K o f z J s f x X T r w b n 6 Q u D h 9 7 w v e f P E i y m R y r a / D H N r e 2 d 3 r 7 R f P j i s H F V r 9 e O + 5 I l A u I c 4 5 W L o Q Y k p i X B P E U X x M B Y Y M o / i g T d 7 y P z B G x a S 8 O h F z W P s M j i N S E A Q V B q 9 j p / J l M F m e D m p N e y W v S x r U z i 5 a I C 8 u p O 6 c T f 2O U o Y j h S i U M q R Y 8 f K T a F Q B F G 8 K I 8 T i W O I Z n C K R 1 p G k G H p p s u N F 9 a F J r 4 V c K F P p K w l / T 2 R e h 7 T X Q y q U K 5 x y K S c M + 8 / L 4 O F n i + z i w o t q b c P s V / o Z U R x T m W R O 0 p U c O u m J I o T h S O 0 e l W Q U E t x K 8 v a 8 o n A S N G 5 F h A J o o O x U A g F R E r / S F l H 7 vw N e F P 0 2 y 3 n q t V + u m 5 0 7 v P w S + A M n I M m c M A N 6 I B H 0 A U 9 g A A D 7 + A D f B r f 5 q 5 Z M a u r V t P I Z 0 7 A W p m n P 5 K i t 1 k = < / l a t e x i t > AB 2 < l a t e x i t s h a 1 _ b a s e 6 4 = " N n r 9 i m 8 i 9 v w G v i 1 6 z 4 V 0 2 m o 9 X 9 V Y 7 D 7 8 M T s A Z u A A e u A Y t 8 A A 6 o A s Q e A X v 4 A N 8 l r 4 t 2 6 p a t W W r V c p n a m C l r N M f p s S 5 R A = = < / l a t e x i t > B1 < l a t e x i t s h a 1 _ b a s e 6 4 = " M w o R b N 6 r F t 4 d S m R O i P V N / b y S e F 5 i p L K l e 0 m m g 9 T T w / v M y c a 3 n 6 + y i t Z Y 2 6 c f g r / U y B a U U e p 3 b j X F 4 3 0 9 4 G M U I I Z u / a h g L G 6 W d N W / 7 X A F D M T W E M s V N M T Y b U 0 U Z m v 8 p m M r d v w W v k l a t 6 t 5 U a y + 3 l X p j U X 6 e l M k F u S Y u u S N 1 8 k S e S Z M w k p B 3 8 k E + c 9 / W q V W 2 z u e j V m 6 x c 0 K W Y F 3 9 A J d j v W I = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " A C e k 7 G x V 9 w S r B / M I Z X O i D Q z 8 X V c = " > A A A C V n i c d Z H N S s N A F I U n 0 d p a / 1 p d u g k W w V V J q q A 7 i 2 5 c t m B / o A 1 l M r l p h 0 4 y Y W Y i 1 N A n c K v v p i 8 j T t o s r I 0 X B g 7 f u Z e 5 c 8 a L G Z X K t r 8 M c 2 e 3 t F e u 7 F c P D o + O T 2 r 1 0 7 7 k i S D Q I 5 x x M f S w B E Y j 6 C m q G A x j A T j 0 G A y 8 + W P m D 1 5 A S M q j

Figure 17 :
Figure 17: Computation of the operator state norm N = ρAB|ρAB .Left: Euclidean path-integral on the doubly-branched manifold in the boundary CFT; Right: bulk saddle point contribution B = B1 ∪ B2 to N .The two wedges intersect at the fixed Cauchy surface Σ(h) = B1 ∩ B2 under reflection, which contains the cosmic-brane γ AB 2 .Due to the brane backreaction, the bulk portions matching to AB and A B towards the boundary only extend an angular range of π around γ AB 2 in the bulk.

2 <Figure 18 :
Figure 18: Prescription for obtaining the bulk geometry dual of operator state |ρAB from that of |ψ .Left: start from the bulk dual of |ψ ; Center: insert a cosmic brane γ AB 2 , let it back-react on the dual geometry and settle down; Right: spatial geometry dual to |ρAB is obtained by identifying two copies of the relevant bulk regions across γ AB 2 .
t e x i t s h a 1 _ b a s e 6 4 = " L J h S d J f 7 i Z r 8 O G 3 x 5 p / o D e m p a h o = " > A A A C Y n i c d V H N S 8 M w H E 3 r 1 6 x f m z v q I T g E Q R j t F P T m 0 I s H D x P c B 2 x 1 p G n m w p K m J q k w y v 4 O r / p n e f c P M d 1 6 c K 4 + C D z e + / 3 I y 0 s Q M 6 q 0 6 3 5 Z 9 t r 6 x u Z W a d v Z 2 d 3b P y h X D j t K J B K T N h Z M y F 6 A F G E 0 I m 1 N N S O 9 W B L E A 0 a 6 w e Q u 8 7 t v R C o q o i c 9 j Y n P 0 U t E R x Q j b S R / w J E e Y 8 T S h 9 n z + b B c c + v u H H C V e D m p g R y t Y c W 6 G Y Q C J 5 x E G j O k V N 9 z Y + 2 n S G q K G Z k 5 g 0 S R G O E J e i F 9 Q y P E i f L T e e o Z P D V K C E d C m h N p O F d / b 6 R B w M 1 U l l A t 6 Y g r N e X B f 1 4 m F n q h y i 4 q t J R J P y Z h o Z c p W g i m i t x + o k f X f k q j O N E k w o t X j R I G t Y B Z 3 z C k k m D N p o Y g L K k p B u I x k g hr 8 y u O q d z 7 W / A q 6 T T q 3 k W 9 8 X h Z a 9 7 m 5 Z f A E T g B Z 8 A D V 6 A J 7 k E L t A E G r + A d f I B P 6 9 t 2 7 I p d X Y z a V r 5 T B U u w j 3 8 A G N q 5 f g = = < / l a t e x i t > L < l a t e x i t s h a 1 _ b a s e 6 4 = " C A D M M y C l B e Z f 8 e c D L d 0 s F 7 Z 2

1 <
4 y d L p 7 v v k 8 z m I G V X a d b 8 s e 2 1 9 Y 3 O r t O 3 s 7 O 7 t H 5 Q r h x 0 l E o l J G w s m Z C 9 A i j A a k b a m m p F e L A n i A S P d Y H K X + d 0 3 I h U V 0 Z O e x s T n 6 C W i I 4 q R N p I / 4 E i P M W L p w + z 5 fF i u u X V 3 D r h K v J z U Q I 7 W s G L d D E K B E 0 4 i j R l S q u + 5 s f Z T J D X F j M y c Q a J I j P A E v Z C + o R H i R P n p P P U M n h o l h C M h z Y k 0 n K u / N 9 I g 4 G Y q S 6 i W d M S V m v L g P y 8 T C 7 1 Q Z R c V W s q k H 5 O w 0 M s U L Q R T R W 4 / 0 a N r P 6 V R n G g S 4 c W r R g m D W s C s b x h S S b B m U 0 M Q l t Q U A / E Y S Y S 1 + R X H V O 7 9 L X i V d Bp 1 7 6 L e e L y s N W / z 8 k v g C J y A M + C B K 9 A E 9 6 A F 2 g C D V / A O P s C n 9 W 0 7 d s W u L k Z t K 9 + p g i X Y x z 8 c v r m A < / l a t e x i t > t = t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " 9 8 1 h y c D t a 3 D X k M m m 5 C L Q T h C 4 w V U = " > A A A C W n i c d Z H L S s N A F I Y n 8 d a L l 1 b d u Q k W w V V J V N C N W H T j s o J p C 2 0 o k 8 m k H T r J h J k T o Y Q + g 1 t 9 N M G H c d J m Y W 0 8 M P D z / e c w Z / 7 x E 8 4 U 2 P a X Y W 5 t 7 + z u V a q 1 + v 7 B 4 V G j e d x T I p W E u k R w I Q c + V p S z m L r A g N N B I i m O f E 7 7 / u w p 9 / t v V C o m 4 l e Y J 9 S L 8 C R m I S M Y N H L h H s b 2 u N G y 2 / a y r E 3 h F K K F i u q O m 8 b D K B A k j W g M h G O l h o 6 d g J d h C Y x w u q i N U k U T T G Z 4 Q o d a x j i i y s u W 2 y 6 s C 0 0 C K x R S n x i s J f 0 9 k f l + p L s i D F O 1 x n G k 1 D z y / / N y W O o F K r + o 1 F J 6 + y k N S r 2 c g B B c l b n D F M I 7 L 2 N x k g K N y e p V Y c o t E F a e s x U w S Q n w u R a Y S K a D s c g U S 0x A / 0 Z N R + 7 8 D X h T 9 K 7 a z n X 7 6 u W m 1 X k s w q + g M 3 S O L p G D b l E H P a M u c h F B D L 2 j D / R p f J u m W T X r q 1 b T K G Z O 0 F q Z p z 9 i H r X f < / l a t e x i t > l a t e x i t s h a 1 _ b a s e 6 4 = " + d 1 K i B w q y w Z q B k Q 3 8 g 5

2 <
s c c a V m P P j P y 2 G p F 6 r 8 o l J L m e 1 j E p Z 6 O d F C M F X m j j I d 3 f p z m q S Z J g l e v S r K G N Q C 5 k H D k E q C N Z s Z g b C k J h i I Y y Q R 1 u Y 7 H B O 5 9 z f g T d F v t 7 y r V v v p u t m 5 L 8 K v g j N w D i 6 B B 2 5 A B z y C L u g B D F 7 A O / g A n 9 a 3 X b E d u 7 Z q t a 1 i 5 g S s l X 3 6 A 3 c C t m M = < / l a t e x i t > t = − ∞ A'B' L + < l a t e x i t s h a 1 _ b a s e 6 4 = " t O F u Z h 2 t H S k d x z p M w z k A N U R I C 0 E = " > A A A C Z n i c d V H d S s M w G E 3 r 3 5 y 6 H 0 W 8 8 C Y 4 B E E Y 7 R T 0 z o F e e O H F B P c D W y 1 p m m 5 h a V O S V B h l T + K t P p R v 4 G O Y b r 1 w r h 4 I H M 7 5 P n J y 4 s W M S m V Z X 4 a 5 s b m 1 v V P a L e / t H 1 S q tf p h T / J E Y N L F n H E x 8 J A k j E a k q 6 h i Z B A L g k K P k b 4 3 v c / 8 / h s R k v L o R c 1 i 4 o R o H N G A Y q S 0 5 N a q o x C p C U Y s f Z q / X r o t t 9 a w m t Y C c J 3 Y O W m A H B 2 3 b t y N f I 6 T k E Q K M y T l 0 L Z i 5 a R I K I o Z m Z d H i S Q x w l M 0 J k N N I x Q S 6 a S L 5 H N 4 r h U f B l z o E y m 4 U H 9 v p J 4 X 6 q k s o 1 z R U S j l L P T + 8 z K x 0 P N l d l G h J X X 6 C f E L v U x R n D N Z 5 A 4 T F d w 6 K Y 3 i R J E I L 1 8 V J A w q D r P O o U 8 F w Y r N N E F Y U F 0 M x B M k E F b 6 Z 8 q 6 c v t v w e u k 1 2 r a V 8 3 W 8 3 W j / Z C X X w K n 4 A x c A B v c g D Z 4 B B 3 Q B R g k 4 B 1 8 g E / j 2 6 y Y x + b J c t Q 0 8 p 0 j s A I T / g A N 5 r p W < / l a t e x i t > L l a t e x i t s h a 1 _ b a s e 6 4 = " 1 j 7 l W v O 4 X k v Z P G h q k P H N q Q X p U Z Y = " > A A A C Z n i c d V H d S s M w G E 3r 3 5 y 6 H 0 W 8 8 C Y 4 B G 8 c 7 R T 0 z o F e e O H F B P c D W y 1 p m m 5 h a V O S V B h l T + K t P p R v 4 G O Y b r 1 w r h 4 I H M 7 5 P n J y 4 s W M S m V Z X 4 a 5 s b m 1 v V P a L e / t H 1 S q tf p h T / J E Y N L F n H E x 8 J A k j E a k q 6 h i Z B A L g k K P k b 4 3 v c / 8 / h s R k v L o R c 1 i 4 o R o H N G A Y q S 0 5 N a q o x C p C U Y s f Z q / X r o t t 9 a w m t Y C c J 3 Y O W m A H B 2 3 b t y N f I 6 T k E Q K M y T l 0 L Z i 5 a R I K I o Z m Z d H i S Q x w l M 0 J k N N I x Q S 6 a S L 5 H N 4 r h U f B l z o E y m 4 U H 9 v p J 4 X 6 q k s o 1 z R U S j l L P T + 8 z K x 0 P N l d l G h J X X 6 C f E L v U x R n D N Z 5 A 4 T F d w 6 K Y 3 i R J E I L 1 8 V J A w q D r P O o U 8 F w Y r N N E F Y U F 0 M x B M k E F b 6 Z 8 q 6 c v t v w e u k 1 2 r a V 8 3 W 8 3 W j / Z C X X w K n 4 A x c A B v c g D Z 4 B B 3 Q B Rg k 4 B 1 8 g E / j 2 6 y Y x + b J c t Q 0 8 p 0 j s A I T / g A R z r p Y < / l a t e x i t > t = t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " 9 8 1 h y c D t a 3 D X k M m m 5 C LQ T h C 4 w V U = " > A A A C W n i c d Z H L S s N A F I Y n 8 d a L l 1 b d u Q k W w V V J V N C N W H T j s o J p C 2 0 o k 8 m k H T r J h J k T o Y Q + g 1 t 9 N M G H c d J m Y W0 8 M P D z / e c w Z / 7 x E 8 4 U 2 P a X Y W 5 t 7 + z u V a q 1 + v 7 B 4 V G j e d x T I p W E u k R w I Q c + V p S z m L r A g N N B I i m O f E 7 7 / u w p 9 / t v V C o m 4 l e Y J 9 S L 8 C R m I S M Y N H L h H s b 2 u N G y 2 / a y r E 3 h F K K F i u q O m 8 b D K B A k j W g M h G O l h o 6 d g J d h C Y x w u q i N U k U T T G Z 4 Q o d a x j i i y s u W 2 y 6 s C 0 0 C K x R S n x i s J f 0 9 k f l + p L s i D F O 1 x n G k 1 D z y / / N y W O o F K r + o 1 F J 6 + y k N S r 2 c g B B c l b n D F M I 7 L 2 N x k g K N y e p V Y c o t E F a e s x U w S Q n w u R a Y S K a D s c g U S 0 x A / 0 Z N R + 7 8 D X h T 9 K 7 a z n X 7 6 u W m 1 X k s w q + g M 3 S O L p G D b l E H P a M u c h F B D L 2 j D / R p f J u m W T X r q 1 b T K G Z O 0 F q Z p z 9 i H r X f < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " p U F k 5 r 9 x A 5 / V O L t 8 P 6 D r U G p S p 0 A = " > A A A C Z n i c d V H d S s M w G E 3 r 3 5 y 6 H 0 W 8 8 C Y 4 B E E Y 7 R T 0 z o E I X n g x w f 3 A V k u a p l t Y 2 p Q k F U b Z k 3 i r D + U b + B i m W y + c q w c C h 3 O + j 5 y c e D G j U l n W l 2 F u b G 5 t 7 5 R 2 y 3

1 <
B e I I E w k r / T F l X b v 8 t e J 3 0 W k 3 7 q t l 6 v m 6 0 H / L y S + A U n I E L Y I M b 0 A a P o A O 6 A I M E v I M P 8 G l 8 m x X z 2 D x Z j p p G v n M E V m D C H w x B u l Y = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " 9 8 1 h y c D t a 3 D X k M m m 5 C LQ T h C 4 w V U = " > A A A C W n i c d Z H L S s N A F I Y n 8 d a L l 1 b d u Q k W w V V J V N C N W H T j s o J p C 2 0 o k 8 m k H T r J h J k T o Y Q + g 1 t 9 N M G H c d J m Y W 0 8 M P D z / e c w Z / 7 x E 8 4 U 2 P a X Y W 5 t 7 + z u V a q 1 + v 7 B 4 V G j e d x T I p W E u k R w I Q c + V p S z m L r A g N N B I i m O f E 7 7 / u w p 9 / t v V C o m 4 l e Y J 9 S L 8 C R m I S M Y N H L h H s b 2 u N G y 2 / a y r E 3 h F K K F i u q O m 8 b D K B A k j W g M h G O l h o 6 d g J d h C Y x w u q i N U k U T T G Z 4 Q o da x j i i y s u W 2 y 6 s C 0 0 C K x R S n x i s J f 0 9 k f l + p L s i D F O 1 x n G k 1 D z y / / N y W O o F K r + o 1 F J 6 + y k N S r 2 c g B B c l b n D F M I 7 L 2 N x k g K N y e p V Y c o t E F a e s x U w S Q n w u R a Y S K a D s c g U S 0 x A / 0 Z N R + 7 8 D X h T 9 K 7 a z n X 7 6 u W m 1 X k s w q + g M 3 S O L p G D b l E H P a M u c h F B D L 2 j D / R p f J u m W T X r q 1 b T K G Z O 0 F q Z p z 9 i H r X f < / l a t e x i t > L l a t e x i t s h a 1 _ b a s e 6 4 = " a B n 4 C E B N 1 x t 9 c p 6 C k z N W z o E c G n 4 = " > A A A C Z n i c d V H d S s M w G E 3 r 3 5 y 6 H 0 W 8 8 C Y 4 B G 8 c 7 R T 0 z o F e e O H F B P c D W y 1 p m m 5 h a V O S V B h l T + K t P p R v 4 G O Y b r 1 w r h 4 I H M 7 5 P n J y 4 s W M S m V Z X 4 a 5 s b m 1 v

4 t
7 a c 4 L P Q y o o W g q s g d J j q 6 8 V P C 4 0 R j j p a v i h L q a O F k e T s h k R h p O j M C I k l M M A 6 a Q g m R N r 9 i m 8 i 9 v w G v i 1 6 z 4 V 0 2 m o 9 X 9 V Y 7 D 7 8 M T s A Z u A A e u A Y t 8 A A 6 o A s Q e A X v 4 A N 8 l r 4 t 2 6 p a t W W r V c p n a m C l r N M f p s S 5 R A = = < / l a t e x i t > AB 0 = t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " 9 81 h y c D t a 3 D X k M m m 5 C L Q T h C 4 w V U = " > A A A C W n i c d Z H L S s N A F I Y n 8 d a L l 1 b d u Q k W w V V J V N C N W H T j s o J p C 2 0 o k 8 m k H T r J h J k T o Y Q + g 1 t 9 N M GH c d J m Y W 0 8 M P D z / e c w Z / 7 x E 8 4 U 2 P a X Y W 5 t 7 + z u V a q 1 + v 7 B 4 V G j e d x T I p W E u k R w I Q c + V p S z m L r A g N N B I i m O f E 7 7 / u w p 9 / t v V C o m 4 l e Y J 9 S L 8 C R m I S M Y N H L h H s b 2 u N G y 2 / a y r E 3 h F K K F i u q O m 8 b D K B A k j W g M h G O l h o 6 d g J d h C Y x w u q i N U k U T T G Z 4 Q o d a x j i i y s u W 2 y 6 s C 0 0 C K x R S n x i s J f 0 9 k f l + p L s i D F O 1 x n G k 1 D z y / / N y W O o F K r + o 1 F J 6 + y k N S r 2 c g B B c l b n D F M I 7 L 2 N x k g K N y e p V Y c o t E F a e s x U w S Q n w u R a Y S K a D s c g U S 0 x A / 0 Z N R + 7 8 D X h T 9 K 7 a z n X 7 6 u W m 1 X k s w q + g M 3 S O L p G D b l E H < l a t e x i t s h a 1 _ b a s e 6 4 = " O 8I H u O L x 1 N Q A V B p O j r H S s Y I 4 z J M = " > A A A C Y H i c d Z H L T g I x F I Y 7 4 w 3 x A u h O N 4 3 E B D d k B k 1 0 J 9 G N S 4 x y S Y C Q T q d A Q z u d t G d M C O E x 3 O p z u f V J 7 A A L k f E k T f5 8 / z n p 6 d 8 g F t y A 5 3 0 5 7 t b 2 z u 5 e b j 9 / c H h 0 X C i W T l p G J Z q y J l V C 6 U 5 A D B M 8 Y k 3 g I F g n 1 o z I Q L B 2 M H l M / f Y b 0 4 a r 6 B W m M e t L M o r 4 k F r l m F M T U C k I 1 t 8 F g O i a a U L B / k r e R + 3 8 D 3 h S t W t W / r t a e b 8 r 1 h 1 X 4 O X S O L l A F + e g W 1 d E T a q A m o k i h d / S B P p 1 v N + c W 3 N K y 1 X V W M 6 d o r d y z H x H Y u A g = < / l a t e x i t >

2 6 B+ 1 (
t 0 ) < l a t e x i t s h a 1 _ b a s e 6 4 = "d g p m C v a F R t j / M N Q 9 K d + A y Q 2 e H b 8 = " > A A A C c 3 i c d V F N T 8 J A E F 3 q F + I X 6 B E P j W i C M Z I W T f Q m w Y t H T U R I A J v t d o A N 2 2 6 z O z U h D R d / j V f 9 N / 4 Q 7 2 6 Rg x / 4 k k 1 e 3 p v J z r z x Y 8 E 1 O s 5 7 z l p a X l l d y 6 8 X N j a 3 t n e K p d 0 H L R P F o M W k k K r j U w 2 C R 9 B C j g I 6 s Q I a + g L a / v g 6 8 9 t P o D S X 0 T 1 O Y u i H d B j x A W c U j e Q V 9 3 v I R Q B p L 6 Q 4 Y l S k z e n 0 8 c R z q + g 5 x 1 6 x 4 t S c G e y / x J 2 T C p n j 1 i v l r n q B Z E k I E T J B t e 6 6 T o z 9 l C r k T 4 A N 8 l r 4 t 2 6 p a t W W r V c p n a m C l r N M f p s S 5 R A = = < / l a t e x i t > t = t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " 9 8 1 h y c D t a 3 D X k M m m 5 C L 5 8 / z n p 6 d 8 g F t y A 5 3 0 5 7 t b 2 z u 5 e b j 9 / c H h 0 X C i W T l p G J Z q y J l V C 6 U 5 A D B M 8 Y k 3 g I F g n 1 o z I Q L B 2 M H l M / f Y b 0 4 a r 6 B W m M e t L M o r 4 k F

2 4 6 B 2
a a U L B / k r e R + 3 8 D 3 h S t W t W / r t a e b 8 r 1 h 1 X 4 O X S O L l A F + e g W 1 d E T a q A m o k i h d / S B P p 1 v N + c W 3 N K y 1 X V W M 6 d o r d y z H x H Y u A g = < / l a t e x i t > 0 (t 0 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " T h s g k M b s H 9 z o A k I n n + j l C R E v F c o = " > A A A C c 3 i c d V H L S s N A F J 3 G d 3 1 V X d Z F s A p 1 0 Z J U Q X c W 3 b h U s A 9 o a 5 h M b t u h k 0 y Y u R F K y M a v c a t / 4 4 e 4 d 1 K 7 8 F E P D B z O u Z e 5 9 1 w / F l y j 4 7 w X r K X l l d W 1 9 Y 3 i 5 t b 2 z m 5 p b 7 + t Z a I Y t J g U U n V 9 q k H w C F r I U U A 3 V k B D X 0 D H n 9 z k f u c J l O Y y e s B p D I O Q j i I + 5 I y i k b z S Y R + 5 C C D t h x T H j I r 0 O s s e a 1 6 j i p 5 z 6 p U q T t 2 Z w f 5 L 3 D m p k D n u 6 B g x / 4 k k 1 e 3 p v J z r z x Y 8 E 1 O s 5 H z l p a X l l d y 6 8 X N j a 3 t n e K p d 0 H L R P F o M W k k K r j U w 2 C R 9 B C j g I 6 s Q I a + g L a / u g 2 8 9 v P o D S X 0 T 2 O Y + i H 9 C n i A 8 4 o G s k r 7 v e Q i w D S X k h x y K h I b y a T x 1 O v d o y e c + I V K 0 7 V m c L + S 9 w 5 q Z A 5 m l 4 p d 9 0 L

Figure 20 : 2
Figure 20: Bulk Lorentzian saddle in the form of two Schwinger-Keldysh two-folds: B(t0) = ∪ a=± i=1,2 Ba i (t0), glued at Σ(t0).The structure of Rindler wedges near the cosmic brane γ AB 2 are shown, whose indices m indicate the corresponding boundary Rindler wedges they asymptote to in Fig. 19.However, the bulk is smooth in the interior, the imaginary part of the local Rindler time near γ AB 2 in the wedges are Im τ = mπ 4 , corresponding to a standard Rindler temperature T Rindler = 1/2π.As a result, a single two-folds (e.g.B1) has a Rindler temperature of T Rindler = 1/π near γ AB 2 after identifying m = 0 and m = 4.

4. 4
Operator Entanglement of e −β/2H |B Let us now consider the holographic picture for the operator entanglement of the generalized operator state |ρ 1 2 t e x i t s h a 1 _ b a s e 6 4 = " 2 0 Q B 9 9 K m Q ym + i M 1 F K q V w r X C S t 3 U = " > A A A C Z 3 i c d Z H L T g I x F I b L e E O 8 M G h i T N y M E h P c k B k 0 0 Z 2 I G 5 e a y C U B J J 1 O g Y Z 2 O m n P k J A J b + J W 38 l H 8 C 3 s A A s R P E m T P 9 9 / T n r 6 1 4 8 4 0 + C 6 X x l r Y 3 N r e y e 7 m 9 v b P z j M 2 4 W j h p a x I r R O J J e q 5 W N N O Q t p H R h w 2 o o U x c L n t O m P H l O / O a Z K M x m + w i S i X Y E H I e s z g s G g n m 1 3 B l g I / J Y 8 1 K Y l 6 L l X P b v o l t 1 Z O a v C W 4 g i W t R z r 5 C 5 7 w S S x I K G Q D j W u u 2 5 E X Q T r I A R T q e 5 T q x p h M k I D 2 j b y B A L q r v J b P W p c 2 l I 4 P S l M i c E Z 0 Z / T y S + L 0 y X w D D U S x w L r S f C / 8 9 L 4 b b y l s n 1 u m 8 1 c o s Z o 7 R U l n n P y b k u l Y = < / l a t e x i t > t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " X / h 9 F 9 Q c 7 I I 4 E b x 3 6 6 M Y e n b H u Y 4 = " > A A A C W H i c d Z H L S s N A F I Z P 4 6 U X b 2 1 d u g a a x Y T O y I Q N j I x I y P R o s d p 1 a V 8 Y 4 t u B V O Z E a K / o 7 4 m F 5 4 W m K y Q 4 1 R u c h F r P Q + 8 / L 4 W 5 n q / T i 3 I t b b a f M j / X S w l K K X S e O 0 g w u B s t e B Q n y C K 6 f l W Q C B u l n a Z s + 1 wx i m J u B K G K m 2 B s O i W K U D R / U T G R u 3 8 D 3 h b d q 6 Z 7 3 b x 6 v m m 0 H r L w S 3 A G 5 3 A J L t x C C 5 6 g D R 2 g M I F 3 + I D P w r c F V t E q r 1 u t Q j Z z C h t l 1 X 8 A v H m 1 G g = = < / l a t e x i t > A (t 0 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " w O i E q k B c C r L E 3 q + I O 4 4 u t U 1 w E D c = " > A A A C Z n i c d Z H L T g I x F I b L e E O 8 A B r j w k 0 j M c E N m U E T 3 Y m 6 c Y m J X B J A 0 u k U a G i n k / a M C Z n w J G 7 1 o X w D H 8 M O s B D B k z T 5 8 / 3 n p K d / / U h w A 6 7 7 l X E 2 N r e 2 d 7 K 7 u b 3 9 g 8 N 8 o X j U N C r W l D W o E k q 3 f W K Y 4 C F r A A f B 2 p F m R P q C t f z x Y + q 3 3 p g 2 X I U v M I l Y T 5 J h y A e c E r C o X 8 h 3 h 0 R K 8 p r c T 8 v Q d y / 7 h Z J b c W e F V 4 W 3 E C W 0 q H q / m L n r B o r G k o V A B T G m 4 7 k R 9 B K i g V P B p r l u b F h E 6 J g M W c f K k E h m e s l s 8 y m + s C T A A 6 X t C Q H P 6 O + J x P e l 7 Z I E R m a J E 2 n M R P r / e S l c 6 w U m v W i t Z e z 2 I x a s 9 V I C S g m z z u 3 E M L j t J T y M Y m A h n b 9 q E A s M C q e Z 4 4 B r R k F M r C B U c x s M p i O i C Q X 7 M z k b u f c 3 4 F X R r F a 8 q 0 r 1 + b p U e 1 i E n 0 V n 6 B y V k Y d u U A0 9 o T p q I I p i 9 I 4 + 0 G f m 2 z l 0 T p z T e a u T W c w c o 6 V y 8 A 9 7 3 7 o K < / l a t e x i t > EoW brane AB (t0) < l a t e x i t s h a 1 _ b a s e 6 4 = " 2 0 Q B 9 9 K m Q y m + i M 1 F K q V w r X C S t 3 U = " > A A A C Z 3 i c d Z H L T g I x F I b L e E O 8 M G h i T N y M E h P c k B k 0 0 Z 2 I G 5 e a y C U B J J 1 O g Y Z 2 O m n P k J A J b + J W 3 8 l H 8 C 3 s A A s R P E m T P 9 9 / T n r 6 1 4 8 4 0 + C 6 X x l r Y 3 N r e y e 7 m 9 v b P z j M 2 4 W j h p a x I r R O J J e q 5 W N N O Q t p H R h w 2 o o U x c L n t O m P H l O / O a Z K M x m + w i S i X Y E H I e s z g s G g n m 1 3 B l g I / J Y 8 1 K Y l 6 L l X P b v o l t 1 Z O a v C W 4 g i W t R z r 5 C 5 7 w S S x I K G Q D j W u u 2 5 E X Q T r I A R T q e 5 T q x p h M k I D 2 j b y B A L q r v J b P W p c 2 l I 4 P S l M i c E Z 0 Z / T y S + L 0 y X w D D U S x w L r S f C / 8 9 L 4V o v 0 O l F a y 1 t t h / S Y K 2 X E p C S 6 3 V u O 4 b + X T d h Y R Q D D c n 8 V f 2 Y O y C d N H Q n Y I o S 4 B M j M F H M B O O Q I V a Y g P m a n I n c + x v w q m h U y t 5 1 u f J y U 6 z W F u F n 0 R m 6 Q C X k o V t U R U / o G d U R Q W P 0 j j 7 Q Z +b b y l s n 1 u m 8 1 c o s Z o 7 R U l n n P y b k u l Y = < / l a t e x i t > B A t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " X / h 9 F 9 Q c 7 I I 4 E b x 3 6 6 M Y e n b H u Y 4 = " > A A A C W H i c d Z H L S s N A F I Z P 4 6 U X b 2 1 d u g a a x Y T O y I Q N j I x I y P R o s d p 1 a V 8 Y 4 t u B V O Z E a K / o 7 4 m F 5 4 W mK y Q 4 1 R u c h F r P Q + 8 / L 4 W 5 n q / T i 3 I t b b a f M j / X S w l K K X S e O 0 g w u B s t e B Q n y C K 6 f l W Q C B u l n a Z s + 1 w x i m J u B K G K m 2 B s O i W K U D R / U T G R u 3 8 D 3 h b d q 6 Z 7 3 b x 6 v m m 0 H r L w S 3 A G 5 3 A J L t x C C 5 6 g D R 2 g M I F 3 + I D P w r c F V t E q r 1 u t Q j Z z C h t l 1 X 8 A v H m 1 G g = = < / l a t e x i t > EoW brane AB (t0) < l a t e x i t s h a 1 _ b a s e 6 4 = " 2 0 Q B 9 9 K m Q y m + i M 1 F K q V w r X C S t 3 U = " > A A A C Z 3 i c d Z H L T g I x F I b L e E O 8 M G h i T N y M E h P c k B k 0 0 Z 2 I G 5 e a y C U B J J 1 O g Y Z 2 O m n P k J A J b + J W 38 l H 8 C 3 s A A s R P E m T P 9 9 / T n r 6 1 4 8 4 0 + C 6 X x l r Y 3 N r e y e 7 m 9 v b P z j M 2 4 W j h p a x I r R O J J e q 5 W N N O Q t p H R h w 2 o o U x c L n t O m P H l O / O a Z K M x m + w i S i X Y E H I e s z g s G g n m 1 3 B l g I / J Y 8 1 K Y l 6 L l X P b v o l t 1 Z O a v C W 4 g i W t R z r 5 C 5 7 w S S x I K G Q D j W u u 2 5 E X Q T r I A R T q e 5 T q x p h M k I D 2 j b y B A L q r v J b P W p c 2 l I 4 P S l M i c E Z 0 Z / T y S + L 0 y X w D D U S x w L r S f C / 8 9 L 4 b b y l s n 1 u m 8 1 c o s Z o 7 R U l n n P y b k u l Y = < / l a t e x i t > AB (t0) < l a t e x i t s h a 1 _ b a s e 6 4 = " 2 0 Q B 9 9 K m Q y m + i M 1 F K q V w r X C S t 3 U = " > A A A C Z 3 i c d Z H L T g I x F I b L e E O 8 M G h i T N y M E h P c k B k 0 0 Z 2 I G 5 e a y C U B J J 1 O g Y Z 2 O m n P k J A J b + J W 3 8 l H 8 C 3 s A A s R P E m T P 9 9 / T n r 6 1 4 8 4 0 + C 6 X x l r Y 3 N r e y e 7 m 9 v b P z j M 2 4 W j h p a x I r R O J J e q 5 W NN O Q t p H R h w 2 o o U x c L n t O m P H l O / O a Z K M x m + w i S i X Y E H I e s z g s G g n m 1 3 B l g I / J Y 8 1 K Y l 6 L l X P b v o l t 1 Z O a v C W 4 g i W t R z r 5 C 5 7 w S S x I K G Q D j W u u 2 5 E X Q T r I A R T q e 5 T q x p h M k I D 2 j b y B A L q r v J b P W p c 2 l I 4 P S l M i c E Z 0 Z / T y S + L 0 y X w D D U S x w L r S f C / 8 9 L 4 V o v 0 O l F a y 1 t t h / S Y K 2 X E p C S 6 3 V u O 4 b + X T d h Y R Q D D c n 8 V f 2 Y O y C d N H Q n Y I o S 4 B M j M F H M B O O Q I V a Y g P m a n I n c + x v w q m h U y t 5 1 u f J y U 6 z W F u F n 0 R m 6 Q C X k o V t U R U / o G d U R Q W P 0 j j 7 Q Z + b b y l s n 1 u m 8 1 c o s Z o 7 R U l n n P y b k u l Y = < / l a t e x i t > A (t 0 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " w O i E q k B c C r L E 3 q + I O 4 4 u t U 1 w E D c = " > A A A C Z n i c d Z H L T g I x F I b L e E O 8 A B r j w k 0 j M c E N m U E T 3 Y m 6 c Y m J X B J A 0 u k U a G i n k / a M C Z n w J G 7 1 o X w D H 8 M O s B D B k z T 5 8 / 3 n p K d / / U h w A 6 7 7 l X E 2 N r e 2 d 7 K 7 u b 3 9 g 8 N 8 o X j U N C r W l D W o E k q 3 f W K Y 4 C F r A A f B 2 p F m R P q C t f z x Y + q 3 3 p g 2 X I U v M I l Y T 5 J h y A e c E r C o X 8 h 3 h 0 R K 8 p r c T 8 v Q d y / 7 h Z J b c W e F V 4 W 3 E C W 0 q H q / m L n r B o r G k o V A B T G m 4 7 k R 9 B K i g V P B p r l u b F h E 6 J g M W c f K k E h me s l s 8 y m + s C T A A 6 X t C Q H P 6 O + J x P e l 7 Z I E R m a J E 2 n M R P r / e S l c 6 w U m v W i t Z e z 2 I x a s 9 V I C S g m z z u 3 E M L j t J T y M Y m A h n b 9 q E A s M C q e Z 4 4 B r R k F M r C B U c x s M p i O i C Q X 7 M z k b u f c 3 4 F X R r F a 8 q 0 r 1 + b p U e 1 i E n 0 V n 6 B y V k Y d u U A 0 9 o T p q I I p i 9 I 4 + 0 G f m 2 z l 0 T p z T e a u T W c w c o 6 V y 8 A 9 7 3 7 o K < / l a t e x i t > B t 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " X / h 9 F 9 Q c 7 I I 4 E b x 3 6 6 M Y e n b H u Y 4 = " > A A A C W H i c d Z H L S s N A F I Z P 4 6 U X b 2 1 d u g a a x Y T O y I Q N j I x I y P R o s d p 1 a V 8 Y 4 t u B V O Z E a K / o 7 4 m F 5 4 W mK y Q 4 1 R u c h F r P Q + 8 / L 4 W 5 n q / T i 3 I t b b a f M j / X S w l K K X S e O 0 g w u B s t e B Q n y C K 6 f l W Q C B u l n a Z s + 1 w x i m J u B K G K m 2 B s O i W K U D R / U T G R u 3 8 D 3 h b d q 6 Z 7 3 b x 6 v m m 0 H r L w S 3 A G 5 3 A J L t x C C 5 6 g D R 2 g M I F 3 + I D P w r c F V t E q r 1 u t Q j Z z C h t l 1 X 8 A v H m 1 G g = = < / l a t e x i t > A EoW brane A (t 0 ) < l a t e x i t s h a 1 _ b a s e 6 4 = " w O i E q k B c C r L E 3 q + I O 4 4 u t U 1 w E D c = " > A A A C Z n i c d Z H L T g I x F I b L e E O 8 A B r j w k 0 j M c E N m U E T 3 Y m 6 c Y m J X B J A 0 u k U a G i n k / a M C Z n w J G7 1 o X w D H 8 M O s B D B k z T 5 8 / 3 n p K d / / U h w A 6 7 7 l X E 2 N r e 2 d 7 K 7 u b 3 9 g 0 9 o T p q I I p i 9 I 4 + 0 G f m 2 z l 0 T p z T e a u T W c w c o 6 V y 8 A 9 7 3 7 o K < / l a t e x i t > AB (t0)< l a t e x i t s h a 1 _ b a s e 6 4 = " 2 0 Q B 9 9 K m Q y m + i M 1 F K q V w r X C S t 3 U = " > A A A C Z 3 i c d Z H L T g I x F I b L e E O 8 M G h i T N y M E h P c k B k 0 0 Z 2 I G 5 e a y C U B J J 1 O g Y Z 2 O m n P k J A J b + J W 38 l H 8 C 3 s A A s R P E m T P 9 9 / T n r 6 1 4 8 4 0 + C 6 X x l r Y 3 N r e y e 7 m 9 v b P z j M 2 4 W j h p a x I r R O J J e q 5 W N

Figure 21 :
Figure21: Phases for the HRT surfaces γ AB (t) (green) and γ A (t) (red) in the holographic pure Bstate.Left: for v d t < f A, the corresponding entanglement entropy S op A (t) grows linearly in t; Center: for f LA < v d t < /2, after a phase transition for γ A (t) itself, S op A (t) settles down to be volume-law; Right: for v d t > /2, due to the phase transition of the underlying bulk state, i.e. the identification surface γ AB (t), γ A (t) is forced to jump discontinuously to be short-ranged.
s 1 N y 6 u y h n U 3 g r U S O r a g 2 r 1 v 0 g l C y N I E Y m q N Z 9 z 0 3 Q z 6 h C z g T M y 4 N U Q 0 L Z l I 6 h b 2 R M I 9 B + t t h 0 7 l w a E j o j q c y J 0 V n Q 3 x N Z E E S m K 6 I 4 0 W u c R l r P o u A / L 4 d b v V D n F 2 2 1 t N l + A u F W L y c o p d D b 3 H 6 K o z s / 4 3 G S I s R s + a p R K h y U T p 6 x E 3 I F D M X M C M o U N 8 E 4 b E I V Z W h + o m w i 9 / 4 G v C l e G n X v u t 5 o 3 9 S a T 6 v w S + S c X J A r 4 p F b 0 i S P p E U 6 h B E g b + S d f F i f 1 r d d s I v L V t t a z Z y R t b I r P 4 F E t X 0 = < / l a t e x i t > AB (t) < l a t e x i t s h a 1 _ b a s e 6 4 = " X x 5 8 T F u q D G Y I 4 T Y y F y 3 x 0U A U m E s = " > A A A C Y 3 i c d Z H L S g M x F I b T 8 V Z b L 2 1 1 J 8 J g E X R T Z q q g O 6 t u X E k F e 4 H e y G T S N j S Z D M k Z o Q x 9 D 7 f 6 V j 6 A 7 2 G m n Y W 1 9 U D g 5 / v P I S d / v J A z D Y 7 z l b E 2 N r e 2 d 7 K 7 u f z e / s F h o V h q a h k p Q h t E c q n a H t a U s 4 A 2 g A G n 7 V B R L D x O W 9 7 k M f F b b 1 R p J o N X m I a 0 J / A o Y E N G M B j U7 4 6 w E L g f 3 z / M L u B y U C g 7 F W d e 9 q p w U 1 F G a d U H x c x d 1 5 c k E j Q A w r H W H d c J o R d j B Y x w O s t 1 I 0 1 D T C Z 4 R D t G B l h Q 3 Y v n a 8 / s c 0 N 8 e y i V O Q H Y c / p 7 I v Y 8 Y b o E h r F e 4 l h o P R X e f 1 4 C 1 3 q + T i 5 a a 2 m z / Z j 6 a 7 2 E g J R c r 3 M 7 E Q x v e z E L w g h o Q B a v G k b c B m k n g d s + U 5 Q A n x q B i W I m G J u M s c I E z L f k T O T u 3 4 B X R b N a c a 8 q 1 Z f r c u 0 5 D T + L T t A Z u k A u u k E 1 9 I T q q I E I U u g d f a D P z L e V t 0 r W 8 a L V y q Q z R 2 i p r N M f O Z 6 5 j g = = < / l a t e x i t > A (t) < l a t e x i t s h a 1 _ b a s e 6 4 = " n 7 / N V k a U M 9 c u w / g 7 3 d u V J Z K p p O Y = " > A A A C Y n i c d Z H L T g I x F I b L e E O 8 g S x 1 M Z G Y 4 I b M o I n u x L j R H S Z y S W A k n U 6 B h n Y 6 t m d M y I T n c K u P 5 d 4 H s Q O z E M G T N P n z / e e k p 3 / 9 i D M N j v O V s z Y 2 t 7 Z 3 8 r u F v f 2 D w 6 N i 6 b i t Z a w I b R H J p e r 6 W F P O Q t o C B p x 2 I 0 W x 8 D n t + J P 7 1 O + 8 U a W Z D J 9 h G l F P 4 F H I h o x g M M j r j 7 A Q +

ln 2 p 3 <
6 2 m k c + c o L U y T 3 8 A C l G 2 p Q = = < / l a t e x i t > 2⇡ l a t e x i t s h a 1 _ b a s e 6 4 = " b 4E h p Q J a O p 8 5 g w 6 Q a X A q r c g V / W g = " > A A A C g X i c d Z H f a t s w F M Y V d / 2 X r m 2 y X v b G L A y 6 Q o P t F F r o R c t 2 s 8 s O l q Q Q h y D L x 4 2 I L L n S c S E I v 0 O f Z r f t a + x t K i e5 W N f s g O D j 9 x 0 h n e 8 k h e A G g + B P w 9 v 4 s L m 1 v b P b 3 P u 4 f 3 D Y a n 8 a

Figure 23 :
Figure 23: Decomposing γ A (t) into two parts: a part outside the horizon (transparent red); a part behind the horizon (solid red).

Figure 24 :
Figure 24: The growth-plateau-drop pattern in three classes of models.