Holographic Quantum Circuits from Splitting/Joining Local Quenches

We study three different types of local quenches (local operator, splitting and joining) in both the free fermion and holographic CFTs in two dimensions. We show that the computation of a quantity called entanglement density, provides a systematic method to capture essential properties of local quenches. This allows us to clearly understand the differences between the free and holographic CFTs as well as the distinctions between three local quenches. We also analyze holographic geometries of splitting/joining local quenches using the AdS/BCFT prescription. We show that they are essentially described by time evolutions of boundary surfaces in the bulk AdS. We find that the logarithmic time evolution of entanglement entropy arises from the region behind the Poincare horizon as well as the evolutions of boundary surfaces. In the CFT side, our analysis of entanglement density suggests such a logarithmic growth is due to initial non-local quantum entanglement just after the quench. Finally, by combining our results, we propose a new class of gravity duals, which are analogous to quantum circuits or tensor networks such as MERA, based on the AdS/BCFT construction.


Introduction
Conformal field theories (CFTs) in two dimensions are very special among quantum field theories (QFTs) in that infinite dimensional conformal symmetries constrain their properties. Owing to this, we can analytically calculate much more physical quantities than those in ordinary QFTs. Nevertheless, we can observe a variety of qualitative difference in the dynamics of two dimensional CFTs (2d CFTs). We can easily come up with two extreme examples of 2d CFTs. One is free CFTs such as massless free fermion or scalar theories. Another one is the strong coupling limits of 2d CFTs with large central charges c, namely the holographic CFTs, which have dual descriptions via the AdS/CFT correspondence [1,2].
One of the most interesting quantities to characterize the dynamical property of a given quantum state is the entanglement entropy (EE) [3,4,5,6,7]. In this paper we would like to explore how differences of 2d CFTs appear in the time evolution of certain excited states. One useful class of excited states is called local quenches, where we exert a local excitation on a vacuum state. Since it initially modifies the state only locally, they provide clean examples where we can interpret the time evolution easier. We consider three different types of local quenches: (i) local operator quenches, (ii) splitting local quenches, and (iii) joining local quenches. We sketched them in Fig.1.
The first one (i) is simply defined by acting a local operator on a vacuum and we analyze its time evolution. This was first introduced in [8,9] and there have already been many related results both in field theory analysis and holographic analysis. The third one (iii) is defined by joining two semi-infinite lines, as first introduced in [10]. The second one (ii) is triggered by splitting a connected line into two disconnected ones. This is a new setup which the present paper will discuss in detail.
In this paper we will point out that it is very helpful and systematic to analyze not the entanglement entropy itself but its second derivatives, called the entanglement density (ED) [11]. Indeed, this quantity extracts the essential behaviors of entanglement entropy under local quenches in term of clear peaks in its graph. By studying the behavior of entanglement density, we can clearly see both similarities and differences among the above three types of local quenches in free CFTs and holographic CFTs, as we will explain later. In particular, this observation resolves an apparent puzzle on the known logarithmic growing entanglement entropy under holographic operator local quenches.
Another purpose of this paper is to explore the spacetime geometries of holographic local quenches. Since the gravity dual geometry of (i) local operator quenches was already given in [11], we will focus on (ii) and (iii). For (iii) joining local quenches, the construction of gravity dual by employing the AdS/BCFT formulation [12] was given in [13] and we will study more details of the spacetime geometry by using this description in this paper. In addition we will provide the gravity dual geometry for (ii) splitting local quenches. Refer to [14] for other classes of topology changing quantum operations in CFTs such as projections and partial identifications. Our details geometric analysis clearly will explain the two different sources of logarithmic time evolutions of entanglement entropy observed in the local quenches.
Our analysis of holographic geometry of local quenches is also motivated by the conjectured connection between the AdS/CFT and tensor networks [15] such as the MERA [16]. In tensor network descriptions, we normally consider discretized lattice theories whose continuum limits correspond to CFTs. Therefore, it is not directly related to the continuous AdS spacetimes, but to their discretized versions as in [15,17,18]. One way to resolve this problem is to consider a continuous tensor networks as in the continuous MERA [19,20] or path-integral approaches [21,22,23,24,25].
However, it is also intriguing to try to realize a discretized version of AdS in a way that it is naturally derived from the conventional AdS/CFT. The gravity duals of local quenches are useful for this purpose because we can model a class of quantum gates in tensor networks by arranging the splitting and joining procedures in CFTs. As we will discuss in the final part of this paper, we presents a sketch of gravity duals of MERA tensor networks, by combining the holographic local quenches. This indeed qualitatively supports the conjecture.
This paper is organized as follows: In section 2, we review basic methods of calculating entanglement entropy by conformal map and also by AdS/(B)CFT. In section 3, after we review the definition of entanglement density, we study its behaviors for global and local operator quenches. In section 4, we analyze the splitting local quenches and its holographic dual. In section 5, we analyze the joining local quenches and its holographic dual. In section 6, we explain how the logarithmic growth of entanglement entropy arises from a geodesic length in each gravity dual. In section 7, we provide qualitative tensor network description of local quenches. In section 8, we combine the results in this paper to provide a new gravity dual of a MERA-like (discretized) tensor network, in the framework of AdS/BCFT. In section 9, we summarize our conclusions and discuss future problems. In appendix A, we show the entanglement density can reproduce the correct entanglement entropy even when the subsystem consists of multiple disconnected intervals in two dimensional massless Dirac fermion CFT. In appendix B, we will present the detailed computations of time evolutions of holographic entanglement entropy under splitting/joining local quenches.

Entanglement Entropy and AdS/(B)CFT
Here we briefly review calculations of entanglement entropy (EE) in two dimensional CFTs, based on conformal mappings and computations of holographic entanglement entropy (HEE) based on the AdS/CFT and AdS/BCFT. The entanglement entropy is defined by the von-Neumann entropy S A = −Tr[ρ A log ρ A ], where ρ A is the reduced density matrix defined by tracing out the original quantum state over all parts of Hilbert space other than A.

EE for CFT Vacuum in Flat Space
Consider a two dimensional CFT (2d CFT) on a plane R 2 , which is described by the complex coordinate (w,w). We define the time and space coordinate (τ, x) as The time τ is Euclidean time and is analytically continued to the real time by A two point function of a primary operator O in a 2d CFT behaves as O(w 1 ,w 1 )O(w 2 ,w 2 ) = 1 |w 1 − w 2 | 2(h+h) , (2.3) where the primary operator O has the chiral/anti-chiral conformal dimension given by (h,h). For the twist operator σ n , we have h =h = c 24 (n − 1/n). When we compute the entanglement entropy (EE) for the subsystem A defined by the interval 0 ≤ x ≤ l at a fixed time τ = 0, the relevant two point function is σ n (l, l)σ n (0, 0) = 1 l c 6 (n−1/n) .

EE for Excited States from Conformal Transformation
To calculate the EE for a special class of excited states in 2d CFTs, we can employ the conformal transformation. We transform the original coordinate system (w,w) into a new one (ξ,ξ): ξ = f (w). (2.8) Also we need to note that the original UV cut off , which is a lattice spacing in the w coordinate, is mapped to the one, called˜ a,b , in the new coordinate (ξ,ξ) as Finally the EE for the excited state is found to be (2.10) In our analysis of joining/splitting local quenches, the EE is computed from a two point function of twist operators in the presence of a conformal boundary, called boundary conformal field theory (BCFT). Since two point functions in BCFT are similar to four point functions in CFTs without boundaries, we do not have any universal expression such as (2.3). However, we can obtain a definite analytical expression for special CFTs such as the Dirac fermion CFT and the holographic CFTs, whose details will be discussed later.

Conformal Map in AdS 3 /CFT 2
The AdS/CFT correspondence argues that gravitational theories on AdS 3 is equivalent to 2d holographic CFTs, which live on the AdS boundary [1]. The physical equivalence is formulated such that the partition function in each side agrees with the other one [2], so called the bulk-boundary correspondence.
Consider the Poincaré metric of AdS 3 , which is dual to the vacuum of 2d CFT (CFT 2 ): where we set the AdS radius to 1 for simplicity. The conformal transformation (2.8) is equivalent to the following coordinate transformation in AdS 3 (see e.g. [26]): (2.12) The metric in the coordinate (w,w, z) reads are the chiral and anti-chiral energy stress tensor.

Holographic Entanglement Entropy in AdS 3 /CFT 2
In Euclidean setups, the holographic entanglement entropy (HEE) [27,28] is given by 15) in terms of the length L of shortest geodesic which connects the two boundary points of the subsystem A. We also need to impose the homology constraint that this geodesic is homologous to the subsystem A in the AdS geometry. If we consider a connected geodesic which connects w a and w b , then its length is given by where note the relation (2.9), following from the above coordinate transformation around η = 0: η |f |z. Then the HEE L ab 4G N = c 6 L ab reproduces the formula (2.10).

