Time Evolution of Entanglement Negativity from Black Hole Interiors

We investigate the time evolution of entanglement negativity following a global quench for mixed state configurations of two disjoint and adjacent intervals in a ($1+1$)-dimensional conformal field theory ($CFT_{1+1}$) dual to the eternal black hole sliced in half by an end of the world brane, through the $AdS_3/CFT_2$ correspondence. To this end we obtain the time evolution of the holographic entanglement negativity for such mixed states from a dual bulk eternal black hole geometry and elucidate the relevant geodesic structures. The holographic entanglement negativity for such mixed states, following a global quench is described by half of the results for the eternal black hole. Significantly, in this case our holographic results match exactly with the corresponding $CFT_{1+1}$ computations.


Introduction
Over the last few years quantum entanglement has emerged as a central issue in the study of diverse physical phenomena ranging from quantum many-body systems in out-of-equilibrium to the process of black hole formation and information loss problem. The entanglement for a bipartite pure state is characterized by the entanglement entropy which is the von Neumann entropy of the reduced density matrix. In (1 + 1)-dimensional conformal field theories (CF T 1+1 ) the entanglement entropy may be computed through a replica technique as described in [1,2]. However entanglement entropy fails to characterize mixed state entanglement as it typically involves correlations irrelevant to the entanglement of the specific mixed state. This subtle issue was addressed in quantum information theory by Vidal and Werner in [3] where the authors proposed a computable measure termed as entanglement negativity which characterizes the upper bound on the distillable entanglement for the mixed state 1 . The entanglement negativity is defined as the logarithm of the trace norm for the partially transposed density matrix with respect to a subsystem. In a series of interesting communications the authors in [5][6][7] described the computation of the entanglement negativity for bipartite mixed states in a CF T 1+1 through a suitable replica technique.
A holographic characterization of the entanglement entropy was proposed by Ryu and Takayanagi [8,9] through the AdS/CF T correspondence. According to their conjecture the universal part of the entanglement entropy of a subsystem in a dual CF T d was given by the area of a co-dimension two static minimal surface in the bulk AdS d+1 geometry homologous to the subsystem. Naturally, these led to intense research activity in various aspects of entanglement related issues in holographic CF T s [10][11][12][13][14][15][16][17] (and references therein). In a series of interesting communications, a proof of this conjecture was developed from the bulk perspective, initially in the context of AdS 3 /CF T 2 and later generalized to the AdS d+1 /CF T d scenario [18][19][20]. A covariant extension of the RT conjecture was proposed for the entanglement entropy of a subsystem in a CF T d dual to a non static bulk AdS d+1 geometry in [21] and the corresponding proof in [22].
The above developments naturally lead to the interesting issue of a holographic characterization of the entanglement negativity for bipartite pure and mixed states in dual CF T s [23,24]. In the recent past two of the present authors (VM and GS) in a collaboration, proposed a holographic construction for the entanglement negativity of bipartite pure and mixed states in a dual CF T 1+1 through AdS 3 /CF T 2 framework and its covariant extension [25,26]. This construction was rigorously substantiated through a large central charge analysis employing the monodromy technique in [27] although a bulk proof for this conjecture along the lines of [28] is an outstanding issue. Subsequently, their proposal was also extended to higher dimensions in the context of the AdS d+1 /CF T d in [29]. It should be mentioned here that a proof of the higher dimensional extension for the holographic entanglement negativity conjecture described above from a bulk perspective along the lines of [20] is still an outstanding open issue. Subsequent to this in [30,31], a holographic characterization for the entanglement negativity of a mixed state of adjacent intervals in CF T 1+1 s dual to bulk AdS 3 geometries, and its covariant generalization was proposed. A higher dimensional generalization of the above holographic construction for such mixed states of adjacent subsystems in a CF T d dual to bulk AdS d+1 geometry, was subsequently advanced in [32,33]. Furthermore very recently a holographic construction for the entanglement negativity of mixed states of two disjoint intervals in a dual CF T 1+1 and its covariant extension have been developed in [34,35].
In a somewhat different context the study of quantum systems out of equilibrium to elucidate their dynamical evolution has seen a surge of interest in recent times [2,[36][37][38][39]. Quantum quench is one of the processes through which a quantum system may be placed in an out of equilibrium configuration. In a global quench scenario a quantum system is initially prepared in the ground state of a translationally invariant Hamiltonian H 0 and is then allowed to evolve unitarily with respect to another Hamiltonian H, where H and H 0 are related to each other through an experimentally tunable parameter. In a CF T 1+1 the resulting state after a global quench may be regarded as a boundary state (B-state) as it is translational and conformal invariant below the energy scale of the global quench. The corresponding time evolution of the entanglement entropy after a global quench for a subsystem maybe computed through the replica technique [2,36,37] in a CF T 1+1 . It was shown that in the spacetime scaling limit 2 the entanglement entropy grows linearly for early times and becomes extensive beyond a certain range. It is possible to describe such a time evolution of the entanglement following a global quench, through a quasi-particle scenario for the entanglement propagation [36,37,43,44]. In this instance the initial state which is at a high energy relative to the ground state of the post-quench Hamiltonian, acts as a source of quasi-particle excitations.
Following the developments described above the corresponding time evolution of the entanglement negativity for mixed states of two adjacent and disjoint intervals in a CF T 1+1 following a global quench was investigated in [40]. Very interestingly they demonstrated that similar to the entanglement entropy, the time evolution of entanglement negativity also admits of a quasiparticle explanation for the mixed state entanglement propagation. However it was observed that the initial increase in the entanglement negativity occurs slightly after the time predicted by the quasi-particle picture as a consequence of the late birth of entanglement in a lattice system which disappears in the continuum limit [40].
On the other hand the time evolution of holographic entanglement entropy after a global quench in a CF T dual to a bulk Vaidya-AdS geometry describing black hole formation has been extensively studied in the literature [21,[45][46][47][48][49][50][51][52][53][54]. Although as discussed in [55], the precise holographic dual of the global quench scenario where the initial state is a time dependent pure boundary state (B-state) in a CF T 1+1 as described in [36] is the eternal black hole geometry sliced in half by an end-of-the-world (ETW) brane or simply a single sided black hole 3 . It is well known from [55,57,58] that the eternal black holes are dual to the thermofield double state in CF T s. These involve two-sided Penrose diagram and are static under time evolution as the time directions are taken to be opposite for the two sides. However, the time dependence may be introduced for this geometry by considering forward time evolution for both the exterior regions of the Penrose diagram.
The computation of the corresponding holographic entanglement entropy for a subsystem A in a CF T for thermofield double state dual to the bulk eternal black hole geometry requires the consideration of the subsystem in question on both sides of the Penrose diagram. The entanglement entropy may then be computed from the area of the extremal surface anchored on both the subsystems in accordance with the RT and the HRT conjectures. It was observed in [55] that the eternal black hole in AdS 3 may be mapped to a BTZ black brane geometry. In this case it was demonstrated that for early times the extremal surface passes through the interior region of the black brane geometry connecting one asymptotic region to the other. Hence the linear growth of the entanglement entropy was related to the growth of the extremal surface along the nice slice 4 in the interior region of the black brane. For late times the extremal surface reduces to the static minimal surface located exterior to the horizon and the entanglement entropy is rendered extensive [55]. It is argued in [55], the holographic entanglement entropy of the configuration in question after a global quench scenario as discussed in [36] is just half of the corresponding holographic entanglement entropy for the eternal black hole geometry.
Quite naturally the developments described above leads us to the interesting issue of a holographic characterization for the time evolution of the entanglement negativity for mixed state configurations following a global quench in a dual CF T 1+1 . As mentioned earlier the correct holographic dual for the global quench scenario is described by an eternal black hole geometry sliced in half by an ETW brane. Hence to this end it is required to obtain the holographic entanglement negativity for such mixed states from the dual eternal black hole geometry. In this context the covariant holographic characterization of the time dependent entanglement negativity for various mixed state configurations in a dual CF T 1+1 [31,35] maybe utilized to investigate this significant issue.
In this article we address the significant issue described above and investigate the time evolution of holographic entanglement negativity for mixed state configurations of disjoint and adjacent intervals following a global quench in a CF T 1+1 . For this purpose it is necessary to consider the bulk eternal black hole geometry with two copies of the mixed state configurations in the dual CF T 1+1 and its thermal double. The holographic entanglement negativity for the above configuration may then be computed from the prescription described in [31,35]. We observe an interesting structure of bulk geodesics which connects both the asymptotic regions which evolves with time to a geodesic structure that is located in the exterior region. The corresponding holographic entanglement negativity for the mixed state configurations following a global quench in a CF T 1+1 is then given by half of that obtained for the eternal black hole geometry. It is observed that the holographic entanglement negativity for the mixed state configurations in question exhibits a linear growth during the initial part of its time evolution. For late times it decreases linearly with an identical slope from a maximum value and goes to zero after a certain value of time which indicates the process of thermalization and black hole formation in the bulk [44,48,49,60]. Remarkably the holographic entanglement negativity of the mixed state configurations in question, computed through our construction exactly reproduces the corresponding CF T 1+1 replica technique results following a golbal quench as described in [40].
This article is organized as follows. In Sec.
(2), we briefly review the time evolution of entanglement entropy for a single interval after a global quench in a CF T 1+1 . The evolution of entanglement negativity for two disjoint and adjacent intervals after a global quench in a CF T 1+1 is reviewed in Sec. (3). In Sec. (4) we review the time evolution of the corresponding holographic entanglement entropy for a single interval from a dual bulk eternal black hole geometry. In the next Sec. (5), we establish our construction for the time evolution of the holographic entanglement negativity for two disjoint and adjacent intervals from the dual bulk eternal black hole geometry. The holographic entanglement negativity for the mixed state configurations following a global quench is then given by half of that obtained for the eternal black hole geometry and exactly matches with the corresponding CF T 1+1 results. We then present a summary and our conclusions in the final Sec. (6).

Entanglement entropy
In this section we briefly review the time evolution of the entanglement entropy following a global quench as described in [2,36,37,40]. In a global quench scenario, the system is initially prepared in a pure state |ψ 0 which is the ground state of a Hamiltonian H 0 at time t = 0. For time t > 0 the same system is allowed to evolve unitarily with respect to another Hamiltonian H, where H and H 0 are related to each other by an experimentally tunable parameter. The unitary evolution of the density matrix ρ 0 = |ψ 0 ψ 0 | for t > 0 is governed by the post-quenched Hamiltonian H as where ρ(t) is the density matrix at a time t which is always in a pure state as a consequence of unitary evolution. The entanglement entropy of a subsystem A is obtained as the von Neumann entropy of the reduced density matrix ρ A (t) = Tr B ρ(t) where B = A c represents the rest of the system, which is given as In a CF T 1+1 for simplicity the initial state may be chosen as |ψ 0 = e −τ 0 H/4 |B where |B is a conformal invariant boundary state (B-state) and τ 0 is the correlation length in the initial state. The entanglement entropy S A (t) may be obtained through a replica technique as follows The quantity Trρ n A for a single interval A of length ℓ as depicted in Fig. 1(a) in the CF T 1+1 , is related to the two point twist correlator on a strip of width 2τ 0 as [2,36,37,40] where w i = u i + iτ are the complex coordinates on the strip (u i ∈ R and 0 < τ < 2τ 0 ). The above twist correlator may be obtained through a conformal map z i = exp(πw i /2τ 0 ) from the strip to the upper half plane (UHP) which is then given as where η i,j are the cross ratios, and the scaling dimension ∆ n of the twist operators T n andT n is given as follows The function F n ({η j,k } in eq. (2.5) is non-universal and depends on the full operator content of the theory. Note that only the limiting behavior of η i,j → 0 and η i,j → 1 are important and in these limits the non-universal function F n ({η j,k } is just a constant [40]. The entanglement entropy S A may now be obtained by utilizing eqs. (2.3) and (2.5) and considering the spacetime scaling limit (t, ℓ ≫ τ 0 ), which is given as follows (see [40] for detailed computation) where the function q(t, ℓ) is defined as follows Note that the entanglement entropy S A increases linearly with time for t < ℓ/2 and saturates to a thermal value beyond the time t ∼ ℓ/2. It may be observed that the entanglement entropy for large times is same as the thermal entropy of a mixed state at a large but finite temperature T = 1 β = 1 4τ 0 . This saturation to a thermal entropy at large times leads to the fascinating phenomena of CF T thermalization as described in [40,44].

Entanglement negativity
In this section we briefly recapitulate the evolution of the entanglement negativity for mixed states of two disjoint and adjacent intervals in a CF T 1+1 after a global quench scenario as described in [40]. As mentioned in the introduction, for a bipartite system A = A 1 ∪ A 2 in a mixed state, the entanglement between the subsystems A 1 and A 2 is characterized by the entanglement negativity, which is defined as [3] Here the partial transposed reduced density matrix ρ T 2 A is obtained as follows where |e (1) i and |e (2) j are the bases for the the Hilbert spaces H A 1 and H A 2 of the subsystems A 1 and A 2 respectively. In a CF T 1+1 the entanglement negativity is obtained through a suitable replica technique through the quantity Tr ρ T 2 A n for even n = n e and the analytic continuation to n e → 1 as follows [5][6][7]

Two disjoint intervals
Now consider a bipartite system A = A 1 ∪ A 2 of two disjoint intervals A 1 and A 2 in a CF T 1+1 of lengths ℓ 1 = |u 2 − u 1 | and ℓ 2 = |u 4 − u 3 | separated by a distance ℓ s = |u 3 − u 2 | as depicted in Fig. 1(b). In this case the quantity Tr(ρ T 2 A ) n is given by a four-point twist correlator on a strip of width 2τ 0 as [40] Tr where ∆ = 4∆ n and z i = e π(w i /(2τ 0 ) with w i = u i + iτ are the complex coordinates on the strip. The scaling dimension ∆ (2) n of the twist operators T 2 n andT 2 n in eq. (3.4) is given as where ∆ n is defined in (2.6).
The time evolution of the entanglement negativity E may then be obtained by utilizing the eqs. (3.3) and (3.4) considering the spacetime scaling limit (t, ℓ ≫ τ 0 ) as [40] where the function q(t, u i − u j ) with (i, j ∈ 1, 2, 3, 4) is defined in the eq. (2.8). Note that the function q(t, u i − u j ) depends non trivially on the relative values of the time t and the lengths of the intervals ℓ 1 , ℓ 2 and ℓ s leading to various interesting limits for the entanglement negativity in this case. We have computed the entanglement negativity in this case for each of these significant limits as which traces the time evolution of this quantity as described below.
In this limit of early time none of the functions q(t, u i − u j )s defined in eq. (2.8) contribute in eq. (3.6) and consequently the entanglement negativity vanishes For this limit only the function q(t, u 3 − u 2 ) contributes and the other q(t, u i − u j )s vanish in eq. (3.6), leading to a non-zero time dependent value of the entanglement negativity which is given as follows In this case the functions q(t, u 3 − u 1 ) and q(t, u 3 − u 2 ) are non-zero in eq. (3.6), and the entanglement negativity saturates to a constant value proportional to the length of the interval A 1 given by In this limit the functions q(t, u 3 − u 1 ), q(t, u 4 − u 2 ) and q(t, u 3 − u 2 ) are non-zero in eq. (3.6), and the entanglement negativity may be expressed as For this limit of late times all of the functions q(t, u i − u j )s in the eq. (3.6) contribute but the entanglement negativity vanishes due to mutual cancellations The above results, obtained by us for the time evolution of the entanglement negativity following a global quench, may be described as follows. The entanglement negativity for early time described in eq. (3.7) is zero and increases linearly during the initial period of the time evolution as described in eq. (3.8). This saturates to a constant value described in eq. (3.9) for an intermediate range of time. As time increases further the entanglement negativity decreases linearly as described in eq. (3.10) and vanishes at late times describing the thermalization of the mixed state under consideration in the CF T 1+1 .

Two adjacent intervals
Having described the case for the mixed state of disjoint intervals above we now proceed to analyse the case for the mixed state of adjacent intervals A 1 and A 2 . For this case the quantity Tr(ρ T 2 A ) n is given by a three-point twist correlator on a strip as follows [40] Tr where ∆ = 2∆ n + ∆ n . The entanglement negativity E may then be obtained from eqs. (3.12) and (3.3) by employing the spacetime scaling limit (t, ℓ ≫ τ 0 ) as [40] (3.13) In the above expression as earlier the function q(u i , u j ) in eq. (2.8) non trivially depends on the relative values of the time t and the lengths of the intervals ℓ 1 and ℓ 2 leading to certain significant limits. In what follows, as earlier we compute the entanglement negativity of the mixed state in question for these interesting limits which describes the time evolution of this measure.
For this early time the functions q(t, u i − u j )s in the eq. (3.13) vanish and the entanglement negativity may then be given as E = πc 8τ 0 t , (3.14) 3.2.2 ℓ 1 /2 < t < ℓ 2 /2 < (ℓ 1 + ℓ 2 )/2 In this case only the function q(t, u 2 − u 1 ) is non-zero in the eq. (3.13) and the entanglement negativity is proportional to the length of the subsystem A 1 as 3.2.3 ℓ 1 /2 < ℓ 2 /2 < t < (ℓ 1 + ℓ 2 )/2 In this case the functions q(t, u 2 − u 1 ) and q(t, u 3 − u 2 ) in the eq. (3.13) are non-zero leading to a non-zero time dependent value of the entanglement negativity given as follows In the limit of late time all of the functions q(t, u i − u j )s in the eq. (3.13) are non-zero which leads to the vanishing value of the entanglement negativity The evolution of the entanglement negativity from the above results is described as follows. It is observed that at the initial part of of its time evolution the entanglement negativity increases linearly with time described in eq. (3.14) and saturates to a constant value for an intermediate time as given in eq. (3.15). As time progress further negativity decreases linearly as described in eq. (3.16) and vanishes at late times indicating the phenomena of thermalization of the mixed state in consideration.

Holographic entanglement entropy
In this section we briefly review the holographic computation of the entanglement entropy for a single interval in a CF T 1+1 following a global quench in the AdS 3 /CF T 2 framework. As described in [61] the correct holographic dual of the global quench scenario is given by the single sided eternal black hole obtained by slicing the Penrose diagram of the usual eternal black hole in half with an end of the world (ETW) brane. In this instance the minimal area extremal surfaces (geodesics for AdS 3 /CF T 2 ) terminate on the ETW brane and the corresponding holographic entanglement entropy is just half of that obtained for the usual eternal black hole geometry.
In the above context we begin by reviewing the time evolution of the holographic entanglement entropy from a bulk eternal BTZ black hole geometry described by a two-sided Penrose diagram as depicted in Fig. 2. As described in [55,62] this bulk geometry is dual to the two copies of the CF T (CF T R ⊗ CF T L ), in the thermofield double state 5 . To obtain the holographic entanglement entropy of a single interval A in this scenario, it is required to consider two copies of A in the two CF T s defined on either side of the Penrose diagram as shown in Fig. 3(b). It is well known that the bulk eternal BTZ black hole is a quotient space of AdS 3 and hence the angular coordinate may be unwrapped to map the the two boundaries of the BTZ Penrose diagram to the corresponding Rindler wedges [55,62] as described in Fig. 3(a). The bulk geometry in this case is then an eternal BTZ black string.
It is now useful to describe the different regions of the BTZ black string and relate each coordinate patch to the Poincaré coordinates, The exterior metric of the BTZ black string is given as  where the horizon is located at ρ = 0, and the boundary is at ρ → ∞. The BTZ coordinates (t, x, ρ) covers a part of the Poincaré patch, close to the boundary and are related to the Poincaré coordinates as follows

2)
The coordinates (t, x) cover the Rindler wedge denoted as the region I of the Minkowski diamond in Fig. 3(a). It is possible to reach the other Rindler wedge denoted by the region III in Fig.  3(a) by analytically continuing the time t → t + iπ. This is the continuation through which it is possible to reach from one side of the exterior region of the black string Penrose diagram to the other. The interior metric to the horizon ρ = 0 is obtained by analytically continuing the coordinates ρ = iα, t =t − iπ/2 in eq. (4.2) as follows ds 2 = sin 2 αdt 2 + cos 2 αdx 2 − dα 2 . These coordinates cover some part of the interior region and meet the boundary along the light cone x 2 0 − x 2 1 = 0 as depicted in Fig. 3(a). The corresponding future bulk region of the interior to the BTZ black string is obtained by setting α = π/2 − iρ, x =x − iπ/2 in eq. (4.4), which is given as ds 2 = cosh 2ρ dt 2 − sinh 2ρ dx 2 + dρ 2 . The boundary of this future region in the Poincaré coordinates is given as z ≤ x 2 0 − x 2 1 . We now consider the interval A of length ℓ on both the Rindler wedges I and III of the Minkowski diamond in Fig. 3(b). These two intervals are described as (P 1 , P 3 ) and (P 2 , P 4 ) in the x 1 -z plane of the Poincaré coordinates as depicted in Fig. (4). The coordinates of the end points of the intervals at a time t are given from eq. (4.3) as follows [55] P 1 =(sinh t, − cosh t), P 3 = e −ℓ P 1 , At early time t ≤ ℓ/2, the entanglement entropy receives contribution form the lengths of the geodesics L 12 and L 34 between the points P 1 to P 3 and P 2 to P 4 respectively as depicted in Fig. 4. These geodesics pass through the interior region in eq. (4.4) of the BTZ black string connecting one Rindler wedge to the other (the dashed geodesics in Fig. 4). These two geodesics are semicircles and they have identical regularized lengths. The entanglement entropy at early time t ≤ ℓ/2 is then given as follows [55] S (1) (4.7) Here S div = c 6 ln β 4πǫ is the divergent part of the entanglement entropy, β = 4τ 0 is the inverse temperature and ǫ is the UV cut off. On the other hand for times t > ℓ/2 the entanglement entropy receives contribution from the lengths of the geodesics L 13 and L 42 between the points P 1 to P 3 and P 2 to P 4 respectively. These geodesics are confined to the exterior region of the BTZ black string outside the horizon as described in eq. (4.2) (depicted by the solid lines in the Fig. 4). The entanglement entropy for time t > ℓ/2 is then given as It is observed from the above development that at early time the geodesics pass through the interior region of the BTZ black brane geometry from one Rindler wedge I to the other III and the linear increase of the holographic entanglement entropy is related to the growth of the nice slice in the interior. On the other hand for late times the geodesics are confined to the exterior region sitting outside the horizon and the corresponding entanglement entropy becomes thermal.
As mentioned earlier the bulk dual geometry corresponding to a global quench in a CF T 1+1 where the initial state is the conformal invariant boundary state |B [36] is described by the eternal black hole sliced in half by an ETW brane (single sided black hole). Hence the corresponding holographic entanglement entropy for the interval A in this case is then described by half of that obtained in eqs. (4.7) and (4.8) respectively. Remarkably, these results match exactly with the corresponding CF T 1+1 computation through the replica technique in eq. (2.7).

Holographic entanglement negativity
In this section we present our holographic construction for the time evolution of the entanglement negativity of two disjoint and adjacent intervals in a CF T 1+1 following a global quench as described in [40]. Motivated by the holographic construction for the time evolution of the entanglement entropy in [61] reviewed above, we begin by computing the holographic entanglement negativity for the mixed states under consideration from the bulk dual eternal black hole geometry. The corresponding holographic entanglement negativity following a global quench may then be obtained by considering half of the corresponding result for the bulk dual eternal black hole geometry.

Two disjoint intervals
We first consider the case of two disjoint intervals A 1 and A 2 of lengths ℓ 1 and ℓ 2 with an interval A s of length ℓ s separating them. As earlier for the dual eternal black hole geometry it is required to consider the intervals A 1 , A 2 and A s in both the Rindler wedges I and III of the Minkowski diamond in Fig. 5(a). The end points of the intervals A 1 , A 2 and A s in the Rindler wedges I and III projected on the x 1 -z plane of the Poincaré coordinates are denoted as (a, b), (c, d), (b, c) and (a ′ , b ′ ), (c ′ , d ′ ), (b ′ , c ′ ) respectively which is depicted in Fig. 5(b). The holographic entanglement negativity for this configuration may then be obtained through the construction described in [34,35]. From this prescription the holographic entanglement negativity for two disjoint intervals in dual CF T 1+1 s is given as follows where S A i is the holographic entanglement entropy of a subsystem A i and I (A i , A j ) is the holographic mutual information 6 between A i and A j (i, j ∈ 1, 2, s). It is extremely interesting to note that substituting the expression for the entanglement entropy of a single interval eq. (2.7) in eq. (3.6) for the entanglement negativity of two disjoint intervals leads us to the identical expression described in eq. (5.1) as follows Note that the entanglement negativity in eq. (5.2) is obtained through a replica technique and considering a space time scaling limit whereas that in eq. (5.1) is purely holographic. This is an extremely strong consistency check for our construction and the associated geodesic structure to be described subsequently, which indicates that the holographic entanglement negativity obtained from our construction in [34,35] is indeed universal.
It is now possible to describe the time evolution of the holographic entanglement negativity for the mixed state configuration of disjoint intervals from the bulk dual eternal black hole geometry by utilizing eqs. (5.1), (4.7) and (4.8). It is observed that distinct sets of geodesics contribute to the holographic entanglement negativity for various values of time t relative to the lengths of the intervals ℓ 1 , ℓ 2 and ℓ s . In what follows we consider these significant limits and describe the relevant geodesic structures leading to the holographic entanglement negativity for each of these scenarios. 5.1.1 t < ℓ s /2 < ℓ 1 /2 < ℓ 2 /2 In the limit of early times all the relevant geodesics pass through the interior region of the BTZ black string connecting one Rindler wedge I to the other III, depicted by dashed curves in the Fig. 6(a). The corresponding holographic entanglement negativity is obtained by substituting the holographic entanglement entropy for the early time limit given in eq. (4.7) in the eq. (5.1), which vanishes due to mutual cancellations as follows In this limit the geodesics (solid curves) that joins c ′ and b ′ and c and b are confined to the exterior region, and rest of the geodesics (dashed curves) pass through the interior region of the BTZ black string as depicted in Fig. 6(b). The holographic entanglement negativity in this case may then be obtained by utilizing both the eqs. (4.7) and (4.8) in the eq. (5.1), which leads to the following expression 5.1.3 ℓ 2 /2 < (ℓ 1 + ℓ s )/2 < t < (ℓ 2 + ℓ s )/2 In this case the geodesics (solid curves) between the points d ′ to b ′ , d to b, c ′ to b ′ and c to b are confined in the exterior region while the other geodesics (dashed curves) between the points d ′ to d and c ′ to c pass through the interior region of the BTZ black string as depicted in Fig.  6(c). The holographic entanglement negativity in this limit is obtained by utilizing eqs. (4.7) and (4.8) in the eq. (5.1) leading to the following expression In this limit the geodesics (dashed curves) between the points d ′ to d and a ′ to a pass through the interior region whereas rest of the geodesics (solid curves) are confined to the exterior region of the BTZ black string. The holographic entanglement negativity in this case may be obtained by utilizing eqs. (4.7) and (4.8) in the eq. (5.1), as follows In the late time limit all the geodesics (solid curves) are confined to the exterior region of the BTZ black string as shown in Fig. 6(e). Consequently, the holographic entanglement negativity is obtained by utilizing eqs. (4.7) and (4.8) in the eq. (5.1) which vanishes due to mutual cancellations as follows =0. (5.7) It is now possible to present a holographic description for the time evolution of the entanglement negativity for mixed state configurations of two disjoint intervals in a CF T 1+1 following a global quench. As described earlier the corresponding bulk dual geometry is the single sided black hole obtained by slicing the eternal black hole geometry in half with an ETW brane. So the time evolution of the entanglement negativity for the mixed state in question following a global quench may be obtained by considering half of the corresponding results for the eternal black hole geometry computed above. Quite significantly, the time evolution of the holographic entanglement negativity obtained from our construction in this case matches exactly with the corresponding replica technique results for a CF T 1+1 as described in section (3.1). Once again this is a robust consistency check for our holographic construction.
In Fig. 7 we have plotted the time evolution of the holographic entanglement negativity for the mixed state in question in a CF T 1+1 following a global quench. It is observed from Fig. 7 that the holographic entanglement negativity for early times is zero and grows linearly to attain a maximum value at an intermediate time after which remains constant upto certain value of time. As time increases further it decreases linearly and vanishes for late times. Note that this behavior for the entanglement negativity of the mixed state in question consistently matches with that described for specific CF T 1+1 s in [40].

Two adjacent intervals
We now consider the time evolution of the entanglement negativity for the case of two adjacent intervals A 1 and A 2 of lengths ℓ 1 and ℓ 2 respectively. Proceeding similarly as the case for two disjoint intervals described above we consider the two adjacent intervals A 1 and A 2 in both the Rindler wedges I and III of the Minkowski diamond as shown in Fig. 8(a). The end points of the intervals A 1 and A 2 in the Rindler wedges projected on the x 1 − z plane of the Poincaré coordinates are denoted as (a, b), (b, c) and (a ′ , b ′ ), (b ′ , c ′ ) respectively which is depicted in Fig. 8(b). The holographic entanglement negativity for this configuration may then be obtained through the construction described in [30,31]. From this prescription the holographic entanglement negativity for the mixed state of two adjacent intervals in a dual CF T 1+1 is given as follows It is very significant that in this case also, substituting the expression for the entanglement entropy of a single interval eq. (2.7) in eq. (3.13) for the entanglement negativity of two adjacent intervals leads us to the identical expression described in eq. (5.8) as follows Once again the entanglement negativity in eq. (5.9) is obtained through a replica technique considering the space time scaling limit whereas the eq. (5.8) is purely holographic which constitutes a strong consistency check for our construction.
Following the analysis described earlier in the last subsection for the mixed state of disjoint intervals it is also possible to obtain the time evolution of the entanglement negativity for the case of two adjacent intervals under consideration here. Once again for this purpose it is required to consider the dual bulk eternal black hole geometry and utilize eqs. (5.8), (4.7) and (4.8) to obtain the corresponding holographic entanglement negativity. In this case also distinct sets of geodesics contribute to the holographic entanglement negativity for various values of the time t relative to the lengths of the intervals ℓ 1 , ℓ 2 and ℓ s . As earlier we consider these interesting limits and elucidate the corresponding geodesic structures leading to the holographic entanglement negativity for each of these cases.
5.2.1 t < ℓ 1 /2 < ℓ 2 /2 < (ℓ 1 + ℓ 2 )/2 In the limit of early times all of the geodesics pass through the interior region of the BTZ black string connecting the one Rindler wedge I to the other III, depicted by dashed curves in the Fig. 9(a). The corresponding holographic entanglement negativity is obtained by employing eq. (4.7) in the eq. (5.8), which is given as follows (5.10) 5.2.2 ℓ 1 /2 < t < ℓ 2 /2 < (ℓ 1 + ℓ 2 )/2 In this limit the geodesics (solid curves) between the points c ′ to b ′ and c to b are confined in the exterior region, rest of the geodesics pass through the interior region of the BTZ black string connecting the Rindler wedge I to the other III as shown in Fig. 9(b). The holographic entanglement negativity in this limits is given as In this case only the geodesics (dashed curves) between the points c ′ to c and c to b pass through the interior region and rest of the geodesics remain in the exterior region of the BTZ black string as depicted in Fig. 9(c). The holographic entanglement negativity in this case is given as follows

ℓ
In this limit all the geodesics (solid curves) are confined in the exterior region of the BTZ black string as depicted in Fig. 9(d). The resulting holographic entanglement negativity vanishes due to mutual cancellations as follows As mentioned earlier the time evolution holographic entanglement negativity of two adjacent intervals after a global quench in a CF T 1+1 dual to a single sided black hole is obtained by taking half of the corresponding results for the eternal black hole case. Interestingly, it is observed that the holographic entanglement negativity in question following a global quench match exactly with the corresponding CF T 1+1 results computed through the replica technique (3.2).
The time evolution of the holographic entanglement negativity is shown in Fig. 10. The holographic entanglement negativity grows linearly for the initial part of its evolution and reaches to its maximum value at an intermediate time after which it becomes a constant upto certain range of time. As time increases further it decreases linearly and becomes zero for late times. Note that this behavior of the holographic entanglement negativity is consistence with that described in [40].

Summary and discussion
To summarize, in this article we have investigated the time evolution of the entanglement negativity for mixed states of two disjoint and adjacent intervals following a global quench in a holographic CF T 1+1 through the AdS 3 /CF T 2 correspondence. In this context we have considered the holographic dual to the global quench scenario as the two sided eternal black hole sliced in half by an end of the world (ETW) brane, where the initial state is the conformal invariant boundary state (B-state). The corresponding holographic entanglement negativity is given by half of the result for the eternal black hole. To this end we obtained the holographic entanglement negativity for the mixed states in question by considering the bulk dual eternal black hole geometry. In this case these intervals were considered on both the Rindler wedges of the Minkowski diamond and projected on the x 1 − z plane of the Poincaré coordinates.
To describe the time evolution of the entanglement negativity for the above mixed state configurations in the CF T dual to the bulk eternal black hole geometry, we utilize a covariant holographic construction proposed earlier in [30,31,34,35]. This construction involves a specific algebraic sum of the lengths of bulk geodesics anchored on appropriate intervals. For our case the bulk configuration is described by the BTZ black string and it is required to consider geodesics located exterior to the horizon and those passing through the interior region connecting the two Rindler wedges of the Minkowski diamond. Interestingly it is observed that the holographic entanglement negativity at early times receive contributions from the lengths of the geodesics passing through the interior region of the bulk BTZ black string. For late times on the other hand the contribution to the holographic entanglement negativity arise from the geodesics restricted to the region exterior to the horizon in both the Rindler wedges.
The corresponding time evolution of the entanglement negativity following a global quench for the disjoint and the adjacent intervals, obtained from half the result for that of the bulk eternal black hole, exhibits a similar behavior. For early times the entanglement negativity is zero and increases linearly to reach a saturation at an intermediate value and remains constant upto a certain range of time. Beyond this range the entanglement negativity decreases linearly with time and vanishes for late times which corresponds to phenomena of the CF T thermalization indicating the formation of a black hole in the bulk. Interestingly, the holographic entanglement negativity obtained by us match exactly with the corresponding CF T 1+1 replica technique results providing a robust consistency check for our construction. Our analysis described in this article provides interesting insight into the structure and evolution of entanglement negativity in conformal field theories following a global quench from a holographic perspective and indicates a possible higher dimensional extension to a generic AdS d+1 /CF T d scenario although that will require certain substantiation and relevant consistency checks. We hope to return to these and other connected issues in the near future.