EE from AdS/BCFT
Since our coming analysis of joining/splitting local quenches require us to consider a CFT, let us extend the AdS/CFT to a setup where a holographic CFT is defined on a manifold M with boundaries ∂M . In particular, when a linear combination of conformal symmetry is preserved on ∂M , we call it a boundary conformal field theory (BCFT). The holographic dual of a BCFT, called AdS/BCFT, can be constructed in the following way [12] (see also [29] for an earlier argument). For a typical AdS/BCFT setup, refer to the left picture of Fig.2. Consider a surface Q which ends on ∂M and extends into the bulk. We impose the following condition on Q where h µν is the induced metric on Q and K µν is the extrinsic curvature on Q; K is the trace h µν K µν . The constant T BCF T describes the tension of the 'brane' Q and can take both positive and negative values in general. This boundary condition (2.17) arises naturally in the AdS/BCFT setup as follows (for more detail refer to [12]). If we consider the standard gravity action given by the Einstein-Hilbert action plus the Gibbons-Hawking boundary term, it is well-known that there are two boundaries conditions: Dirichlet and Neumann. For AdS/BCFT, we choose the Neumann boundary condition as we want to keep the boundary Q dynamical as is so in the BCFT boundary. The Neumann boundary condition is given by K µν − Kh µν = 0. To generalize this boundary condition, we add the tension term of the surface Q to the bulk action, given by T BCF T Q √ h. This modified the Neumann boundary condition into the form (2.17). It is also useful to note that in explicit examples, we can confirm that this boundary condition (2.17) preserves the boundary conformal symmetry.
The gravity dual of a CFT on M is given by the AdS gravity solution restricted on the space N , defined by the bulk region surrounded by M and Q. To find such a solution, we need to solve the Einstein equation with the boundary condition (2.17), where the presence of Q gives back-reactions and modifies the bulk metric [12,30,31]. In our examples which we will discuss later, we can analytically find gravity duals of two dimensional BCFTs by using the bulk extension of the conformal map (2.12).
The region N gets larger as the tension T BCF T increases and this suggests that T BCF T estimates the degrees of freedom on the boundary ∂M . Indeed, as shown in [12] in the AdS 3 case, the tension is monotonically related to the boundary entropy S bdy introduced in [32] as follows: S bdy = c 6 arctanh(T BCF T ). (2.18) For notational simplicity, we introduce a positive parameter k(> 0) by A larger k means a larger boundary entropy or tension T BCF T . In particular, we have k = 1 for T BCF T = S bdy = 0. Refer to e.g. [33,34] for studies of HEE in higher dimensional setups. Later we will employ the AdS/BCFT to calculate the holographic entanglement entropy (HEE). Consider the holographic entanglement entropy S A for an interval A in the AdS 3 /BCFT 2 setup. As in the right picture, there are two possibility: connected one S con A and the disconnected one S dis A . The latter arises because the geodesic which connects the two end points of A can end on the boundary surface Q in the middle as depicted in the right picture of Fig.2. The correct holographic entanglement entropy is given by the one with a smaller length.
The connected HEE S con A is not affected by the boundary surface Q. Therefore, this simply agrees with the estimation (2.10) of the CFT without boundary. On the other hand, the disconnected HEE S dis A is highly affected by the boundary Q as the geodesic ends on Q. In the Poincare AdS setup, S dis A in Fig.2 is computed as (refer to [12] for the derivation): where s a,b are the distance between the surface Q and the two end points of the interval A;˜ a,b are the UV cut off in the Poincare AdS at the two boundary points. Intuitively, as the surface Q gets further from the subsystem A, then geodesic length (i.e. HEE) gets larger. Thus as the tensor T BCF T gets larger, the HEE gets larger. This effect is universally described by the constant term proportional to S bdy in (2.20) as follows form the result in [12]. Note that this is the sum of two disconnected geodesics and therefore has the doubled contribution of S bdy . We can also understand this from the CFT viewpoints. In the holographic CFTs (or large central charge CFTs with sparse spectrum). As in the standard large c arguments [68], we can approximate a correlation function by a semiclassical saddle point. In our setup with a boundary, there are two saddles. One of them is obtained by contracting two twist operators, which give the connected EE S con A . The other is given by contracting each of them with its mirror operator across the boundary, which leads to the disconnected one S dis A . Both of them agree with the results from AdS/BCFT. The S bdy dependence in (2.20) occurs because the branch cut which extends from a twist operator ends on the boundary as is already known in the standard calculation of EE in BCFT [6].
More generally, including higher dimensional setups, we can calculate the holographic entanglement entropy in the following way. In the standard holographic entanglement entropy without any boundaries, S A is given by the area of minimal or extremal surface Γ A which ends on ∂A and which is homologous to A. In the presence of boundary surfaces Q, we impose the homology condition by regarding the surfaces Q as trivial spaces with the zero size.

Entanglement Density
Here we study a quantity called entanglement density (ED) introduced in [11]. First we would like to note that this is not a new quantity as it follows from the values of the entanglement entropy (EE). However, the ED provides a helpful way to display the behavior of EE in complicated systems, because we can capture the essence of time evolution from the behavior of ED as we explain later.
The ED is defined from the data of EE in a 2d CFT when the subsystem A is an interval. Since we have in mind generic excited states, we do not require the translational symmetry. When A is given by an interval a ≤ x ≤ b at a fixed time t, the EE is written as S A (a, b, t) (let us assume a < b). The entanglement density n(a, b, t) is defined by n(a, b, t) = 1 2 where we introduced the center of the interval 1 ξ and length l such that For example, the ED for a CFT vacuum, denoted by n 0 (ξ, l, t), takes the universal form: By definition, the EE is represented as a double integral of ED: Therefore, if we assume all entanglement comes from the bipartite quantum correlation, the ED measures the number of EPR pairs between x = a and x = b. However we should note that this interpretation is too naive as in general we have multi-partite entanglement in QFTs. This issues becomes serious in the later time behavior of holographic local operator quench [8], as pointed out in [35]. On the other hand, for the massless Dirac fermion CFT, this interpretation works so well that we can reproduce correct results from the universal ED (3.3) even when A consists of multiple disconnected intervals as we will explain in the appendix A. One of the purposes of the present paper is to emphasize that the entanglement density at least provides a useful and simple way to extract the essential behaviors of the entanglement entropy in a systematical way. We will see this in many examples.
We are especially interested in the difference ∆n(ξ, l, t) = n(ξ, l, t) − n 0 (ξ, l, t), (3.5) where n 0 (ξ, l, t) is given by (3.3). In other words, ∆n(ξ, l, t) is the entanglement density for the entanglement entropy growth where S 0A is the entanglement entropy for the ground state.

Properties
As already noted in [11], the growth of ED ∆n(ξ, l, t) enjoys several interesting properties. First of all, owing to the first law of entanglement entropy [36,37,38], in the small size limit l → 0, we have ∆n(ξ, 0, t) = − π 3 T tt (ξ, t), where T tt is the energy density. Moreover, the growth of ED satisfies a sum rule: da db∆n(a, b, t) = 0 if the total state is pure. In [11], this was proved when we impose the periodic boundary condition, which identifies x = a with x = b. More generally, we can prove this even when the total space is an interval 0 ≤ x ≤ L as we will explain below. The precise statement of the sum rule is at any time t. To prove this, we can rewrite it as follows First of all, the quantity ∂ b ∆S A | a=b is vanishing because the first law (3.7) tells us the behavior ∆S A ∝ (b − a) 2 when |b − a| is very small. The second term in the right hand side also follows because where we employed the first law behavior and the pure state property S A = S A C (A C is the complement of A). It is also intriguing to note that we can also define the entanglement density for Rényi entanglement entropy (2.6). Since this quantity also has the property ∆S (n) when |b − a| is very small, the Rényi entanglement density ∆n (n) (a, b, t) also satisfies the sum rule (3.8).

Example1: Global Quenches
One of the simplest but non-trivial examples of homogeneous excited states in CFTs is the global quenches. This is triggered by a sudden change of the Hamiltonian from a gapped one to a critical one H at a specific time t = 0. As argued in [39], we can model this process by approximating the state just after the quench by a regularized boundary state |Ψ(t = 0) = e −αH |B , where α is an infinitesimally small parameter of the regularization and |B is a boundary state [40,41]. In general, The EE of a subsystem with length l shows the linear growth and the saturation [39] This linear growth can be explained by the entangled pair creations and their relativistic propagations [39]. In this case the ED looks like ∆n(ξ, l, t) = πc 24α δ(l − 2t) + ∆n(ξ, l, t) U V . (3.12) The UV contribution ∆n(ξ, l, t) U V ∼ O(T tt ) is localized at short distances l ≤ α such that it obeys the sum rule (3.8).
If we consider the massless free Dirac fermion CFT in two dimension (c = 1) as a solvable example, we explicitly find the following expression of the EE 2 Thus, the ED reads . (3.14) Indeed we can confirm the behavior (3.12) in this explicit example. We numerically plotted the profile of ∆n(l, t) in Fig.3. We can also confirm that the sum rule ∞ 0 dl∆n(l, t) = 0, explicitly. Similar arguments can also be applied to holographic calculations of global quenches [43,44,45,46].

Example2: Local Operator Quenches
Next we consider the time evolution of an excited state produced by a local insertion of a primary operator O at x = 0 at the time t = 0 [8,9]: where α is again an infinitesimally small regularization parameter; N O is the normalization factor so that the state has the unit norm. We call its time evolution a local operator quench (refer to the left picture of Fig.1). We will discuss other types of local quenches (splitting and joining ones) in later sections, which are the main setups we consider in this paper. However, here we briefly discuss the operator quenches because they are instructive for our later arguments and the behaviors of their ED have not been studied well before. Let us focus on two dimensional free CFTs or more generally rational CFTs (RCFTs) such as free scalar fields or minimal models, for simplicity. 3 We again choose the subsystem A to be the interval a ≤ x ≤ b. As found in [8,9,47], the evolution of (Rényi) EE looks like (in the limit α → 0) ∆S (n) for any n, including the EE n = 1. Here we assumed b > |a|. The quantity d O (≥ 1) is called the quantum dimension, which is a finite constant. This behavior is obvious from the picture that the operator insertion creates an entangled pair at x = 0 and then it splits into the left and right moving quanta, which are entangled with each other as argued in [8,9,47]. The entanglement density (ED) is found to be 4 It is also straightforward to reproduce the behavior (3.16) by integrating this entanglement density over a suitable region of (ξ, l) as we explain in the left and middle pictures of Fig.4.
Let us confirm this in an explicit example of a massless free scalar field φ. To obtain an analytical expression, we calculate the growth of 2nd Rényi entanglement entropy ∆S (2) A . We choose the operator O to be O = e iφ + e −iφ . In this case we have d O = 2 and ∆S A = log 2 for a < t < b. This is simply because this operator creates the genuine Bell pair [8,9,47]. The explicit form of ∆S (2) A was obtained in [47] and is given by is the cross ratio of the four locations: Indeed, in the limit α → 0, we can confirm (z,z) → (0, 0) for the regions 0 < t < |a| and t > b , while we have (z,z) → (1, 0) when |a| < t < b. This leads to the result (3.16). We plotted the ED obtained from this Rényi entropy in Fig.5 for α = 0.1. Indeed, we can numerically confirm the behavior (3.18) as well as the sum rule. It is also useful to note that the ED is vanishing completely at t = 0 as the positive and negative delta-functional peaks coincide and cancel each other.

Holographic Case
4 To see this, note that for any a, b, we can write The brown shaded region in the right picture describe the range of (ξ, l) where the local excitation at x = t = 0 can make any physical influence assuming causality. In all these three pictures, the purple dot at (ξ, l) = (0, 2t) represents the peak of ED due to the local quench. In the right picture, the blue curve near this dot (given by l = 2 ξ 2 + t 2 ) describes the delta functional peaks which are peculiar to holographic local quenches. Note that they are outside of the causality region. The calculations of S A for a t b correspond to integrating the red region in the right picture. We note the left part of the blue curve gradually gets into the red region, which gives the logarithmic growth of EE. One may wonder if a similar result is true for local operator quenches in holographic CFTs. A holographic calculation for the operator local quench was given in [11] (see also [67]) and its entanglement entropy was reproduced from CFT analysis in [68]. Especially, the resulting entanglement entropy shows a logarithmic growth when |a| t b. We will explain the geometric explanation in gravity dual in section 6. This logarithmic growth is clearly different from the previous RCFT case (3.16). This looks puzzling because the relativistic propagation picture leads to the behavior (3.16) and there seems to be no way to explain the slow logarithmic growth (3.21) in relativistic theories.
Interestingly, we can resolve this puzzle by looking at the ED instead of EE. Our argument here will be very brief because we will study this more closely in the joining local quench model, which has the similarity on this aspect. The plot of ED shows that the peak at (ξ, l) = (2t, 0) is not localized as opposed to the Dirac fermion case, but it is continuously distributed on a curve, which is depicted as the blue curve in the right picture of Fig.4. We observe that the left part of peak curve (blue) enters into the integration region (red) as time evolves. Thus, we find that S A increases gradually under the time evolution. This qualitatively explains the behavior (3.21).
Note that the continuous peaks are outside of the region which we expect from the causality propagation of the local excitation inserted at x = t = 0. The holographic analysis in [11] shows they are on the curve l = 2 ξ 2 + t 2 . In particular, at t = 0, the peaks are distributed on l = 2|ξ|. This is clear different from the RCFT examples, where the ED vanishes everywhere at t = 0. From this, we learn that in holographic CFTs, the operator local quench (3.15) initially generates highly non-local entanglement, as opposed to what we expect from the naive causality argument. The relativistic propagation of this non-local entanglement is the reason why we find the logarithmic growth (3.21). This property is also expected from the fact that the insertion of local operator is not a local unitary transformation [56]. Note that this initial non-local entanglement can be negligible in the operator local quench in RCFTs, whose ED only has the localized peak at (ξ, l) = (2t, 0).

Splitting Local Quench
Now we move on to the main analysis of this paper. Consider the splitting process in a 2d CFT. Namely, we start with a 2d CFT on a connected line −∞ < x < ∞ and at t = 0 we cut it into two halves at x = 0, as sketched in the middle picture of Fig.1. The preparation of the quantum state at t = 0 just after the splitting process can be done by considering the Euclidean path-integral as in Fig.6. The parameter α corresponds to the regularization of local quench. This setup is described by the conformal transformation:

Splitting Local Quench in Dirac Free Fermion CFT
One example where we can calculate the time evolution of EE under the splitting quench is the massless Dirac fermion CFT in two dimensions (c = 1). We choose the subsystem A to be an interval a ≤ x ≤ b, or equally l = b − a and ξ = (a + b)/2. Due to the parity symmetry x → −x, we can assume without losing generality. In the rest of this paper, we always focus on the entanglement entropy S A for this interval.
In the Dirac free fermion CFT, the twist operator σ n is described by a bosonization σ n = e i n X , where X is the free massless scalar field dual to the Dirac fermion ψ via the standard relation ψ = e iX . Since we know the analytical form of correlation functions of free scalar on an upper half plane, we can calculate the two point function of the twist operators even in the presence of the boundary as follows (this is identical to eq.(2.31) of [14] see also [7,69,42]) : where F is given by . The analytical expressions of S A when α is infinitesimally small are given as follows. In the early time period, 0 < t < |a|, we have When |a| < t < b, we have At late time, t > b, we have First, from (4.5),(4.6) and (4.7), we find that when the subsystem A is symmetric around the origin i.e. ξ = 0, the EE coincides with that for the ground state S (0) A : We also plotted the numerical values of EE in Fig.7. When ξ = 0, we find S A − S (0) A gets non-trivially positive only for the period |a| < t < b. This agrees with the relativistic particle propagation picture of free massless Dirac fermion CFT. For 0 < t < |a|, the EE is the same as that for the vacuum. At the late time region t > b, it approaches to the EE for the ground state of the separated system. This is manifest in the left picture in Fig.7. The middle picture in Fig.7 shows that when we take |a| b, the EE is gradually decreasing just after the sudden initial rise. Indeed, from (4.6), we find the analytical profile for |a| t < b: The right picture in Fig.7 shows that the EE is reduced when either end point of A gets closer to x = 0, i.e. at the splitting point. This is because in this case, the EE only comes from that for a half-line.  Figure 8: The plots of the entanglement density ∆n(ξ, l, t) for α = 0.1 at t = 0 (left two graphs) and t = 1 (right two graphs). The upper and lower graphs are the contour and 3D plots. The horizontal and depth coordinate corresponds to ξ and l, respectively.
The entanglement density (ED) can also be computed from the second differentiation of the EE (3.1). The results are plotted in Fig.8. We can clearly see the positive peak at (l, ξ) = (2t, 0) as well as the negative peaks at (l, ξ) = (0, ±t), similarly to the local operator quench in RCFTs (3.18). The former agrees with the picture of propagating relativistic particles. The latter coincides the peaks of the energy density predicted from the first law of EE. In addition, there are also peaks near the origin ξ = l = 0 and this is due to the cut along x = 0 where the space is divided into the left and right part, which is special to the splitting local quenches.
It is also interesting to note that a negative value region ∆n < 0 is expanding between these three peaks under the time evolution. This is responsible for the decreasing behavior of the EE (4.9) when we integrate the ED over the region given by the middle picture in Fig.4.

Splitting Local Quench in Holographic CFT
Another example which allows us to calculate the EE analytically is 2d holographic CFTs. For the holographic CFT, we can apply the coordinate transformation (2.12) to calculate the EE. The dual Euclidean geometry is given by Im[ξ] ≥ 0 in the Poincaré AdS coordinate (2.11). The Lorentzian time evolution is obtained by performing the analytic continuation.
The holographic entanglement entropy (HEE) for the connected geodesic S con A of the subsystem A(= [w a , w b ]) is computed by using (2.16) as follows: where f (w) is given by (4.1).
On the other hand, the HEE for the disconnected geodesics is found by applying the AdS/BCFT as we explained in eq.(2.20). The final result reads Note that the boundary entropy S bdy contribution arises because the location of boundary surface gets tilted according to the value of the tension T BCF T following (2.18), so that it solves the boundary condition (2.17). The boundary entropy S bdy is parameterized by k as in (2.19). Since S dis is the sum of two disconnected geodesic, we have the doubled contribution of S bdy in (4.11). Refer to the right picture of Fig.2.
The final HEE is given by the smaller of the two: The behaviors of EE are plotted in Fig.9 by choose the vanishing tension T BCF T = 0 or equally vanishing boundary entropy S bdy = 0.
When α is infinitesimally small, we obtain the following analytical expressions (we take the boundary entropy arbitrary) as explained in the appendix B: In the early time period, 0 < t < |a|, we have When |a| < t < b, we have (4.14) At late time, t > b, we have It is clear that at the late time limit t → ∞, the EE approaches to a constant value which agrees with the expected result for the separated two half lines. The behavior of entanglement entropy is also numerically plotted in Fig.9 by choosing the vanishing tension T BCF T = 0 or equally vanishing boundary entropy S bdy = 0. The left graph looks very similar to the one in the Dirac fermion case, which is interpreted by the propagation of relativistic particles. However, it is intriguing to note that if the boundary entropy S bdy is positive and very large, then S con A can dominate in some region and gives a qualitative discrepancy from the Dirac fermion result.
It is interesting to ask the time evolution of EE when we choose the subsystem A to be almost a half of the total system i.e. a t b. If we choose the boundary entropy S bdy is very large such that k 1 (remember the definition of k (2.19)), then S con A is favored for the period a t (2k + 1)a, which results in the logarithmic growth This is peculiar to the holographic CFTs. However, if k is order one, this logarithmic growth is missing and the EE monotonically decreases as in the previous Dirac fermion example.
The right picture (assuming S bdy = 0) in Fig.9 is qualitatively similar to what we obtained for the free Dirac fermion CFT. For finite α > 0, one may notice that the entanglement entropy for the connected curve S con A is not continuous as in the third graph in Fig.9. However, this occurs when one of the endpoints of the subsystem A coincides with x = 0 and thus this discontinuity should happen because the left and right at x = 0 are disconnected in the splitting quench. This can be easily understood if we go to the Poincaré coordinate (2.11).
The entanglement density is plotted in Fig.10. In the upper graphs, we smeared the derivatives in (3.1) by replacing it with the finite difference: e.g. ∂ a S(a, b, t) → (S(a+δ/2, b, t)−S(a−δ/2, b, t))/δ. We did so because in this holographic case, we expect genuine delta-functional behaviors due to the phase transitions of HEE. In addition to the expected positive peak at (ξ, l) = (0, 2t) and the two negative ones at (ξ, l) = (±t, 0), A at α = 0.1. We assume the tension of the boundary surface Q is vanishing T BCF T = 0. The blue/red graph describes the connected/disconnected contribution. The left graph describes the time evolution when we take (a, b) = (15,20). The middle one is the time evolution for (a, b) = (0. 1,20). The right one is the plot for various ξ when we fix l = 2 and t = 0.
we observe continuous peaks along the following two curves where we used the positive parameter k defined in (2.19). These two curves are located at phase transition points between S con and S dis as we can see from (4.14), where (4.18) and (4.19) are situated in a > 0 and a < 0, respectively. Their profile is plotted in the lower left picture of Fig.10. Both of them coincides with l = 2|ξ| at t = 0. In the limits t → ∞ and |ξ| → ∞, the curves are approximated by The presence of these peak curves is peculiar to holographic CFTs and is missing in the free Dirac fermion CFTs or more generally RCFTs. They are outside of the causality zone and therefore the initial state already has highly non-local entanglement. In the lower right picture of Fig.10, we sketched the calculation of EE, which is given by the integration of ED over the red region. For this, we notice that the blue curve (4.19) gets into the red region, while the green curve (4.18) goes out. Thus, we can explain the absence of the logarithmic growth of EE for the splitting local quenches due to the cancellations between these two opposite effects. 5 Also we will provide the geometric explanation of the log t behaviors in section 6.  Figure 10: The behaviors of a smeared entanglement density (ED) for the holographic splitting local quench. The upper left and right graph are the plots of ∆n(ξ, l, t) at t = 0 and t = 2, respectively with α = 0.1 (again we set T BCF T = 0). The horizontal and depth coordinate corresponds to ξ and l, respectively. Here we smeared the originally delta functional behavior of the ED by replacing the second derivatives w.r.t l and ξ with finite-differences δ = 0.2. We observe delta functional behaviors on the two curves (4.18) and (4.19) at k = 1. The lower left picture shows these two curves at t = 2 for k = 0.5, 1, 2.
The lower right picture describes the calculation of S A for a t b by integrating the ED over the red region. The green and blue curves represent curve 1: (4.18) and curve 2: (4.19). We find that the blue curve gets into the red region, while the green curve goes out. The cancellations between them leads to the absence of logarithmic growth.

Holographic Geometry of Splitting Local Quench
Now let us study the holographic geometry dual to the splitting local quench. The metric of the gravity dual is given by (2.13) with where w = x + iτ andw = x − iτ . For its Lorentzian extension, we can set w = x − t andw = x + t. The important ingredient of AdS/BCFT as reviewed in section 2.5 is the boundary surface Q which extends from the AdS boundary z = 0 to the bulk AdS. At z = 0, Q coincides with the cut in Fig.6, which describes the splitting process in the 2d CFT. Therefore, we have to be careful in its global geometry i.e. which part of the geometry we should pick up and where the boundary surface Q in the AdS/BCFT prescription is located. Let us consider the Euclidean geometry dual to the left CFT picture of Fig.6. It is important to note that the gravity dual is precisely defined by mapping an upper half Im[ξ] ≥ 0 of the Poincaré AdS 3 (2.11) using the transformation (2.12) with (4.1), assuming that the tension is vanishing T BCF T = 0. The boundary Q in the AdS/BCFT, given by Im[ξ] = 0 in the latter setup, is mapped to a region on x = 0.
First we focus on the time slice τ = 0 (i.e. Im[w] = 0) of the gravity dual. The metric on this time slice is given by (4.21) We can confirm that this slice is mapped into a quarter of the sphere given by |ξ| 2 +η 2 = 1 with η ≥ 0 and Im[ξ] ≥ 0. However, to realize one to one map, we need to remove the region and identify x with −x for any x along the curve z = 2(x 2 +α 2 ) α . A sketch of the map (2.12) at t = 0 is depicted in Fig.11. The boundary surface Q extends from the AdS boundary z = 0 toward the IR region but it ends at z = 2α due to the identification. This shows that the two half lines x > 0 and x < 0 at the AdS boundary z = 0, is connected in the bulk through the horizon given by the identified curve z = 2(x 2 +α 2 ) α . Therefore the CFTs on these two half lines are entangled and the entanglement entropy is estimated as the length of this curve: where z ∞ = 2x 2 ∞ /α is the IR cut off. This agrees with our expectation because the state before the splitting of the system was the ground state of the CFT on a line. Note that in the Poincaré coordinate, this is equal to a quarter of the largest circle in the sphere |ξ| 2 + η 2 = 1, depicted as the red curve in the left picture of Fig.11.
Next we would like to examine how the boundary Q in the AdS/BCFT looks like in our gravity dual. To see this we analyze which region is mapped into Im[ξ] = 0. It is straightforward to see that the boundary Q should be on the slice x = 0 (i.e. Re[w] = 0). Notice that we need to distinguish the two segments Q + and Q − of the boundary surface Q at x = +δ and x = −δ (δ > 0 is infinitesimal), which corresponds to Re[ξ] < 0 and Re[ξ] > 0 in the Poincaré AdS. We focus on the time period −πα < τ < πα. In Fig.12, we sketched how the boundary surface Q + looks like in the (w,w, z) coordinate and how it is mapped into the boundary Im[ξ] = 0 and Re[ξ] < 0 in the Poincaré coordinate for x = +δ.
The tip z = 2α of the identification region (4.22) at the time slice τ = 0 is extended along for − α 2 < τ < α 2 . This is the red curve in Fig.12. Note that this tip is the end point of the boundary surface Q in the case of T BCF T = 0. Therefore the two sheets (x = +δ and x = −δ) are identified along this curve (4.24) with each other.
There is one more non-trivial issue. To secure one to one mapping into the half of Poincaré AdS, we need to remove the white region in the left picture in Fig.12 and join the two blue curves with each other such that it agrees with the right picture i.e. Im Note that there are one more boundary surface corresponding to x = −δ → 0, which is given by the same structure and is mapped to the Re[ξ] > 0 region. The white region in − α 2 < τ < α 2 should also be removed to secure the map is one to one.
in the Poincaré coordinate. Finally we can confirm that the colored regions in the left picture are mapped into those in the right one in Fig.12. The full boundary surface Q is given by joining two copies of such a space. We will not get into details of the geometry in other parts as we can understand how the Lorentzian space looks like from the above observations. Nevertheless, it is helpful to point out that the identification surface (4.24) at non-zero x is found to be In the Lorentzian signature, this surface looks like Thus the Lorentzian geometry for t > 0 in the (w,w, z) is given by removing the part z > G L (x, t) and by identifying each two points on (4.26) as (τ, x, z) ∼ (τ, −x, z). The boundary surface Q consist of identical two surfaces Q + and Q − which extend in z < 2(α 2 +t 2 ) √ α 2 +4t 2 and are localized at x = δ and x = −δ, respectively. The final spacetime of the gravity dual is depicted in Fig.13.
At late time we can approximate this as z t. (4.27) This behavior z ∼ t at x = 0 clearly shows that the two half lines at the boundary are not causally connected in a marginal sense. Therefore the region x > 0 and x < 0, which are separated at the AdS boundary z = 0 are connected in a space-like way (i.e. it is a non-traversable wormhole) through this horizon in the bulk. This is consistent with the fact that in the CFT side there is no direct interaction between these two CFTs for t > 0, though they are entangled. As in our arguments in this subsection we assumed T BCF T = 0 for simplicity, we would like to briefly mention how the gravity dual is modified when T BCF T = 0. When T BCF T is non-zero, the boundary surface Q in (ξ,ξ, η) will be tilted as When T BCF T < 0, the boundary surfaces Q + and Q − tilted in the bulk such that they still coincide at z = 0 and the angle between them takes a fixed positive value. Therefore, the bulk region gets squeezed and the entanglement entropy S A (4.23) is decreased. When T BCF T > 0, the boundary surfaces are modified in an opposite way and the bulk region expands. The entanglement entropy S A gets increased.
In this subsection we studied the Lorentzian motion of the boundary surface Q. For the holographic entanglement entropy, we actually need the profile of geodesic in Lorentzian geometry, though in the previous subsection we employed the formal Wick rotation arguments to avoid this analysis. We will explicitly study how the Lorentzian geodesic looks like in section 6.

Splitting Local Quench in A Simple Spin System
Before we go on, we would like to present numerical results for a splitting local quench in a simple spin model for a reference. We start with a spin-1/2 transverse Ising chain with L sites and free boundary condition. Cut the interaction between site j and j + 1 to separate it into two independent chains at t = 0. In this case Hamiltonian before quench H bef ore and that after quench H af ter turns out to be Then let us start with the ground state of H bef ore and investigate its dynamics under H af ter at t ≥ 0. We presented our numerical results in Fig.14. We can see the behavior of its EE at an early time is very similar to the result in free fermion (4.6) and the disconnected contribution in holographic CFTs (4.14), qualitatively. At late time we observe oscillatory behavior because our spin system has a finite size.

Joining Local Quench
Consider a 2d CFT on two semi-infinite lines x > 0 and x < 0. We join each endpoint x = 0 at t = 0 as in the right picture of Fig.1. This setup is described by the path-integral  Figure 15: The geometries which realize the joining local quenches. The left figure describes the space on which we perform the Euclidean path-integral. This is mapped into an upper half plane by the map (5.1). The right picture describes the path-integral realization of the time evolution after the joining process happened at t = 0, where the Euclidean path-integral for τ < 0 creates the state just after this process. depicted in Fig.15. We can map the Euclidean space (the left picture) into an upper half plane by the conformal map: The parameter α again corresponds to the regularization of local quench. We choose the subsystem A to be a ≤ x ≤ b at time t as before. This corresponds to where we performed the analytical continuation of the Euclidean time τ = it.

Joining Local Quench in Dirac Free Fermion CFT
Consider a massless Dirac fermion CFT in this joining quench. The time evolution of EE can be computed from the formula (4.3) and (4.4) as in the previous analysis of splitting quench. The analytical expressions of S A at the time t when α is infinitesimally small, are given as follows (we choose A to be the interval a ≤ x ≤ b and we can again assume (4.2) without losing generality). In the early time period, 0 < t < |a|, we have When |a| < t < b, we have It is useful to note that when b = −a > 0, S A takes a constant value for any t: A , as we found also for the splitting quench. The behavior of entanglement entropy is also numerically plotted in Fig.16. The left graph shows that the EE gets larger during only the time period where one of the entangled pair is included in the subsystem A, while the EE is the same as its vacuum value for the other time period. This fact is also true in the above analytical results as well as in all of our numerical plots. This property is consistent with the relativistic particle propagation picture for the massless Dirac fermion CFT.
The middle graph shows a logarithmic growth at an early time. Indeed, from (5.4), we find that when |a| t b, the EE shows the logarithmic behavior: which agrees with the known local quench behavior [10]. The entanglement density (ED) can also be computed by taking differentiations of this entropy. The results are plotted in Fig.17. We can clearly see the positive peak at (l, ξ) = (2t, 0) as well as the negative peaks at (l, ξ) = (0, ±t). The former agrees with the relativistic particle picture. The latter coincides the peaks of the energy density. Even though in the previous splitting case of the Dirac fermion CFT, the negative region ∆n < 0 around (l, ξ) = (t, 0) was expanding, in the present joining case, the expanding region at the same place is positive ∆n > 0. This indeed violates the naive causality constraint (refer to the right picture in Fig.4) and we can again interpret this as the presence of non-local entanglement in the initial state just after the quench. This leads to the logarithmic growth at late time in the joining case (5.6) as opposed to the splitting case. Figure 17: The plots of the entanglement density ∆n(ξ, l, t) for α = 0.1 at t = 0 (left) and t = 1 (right). The horizontal and depth coordinate corresponds to ξ and l, respectively.

Joining Local Quench in Holographic CFT
For the holographic CFT, we can substitute the coordinate transformation (5.1) to (2.12) in order to calculate the EE. The dual geometry is given by Im[ξ] ≥ 0 in the Poincaré AdS coordinate (2.11). The HEE can be computed as (4.10) and (4.11) using the AdS/BCFT, as already analyzed in the earlier paper [13].
The analytical expressions of S con A and S dis A at the time t when α is infinitesimally small, are given as follows (we choose A to be the interval a ≤ x ≤ b with (4.2)). In the early time period, 0 < t < |a|, we have When |a| < t < b, we have In the late time limit t b, we always find S con A = S (0) A dominates as we expect. The behavior of entanglement entropy is also numerically plotted in Fig.18. By comparing the holographic results (Fig.18) for vanishing boundary entropy S bdy = 0 with those for the Dirac fermion's one (Fig.16), we find they agree with each other qualitatively, if we ignore the phase transition behavior in the former. If we take |a| t b, then we find the logarithmic growth: which indeed agrees with the Dirac fermion result (5.6). However if the boundary entropy is large enough for S con A to be favored, then we find a slower logarithmic growth for |a| t b: The entanglement density is plotted in Fig.19. In the upper two graphs, we again smeared the derivatives in (3.1) by replacing it with the finite difference. We observe that in addition to the expected positive peak at (ξ, l) = (0, 2t) and the two negative ones at A for α = 0.1 for the holographic joining local quench. We set the tension vanishing T BCF T = 0. The blue and red curve describes S con A and S dis A , respectively. The left graph describes the time evolution when we take (a, b) = (15,20). The middle one is the time evolution for (a, b) = (0. 1,20). The right one is the plot for various ξ when we fix l = 2 and t = 0.
(ξ, l) = (±t, 0), we find a continuous peak along the following curve (we again defined k as in (2.19)): The profile of this curve is plotted in the lower left picture of Fig.19. This continuous peak on (4.19) arises since there is a phase transition from S dis to S con as we can see from (5.8). At t = 0, this curve is reduced to the line l = 2k k+1 |ξ|. In the limit t → ∞ and |ξ| → ∞ we have The presence of the codimension one peak is again peculiar to holographic CFTs and is missing in the free Dirac fermion CFTs or more generally RCFTs. Its presence shows that initial entanglement gets modified in highly non-local way. Moreover, this time, the curve (5.12) extends to the region outside of the region l > 2|ξ| (i.e. the red region). Therefore the blue curve enters into the red region from both the left and right side as in the lower right picture of Fig.19. This explains the logarithmic growth of S A (5.10) with the doubled coefficient compared with that for the local operator quench (3.21). Moreover, if k gets very large, the curve (5.12) approaches to l 2|ξ| in the limit |ξ| → ∞. Therefore the situation gets very similar to the local operator quench and therefore this explains the behavior (5.11). Also we will provide the geometric explanation of the log t behaviors in section 6.

Holographic Geometry of Joining Local Quench
Now let us study the geometry of the gravity dual for the joining local quench. Since this has many similarities with that for the splitting local quench, our explanation will be  Figure 19: The behaviors of a smeared entanglement density for the holographic joining local quench. The upper left and right graph are the plots of ∆n(ξ, l, t) at t = 0 and t = 2, respectively, with α = 0.1 (again we set T BCF T = 0). The horizontal and depth coordinate corresponds to ξ and l, respectively. Here we smeared the originally delta functional behavior of the entanglement density by replacing the second derivatives w.r.t l and ξ with finite-differences δ = 0.1. We observe a delta functional behavior on the curve (5.12) at k = 1. The lower left picture describes the curve (5.12) for k = 0.5 (blue), k = 1 (orange) and k = 2 (green) at t = 2 as well as the line l = 2|ξ| (red). The lower right picture shows the computation of S A for a t b by integrating the ED over the red region. The blue curve represents the peaks on (4.19). As the time evolves, the blue curve enters into the red region from both the left and right side. This explains the logarithmic growth of EE with the doubled coefficient.
brief. Indeed since the energy stress tensor (4.20) remains the same, the metric of gravity dual for the joining quench is the same as that for the splitting one. The difference comes from the location of the boundary surface Q. The time slice t = 0 in the (w,w, z) coordinate is again mapped into the quarter sphere |ξ| 2 + η 2 = 1 in the half Poincaré AdS: Im[ξ] > 0. However the details of the mapping is different from the previous one as sketched in Fig.20. In the present case, the region (4.22) corresponds to Im[ξ] < 0 and thus should be removed. Its boundary (the red curve) given by z = 2(x 2 + α 2 )/α now represents the boundary surface Q. The region −iα < τ < iα and x = 0 is now mapped into a region Re[ξ] = 0 in the half Poincaré AdS 3 . Since this map is identical to that for the splitting quench i.e. Fig.12, we will omit its detail.
Finally we obtain the spacetime geometry sketched in Fig.21. The boundary surface Q, expressed by the red surface, is given by the same expression (4.25) and (4.26) in the Euclidean and Lorentzian case, respectively. In the Euclidean geometry, there is the green surface Q which extends from the CFT boundary and this is actually identified with the red surface Q, whose details we will omit as they are involved with complicated numerics. In some sense, in this geometry, the location of the boundary surface Q is opposite to that for the splitting local quenches.
In summary, the time evolution after the joining local quench is described by the spacetime with the boundary surface Q which moves toward the horizon at the speed of light. 6 This intuitively agrees with our expectation that after the joining, the CFT state gradually approaches to the vacuum state on a connected line. It is also intriguing to note that in the Lorentzian description, the boundary surface Q does not end on the AdS boundary but is localized in the bulk.
So far we implicitly assumed T BCF T = 0. For more general cases T BCF T = 0, the AdS/BCFT setup in (ξ,ξ, η) coordinate is given by the boundary surface (4.28). Then we can map Q into the original coordinate (w,w, z) as in Fig.20. When T BCF T > 0 (T BCF T < 0), the boundary surfaces Q is moved in the larger (smaller) z direction. In other words, the physical region of the gravity dual, surrounded by the boundary surface Q and the AdS boundary z = 0, expands as T BCF T gets larger as expected.
In this subsection we studied the Lorentzian motion of the boundary surface Q. For the holographic entanglement entropy, we actually need the profile of geodesic in Lorentzian geometry, though in the previous subsection we employed the formal Wick rotation arguments to avoid this analysis. We will explicitly study how the Lorentzian geodesic looks like in section 6.

Log t Geodesic Length in Holographic Lorentzian Geometry
In our previous calculations of holographic local quenches we mainly worked in Euclidean setups and we took the analytic continuation to Lorentzian time evolutions only for the  final results of entanglement entropy. This procedure went straightforwardly or mechanically due to the simplicity of Euclidean gravity duals, while this raises the question whether we really have sensible Lorentzian gravity duals whose geodesic lengths give the correct holographic entanglement entropy. We will confirm this by explicitly finding out the Lorentzian geodesics in this section. Note that in the true Lorentzian gravity dual, coordinates (t, x, z) should take real values, while the coordinates in Poincare coordinate (ξ,ξ, η) take complex values in general due to the coordinate transformation ( A (5.7)-(5.9) are characterized by two logarithmic behavior log t/ and log t/α. Moreover, we can find the same expression of S A for these different types of local quenches. Our analysis of Lorentzian geodesics in this section, will also explain such logarithmic behaviors in a systematical way.

Late Time Behavior
First note that in the late time limit t 0 b > a, we find First we would like to point out that both are identified with the twice of the geodesic length between the boundary point (t, x, z) = (t, 0, 0) and the tip of the surface Q, given by When t α, we have These geodesics are sketch as Γ ± in Fig.13 and Γ in Fig.21. In the former (splitting quench) case, we consider a connected geodesic and therefore should avoid touching on the boundary surface Q. On the other hand, in the latter (joining) case, the disconnected geodesic should end on the boundary surface Q. Therefore these two should approximately coincide in the late time limit, by setting a, b → 0.
Next we would like to explain why the geodesic length looks like (6.1). Instead of finding a geodesic directly in (t, x, z) coordinate, it is easier to start with the transformed coordinate (Re[ξ], Im[ξ], η). However, we have to note that since we are considering an analytical continuation into the Lorentzian geometry, the coordinates (Re[ξ], Im[ξ], η) take complex values in general, though (t, x, z) should take real values. By remembering the conformal map (5.1) for the joining quench and the holographic transformation (2.12), it is straightforward to find that the correct geodesic is given by 7 ξξ + η 2 = e −iβ , (6.4) where the real parameter β(> 0) is related to the boundary time t 0 at z = 0 via Therefore for t 0 α we have In terms of the original coordinate, (6.4) is equivalent to which is plotted in Fig.22.
At late time, we have z(t) t − t 0 2 + O(1/t) and thus the geodesic extends to infinity almost in a light-like way. In the limit t → t 0 α we have Indeed we can confirm that this is a solution to the geodesic equation in the holographic spacetime, obtained from (4.20) and (2.13), given by (on x = 0): We can estimate the geodesic length of (6.7) as The twice of this length can explain only one of the logarithmic term in (6.1). The missing contribution actually comes from another part of the geodesic which is not covered by the coordinate (6.9). Indeed, if we parameterize the geodesic (6.4) as ξ = e −iβ/2 cos θ, η = e −iβ/2 sin θ, (0 ≤ θ < π 2 ). (6.12) Note that θ = 0 corresponds to the AdS boundary η = 0 and θ = π/2 corresponds to the boundary surface Q on Im[ξ] = 0. Actually from the map (2.12), we find that the trajectory (6.7) for t 0 ≤ t ≤ ∞ only covers a part of the full geodesic (6.12) given by 0 ≤ θ ≤ β/2. The geodesic length for this part is computed as where η =ˆ is the cut off corresponding to the original one z = . The first term in the right hand side is the full contribution which agrees with a half of (6.1). The second term is estimated when t 0 α as log cos(β/2) 1 + sin(β/2) − log t 0 α . (6.15) This shows that the part of geodesic length in the Poincaré-like coordinate patch (6.9) gives (identical to (6.11)) To understand this outside contribution well 8 , let us ignore the back reaction and focus on the Poincaré AdS (restricted on x = 0) by setting f (t, z) = 0 in (6.9). This is embedded in the global AdS 3 (refer to Fig.22) at θ = 0, via the map (refer to e.g. [8]) Here the parameter α is introduced such that the static trajectory at r = 0 is mapped into that for an accelerated particle z = √ t 2 + α 2 in the Poincaré AdS, which imitates the behavior (6.3). At the AdS boundary, the late time limit t 0 α in the Poincaré AdS corresponds to the time T π − α/t 0 and the cut off r ∞ = t 2 0 /(α ) in the global AdS. We are considering a geodesic which starts from this boundary point and ends on the tip of the surface Q, namely the straight line r = 0 in the global AdS 3 . It is clear that this geodesic is given by T =constant and the radial coordinate r runs from r = r ∞ to r = 0. Along this geodesic, the part outside of Poincaré patch is 0 < r < r * , where r * is the solution to √ r 2 + 1 cos T + r = 0, from which we find r * t 0 /α 1. Thus, in this simplified model, the geodesic length in the Poincaré patch and that for the outside read (refer to Fig.22): These reproduce the previous results from the exact geodesic (6.16) and (6.17).

Middle Time Period
Now we move on to the middle time period and we focus on the region a t b to see the logarithmic behavior clearly. Then we find the following coincidences between the three different types of local quenches 9 : On the other hand, the disconnected geodesic length in joining quench behaves differently: The coincidence (6.21) can be easily understood because they are described by a geodesic which connects the two end points, where the existence of slits is not a crucial obstruction, and because the background metric is the same (refer to the brown curve in Fig.22). We can regard the first term (c/6) log(t/α) as the geodesic length outside of the Poincaré patch as in (6.17). The second term (c/3) log(b/ ) is simply interpreted as the familiar geodesic length in the Poincaré AdS 3 .
On the other hand, the geodesic length (6.22) is divided into two contributions: the one which connects x = a and Q and the other one which connects x = b and Q. The former one is clearly equal to (c/6)(L poincaré + L outside ) in (6.16) and (6.17). The latter one is simply given by the standard result of AdS/BCFT: (c/6) log(b/ ). These explain the behavior (6.22).

Tensor Network Interpretation of Local Quenches
Here we would like to present a sketch of tensor network description of time evolution of our local quench states and their entanglement entropies so that they matches with our gravity duals. Our argument in this section will be heuristic and qualitative.
We employ the MERA (Multiscale Entanglement Renormalization Ansatz) [16] for our qualitative interpretation, though our arguments can equally hold for other holographic tensor networks. MERA is a class of tensor networks which are constructed by two kinds of tensors: isometry and disentangler (Fig.23). This is manifestly scale invariant and is known to describe ground states of critical quantum systems with high accuracy. In tensor networks, by contracting all interior tensors, a desired quantum state is realized at the boundary. We regard the physical system as the boundary, and the tensor network structure as the bulk in AdS/CFT [15].
Consider a curve γ in the bulk, whose length can be defined by the number of tensors it crosses 10 . For a subsystem A on the boundary, its EE S A is upper bounded by where γ min A is the minimal curve, whose length is proportional to log L A , with L A the size (i.e. the number of spins) of subsystem A.
A qualitative correspondence between MERA and AdS is as follows [15]. If we set the coordinate x parallel to the layers, and the coordinate w = log 2 (z/ ) vertical to them (the right picture in Fig.23), then the geometry of MERA is interpreted as a discrete version of the time slice of AdS (i.e. hyperbolic space), whose metric is where is introduced as an UV cutoff. Refer to [25] for a possible interpretations of various slices in AdS spacetime as tensor networks. Now we would like to consider an tensor network interpretation of our three types of local quenches (local operator/splitting/joining). We will take a heuristic strategy: first we assume specific structures of tensor networks under the time evolutions and then confirm if they reproduce the correct time evolutions of entanglement entropy. In our construction Figure 24: A tensor network description of local operator quench at t = 4 . The blue dot is the large modification of tensor dual to the falling particle. The blue curves describe the image of shock wave deformation of tensors. The parameters of tensors are expected to be modified, and tensors close to the blue curve are expected to accept the strongest modification.
of tensor networks corresponding to local quenches, we assume that an appropriate time slice t = const. of a gravity dual spacetime after a local quench, is approximately a tensor network at the time t. We note that the information of tensor network at each time can be separated in to two ingredients: one is the geometric structure of tensor network itself, and another one is the parameters of each tensor. 11 One essential observation we will make below is that we can divide the growth of EE into a "shock wave contribution" and "splitting/joining contribution", which can also be suggested from the holographic analysis in the previous section. In the discussion below, we will regard the "shock wave contribution" as modifications of the parameters and "splitting/joining contribution" as modifications of the network structure. We will see that under this assumption, by choosing appropriate tensor network structure and neglecting modifications on the parameters, we can realize qualitative features of the entanglement structure in the AdS/CFT setup.

Local Operator Quench
A local operator quench in a holographic CFT can be regarded as a particle falling towards the AdS black hole horizon on the gravity side [11]. This particle trajectory is analogous to the tip of the boundary surface Q in the splitting/joining local quench. In this picture, its back-reaction on spacetime spreads as a shock wave [71], and the EE of a subregion on the boundary will change when the shock wave crosses its minimal surface.
Accordingly in the tensor network, we expect the local operator quench causes a modification of tensor at x = 0 and w = 0 due to the initial operator insertion, though the structure of tensor network does not change. As the time evolves, this modification propagates into the interior as w ∼ log t. At the same time, tensors are modified along the shock wave. These modifications increases quantum entanglement locally. These are depicted in Fig.24.
The behavior of the EE is consistent with the causal behavior that the EE gets nontrivial only for the time interval a < t < b, as we found both for the RCFT and holographic CFT. In the holographic case (3.21), however, we observed a logarithmic growth Our entanglement density analysis showed such a logarithmic contribution is due to the presence of non-local entanglement at the initial state. Our geometric analysis in the gravity dual shows this contribution comes from outside of the Poincaré patch (6.17). In tensor network, this hidden contribution is expected to come from modifications of tensor parameters, though we cannot figure it out precisely. Since this is owing to the huge back-reactions, we would like to still call this a "shock wave contribution". To find out the details of this contribution, we need to make the conjectured correspondence between the AdS/CFT and tensor networks more precisely beyond qualitative arguments, which is not available at present.

Splitting Quench
Now we moved on to the splitting quench. We split the disentangler in the MERA at x = 0, z = (w = 0), initially (refer to [72] for introducing boundaries in MERA). This is motivated by the gravity dual where the boundary surface extends to the bulk as in Fig.13. We identify the boundary in the gravity dual as the termination of tensor network (i.e. cutting the disentanglers or unitaries). This elimination of disentanglers propagates under the time evolutions as in the gravity dual. Therefore at the time t = 2 n , the disentangler at x = 0, w = n is eliminated. At the same time, from the gravity dual, we expect there are shock wave propagations as in the local operator quench case. Thus we also expect that in our tensor network description a shock wave spreads from w = n at the same time t = 2 n . Both of them will contribute to the change of EE. We call them "shock wave contribution" and "splitting contribution" respectively.
From the CFT viewpoint, we can heuristically understand this prescription as follows. When we triggered the splitting quench, we eliminate the nearest neighbor interactions between two adjacent lattice sites. This is equivalent to cutting the disentangler in the very UV. After the time evolution this cutting of entanglement propagates towards the IR.
We show in Fig.25 how the network changes after splitting quench. If we focus on the subsystem A = {x|x ∈ (0, l )} and neglect modifications of tensor parameters (equally neglect "shock wave contribution"), we will find S A is given by the length of minimal disconnected curve, which decreases logarithmically in time at first, and then stays as a constant. In other words, splitting contribution gives a negative logarithmic growth ∝ − log t/ to S dis A . On the other hand, as we have mentioned in the case of local operator quench (7.4), the shock wave contribution should give a positive logarithmic growth ∝ log t/α to EE. In this way, we reach the following estimation: log t splitting +..., (7.5) where we determined the coefficient of the second logarithmic term so that we do not totally have any logarithmic growth as in our previous calculation of HEE (4.14).

Joining Quench
For the joining local quench, we start with two copies of semi infinite MERA, which are disconnected initially. The quench is triggered by adding a new disentangler on x = 0, z = (w = 0), at t = . Then, besides a shock wave spreading at w = n in the bulk, a new disentangler on x = 0 comes out to connect the nth layer at t = 2 n . Both of them will contributes to the EE evolution. We call them "shock wave contribution" and "joining contribution" respectively. We show in Fig.26 how the network will change after joining quench. If we focus on block A = {x|x ∈ (0, l )} and neglect modifications of tensor parameters (equally neglect "shock wave contribution"), we will find S A is firstly given by the length of minimal disconnected curve, which increases logarithmically in time. In other words, the joining contribution gives a positive logarithmic growth to S dis A , which is expected to be equal to the absolute value of splitting contribution, since they should cancel with each other if they occur at the same time. On the other hand, since the shock wave contribution should give a positive logarithmic growth that cancels with split contribution, it is proper to think that the two contributions together give a double logarithmic growth, which is exactly consistent with the results in Section 5.2. That is, In this way, we can explain the logarithmic growth we saw in (5.8) Figure 26: Time evolution of the tensor network after joining quench. The left, middle and right picture correspond to t = , 2 and 4 , respectively. The red curve shows the minimum curve in the bulk. The blue dot and green curve are dual to the falling particle and boundary surface Q. Here, A = {x|x ∈ (0, 4 )}.

Holographic Quantum Circuits from AdS/BCFT
We can simply summarize gravity dual spacetimes of the splitting and joining local quench as in the left and middle picture of Fig.27 in the limit α → 0. If we combine them such that we join two CFTs at t = t 1 and split them again at t = t 2 , then the expected holographic dual is given by the right picture in Fig.27. The important lesson from this example is that we can create a slit (S in that picture) which looks "floating" in the bulk Figure 27: A simplified sketch of the gravity dual spacetimes for the splitting (left), the joining (middle) local quench and their combination (right). The brown surfaces describe the boundary surfaces Q. In the right picture, the boundary surface Q intersects with the blue horizontal plane which describes the time slice t = 0. Their intersection S makes a slit on the time slice.
on a time slice of gravity dual. Of course, the location of such a floating slit moves under the time evolution following the condition (2.17).
As an ambitious attempt which finalizes this paper, we would like to propose a holographic counterpart of tensor network states in the AdS/BCFT setup by repeating such procedures. 12 Consider the MERA tensor network [16] and try to realize an analogous geometry on a time slice of a gravity dual. As we explained in section 7, the MERA network consists of the disentangling and coarse-graining operations, which cut the entanglement in the ground state of a critical spin chain. We depicted a simple setup to achieve this in Fig.28. In this analogy, each "spin" in MERA corresponds to the (blue) region between two (brown) slits. A coarse-graining and disentangler correspond to terminating and creating a hole in a (brown) slit, respectively. The standard estimation of EE S A in the MERA nicely agrees with the HEE calculation in the AdS/BCFT as we show in Fig.28, where we simply ignored the back-reactions. The entanglement entropy reduced by the coarse-graining/disentangler operation correspond to the green/red dotted line in the pictures.
Notice that in this section we consider a different interpretation of the AdS spacetime in terms of tensor networks, as compared with the arguments in section 7. Here we approximately regard the narrow strip in AdS space as a link in a tensor network and fill the gap between the discretized lattice network (such as MERA) between the continuous AdS spacetime by cutting the latter spacetime into lattices by inserting various boundary surfaces of AdS/BCFT. This idea opens up a new approach to the conjectured connection between the AdS/CFT and tensor networks.
Next we discuss how to realize a gravity dual whose time slice is given by the geometry in Fig.28. One way is to perform an Euclidean path-integration on the manifold given by the left picture in Fig.28 with z interpreted as the Euclidean time −τ and consider its gravity dual. Note that the initial state z = −τ = ∞ is expected to be a completely , and the estimation of EE S A in the AdS/BCFT and tensor network (right two pictures). The blue region in each picture is the bulk time slice and we added its tensor network description in its next right. The brown thick lines are the slits created by the boundary surface Q in the AdS/BCFT. The red lines and green triangles in the tensor network describe the disentanglers and the coarse-graining tensors. The red and green dotted lines are the geodesics whose length compute the holographic entanglement entropy, corresponding to the entanglement removed by the disentanglers or coarse-graining tensors. We ignore the back-reactions just for simplicity and choose the locations of the slits such that the values of the EE for each red/green dotted lines are the same, which leads to the scale invariance. By comparing the right two pictures, we find that the calculation of HEE in our AdS/BCFT setup is equivalent to that for the tensor network. disentangled state. Therefore we choose the state at τ = −∞ to be the (regularized) boundary state |B [74,75]. Indeed, we expect that the time slice τ = 0 of such a Euclidean gravity dual is given by what we want, according to what we learn from the holographic local quenches summarized in Fig.27.
Another idea is based on a Lorentzian path-integration in the dual CFT. Now we regard z as the real time −t in Fig.28 and again choose the initial state at t = −∞ to be the boundary state |B . Since each end points of the slits moves toward the horizon at the speed of light, we again expect that the gravity dual of this time-dependent setup is approximately given by the one in Fig.28.

Conclusions
In this paper, we studied the time evolutions of entanglement entropy (EE) under three types of local quenches (local operator, splitting and joining) in two dimensional CFTs.
Our main examples of CFTs are the free massless Dirac fermion CFT and holographic CFTs in two dimensions. The typical behavior of EE for each of the three types of local quenches in holographic CFTs is the logarithmic growth for a large subsystem. On the other hand, for operator local quenches in free or RCFTs, we observe simple step functional behaviors. However, for splitting and joining local quenches, we still have logarithmic time evolutions even for free or RCFTs. We also noted that there are two types of logarithmic growth: (c/6) log(t/α) and (c/6) log(t/ ), where α is a regularization parameter of the local quench and is the lattice spacing of the CFT. In this paper, we got systematic understandings of these differences from various points of view, including entanglement density (ED), holographic geometry, and tensor network descriptions. We also calculated evolution of EE in spin systems for splitting quenches.
In the holographic CFT case, the EE is computed as the length of geodesics in the dual AdS 3 spacetime based on the holographic entanglement entropy (HEE). Furthermore, in the splitting/joining quench case, these geodesics can end on the boundary surfaces Q, following the AdS/BCFT prescription. Therefore there are two kinds of geodesics: connected and disconnected. Thus the HEE is given by min{S con A , S dis A }. One of them, S dis A , depends on the tension parameter T BCF T in AdS/BCFT. In the splitting and joining quench, both gravity duals share the same metric but have different boundary surfaces Q. In order to see the logarithmic growth of HEE, we have to consider a region which is not covered by the Poincaré patch because the geodesic penetrates the Poincaré horizon. We found that the growth (c/6) log(t/ ) comes from the geodesic length in the Poincaré patch, while the other one (c/6) log(t/α) comes from that near the surface Q, which is hidden inside the Poincaré horizon.
In the splitting quench, the boundary surface Q expands from the AdS boundary toward the bulk at the speed of light. We can regard the two semi-infinite lines at each time are connected through an expanding horizon or equally (non-traversal) wormhole. Note that since Q is time-like, there is no causal influence between them. In the joining quench, the boundary surface Q can be regarded as a "falling string" towards the horizon, which is analogous to the falling particle in the holographic local operator quench. In both examples, the precise location of Q depends on the tension T BCF T in a way that when T BCF T > 0 (or T BCF T < 0), the region of gravity dual expands (or shrinks).
Entanglement density (ED) turns out to be a good tool to systematically probe entanglement structures between two space points. The analysis with the ED helped us to figure out that the initial entanglements created by our "local" quenches is highly nonlocal, except for the local operator quench in RCFTs. In general, the "local" quenches are not described as local unitary transformations. In holographic case, the ED can straightforwardly explain logarithmic time evolutions under different local quenches. However, the two evolutions (c/6) log(t/α) and (c/6) log(t/ ) cannot be easily distinguished from this perspective.
The gravity duals of local quenches also have a qualitative interpretation in terms of MERA-like tensor networks, where ∼ log(t/ ) behavior can be understood as splitting/joining disentanglers. However, we could not find a clear realization of the ∼ log(t/α) behavior, called "shock wave" contribution. One possibility is to explain this by inserting extra tensors into the networks and we would like to leave more details for a future problem.
In the final part of this paper, we presented an analogue of quantum circuits using both of splitting and joining quenches and gave their gravity duals based on the AdS/BCFT. This offers a method of discretizing gravity dual spacetimes. It is expected that this can further lead to a deeper understanding of the conjectured AdS/tensor network correspondence.

A Vacuum Entanglement Density in Massless Dirac
Fermion CFT Here we would like to show that the simple profile of entanglement density (we set c = 1 in (3.3)) n 0 (x, y) = 1 6(x − y) 2 , (A.1) reproduces arbitrary entanglement entropy in the massless Dirac fermion CFT, including cases where the subsystem A consists of multiple disconnected intervals: A = I 1 ∪ I 2 ∪ · · · ∪ I n . (A.2) We choose the n intervals are parameterized by I i = [a i , b i ] for i = 1, 2, · · ·, n. For convenience, we introduce b 0 = −∞ and a n+1 = ∞. Then, by summing all bipartite entanglement, using the entanglement density (A.1), S A is estimated as follows: This indeed agrees with the known expression of entanglement entropy in the massless Dirac fermion CFT [7]. A closely related property is that the multi partite mutual information is vanishing [76].

B Detailed Computations of Evolutions of Holographic Entanglement Entropy
Here we present the details of calculations of HEE under splitting/joining local quenches. We define the Euclidean coordinate (τ, x) and complex coordinate (w,w) as Both splitting/joining quench setup are described by almost the same conformal maps: (splitting) f + (w) = i ζ(w) arg(ζ) ∈ (−π, π] (B.2) (joining) f − (w) = i −ζ(w) arg(−ζ) ∈ (−π, π]. (B.3) Here we define ζ(w) by The HEE for the connected/disconnected geodesics S con A,± , S dis A,± for the subsystem A(= [w a , w b ]) is computed by (when the tension T BCF T is vanishing): S dis A,± = c 12 log 16(Imf ± (w a )) 2 (Imf ± (w b )) 2 4 |f ± (w a )| 2 |f ± (w b )| 2 . (B.6) The plus and minus subscript in S A,± corresponds to splitting and joining quench, respectively. When computing these quantities, it will be helpful to use the results below.