This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

{NoHyper}

Holographic entanglement negativity for a single subsystem in conformal field theories with a conserved charge
 

Sayid Mondal E-mail:  [email protected] Center for High Energy Physics and Department of Physics Chung-Yuan Christian University, Chung-Li 320, Taiwan Boudhayan Paul E-mail:  [email protected] Department of Physics Indian Institute of Technology Kanpur 208016, India Gautam Sengupta E-mail:  [email protected] Department of Physics Indian Institute of Technology Kanpur 208016, India Punit Sharma E-mail:  [email protected] Department of Physics and Astronomy The University of Iowa Iowa City, IA 52242, USA
Abstract

We utilize a holographic construction to compute the entanglement negativity for bipartite mixed state configurations of a single subsystem in CFTdCFT_{d}s with a conserved charge dual to bulk AdSd+1AdS_{d+1} geometries. In this context, we obtain the holographic entanglement negativity for single subsystems with long rectangular strip geometries in CFTdCFT_{d}s with a conserved charge dual to bulk non extremal and extremal Reissner-Nordström (RN)-AdSd+1AdS_{d+1} black holes at the leading order in a perturbation theory. We demonstrate that the holographic entanglement negativity computed involves the elimination of thermal contributions at the leading order confirming earlier results in the literature. This also conforms to quantum information theory expectations and constitutes further consistency checks for the holographic construction.

 


 


Introduction

The recent past has witnessed the emergence of quantum entanglement as a central theme in the study of many body systems relevant to diverse fields from condensed matter physics to quantum gravity and black holes. For bipartite pure states the entanglement entropy serves as an appropriate entanglement measure. Remarkably for (1+1)(1+1)-dimensional conformal field theories (CFT1+1CFT_{1+1}s) an explicit computation for the entanglement entropy utilizing the replica technique was established in [1, 2, 3, 4, 5]. For bipartite mixed states however the entanglement entropy incorporates irrelevant correlations and is hence unsuitable for characterization of the entanglement. In quantum information theory, several measures have been proposed for the characterization of mixed state entanglement which are in general not easily computable. Vidal and Werner in a significant work [6] proposed a computable measure termed entanglement negativity (logarithmic negativity) for bipartite mixed state entanglement which characterizes an upper bound on the distillable entanglement for the mixed state.111Distillable entanglement is the amount of pure entanglement that can be extracted from a given state by means of local operations and classical communication (LOCC). This measure was defined as the logarithm of the trace norm for the partially transposed reduced density matrix with respect to one of the subsystems and shown to be an entanglement monotone [7].222For an excellent review see [8]. Interestingly in a series of articles [9, 10, 11] the entanglement negativity for various bipartite pure and mixed states in a CFT1+1CFT_{1+1} could be obtained through an appropriate replica technique.

In a different context a holographic characterization of the entanglement entropy for dd-dimensional conformal field theories (CFTdCFT_{d}s) dual to bulk AdSd+1AdS_{d+1} geometries was conjectured in a seminal work by Ryu and Takayanagi (RT) [12, 13]. Their proposal described the holographic entanglement entropy of a subsystem in the dual CFTdCFT_{d} in terms of the area of a bulk codimension two static minimal surface homologous to the subsystem (RT Surface). The RT conjecture and its covariant generalization by Hubeny, Rangamani and Takayanagi (HRT) [14] were proved in a series of subsequent works [15, 16, 17, 18, 19, 20]. The RT/HRT conjectures inspired the remarkable development of the field of holographic quantum entanglement in dual CFTdCFT_{d}s summarized in [21, 22, 23, 24, 25, 26] and references therein.333Note that the field is extremely diverse and the list of references is not meant to be exhaustive.

The above developments naturally brought the critical issue of a holographic characterization for the entanglement negativity in dual CFTdCFT_{d}s into sharp focus. This significant issue was initially partially addressed for pure states in [8] although a clear holographic prescription for generic bipartite mixed states in dual CFTdCFT_{d}s remained an outstanding issue. In an interesting communication the authors in [27] proposed a holographic entanglement negativity conjecture for the bipartite states of a single interval in CFT1+1CFT_{1+1}s dual to bulk pure AdS3AdS_{3} geometries and Euclidean BTZ black holes. This proposal involved a certain algebraic sum of the lengths of bulk geodesics homologous to appropriate combinations of the single interval and other auxiliary intervals. Remarkably the results for the holographic entanglement negativity computed from their proposal were in exact agreement with the corresponding CFT1+1CFT_{1+1} replica technique results in the large central charge limit described in [9, 10, 11]. This construction was further substantiated through a rigorous large central charge analysis utilizing the monodromy technique in [28] for the AdS3/CFT2AdS_{3}/CFT_{2} scenario. Subsequently in a series of articles this holographic construction and its covariant generalization were utilized to obtain the entanglement negativity for diverse mixed state configurations in the AdS3/CFT2AdS_{3}/CFT_{2} scenario [29, 27, 30, 31, 32, 33].

A higher dimensional generalization of the above holographic conjecture for generic AdSd+1/CFTdAdS_{d+1}/CFT_{d} scenarios involving specific algebraic sums of the areas of the RT surfaces for certain combinations of subsystems, was subsequently proposed and utilized to compute the entanglement negativity of bipartite states of subsystems with long rectangular strip geometries in dual CFTdCFT_{d}s [34, 35, 36, 37]. The results conformed to quantum information theory expectations and reproduced certain significant features for the entanglement negativity also observed for the AdS3/CFT2AdS_{3}/CFT_{2} scenario which served as crucial consistency checks. Very recently a plausible bulk proof of these holographic entanglement negativity conjectures from a gravitational path integral perspective, along the lines of [18, 19] for spherical subsystem geometries based on [38] was described in [39].444A review of this plausible proof for single subsystems with spherical geometries has been reviewed in appendix A. However a generalization of this for generic subsystem geometries is still a non trivial outstanding issue.

We should mention here that subsequent to the developments described above, the authors in [40, 41] proposed an alternate holographic entanglement negativity conjecture for bipartite states in dual CFTdCFT_{d}s involving the minimal entanglement wedge cross section (EWCS) [42] backreacted by the cosmic brane for the conical defect of the bulk replicated geometry. For spherical entangling surfaces this backreaction could be computed and resulted in a dimension dependent constant numerical factor for the entanglement negativity. For the AdS3/CFT2AdS_{3}/CFT_{2} scenario the holographic results from their proposal matched with the corresponding field theory replica technique results in the large central charge limit modulo certain constants arising from the Markov gap [43]. The above holographic construction was further refined in [44] to resolve an outstanding issue and demonstrated the equivalence of the two proposals modulo the constants mentioned above.555Note that the constant backreaction factor χd=32\chi_{d}=\frac{3}{2} for the AdS3/CFT2AdS_{3}/CFT_{2} scenario was incorporated in the proposal based on the algebraic sum of the lengths of bulk geodesics homologous to appropriate combinations of intervals described in [27]. However this alternative proposal requires further substantiation for the higher dimensional AdSd+1/CFTdAdS_{d+1}/CFT_{d} scenario.666See also the discussions in [45].

In this article we extend the higher dimensional holographic construction to compute the entanglement negativity for bipartite states of single subsystems described in [27] to CFTdCFT_{d}s with a conserved charge dual to bulk non extremal and extremal AdSd+1AdS_{d+1}-RN black holes. Our results provide further substantiation for the holographic conjecture for the entanglement negativity of bipartite states in higher dimensional CFTdCFT_{d}s. In this context following [27] we also consider subsystems with long rectangular strip geometries and utilize a perturbative technique described in [24, 23, 46, 47], to compute the areas of RT surfaces involving non trivial limits of certain parameters relevant for the bulk AdSd+1AdS_{d+1}-RN geometry.777 Note that an exact closed form analytic evaluation of the area has been performed through the use of Meijer G-functions in [48]. However the leading order behavior necessary for our purpose is correctly captured through the perturbative techniques utilized here. Interestingly as described in the literature, through our construction we reproduce certain significant features (such as elimination of the thermal contributions) of the holographic entanglement negativity for bipartite states in CFT2CFT_{2}s and CFTdCFT_{d}s described in [27, 34]. This conforms to quantum information theory expectations of the negativity being an upper bound on the distillable entanglement.

Note that in the context of the developments described in [40, 41, 44] our results also involve backreaction effects due to the cosmic brane arising from the conical defect of the replicated bulk geometry in the replica limit. However the explicit determination of these effects remains a non trivial outstanding issue for arbitrary subsystem geometries, in particular the long rectangular strip geometries investigated in this article.

In this connection we state here that in [49] the authors provided another construction to characterize the modification of the RT surfaces due to the presence of a charged bulk geometry. In this article they considered a charged topological black hole with a hyperbolic horizon coupled to a dual Maxwell gauge field to compute charged holographic Rényi entropies and the holographic entanglement entropy for spherical entangling surfaces in the dual CFTdCFT_{d}. Interestingly the corresponding charged Rényi entropies and the entanglement entropy for such spherical domains in the CFTdCFT_{d} could be computed using a generalized replica technique for a thermal path integral via a background gauge field with a non trivial Wilson line on the Euclidean time circle following an appropriate conformal transformation.

This article is organized as follows. In section 2 we describe the holographic entanglement negativity construction for the bipartite configuration of a single subsystem characterized by long rectangular strip geometry in CFTdCFT_{d}s dual to bulk AdSd+1AdS_{d+1} geometries. Subsequently in section 3 the construction is employed to compute the holographic entanglement negativity for mixed state configuration of a single subsystem in CFTdCFT_{d}s with a conserved charge dual to bulk non extremal and extremal AdSd+1AdS_{d+1}-RN black holes. Finally in section 4 we summarize our results and present our conclusions. We also provide some involved and lengthy results in the appendices.

Holographic entanglement negativity

In this section we present a brief review of the holographic entanglement negativity conjecture for the bipartite state of a single subsystem in a CFTdCFT_{d} at a finite temperature dual to bulk AdSd+1AdS_{d+1} geometries as described in [9, 10, 11, 34, 27]. In this context we begin with recapitulating the corresponding issue in the AdS3/CFT2AdS_{3}/CFT_{2} scenario [27]. For this purpose we consider a bipartite system (AAcA\cup A^{c}) which is described by an interval AA of length ll and its complement AcA^{c} in the dual CFT1+1CFT_{1+1}. To describe this construction it is also required to consider two auxiliary large but finite intervals B1B_{1} and B2B_{2} each of length bb on either side of AA and adjacent to it such that B=(B1B2)AcB=(B_{1}\cup B_{2})\subset A^{c}. The intervals may be represented as B1=(bl/2,l/2),A=(l/2,l/2),B2=(l/2,b+l/2)B_{1}=(-b-l/2,-l/2),A=(-l/2,l/2),B_{2}=(l/2,b+l/2). The entanglement negativity for the mixed state of the single interval at a finite temperature in the dual CFT1+1CFT_{1+1} may then be defined through a replica technique as follows [9]

=limBAclimne1lnTr(ρABTB)ne,\mathcal{E}=\lim_{B\to A^{c}}\lim_{n_{e}\to 1}\ln\mathrm{Tr}\left(\rho_{AB}^{T_{B}}\right)^{n_{e}}, (1)

where ρABTB\rho_{AB}^{T_{B}} is the partial transpose of the reduced density matrix ρAB\rho_{AB} with respect to the subsystem BB, and the replica limit ne1n_{e}\to 1 is an analytic continuation of even sequences of nen_{e} to ne=1n_{e}=1. For two arbitrary subsystems XX and YY, the partial transpose ρXYTY\rho_{XY}^{T_{Y}} of the reduced density matrix ρXY\rho_{XY}, with respect to the subsystem YY, may be defined as

ei(X)ej(Y)|ρXYTY|ek(X)el(Y)=ei(X)el(Y)|ρXY|ek(X)ej(Y),\matrixelement{e^{(X)}_{i}e^{(Y)}_{j}}{\rho_{XY}^{T_{Y}}}{e^{(X)}_{k}e^{(Y)}_{l}}=\matrixelement{e^{(X)}_{i}e^{(Y)}_{l}}{\rho_{XY}}{e^{(X)}_{k}e^{(Y)}_{j}}, (2)

where |ei(X){|e^{(X)}_{i}\rangle} and |ej(Y)|e^{(Y)}_{j}\rangle describe orthonormal bases for the Hilbert spaces X\mathcal{H}_{X} and Y\mathcal{H}_{Y} respectively. Note that the quantity Tr(ρABTB)ne\mathrm{Tr}(\rho_{AB}^{T_{B}})^{n_{e}} may be expressed as a twist correlator in the CFT1+1CFT_{1+1} appropriate to the bipartite configuration. For the case of the bipartite state of a single interval in the dual CFT1+1CFT_{1+1} this is given by a four point twist field correlator on an infinitely long cylinder and the corresponding entanglement negativity is given as follows [11]

=liml/b0limne1ln[𝒯ne(bl/2)𝒯¯ne2(l/2)𝒯ne2(l/2)𝒯¯ne(b+l/2)cyl],\mathcal{E}=\lim_{l/b\to 0}~{}\lim_{n_{e}\to 1}\ln\left[\left\langle\mathcal{T}_{n_{e}}(-b-l/2)\overline{\mathcal{T}}^{2}_{n_{e}}(-l/2)\mathcal{T}^{2}_{n_{e}}(l/2)\overline{\mathcal{T}}_{n_{e}}(b+l/2)\right\rangle_{cyl}\right], (3)

where l/b0l/b\to 0 is the bipartite limit B=(B1B2)AcB=(B_{1}\cup B_{2})\to A^{c}. It is important to note that the order of the limits in eq. (3) is crucial where the bipartite limit (l/b0l/b\to 0) should be taken only after the replica limit (ne1)(n_{e}\to 1).

The twist fields 𝒯\mathcal{T}s, 𝒯¯\overline{\mathcal{T}}s above provide the boundary conditions between the nen_{e} replicas and reduce to primary fields on the complex plane in the replica limit ne1n_{e}\to 1. As described in [28] the universal part of the four point twist field correlator described above may be obtained in the large central charge limit through a rigorous monodromy analysis as follows

𝒯ne(z1)𝒯¯ne2(z2)𝒯ne2(z3)𝒯¯ne(z4)1xΔne(2)z232Δne(2),\left\langle\mathcal{T}_{n_{e}}(z_{1})\overline{\mathcal{T}}^{2}_{n_{e}}(z_{2})\mathcal{T}^{2}_{n_{e}}(z_{3})\overline{\mathcal{T}}_{n_{e}}(z_{4})\right\rangle_{\mathbb{C}}\approx\frac{1}{x^{\Delta^{(2)}_{n_{e}}}z_{23}^{2\Delta^{(2)}_{n_{e}}}}\,, (4)

where the cross ratio xz12z34z13z24x\equiv\frac{z_{12}z_{34}}{z_{13}z_{24}} and zij=zizjz_{ij}=z_{i}-z_{j} and the scaling dimensions Δne\Delta_{n_{e}} and Δne(2)\Delta_{n_{e}}^{(2)} of the operator 𝒯ne\mathcal{T}_{n_{e}} and 𝒯ne2\mathcal{T}_{n_{e}}^{2} are related to each other as follows

Δne=c12(ne1ne),Δne(2)=2Δne2=c6(ne22ne).\Delta_{n_{e}}=\frac{c}{12}\left(n_{e}-\frac{1}{n_{e}}\right),\qquad\qquad\Delta_{n_{e}}^{(2)}=2\Delta_{\frac{n_{e}}{2}}=\frac{c}{6}\left(\frac{n_{e}}{2}-\frac{2}{n_{e}}\right). (5)

From the two point correlator of the twist fields in a CFT1+1CFT_{1+1}, given as

𝒯ne(zi)𝒯¯ne(zj)|zij|2Δ𝒯ne,\left\langle\mathcal{T}_{n_{e}}(z_{i})\overline{\mathcal{T}}_{n_{e}}(z_{j})\right\rangle_{\mathbb{C}}~{}\sim~{}\left\lvert z_{ij}\right\rvert^{-2\Delta_{\mathcal{T}_{n_{e}}}}, (6)

it is observed that the large central charge form of the four point twist correlator in eq. (4) may be expressed in terms of the products of certain two point twist field correlators. Using the AdS3/CFT2AdS_{3}/CFT_{2} dictionary these two point twist correlators may in turn be expressed in terms of the lengths of space like geodesics in the bulk AdS3AdS_{3} geometry (with length scale RR) homologous to certain combinations of the intervals as follows

𝒯ne(zi)𝒯¯ne(zj)eΔneijR,\left\langle\mathcal{T}_{n_{e}}(z_{i})\overline{\mathcal{T}}_{n_{e}}(z_{j})\right\rangle_{\mathbb{C}}~{}\sim~{}e^{-\frac{\Delta_{n_{e}}\mathcal{L}_{ij}}{R}}\,,
𝒯ne2(zk)𝒯¯ne2(zl)eΔne2klR.\left\langle\mathcal{T}_{\frac{n_{e}}{2}}(z_{k})\overline{\mathcal{T}}_{\frac{n_{e}}{2}}(z_{l})\right\rangle_{\mathbb{C}}~{}\sim~{}e^{-\frac{\Delta_{\frac{n_{e}}{2}}\mathcal{L}_{kl}}{R}}\,. (7)

In the above expressions ij\mathcal{L}_{ij} is the length of the geodesic homologous to the interval described by the endpoints ziz_{i} and zjz_{j}. From eqs. (3), (4) and (7), the holographic entanglement negativity may then be expressed as

=limBAc316GN(3)(2A+B1+B2AB1AB2),\mathcal{E}=\lim_{B\to A^{c}}~{}\left.\frac{3}{16G_{N}^{(3)}}\middle(2\mathcal{L}_{A}+\mathcal{L}_{B_{1}}+\mathcal{L}_{B_{2}}-\mathcal{L}_{A\cup B_{1}}-\mathcal{L}_{A\cup B_{2}}\right), (8)

where we have utilized the Brown-Henneaux formula c=3R2GN(3)c=\frac{3R}{2G_{N}^{(3)}} [50]. As discussed in the Introduction, the expression for the holographic entanglement negativity above contains an overall constant numerical factor due to the backreaction of the cosmic brane for the bulk conical defect in the replica limit.888For the AdS3/CFT2AdS_{3}/CFT_{2} scenario where the entangling surfaces are just points this overall dimension dependent backreaction factor χd\chi_{d} may be explicitly determined as χ2=32\chi_{2}=\frac{3}{2} as described in [51, 40]. Note that using the RT prescription for the AdS3/CFT2AdS_{3}/CFT_{2} scenario the above equation may be expressed as follows

=limBAc34(2SA+SB1+SB2SAB1SAB2).\mathcal{E}=\lim_{B\to A^{c}}\left.\frac{3}{4}\middle(2S_{A}+S_{B_{1}}+S_{B_{2}}-S_{A\cup B_{1}}-S_{A\cup B_{2}}\right). (9)

Given the AdS3/CFT2AdS_{3}/CFT_{2} scenario discussed above, it was conjectured in [34] that in the context of the AdSd+1/CFTdAdS_{d+1}/CFT_{d} framework the corresponding holographic entanglement negativity for the mixed state of a single subsystem AA may be expressed in the bipartite limit (BAc)(B\to A^{c}) as

A=limBAc316GN(d+1)(2𝒜A+𝒜B1+𝒜B2𝒜AB1𝒜AB2),\mathcal{E}_{A}=\lim_{B\to A^{c}}\left.\frac{3}{16G_{N}^{(d+1)}}\middle(2\mathcal{A}_{A}+\mathcal{A}_{B_{1}}+\mathcal{A}_{B_{2}}-\mathcal{A}_{A\cup B_{1}}-\mathcal{A}_{A\cup B_{2}}\right), (10)

where 𝒜γ\mathcal{A}_{\gamma} denotes the area of the RT surface homologous to the subsystem γ\gamma, and GN(d+1)G_{N}^{(d+1)} is the gravitational constant in (d+1)(d+1) dimensions.999As discussed in the Introduction a plausible proof of the holographic conjectures expressed in eqs. (8) and (10) for spherical entangling surfaces discussed in [39] from replica symmetry breaking saddles for the bulk gravitational path integral [38] has been described in appendix A.

For the case when the subsystems B1B_{1} and B2B_{2} are symmetric with respect to the subsystem AA, eq. (10) simplifies to

A=limBAc38GN(d+1)(𝒜A+𝒜B1𝒜AB1).\mathcal{E}_{A}=\lim_{B\to A^{c}}\left.\frac{3}{8G_{N}^{(d+1)}}\middle(\mathcal{A}_{A}+\mathcal{A}_{B_{1}}-\mathcal{A}_{A\cup B_{1}}\right). (11)

Using the standard Ryu-Takayanagi (RT) formula given as

𝒮A=𝒜A4GN(d+1),\mathcal{S}_{A}=\frac{\mathcal{A}_{A}}{4G_{N}^{(d+1)}}\;, (12)

the above conjecture for the holographic entanglement negativity for the bipartite state of a single subsystem AA may be expressed as101010Note that this is a purely holographic result valid for the corresponding large cc limit of the dual CFTCFT. This is explicit for the AdS3/CFT2AdS_{3}/CFT_{2} scenario where exact computations are possible for both the bulk and the boundary. To check its validity for the full quantum information theory measure, it would be required to compute the non universal contributions to the CFT2CFT_{2} correlation functions, which is extremely non trivial.

A=limBAc32(SA+SB1SAB1).\mathcal{E}_{A}=\lim_{B\to A^{c}}\left.\frac{3}{2}\middle(S_{A}+S_{B_{1}}-S_{A\cup B_{1}}\right). (13)

Note that as discussed earlier for the AdS3/CFT2AdS_{3}/CFT_{2} scenario, the expressions in eqs. (9) and (10) for the AdSd+1/CFTdAdS_{d+1}/CFT_{d} framework will also involve effects due to the backreaction of the cosmic brane for the conical defect of the replicated bulk geometry. For arbitrary subsystem geometries the determination of the explicit form of these effects is still a non trivial open issue in particular for the long rectangular strip geometries considered in this article. However for spherical entangling surfaces the explicit form of the backreaction factor is described in [51, 40].

Holographic entanglement negativity for CFTdCFT_{d} dual to AdSd+𝟏AdS_{d+1}-RN

We now utilize the construction reviewed above to compute the holographic entanglement negativity for bipartite states of single subsystems in CFTdCFT_{d}s with a conserved charge, dual to bulk AdSd+1AdS_{d+1}-RN geometries. We employ a perturbative expansion of the areas of the corresponding bulk RT surfaces, in terms of the temperature TT and the chemical potential μ\mu conjugate to the charge QQ of the AdSd+1AdS_{d+1}-RN black hole. It is convenient to express the holographic entanglement negativity in terms of the total energy ε\varepsilon of the dual CFTdCFT_{d}, and an effective temperature TeffT_{\textit{ef{}f}}, both of which are functions of the temperature and the chemical potential. Note that the parameter ε\varepsilon is related to the expectation value of the component T00T_{00} of the energy momentum tensor as described in [46]. To this end we first review the perturbative calculation of the area of the RT surfaces required for the computation of the entanglement negativity in subsection 3.1 before proceeding to the non extremal bulk AdSd+1AdS_{d+1}-RN black holes in subsection 3.2 and finally describing the corresponding bulk extremal AdSd+1AdS_{d+1}-RN black holes in subsection 3.3.

Area of RT surfaces in AdSd+𝟏AdS_{d+1}-RN

The metric for the AdSd+1AdS_{d+1}-RN black hole (with d3d\geq 3) may be expressed as

ds2=1z2(f(z)dt2+dz2f(z)+dx2),ds^{2}=\frac{1}{z^{2}}\left(-f(z)dt^{2}+\frac{dz^{2}}{f(z)}+d\vec{x}^{2}\right), (14)

f(z)=1Mzd+(d2)Q2(d1)z2(d1),At(z)=Q(zHd2zd2),f(z)=1-Mz^{d}+\frac{(d-2)Q^{2}}{(d-1)}z^{2(d-1)}\;,\qquad\qquad A_{t}(z)=Q\left(z_{H}^{d-2}-z^{d-2}\right),

where the AdSAdS length scale has been set to unity, and QQ and MM are the mass and charge of the black hole respectively, and ff is the lapse function. The horizon zHz_{H} is determined by the smallest real root of the equation f(z)=0f(z)=0. The chemical potential μ\mu, conjugate to the charge QQ, is then defined as follows

μlimz0At(z)=QzHd2,\mu\equiv\lim_{z\to 0}A_{t}(z)=Qz_{H}^{d-2}, (15)

while the Hawking temperature is

T=14πddzf(z)|z=zH=d4πzH(1(d2)2Q2zH2(d1)d(d1)).T=-\left.\frac{1}{4\pi}\frac{d}{dz}f(z)\right\rvert_{z=z_{H}}=\frac{d}{4\pi z_{H}}\left(1-\frac{(d-2)^{2}Q^{2}z_{H}^{2(d-1)}}{d(d-1)}\right). (16)

The total energy ε\varepsilon of the system, which is a dimensionless quantity with limits 2(d1)d2ε1\frac{2(d-1)}{d-2}\leq\varepsilon\leq 1, is given as follows [46]

ε(T,μ)=b02n1+1+d22π2b0b1(μ2T2),\varepsilon(T,\mu)=b_{0}-\frac{2n}{1+\sqrt{1+\frac{d^{2}}{2\pi^{2}b_{0}b_{1}}\left(\frac{\mu^{2}}{T^{2}}\right)}}\;, (17)

where the numerical constants b0b_{0} and b1b_{1} are listed in appendix B in eqs. (85) and (86) respectively.

The lapse function, chemical potential and temperature may then be expressed in terms of ε\varepsilon as follows

f(z)=1ε(zzH)d+(ε1)(zzH)2(d1),f(z)=1-\varepsilon\left(\frac{z}{z_{H}}\right)^{d}+(\varepsilon-1)\left(\frac{z}{z_{H}}\right)^{2(d-1)},

μ=1zH(d1)(d2)(ε1),\mu=\frac{1}{z_{H}}\sqrt{\frac{(d-1)}{(d-2)}(\varepsilon-1)}\;\;,

T=2(d1)(d2)ε4πzH.T=\frac{2(d-1)-(d-2)\varepsilon}{4\pi z_{H}}\;. (18)

The effective temperature TeffT_{\textit{ef{}f}}, which describes the number of microstates for a given temperature and chemical potential, may be defined as [46]

Teff(T,μ)d4πzH=T2[1+1+d22π2b0b1(μ2T2)].T_{\textit{ef{}f}}(T,\mu)\equiv\frac{d}{4\pi z_{H}}=\frac{T}{2}\left[1+\sqrt{1+\frac{d^{2}}{2\pi^{2}b_{0}b_{1}}\left(\frac{\mu^{2}}{T^{2}}\right)}\;\right]. (19)

We now proceed to review the computation the area of a RT surface homologous to a single subsystem AA with long rectangular strip geometry in a CFTdCFT_{d} dual to a bulk AdSd+1AdS_{d+1}-RN black hole. The subsystem AA is specified by

xx1[l2,l2],xi[L2,L2],i=2,,d2;x\equiv x^{1}\in\left[-\frac{l}{2},\;\frac{l}{2}\right],\qquad x^{i}\in\left[-\frac{L}{2},\;\frac{L}{2}\right],\qquad i=2,\dots,d-2\,; (20)

with lLl\ll L. The area 𝒜A\mathcal{A}_{A} of the RT surface homologous to this subsystem may then be expressed as

𝒜A=2Ld2zd10l/2dxz(x)2(d1)=2Ld2zd1azdzzd1f(z)[z2(d1)z2(d1)],\mathcal{A}_{A}=2L^{d-2}z_{*}^{d-1}\int_{0}^{l/2}\frac{dx}{z(x)^{2(d-1)}}=2L^{d-2}z_{*}^{d-1}\int_{a}^{z_{*}}\frac{dz}{z^{d-1}\sqrt{f(z)\left[z_{*}^{2(d-1)}-z^{2(d-1)}\right]}}\,, (21)

where zz_{*} is the turning point of the RT surface in the xx direction, and aa is the UV cut off of the CFTdCFT_{d}. The length ll of the rectangular strip in the xx direction is related to zz_{*} as

l=20zdzf(z)[(zz)2(d1)1].l=2\int_{0}^{z_{*}}\frac{dz}{\sqrt{f(z)\left[\left(\frac{z_{*}}{z}\right)^{2(d-1)}-1\right]}}\,. (22)

The integral in eq. (22) may be expressed as the following double sum [46]

l=zd1n=0k=0nΓ[12+n]Γ[d(n+k+1)2k2(d1)]εnk(1ε)kΓ[1+nk]Γ[k+1]Γ[d(n+k+2)2k12(d1)](zzH)nd+k(d2).l=\frac{z_{*}}{d-1}\sum_{n=0}^{\infty}\sum_{k=0}^{n}\frac{\Gamma\left[\frac{1}{2}+n\right]\Gamma\left[\frac{d(n+k+1)-2k}{2(d-1)}\right]\varepsilon^{n-k}(1-\varepsilon)^{k}}{\Gamma[1+n-k]\Gamma[k+1]\Gamma\left[\frac{d(n+k+2)-2k-1}{2(d-1)}\right]}\left(\frac{z_{*}}{z_{H}}\right)^{nd+k(d-2)}. (23)

The area of the RT surface in eq. (21) may likewise be described as a double sum as follows [46]

𝒜A\displaystyle\mathcal{A}_{A} =2d2(La)d2+2Ld2zd2[πΓ[d22(d1)]2(d1)Γ[12(d1)]]\displaystyle=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+2\frac{L^{d-2}}{z_{*}^{d-2}}\left[\frac{\sqrt{\pi}\Gamma\left[-\frac{d-2}{2(d-1)}\right]}{2(d-1)\Gamma\left[\frac{1}{2(d-1)}\right]}\right] (24)
+Ld2(d1)zd2[n=1k=0nΓ[12+n]Γ[d(n+k1)2k+22(d1)]εnk(1ε)kΓ[1+nk]Γ[k+1]Γ[d(n+k)2k+12(d1)](zzH)nd+k(d2)].\displaystyle+\frac{L^{d-2}}{(d-1)z_{*}^{d-2}}\left[\sum_{n=1}^{\infty}\sum_{k=0}^{n}\frac{\Gamma\left[\frac{1}{2}+n\right]\Gamma\left[\frac{d(n+k-1)-2k+2}{2(d-1)}\right]\varepsilon^{n-k}(1-\varepsilon)^{k}}{\Gamma[1+n-k]\Gamma[k+1]\Gamma\left[\frac{d(n+k)-2k+1}{2(d-1)}\right]}\left(\frac{z_{*}}{z_{H}}\right)^{nd+k(d-2)}\right].

In what follows the area and the turning point of the corresponding bulk RT surface are expressed as perturbative expansions in terms of the specified parameters ε\varepsilon and TeffT_{\textit{ef{}f}}, for various limits of the temperature TT and the chemical potential μ\mu of the dual CFTdCFT_{d}. The results are then utilized to compute the holographic entanglement negativity for various mixed state configurations of the single subsystem in the CFTdCFT_{d}s dual to bulk AdSd+1AdS_{d+1}-RN non extremal and extremal black holes, from the conjecture as described in eq. (10).

Here it is required to consider two auxiliary subsystems B1B_{1} and B2B_{2} on either side of the subsystem AA along the direction of the x1x^{1} axis (henceforth referred to as the partitioning direction), specified by

x1[(b+l2),l2],xi[L2,L2];x^{1}\in\left[-\left(b+\frac{l}{2}\right),\;-\frac{l}{2}\right],\qquad x^{i}\in\left[-\frac{L}{2},\;\frac{L}{2}\right]; (25)

and

x1[l2,(b+l2)],xi[L2,L2];x^{1}\in\left[\frac{l}{2},\;\left(b+\frac{l}{2}\right)\right],\qquad x^{i}\in\left[-\frac{L}{2},\;\frac{L}{2}\right]; (26)

respectively, where i=2,,d2i=2,\dots,d-2, with lbl\ll b and bLb\ll L. Here LL denotes the length of the strip in the remaining (d2)(d-2) spatial directions. Note that we have chosen the subsystems B1B_{1} and B2B_{2} to be symmetric along the partitioning direction. This leads to the equality of the areas 𝒜B1=𝒜B2\mathcal{A}_{B_{1}}=\mathcal{A}_{B_{2}} and 𝒜AB1=𝒜AB2\mathcal{A}_{A\cup B_{1}}=\mathcal{A}_{A\cup B_{2}}. With this identification, the holographic entanglement negativity for the above configuration may be computed from eq. (11).

Non extremal AdSd+𝟏AdS_{d+1}-RN black holes

In this subsection we utilize the area results for RT surfaces desribed above to compute the holographic entanglement negativity for finite temperature mixed states of single subsystems with a long rectangular strip geometry in CFTdCFT_{d}s dual to non extremal AdSd+1AdS_{d+1}-RN black holes for various non trivial limits of the temperature TT and the chemical potential μ\mu describing different regimes. This involves the perturbative expansion of the areas of the bulk RT surfaces for appropriate combinations of the subsystem and the auxiliary subsystems, as mentioned in subsection 3.1. We begin with the regime of small chemical potential and low temperature in subsection 3.2.1, proceeding to the case of small chemical potential and low temperature in subsection 3.2.2 and the large chemical potential and low temperature limit in subsection 3.2.3.

Small chemical potential - low temperature

The regime of small chemical potential and low temperature is described by the conditions μl1\mu l\ll 1 and Tl1Tl\ll 1. Depending on whether TμT\ll\mu or TμT\gg\mu [46], this gives rise to two different situations as described below.

(i) Tlμl𝟏Tl\ll\mu l\ll 1

We first consider the case where TμT\ll\mu, which, together with μl1\mu l\ll 1 and Tl1Tl\ll 1, may be recast as Tlμl1Tl\ll\mu l\ll 1. In this limit, the parameters ε(T,μ)\varepsilon(T,\mu) and Teff(T,μ)T_{\textit{ef{}f}}(T,\mu), given by eqs. (17) and eq. (19), may be expanded about Tμ=0\frac{T}{\mu}=0 to the leading order as [46]

εb02nπ2b0b1d(Tμ),\varepsilon\approx b_{0}-\frac{2n\pi\sqrt{2b_{0}b_{1}}}{d}\left(\frac{T}{\mu}\right), (27)

Teff12(μdπ2b0b1+T).T_{\textit{ef{}f}}\approx\frac{1}{2}\left(\frac{\mu d}{\pi\sqrt{2b_{0}b_{1}}}+T\right). (28)

In this regime the turning point of the bulk RT surface remains far away from the black hole horizon (zzHz_{*}\ll z_{H}). Thus eq. (23) may be utilized to expand zz_{*} in terms of (lzH)d\left(\frac{l}{z_{H}}\right)^{d} as

z=lΓ[12(d1)]2πΓ[d2(d1)][112(d+1)21d1dΓ[1+12(d1)]Γ[12(d1)]d+1πd+12Γ[12+1d1]Γ[d2(d1)]dε(lzH)d]\displaystyle z_{*}=\frac{l~{}\Gamma\left[\frac{1}{2(d-1)}\right]}{2\sqrt{\pi}~{}\Gamma\left[\frac{d}{2(d-1)}\right]}\left[1-\frac{1}{2(d+1)}\frac{2^{\frac{1}{d-1}-d}\Gamma\left[1+\frac{1}{2(d-1)}\right]\Gamma\left[\frac{1}{2(d-1)}\right]^{d+1}}{\pi^{\frac{d+1}{2}}\Gamma\left[\frac{1}{2}+\frac{1}{d-1}\right]\Gamma\left[\frac{d}{2(d-1)}\right]^{d}}\varepsilon~{}\left(\frac{l}{z_{H}}\right)^{d}\right] (29)
+𝒪[(lzH)2(d1)].\displaystyle+\mathcal{O}\left[\left(\frac{l}{z_{H}}\right)^{2(d-1)}\right].

The area of the RT surface homologous to a generic subsystem Σ\Sigma may likewise be computed from eq. (24) and re-expressed in terms of the parameters ε\varepsilon and TeffT_{\textit{ef{}f}} as follows [46]

𝒜Σ=2d2(La)d2+𝒮0(Ll)d2+ε𝒮0𝒮1(4πTeffd)dLd2l2+𝒪[(Teffl)2(d1)].\mathcal{A}_{\Sigma}=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\varepsilon\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d}L^{d-2}l^{2}+\mathcal{O}\left[(T_{\textit{ef{}f}}l)^{2(d-1)}\right]. (30)

where ll is the length of the subsystem Σ\Sigma in the partitioning direction, and the numerical constants 𝒮0\mathcal{S}_{0} and 𝒮1\mathcal{S}_{1} are listed in appendix B in eqs. (87) and (88) respectively.

As discussed in [34], the subsystems B1B_{1} and AB1A\cup B_{1} in the CFTdCFT_{d} become infinite along the direction of partition in the bipartite limit BAcB\to A^{c}. Thus the RT surfaces homologous to these subsystems penetrate deep into the bulk to approach the black hole horizon (zzHz_{*}\approx z_{H}). For these surfaces the appropriate Taylor expansion gives rise to the following expression for the area [46]

𝒜Σ=2d2(La)d2\displaystyle\mathcal{A}_{\Sigma}=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2} +(4πTeffd)d1Ld2l\displaystyle+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l (31)
+(4πTeffd)d2Ld2(N0+N1(b0ε))+𝒪[Tμ].\displaystyle+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-2}\left.L^{d-2}\middle(N_{0}+N_{1}(b_{0}-\varepsilon)\right)+\mathcal{O}\left[\frac{T}{\mu}\right].

The numerical constants N0N_{0} and N1N_{1} in eq. (31) are listed in appendix B in eqs. (89) and (90) respectively.

Utilizing eqs. (30) and (31), the holographic entanglement negativity in the limit of small chemical potential and low temperature (with TμT\ll\mu), for the mixed state configuration of the subsystem AA may now be obtained from the construction in [27] by employing eq. (11) as follows111111Eq. (30) is employed to calculate 𝒜A\mathcal{A}_{A}, while 𝒜B1\mathcal{A}_{B_{1}} and 𝒜AB1\mathcal{A}_{A\cup B_{1}} are obtained utilizing eq. (31), followed by implementation of the bipartite limit l/b0l/b\to 0.

A=38GN(d+1)[2d2(La)d2\displaystyle\mathcal{E}_{A}=\frac{3}{8G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}\right. +𝒮0(Ll)d2+ε𝒮0𝒮1(4πTeffd)dLd2l2\displaystyle+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\varepsilon\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d}L^{d-2}l^{2} (32)
(4πTeffd)d1Ld2l]+𝒪[(Teffl)2(d1)].\displaystyle\left.-\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l\right]+\mathcal{O}\left[(T_{\textit{ef{}f}}l)^{2(d-1)}\right].

The corresponding holographic entanglement entropy is given from eq. (12) through the use of eq. (30) as follows

𝒮A=14GN(d+1)[2d2(La)d2+𝒮0(Ll)d2+ε𝒮0𝒮1(4πTeffd)dLd2l2]\displaystyle\mathcal{S}_{A}=\frac{1}{4G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\varepsilon\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d}L^{d-2}l^{2}\right] (33)
+𝒪[(Teffl)2(d1)].\displaystyle+\mathcal{O}\left[(T_{\textit{ef{}f}}l)^{2(d-1)}\right].

Comparing eq. (32) with eq. (33) we arrive at

A=32(𝒮A𝒮Ath),\mathcal{E}_{A}=\left.\frac{3}{2}\middle(\mathcal{S}_{A}-\mathcal{S}_{A}^{\textit{th}}\right), (34)

where 𝒮Ath=14GN(d+1)(4πTeffd)d1V\mathcal{S}_{A}^{\textit{th}}=\frac{1}{4G_{N}^{(d+1)}}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}V with VLd2lV\equiv L^{d-2}l being the volume of the subsystem under consideration.

The first two terms on the right hand side of eq. (32) for the holographic entanglement negativity correspond to the holographic entanglement negativity for the zero temperature pure state configuration of a single subsystem in a CFTdCFT_{d} dual to bulk pure AdSd+1AdS_{d+1} geometry, while the rest of the terms characterize the perturbative corrections due to the temperature and chemical potential of the black hole. The relation in eq. (34) describes that the holographic entanglement negativity is characterized at the leading order by the elimination of the volume dependent thermal entropy from the holographic entanglement entropy for the finite temperature mixed state of a single subsystem reproducing earlier results in the literature. Note that this is in contrast with the holographic entanglement entropy given in eq. (33), which scales with the volume. This is expected as the entanglement negativity characterizes an upper bound on the distillable entanglement in quantum information theory [6, 7]. The elimination of the thermal contributions to the leading order hence serves as a strong consistency check for the higher dimensional holographic construction described in [27, 34].121212This elimination of the thermal contributions is expected to hold to all orders. Note also that an exact analysis of the area of the RT surfaces is there in the literature as mentioned in footnote 7. Note that this is similar to the corresponding lower dimensional AdS3/CFT2AdS_{3}/CFT_{2} scenario where it is exact in the large central charge limit.


(ii) μlTl𝟏\mu l\ll Tl\ll 1

Next we consider the limit TμT\gg\mu. This condition may be combined with μl1\mu l\ll 1 and Tl1Tl\ll 1 into μlTl1\mu l\ll Tl\ll 1. For this case the parameters ε(T,μ)\varepsilon(T,\mu) and Teff(T,μ)T_{\textit{ef{}f}}(T,\mu) described by eqs. (17) and (19) have the following Taylor expansions around μT=0\frac{\mu}{T}=0 [46]

ε=1+d2(d2)16π2(d1)(μT)2+𝒪[(μT)4],\varepsilon=1+\frac{d^{2}(d-2)}{16\pi^{2}(d-1)}\left(\frac{\mu}{T}\right)^{2}+\mathcal{O}\left[\left(\frac{\mu}{T}\right)^{4}\right], (35)
Teff=T[1+d(d2)216π2(d1)(μT)2+𝒪[(μT)4]].T_{\textit{ef{}f}}=T\left[1+\frac{d(d-2)^{2}}{16\pi^{2}(d-1)}\left(\frac{\mu}{T}\right)^{2}+\mathcal{O}\left[\left(\frac{\mu}{T}\right)^{4}\right]\right]. (36)

Once again the conditions Tl1Tl\ll 1 and μl1\mu l\ll 1 entail that zzHz_{*}\ll z_{H}. The expression for zz_{*} may be obtained by expanding eq. (23) in terms of (lzH)d\left(\frac{l}{z_{H}}\right)^{d}, and is same as the one given in eq. (29). The area of the RT surface homologous to a generic subsystem Σ\Sigma is once again determined by expansion of eq. (24) in terms of (lzH)d\left(\frac{l}{z_{H}}\right)^{d} as [46]

𝒜Σ=2d2(La)d2+𝒮0(Ll)d2+ε𝒮0𝒮1(4πTeffd)dLd2l2+𝒪[(Teffl)2(d1)],\mathcal{A}_{\Sigma}=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\varepsilon\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d}L^{d-2}l^{2}+\mathcal{O}\left[(T_{\textit{ef{}f}}l)^{2(d-1)}\right], (37)

where the numerical constants 𝒮0\mathcal{S}_{0} and 𝒮1\mathcal{S}_{1} are the same as in eq. (30), and are given in appendix B in eqs. (87) and (88) respectively.

We note here that although the areas of the bulk RT surface given in eqs. (30) and (37) have identical forms for both Tlμl1Tl\ll\mu l\ll 1 and μlTl1\mu l\ll Tl\ll 1, the underlying parameters ε\varepsilon and TeffT_{\textit{ef{}f}} are different for the two cases as described in eqs. (27) and (28), and eqs. (35) and (36) respectively. As outlined in the discussion just preceding eq. (31), the subsystems B1B_{1} and AB1A\cup B_{1} become infinite along the partitioning direction in the limit BAcB\to A^{c}, and consequently the corresponding RT surfaces extend deep into the bulk so that the turning points of these surfaces are close to the horizon (zzHz_{*}\approx z_{H}). The relevant expression for the area of these surfaces is given as follows [46]

𝒜Σ=2d2(La)d2+(4πTeffd)d1Ld2l+(4πTeffd)d2Ld2γd(μT).\mathcal{A}_{\Sigma}=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-2}L^{d-2}\gamma_{d}\left(\frac{\mu}{T}\right). (38)

The function γd(μT)\gamma_{d}\left(\frac{\mu}{T}\right) in eq. (38) admits a perturbative expansion in μT\frac{\mu}{T}, which has been described in appendix B in eq. (91).

The holographic entanglement negativity for the required finite temperature mixed state configuration of the single subsystem in question with long rectangular strip geometry, in the limit of small chemical potential and low temperature (with TμT\gg\mu) may then be computed131313Once again 𝒜A\mathcal{A}_{A} is computed by applying eq. (37), whereas eq. (38) is utilized to obtain 𝒜B1\mathcal{A}_{B_{1}} and 𝒜AB1\mathcal{A}_{A\cup B_{1}}, and afterwards we take the bipartite limit l/b0l/b\to 0. from eqs. (37) and (38), from the conjecture described by eq. (11), and has a form identical to that in eq. (32). Relations described in eqs. (33) and (34), and the discussions after eq. (34) also follow identically.

Small chemical potential - high temperature

In this subsection we proceed to investigate the holographic entanglement negativity in the limit of small chemical potential and high temperature. As detailed in [46], this regime is given by the conditions μT\mu\ll T and Tl1Tl\gg 1. The parameters ε\varepsilon and TeffT_{\textit{ef{}f}} may be given in the limit μT\mu\ll T (refer to case (ii), subsection 3.2.1 for details) by eqs. (35) and (36). The other condition Tl1Tl\gg 1 implies zzHz_{*}\approx z_{H}. Thus the area of the bulk extremal surface may be obtained by perturbative expansion of eq. (24) around zzH=1\frac{z_{*}}{z_{H}}=1, and is given by eq. (38).

Substituting eq. (38) into eq. (11), the holographic entanglement negativity in the limit of small chemical potential and high temperature, for the mixed state configuration of the single subsystem under consideration may be determined as follows

A=38GN(d+1)[2d2(La)d2+(4πTeffd)d2Ld2γd(μT)].\mathcal{E}_{A}=\frac{3}{8G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-2}L^{d-2}\gamma_{d}\left(\frac{\mu}{T}\right)\right]. (39)

Next we employ eq. (38) to compute the corresponding holographic entanglement entropy from eq. (12) as follows

𝒮A=14GN(d+1)[2d2(La)d2+(4πTeffd)d1Ld2l+(4πTeffd)d2Ld2γd(μT)].\mathcal{S}_{A}=\frac{1}{4G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-2}L^{d-2}\gamma_{d}\left(\frac{\mu}{T}\right)\right]. (40)

Comparing eqs. (39) and (40) we again obtain

A=32(𝒮A𝒮Ath),\mathcal{E}_{A}=\left.\frac{3}{2}\middle(\mathcal{S}_{A}-\mathcal{S}_{A}^{\textit{th}}\right), (41)

where 𝒮Ath=14GN(d+1)(4πTeffd)d1Ld2l\mathcal{S}_{A}^{\textit{th}}=\frac{1}{4G_{N}^{(d+1)}}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l as before.

Once again we observe that the expression for the holographic entanglement negativity to the leading order in eq. (39) for the mixed state in question is purely dependent on the area Ld2L^{d-2} of the entangling surface between the single subsystem with the rest of the system. As earlier this indicates the vanishing of the volume dependent thermal contribution, as may be seen from eq. (41).

Large chemical potential - low temperature

The case for large chemical potential and low temperature regime is described by the conditions μl1\mu l\gg 1 and TμT\ll\mu. As mentioned earlier (case (i), subsection 3.2.1), the condition TμT\ll\mu implies that the parameters ε(T,μ)\varepsilon(T,\mu) and Teff(T,μ)T_{\textit{ef{}f}}(T,\mu) admit of the expansions described by eqs. (27) and (28) respectively. Similar to subsection 3.2.2, the condition μl1\mu l\gg 1 dictates zzHz_{*}\approx z_{H}. Hence, the area of the bulk RT surface may once again be obtained through the perturbative expansion of eq. (24) around zzH=1\frac{z_{*}}{z_{H}}=1 [46], and is given by eq. (31).

Utilizing eq. (31), we may now obtain the holographic entanglement negativity in the limit of large chemical potential and low temperature, for the mixed state configuration of the subsystem AA, from the construction advanced in [27] by employing eq. (11) as follows

A=38GN(d+1)[2d2(La)d2+(4πTeffd)d2Ld2(N0+N1(b0ε))]+𝒪[Tμ].\mathcal{E}_{A}=\frac{3}{8G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-2}\left.L^{d-2}\middle(N_{0}+N_{1}(b_{0}-\varepsilon)\right)\right]+\mathcal{O}\left[\frac{T}{\mu}\right]. (42)

Similarly the entanglement entropy for this regime may also be determined from eq. (12) using eq. (31) as follows

𝒮A=14GN(d+1)[2d2(La)d2\displaystyle\mathcal{S}_{A}=\frac{1}{4G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}\right. +(4πTeffd)d1Ld2l\displaystyle+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l (43)
+(4πTeffd)d2Ld2(N0+N1(b0ε))]+𝒪[Tμ].\displaystyle\left.+\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-2}\left.L^{d-2}\middle(N_{0}+N_{1}\left(b_{0}-\varepsilon\right)\right)\right]+\mathcal{O}\left[\frac{T}{\mu}\right].

Once again comparison of eqs. (42) and (43) leads us to

A=32(𝒮A𝒮Ath),\mathcal{E}_{A}=\left.\frac{3}{2}\middle(\mathcal{S}_{A}-\mathcal{S}_{A}^{\textit{th}}\right), (44)

where as earlier 𝒮Ath=14GN(d+1)(4πTeffd)d1Ld2l\mathcal{S}_{A}^{\textit{th}}=\frac{1}{4G_{N}^{(d+1)}}\left(\frac{4\pi T_{\textit{ef{}f}}}{d}\right)^{d-1}L^{d-2}l.

Again the holographic entanglement negativity to the leading order in eq. (42), obtained in the limit of large chemical potential and high temperature, is purely dependent on the area of the relevant entangling surface. As earlier this is due to the precise elimination of the volume dependent thermal term to the leading order described by eq. (44), in conformity with the standard quantum information theory expectations, and serves to further substantiate the holographic conjecture proposed in [27] and our construction.

Extremal AdSd+𝟏AdS_{d+1}-RN black holes

In this subsection we outline the procedure to obtain the holographic entanglement negativity for the zero temperature mixed state configurations of single subsystems in CFTdCFT_{d}s dual to bulk extremal AdSd+1AdS_{d+1}-RN black holes. The relevant parameters for these bulk extremal black hole geometries are given as [46]

Q2=d(d1)L2(d2)2zH2(d1),Q^{2}=\frac{d(d-1)L^{2}}{(d-2)^{2}z_{H}^{2(d-1)}}\;, (45)
ε=b1,\varepsilon=b_{1}, (46)
μ=1zHb0b12=1zHd(d1)(d2)2,\mu=\frac{1}{z_{H}}\sqrt{\frac{b_{0}b_{1}}{2}}=\frac{1}{z_{H}}\sqrt{\frac{d(d-1)}{(d-2)^{2}}}\;, (47)
Teff=μd2π2b0b1.T_{\textit{ef{}f}}=\frac{\mu d}{2\pi\sqrt{2b_{0}b_{1}}}\;. (48)

In the above equations the parameter QQ describes the charge of the extremal AdSd+1AdS_{d+1}-RN black hole, while TeffT_{\textit{ef{}f}} as earlier denotes the effective temperature.

In the next two subsections we proceed to obtain perturbative expansions of the area of the RT surface homologous to the single subsystem with long rectangular strip geometry, for two different limits of the charge QQ. The expression obtained will then be utilized to compute the holographic entanglement negativity for the mixed state configuration under consideration from the conjecture proposed in [27]. In subsection 3.3.1 we study the case involving a small chemical potential, while subsection 3.3.2 describes the large chemical potential regime.

Small chemical potential

The limit of small chemical potential is defined by the condition μl1\mu l\ll 1, which entails zzHz_{*}\ll z_{H}. The solution of eq. (22) in this limit once again leads to eq. (29) [46]. The area of the RT surface homologous to a generic subsystem Σ\Sigma may be similarly obtained from eq. (24) and re-expressed in terms μ\mu as follows [46]

𝒜Σ=2d2(La)d2+𝒮0(Ll)d2+2(d1)d2𝒮0𝒮1(μ(d2)d(d1))dLd2l2\displaystyle\mathcal{A}_{\Sigma}=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\frac{2(d-1)}{d-2}\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d}L^{d-2}l^{2} (49)
+𝒪[(μl)2(d1)],\displaystyle+\mathcal{O}\left[(\mu l)^{2(d-1)}\right],

where the constants 𝒮0\mathcal{S}_{0} and 𝒮1\mathcal{S}_{1} are identical to those in the earlier cases described in subsection 3.2.1, and are given in appendix B in eqs. (87) and (88) respectively.

As discussed in earlier subsections, the turning points of the RT surfaces corresponding to the infinite auxiliary subsystems approach the black hole horizon (zzHz_{*}\approx z_{H}). The appropriate area expansion for these surfaces are given as [46]

𝒜Σ=2d2(La)d2+(μ(d2)d(d1))d1Ld2l+(μ(d2)d(d1))d2Ld2N(b0),\mathcal{A}_{\Sigma}=\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-1}L^{d-2}l+\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-2}L^{d-2}N(b_{0}), (50)

where N(b0)N(b_{0}) is the value of N(ε)N(\varepsilon) at ε=b0\varepsilon=b_{0}. It is listed in appendix B in eq. (92).

Utilizing eqs. (49) and (50), it is now possible to compute the holographic entanglement negativity in the limit of small chemical potential, for the bipartite configuration of the single subsystem in question by employing eq. (11) as follows

A=38GN(d+1)[2d2(La)d2+𝒮0(Ll)d2+2(d1)d2𝒮0𝒮1(μ(d2)d(d1))dLd2l2\displaystyle\mathcal{E}_{A}=\frac{3}{8G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\frac{2(d-1)}{d-2}\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d}L^{d-2}l^{2}\right. (51)
(μ(d2)d(d1))d1Ld2l]+𝒪[(μl)2(d1)].\displaystyle\left.-\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-1}L^{d-2}l\right]+\mathcal{O}\left[(\mu l)^{2(d-1)}\right].

We then employ eq. (49) to obtain the following expression for the holographic entanglement entropy from eq. (12)

𝒮A=14GN(d+1)[2d2(La)d2+𝒮0(Ll)d2+2(d1)d2𝒮0𝒮1(μ(d2)d(d1))dLd2l2]\displaystyle\mathcal{S}_{A}=\frac{1}{4G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\mathcal{S}_{0}\left(\frac{L}{l}\right)^{d-2}+\frac{2(d-1)}{d-2}\mathcal{S}_{0}\mathcal{S}_{1}\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d}L^{d-2}l^{2}\right] (52)
+𝒪[(μl)2(d1)].\displaystyle+\mathcal{O}\left[(\mu l)^{2(d-1)}\right].

Through the comparison of eqs. (51) and (52) we arrive at the following relation

A=32(𝒮A𝒮Ath),\mathcal{E}_{A}=\left.\frac{3}{2}\middle(\mathcal{S}_{A}-\mathcal{S}_{A}^{\textit{th}}\right), (53)

where we have 𝒮Ath=14GN(d+1)(μ(d2)d(d1))d1Ld2l\mathcal{S}_{A}^{\textit{th}}=\frac{1}{4G_{N}^{(d+1)}}\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-1}L^{d-2}l.

In eq. (51) we note that the first two terms describe the holographic entanglement negativity for the zero temperature pure state configuration of a single subsystem with long rectangular strip geometry in a CFTdCFT_{d} dual to bulk pure AdSd+1AdS_{d+1} geometry. The rest of the terms represent the corrections due to the chemical potential of the black hole. Once again we observe the elimination of the volume dependent “thermal” contribution to the leading order141414The thermal contribution characterizes the entropy due to the ground state degeneracy of the black hole, corresponding to the vacuum degeneracy of the dual CFTCFT. in eq. (53) which is in conformity with quantum information theory and substantiates the holographic conjecture in [27].

Large chemical potential

In this subsection we investigate the large chemical potential regime, specified by the condition μl1\mu l\gg 1. With this condition it may be observed from eq. (47) that the horizon radius is large, and the turning point of the bulk RT surface is thus close to the horizon (zzHz_{*}\approx z_{H}). In this regime, the area of the RT surface homologous to a generic subsystem Σ\Sigma with long rectangular strip geometry may be determined by a perturbative computation of the integral in eq. (24) around zzH=1\frac{z_{*}}{z_{H}}=1, and is given by eq. (50).

The holographic entanglement negativity in the limit of large chemical potential, for the bipartite configuration of the single subsystem under consideration may now be derived through the use of eq. (50) by employing eq. (11) as follows

A=38GN(d+1)[2d2(La)d2+(μ(d2)d(d1))d2Ld2N(b0)].\mathcal{E}_{A}=\frac{3}{8G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-2}L^{d-2}N(b_{0})\right]. (54)

The corresponding expression for the holographic entanglement entropy may then be obtained from eq. (50) utilizing eq. (12) as

𝒮A=14GN(d+1)[2d2(La)d2+(μ(d2)d(d1))d1Ld2l+(μ(d2)d(d1))d2Ld2N(b0)].\mathcal{S}_{A}=\frac{1}{4G_{N}^{(d+1)}}\left[\frac{2}{d-2}\left(\frac{L}{a}\right)^{d-2}+\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-1}L^{d-2}l+\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-2}L^{d-2}N(b_{0})\right]. (55)

Once again comparing eqs. (54) and (55) we arrive at

A=32(𝒮A𝒮Ath),\mathcal{E}_{A}=\left.\frac{3}{2}\middle(\mathcal{S}_{A}-\mathcal{S}_{A}^{\textit{th}}\right), (56)

where as before we have identified 𝒮Ath=14GN(d+1)(μ(d2)d(d1))d1Ld2l\mathcal{S}_{A}^{\textit{th}}=\frac{1}{4G_{N}^{(d+1)}}\left(\frac{\mu(d-2)}{\sqrt{d(d-1)}}\right)^{d-1}L^{d-2}l as the “thermal” entropy.

We once again note that the holographic entanglement negativity computed in eq. (54) for the large chemical potential regime is purely dependent on the area of the appropriate entangling surface, due to the elimination of the volume dependent “thermal” contribution to the leading order described in eq. (56), which agrees with quantum information theory expectations and substantiates the holographic construction in [27] as earlier.

Discussion of the results

The leading order results derived in sections 3.2 and 3.3 for the holographic entanglement negativity of a single subsystem in CFTdCFT_{d}s dual to bulk non extremal and extremal AdSd+1AdS_{d+1}-RN black holes for various limits of the relevant parameters are consistent with the quantum information theory expectations. It is observed that in the limit of small chemical potential and low temperature the leading contribution to the holographic entanglement negativity arises from the zero temperature configuration of single subsystem in the CFTdCFT_{d} dual to the bulk pure AdSd+1AdS_{d+1} geometry, while the correction terms depend on the chemical potential μ\mu and the temperature TT. The reason for this is that the bulk RT surface is located near the boundary far away from the black hole horizon, and thus the leading contribution arises from the near boundary pure AdSd+1AdS_{d+1} geometry.

On the other hand for both the regimes of small chemical potential and high temperature, and large chemical potential and low temperature, the leading part is purely dependent on the area of the entangling surface between the subsystem and the rest of the system, while the volume dependent thermal term is eliminated. In this case the dominant contribution arises from the entanglement between the degrees of freedom at the entangling surface. For the case of small chemical potential with extremal AdSd+1AdS_{d+1}-RN black hole, the leading contribution involves the bulk pure AdSd+1AdS_{d+1} geometry term, and correction terms dependent on the chemical potential μ\mu. However for the large charge regime, the leading contribution arises only from the terms dependent purely on the area of the relevant entangling surface, while once again the volume dependent thermal term drops out.

From the above description we note that the holographic entanglement negativity to the leading order is purely a function of the area of the entangling surface in the CFTCFT even in the regimes of large charge or high temperature, whereas the corresponding holographic entanglement entropy involves volume dependent thermal contributions. This ensues from the precise cancellation of the volume dependent thermal terms in the expressions for the holographic entanglement negativity to the leading order and is expected to hold to all orders. As discussed in the previous sections, this elimination reproduces the corresponding phenomena observed also for the AdS3/CFT2AdS_{3}/CFT_{2} scenario where this elimination is exact in the large central charge limit and is expected to hold to all orders in the perturbation theory. Naturally these results provide strong consistency checks, and serve to further substantiate the higher dimensional holographic construction proposed in [27].

Summary and conclusions

To summarize, we have extended the higher dimensional holographic entanglement negativity proposal for bipartite states described by a single subsystem in CFTdCFT_{d}s dual to bulk AdSd+1AdS_{d+1} geometries described in [34] to CFTdCFT_{d}s with a conserved charge dual to bulk AdSd+1AdS_{d+1}-RN black holes. Our construction involved a specific algebraic sum of the areas of RT surfaces for specific combinations of the subsystem and auxiliary subsystems in the dual CFTdCFT_{d}.151515A plausible proof for the holographic construction from a gravitational path integral perspective in the literature has been reviewed in appendix A. In this context we have computed the holographic entanglement negativity for such bipartite states described by single subsystems with long rectangular strip geometries in CFTdCFT_{d}s with a conserved charge dual to bulk non extremal and extremal AdSd+1AdS_{d+1}-RN black holes.

To elucidate our construction we have reviewed the perturbative techniques utilized to compute the areas of the relevant bulk RT surfaces for combinations of subsystems in the dual CFTdCFT_{d} with a conserved charge as mentioned above. To establish the perturbative techniques in this context, we have described various regimes characterized by the total energy and an effective temperature which are functions of the temperature and the chemical potential conjugate to the charge of the bulk RN black hole in the AdSd+1/CFTdAdS_{d+1}/CFT_{d} scenario. We have obtained the holographic entanglement negativity for bipartite states of single subsystems with long rectangular strip geometries for the regimes described above in CFTdCFT_{d}s with a conserved charge dual to bulk non extremal and extremal AdSd+1AdS_{d+1}-RN black holes.

Interestingly for all the cases we have demonstrated the elimination of the thermal contributions for the holographic entanglement negativity at the leading order which is similar to the case of CFTdCFT_{d}s dual to bulk AdSd+1AdS_{d+1}-Schwarzschild black holes described earlier in the literature. This feature for the holographic entanglement negativity is also observed for the corresponding AdS3/CFT2AdS_{3}/CFT_{2} scenario where it is exact in the large central charge limit. The elimination of the thermal contributions to the leading order mentioned above is expected to hold to all orders in the perturbation theory. This also conforms to quantum information theory expectations for the entanglement negativity which characterizes an upper bound on the distillable entanglement of a mixed state. Furthermore, this serves as a strong non trivial consistency check for the higher dimensional holographic construction for the entanglement negativity as proposed and exemplified in the literature.

We emphasize here that there is an alternative holographic entanglement negativity proposal involving the area of a backreacted minimal entanglement wedge cross section due to a cosmic brane for the conical defect of the bulk replicated geometry. For subsystems with spherical entangling surfaces this backreaction is described by an overall dimension dependent numerical constant which may be explicitly computed. The equivalence of this construction modulo certain constants, with the proposal involving the algebraic sums of the lengths of bulk geodesics for combinations of intervals has been explicitly demonstrated in the literature for the AdS3/CFT2AdS_{3}/CFT_{2} scenario.

Keeping the above discussion in perspective, the holographic construction and the results described in the present article are also expected to involve appropriate backreaction effects. However for the rectangular strip subsystem geometries the explicit computation of these effects remains a non trivial open issue. As mentioned in footnote 15 a plausible bulk proof of our construction for spherical entangling surfaces in the literature has been reviewed in appendix A. However the extension of this proof to general subsystem geometries is an outstanding issue.

Our results described in this article serve to further substantiate the holographic entanglement negativity proposal for bipartite states in CFTdCFT_{d}s in the context of a generic higher dimensional AdSd+1/CFTdAdS_{d+1}/CFT_{d} scenario which may find interesting applications to the study of diverse issues involving mixed state entanglement for various phenomena ranging from many body systems to quantum gravity and black holes. We emphasize that our construction involving the algebraic sums of the areas of RT surfaces may be simpler than determining the EWCS in higher dimensions which is extremely non trivial. Furthermore our results are also expected to provide insight into the generalization of the bulk proof from a gravitational path integral obtained for spherical entangling surfaces to subsystems with arbitrary geometries, which is an extremely difficult outstanding open problem. We hope to address these interesting issues in the near future.

Acknowledgments

We would like to thank Vinay Malvimat and Debarshi Basu for crucial suggestions and discussions. The work of Sayid Mondal is supported in part by the Ministry of Science and Technology of Taiwan (under Grant No. 106-2112-M-033-007-MY3 and Grant No. 109-2112-M-033-005-MY3) and by the National Center for Theoretical Sciences (NCTS) of Taiwan.

Appendix A Derivation of holographic entanglement negativity

In this appendix we follow [39] to provide a brief review of a plausible proof of the holographic conjecture for the entanglement negativity of a single interval in a (1+1)(1+1) dimensional conformal field theory (CFT1+1CFT_{1+1}) dual to a bulk AdS3AdS_{3} geometry. As described in section 2 this involves a specific algebraic sum of the lengths (areas) of bulk space like geodesics (RT surfaces). This proof follows from [38] which provided the areas of backreacting cosmic branes through the analysis of replica symmetry breaking saddles for the bulk gravitational path integral.

Holographic Rényi entanglement entropy

We first breifly review the holographic characterization of the Rényi entropy following [52] in this section. In this work the gravitational path integral technique in [19] was employed to establish the relation between the holographic Rényi entropy of order nn of a subsystem in a dual CFT and the area of the corresponding bulk cosmic brane 𝒞n\mathcal{C}_{n} homologous to the subsystem as

n2n(n1nS(n)(A))=Area(𝒞n)4GN,n^{2}\frac{\partial}{\partial n}\left(\frac{n-1}{n}S^{(n)}(A)\right)=\frac{\mathrm{Area}\left(\mathcal{C}_{n}\right)}{4G_{N}}\,, (57)

where GNG_{N} is the Newton constant in appropriate bulk dimensions. Next we define a quantity 𝒜(n)\mathcal{A}^{(n)} in terms of the nthn^{th} holographic Rényi entropy S(n)(A)S^{(n)}(A) as follows

S(n)(A)=𝒜(n)4GN.S^{(n)}(A)=\frac{\mathcal{A}^{(n)}}{4G_{N}}\,. (58)

Substituting eq. (58) into eq. (57), we find the relation between 𝒜(n)\mathcal{A}^{(n)} and the area of the cosmic brane as

n2n(n1n𝒜(n))=Area(𝒞n).n^{2}\frac{\partial}{\partial n}\left(\frac{n-1}{n}\mathcal{A}^{(n)}\right)=\mathrm{Area}\left(\mathcal{C}_{n}\right). (59)

In the following analysis we will refer to 𝒜(n)\mathcal{A}^{(n)} as the area of the backreacting cosmic brane 𝒞n\mathcal{C}_{n} but it is to be understood in the context of eq. (59). The tension TnT_{n} of the cosmic brane 𝒞n\mathcal{C}_{n} is related to the replica index nn as

Tn=n14nGN.T_{n}=\frac{n-1}{4nG_{N}}\,. (60)

Note that the tension TnT_{n} of the cosmic brane 𝒞n\mathcal{C}_{n} in eq. (60) vanishes in the replica limit n1n\to 1, and the holographic construction in eq. (57) reduces to the Ryu-Takayanagi (RT) conjecture. For n1n\neq 1 however the tension in eq. (60) is non zero, which leads to a non trivial backreaction of the brane on the replicated bulk geometry. In this case the corresponding area of the cosmic brane 𝒞n\mathcal{C}_{n} is difficult to compute. Remarkably the backreaction and the correspondng area of the cosmic brane were explicitly computed for subsystems with spherical entangling surfaces in [53] where the authors employed a conformal map from a hyperbolic cylinder to the causal evolution of a subregion, which is enclosed by a spherical entangling surface in flat space. The entanglement entropy of a spherical region in a conformal field theory (CFT) on a flat space may thus be expressed in terms of an integral of the thermal entropy of a CFT on a hyperbolic cylinder. This translates to the computation of the horizon entropy for a certain topological black hole through the Wald formula in the AdS/CFTAdS/CFT framework. The holographic Rényi entropy S(n)(A)S^{(n)}(A) for a subsystem AA with spherical entangling surface in a dd-dimensional conformal field theory (CFTdCFT_{d}) on a flat space is then given by

S(n)(A)=𝒳d(n)S(A),S^{(n)}(A)=\mathcal{X}^{(n)}_{d}\,S(A)\,, (61)
𝒳d(n)=n2(n1)(2xnd2(1+xn2)),xn=1nd(1+12dn2+d2n2).\mathcal{X}^{(n)}_{d}=\frac{n}{2(n-1)}\left(2-x_{n}^{d-2}\left(1+x_{n}^{2}\right)\right),\qquad x_{n}=\frac{1}{nd}\left(1+\sqrt{1-2dn^{2}+d^{2}n^{2}}\right).

Note that the coordinate volume of the (d1)(d-1) dimensional hyperbolic plane was computed in [53] to evaluate the entanglement entropy S(A)S(A) in eq. (61). The holographic Rényi entropy of order half, which is related to the holographic entanglement negativity of a single subsystem, may then be obtained from eq. (61) as

S(1/2)(A)=𝒳dS(A),S^{(1/2)}(A)=\mathcal{X}_{d}\,S(A)\,, (62)
𝒳d=(12xdd2(1+xd2)1),xd=2d(1+1d2+d24).\mathcal{X}_{d}=\left(\frac{1}{2}x_{d}^{d-2}\left(1+x_{d}^{2}\right)-1\right),\qquad x_{d}=\frac{2}{d}\left(1+\sqrt{1-\frac{d}{2}+\frac{d^{2}}{4}}\right).

In eq. (62) and the rest of this appendix, we focus on the n=1/2n=1/2 case and thus we denote 𝒳d(1/2)\mathcal{X}_{d}^{(1/2)} by 𝒳d\mathcal{X}_{d}. We note here that for the AdS3/CFT2AdS_{3}/CFT_{2} scenario, eq. (62) leads to 𝒳2=3/2\mathcal{X}_{2}=3/2, which is also the numerical factor obtained through the replica technique for a single interval.

Observe that the result in eq. (62) is valid only for a single connected subsystem with spherical geometry in the AdSd+1/CFTdAdS_{d+1}/CFT_{d} framework (which reduces to a single spatial interval in the AdS3/CFT2AdS_{3}/CFT_{2} scenario). For configurations involving multiple disjoint subsystems, the corresponding bulk geometry incorporates multiple disjoint branes whose backreaction on each other complicates the computation of the Rényi entanglement entropy [52]. However in the context of the AdS3/CFT2AdS_{3}/CFT_{2} framework involving two disjoint intervals in proximity (cross ratio of the intervals x1x\to 1) the above result continues to be valid. For this specific configuration, the following function involving integral values of the replica index nn vanishes in the limits x0x\to 0 and x1x\to 1 when n>2n>2 as shown in [18]

Jn=I(n)(A:B)12(1+1n)I(A:B).J_{n}=I^{(n)}(A:B)-\frac{1}{2}\left(1+\frac{1}{n}\right)I(A:B)\,. (63)

In eq. (63) I(n)(A:B)I^{(n)}(A:B) and I(A:B)I(A:B) denote the Rényi mutual information of order nn and the mutual information between the subsystems AA and BB respectively. These are defined as follows

I(n)(A:B)=S(n)(A)+S(n)(B)S(n)(AB),I^{(n)}(A:B)=S^{(n)}(A)+S^{(n)}(B)-S^{(n)}(AB)\,, (64)
I(A:B)=S(A)+S(B)S(AB).I(A:B)=S(A)+S(B)-S(AB)\,. (65)

We may estimate the holographic Rényi mutual information from eq. (63) in the proximity limit x1x\to 1 as

I(n)(A:B)12(1+1n)I(A:B).I^{(n)}(A:B)\approx\frac{1}{2}\left(1+\frac{1}{n}\right)I(A:B)\,. (66)

Eq. (66) may be continued analytically to n1/2n\to 1/2 as

I(1/2)(A:B)32I(A:B).I^{(1/2)}(A:B)\approx\frac{3}{2}I(A:B)\,. (67)

The condition in eq. (67) requires that each of the Rényi entanglement entropies of order half for the subsystems AA, BB and ABAB, and their corresponding entanglement entropy in eqs. (64) and (65) are related through the prefactor 𝒳2=3/2\mathcal{X}_{2}=3/2 as shown in eq. (62). We observe that the backreaction of the cosmic branes on each other remains non trivial if the cross ratio xx is not close to 0 or 11, and the relation in eq. (67) may not hold. Thus our discussion of holographic entanglement negativity from now on will be restricted to the scenario when the intervals are in proximity (x1x\to 1) and the analytic continuation n1/2n\to 1/2 employed in eq. (67) holds. Having discussed the holographic generalization of the Rényi entanglement entropy and some of its properties for subsystems with spherical entangling surfaces, in section A.2 we review the derivation of the holographic entanglement negativity.

Holographic entanglement negativity

In this section we review the derivation of the holographic entanglement negativity of two disjoint subsystems in proximity in terms of the areas of the backreacting cosmic branes homologous to the subsystems as described in [39]. This utilizes the holographic derivation of the entanglement negativity as shown in [38]. To this end we may define the Rényi entanglement negativity of order kk following [9, 10] as

𝒩(k)(A:B)=Tr[(ρABTB)k]=Z[kA,B](Z[])k.\mathcal{N}^{(k)}(A:B)=\mathrm{Tr}\left[\left(\rho_{AB}^{T_{B}}\right)^{k}\right]=\frac{Z\left[\mathcal{M}_{k}^{A,B}\right]}{\left(Z\left[\mathcal{M}\right]\right)^{k}}\,. (68)

Here kA,B\mathcal{M}_{k}^{A,B} denotes the kk-sheeted replicated manifold in which the individual copies of the subsystems AA and BB are respectively glued cyclically and anti-cyclically, \mathcal{M} denotes the original manifold, and ZZ denotes the path integral over the corresponding manifold. The entanglement negativity is then defined by the analytic continuation of even Rényi negativities of order kk (k=2nk=2n) as given in eq. (68) to k1k\to 1 (n1/2n\to 1/2) as

(A:B)=limn1/2log𝒩(k=2n)(A:B).\mathcal{E}(A:B)=\lim_{n\to 1/2}\log\mathcal{N}^{(k=2n)}(A:B)\,. (69)

Note that the analytic continuation of the odd Rényi negativities of order kk (k=2n1k=2n-1) to k1k\to 1 (n1n\to 1) reproduces the trace normalization condition TrρABTB=1\mathrm{Tr}\rho_{AB}^{T_{B}}=1. Utilizing the AdS/CFT dictionary, we may express the ratio of the partition functions on the right hand side of eq. (68) for even k=2nk=2n at the leading order in terms of the corresponding gravitational on-shell actions in the bulk through the saddle point approximation as

Z[2nA,B](Z[])2n=Z[2nA,B](Z[])2n=exp(2nIgrav[]Igrav[2n]),\frac{Z\left[\mathcal{M}_{2n}^{A,B}\right]}{\left(Z\left[\mathcal{M}\right]\right)^{2n}}=\frac{Z\left[\mathcal{B}_{2n}^{A,B}\right]}{\left(Z\left[\mathcal{B}\right]\right)^{2n}}=\exp\left(2n\mathrm{I_{grav}}[\mathcal{B}]-\mathrm{I_{grav}}[\mathcal{B}_{2n}]\right)\,, (70)

where 2nA,B\mathcal{B}_{2n}^{A,B} and \mathcal{B} represent the bulk geometries corresponding to the manifolds 2nA,B\mathcal{M}_{2n}^{A,B} and \mathcal{M} respectively. Remarkably it was shown in [38] that unlike the Rényi negativities, the replica symmetric gravitational saddle is the same for both even and odd kk, and the corresponding entanglement negativity vanishes

(sym)(A:B)=limk1log𝒩(k)(sym)=0.\mathcal{E}^{(\mathrm{sym})}(A:B)=\lim_{k\to 1}\log\mathcal{N}_{(k)}^{(\mathrm{sym})}=0. (71)

Quite interestingly the authors demonstrated that the bulk replica symmetry is broken by the dominant saddle corresponding to the entanglement negativity.

We now briefly describe the holographic construction of the replica non symmetric saddle involving 2n2n copies of the bulk manifold. To this end we consider three codimension one homology hypersurfaces ΣA\Sigma_{A}, ΣB\Sigma_{B} and ΣAB¯\Sigma_{\overline{AB}} along which each such manifold has cuts. These hypersurfaces are non-overlapping and obey the homology conditions ΣX=XγX\partial\Sigma_{X}=X\cup\gamma_{X}, γX\gamma_{X} being the codimension two bulk hypersurface homologous to a subsystem XX (X=A,B,AB¯X=A,B,\overline{AB}). These hypersurfaces are then sewed as follows.

  • ΣA\Sigma_{A}: odd copies of bulk manifold are sewed cyclically and even copies are sewed to themselves,

  • ΣB\Sigma_{B}: even copies of bulk manifold are sewed anti-cyclically and odd copies are sewed to themselves,

  • ΣAB¯\Sigma_{\overline{AB}}: all copies are sewed pairwise.

We note that whereas the replica symmetry holds for the above construction in the boundary, it is explicitly broken while going from 2nn\mathbb{Z}_{2n}\to\mathbb{Z}_{n} in the bulk. This observation motivated the authors in [38] to consider instead the n\mathbb{Z}_{n} quotient of the original manifold given as

^2nA,B(nsym)=2nA,B(nsym)/n.\hat{\mathcal{B}}_{2n}^{A,B(\mathrm{nsym)}}=\mathcal{B}_{2n}^{A,B(\mathrm{nsym)}}/\mathbb{Z}_{n}\,. (72)

The quotienting described in eq. (72) modifies the corresponding bulk gravitational on-shell actions as

Igrav[2nA,B(nsym)]nIgrav[^2A,B(nsym)]=nIgrav[2AB,γA1(n),γB2(n)].\mathrm{I_{grav}}\left[\mathcal{B}_{2n}^{A,B(\mathrm{nsym)}}\right]\equiv n\mathrm{I_{grav}}\left[\hat{\mathcal{B}}_{2}^{A,B(\mathrm{nsym)}}\right]=n\mathrm{I_{grav}}\left[\mathcal{M}_{2}^{AB},\gamma_{A_{1}}^{(n)},\gamma_{B_{2}}^{(n)}\right]. (73)

Here Igrav[2AB,γA1(n),γB2(n)]\mathrm{I_{grav}}\left[\mathcal{M}_{2}^{AB},\gamma_{A_{1}}^{(n)},\gamma_{B_{2}}^{(n)}\right] is the on-shell bulk action for the alternative configuration of the quotiented bulk manifold ^2A,B(nsym)\hat{\mathcal{B}}_{2}^{A,B(\mathrm{nsym)}}, where it has conical defects along codimension two surfaces γA1(n)\gamma_{A_{1}}^{(n)} and γB2(n)\gamma_{B_{2}}^{(n)}, and 2AB\mathcal{M}_{2}^{AB} as its asymptotic boundary. Subscripts 1 and 2 denote the copy where the corresponding surface resides. We are now in a position to obtain an expression for the Rényi negativity by substituting eq. (73) into eq. (70) and inserting the result into eq. (68) as

log𝒩(2n)(even,nsym)=n(Igrav[2AB,γA1(n),γB2(n)]2Igrav[]).\log\mathcal{N}_{(2n)}^{(\mathrm{even,nsym})}=-n\left(\mathrm{I_{grav}}\left[\mathcal{M}_{2}^{AB},\gamma_{A_{1}}^{(n)},\gamma_{B_{2}}^{(n)}\right]-2\mathrm{I_{grav}}\left[\mathcal{B}\right]\right). (74)

Note that eq. (74) was derived in [38] on the basis of the assumption that the system comprising the subsystems AA, BB and AB¯\overline{AB} is in a tripartite pure state. A result for the bulk action away from n=1n=1 described in [54, 55] was employed in [39] to express the on-shell bulk action as

Igrav[2AB,γA1(n),γB2(n)]=2Igrav[]+𝒜(1/2)(γAB)4GN+(11n)𝒜(n)(γA)+𝒜(n)(γB)4GN.\mathrm{I_{grav}}\left[\mathcal{M}_{2}^{AB},\gamma_{A_{1}}^{(n)},\gamma_{B_{2}}^{(n)}\right]=2\mathrm{I_{grav}}\left[\mathcal{B}\right]+\frac{\mathcal{A}^{(1/2)}\left(\gamma_{AB}\right)}{4G_{N}}+\left(1-\frac{1}{n}\right)\frac{\mathcal{A}^{(n)}\left(\gamma_{A}\right)+\mathcal{A}^{(n)}\left(\gamma_{B}\right)}{4G_{N}}\,. (75)

Substituting eq. (75) into eq. (74) and inserting the result further into eq. (69) we arrive at the following expression for the holographic entanglement negativity

(A:B)\displaystyle\mathcal{E}(A:B) =18GN[𝒜(1/2)(γA)+𝒜(1/2)(γB)𝒜(1/2)(γAB)]\displaystyle=\frac{1}{8G_{N}}\left[\mathcal{A}^{(1/2)}\left(\gamma_{A}\right)+\mathcal{A}^{(1/2)}\left(\gamma_{B}\right)-\mathcal{A}^{(1/2)}\left(\gamma_{AB}\right)\right]
=12[S(1/2)(A)+S(1/2)(B)S(1/2)(AB)]\displaystyle=\frac{1}{2}\left[S^{(1/2)}(A)+S^{(1/2)}(B)-S^{(1/2)}(AB)\right]
=12I(1/2)(A:B).\displaystyle=\frac{1}{2}I^{(1/2)}(A:B)\,. (76)

Starting from eq. (76), in section A.3 we follow [39] to describe a plausible proof of the holographic construction for the entanglement negativity of a single interval as advanced in [27].

Holographic entanglement negativity for a single interval

In this section we describe a plausible derivation of the holographic proposal for the entanglement negativity of a single interval in the context of AdS3/CFT2AdS_{3}/CFT_{2}, following the work in [39]. To this end we establish an important relation below for the Rényi mutual information of order half. Consider a tripartite system ABCABC (consisting of the subsystems AA, BB and CC) in a pure state. We have

I(1/2)(A:BC)I(1/2)(A:C)\displaystyle I^{(1/2)}(A:BC)-I^{(1/2)}(A:C)
=[S(1/2)(A)+S(1/2)(BC)S(1/2)(ABC)][S(1/2)(A)+S(1/2)(C)S(1/2)(AC)]\displaystyle=\left[S^{(1/2)}(A)+S^{(1/2)}(BC)-S^{(1/2)}(ABC)\right]-\left[S^{(1/2)}(A)+S^{(1/2)}(C)-S^{(1/2)}(AC)\right]
=S(1/2)(BC)+S(1/2)(AC)S(1/2)(C)S(1/2)(ABC).\displaystyle=S^{(1/2)}(BC)+S^{(1/2)}(AC)-S^{(1/2)}(C)-S^{(1/2)}(ABC)\,. (77)

Since the Rényi entropy of a system in a pure state is zero, we have S(1/2)(ABC)=0S^{(1/2)}(ABC)=0. Furthermore, for a bipartite system in a pure state, the entanglement negativity is given by the Rényi entropy of order half, (A:BC)=S(1/2)(A)\mathcal{E}(A:BC)=S^{(1/2)}(A) and (BC:A)=S(1/2)(BC)\mathcal{E}(BC:A)=S^{(1/2)}(BC). Since (A:BC)=(BC:A)\mathcal{E}(A:BC)=\mathcal{E}(BC:A), we have S(1/2)(A)=S(1/2)(BC)S^{(1/2)}(A)=S^{(1/2)}(BC). Relabeling AA, BB and CC, we further obtain S(1/2)(B)=S(1/2)(AC)S^{(1/2)}(B)=S^{(1/2)}(AC) and S(1/2)(C)=S(1/2)(AB)S^{(1/2)}(C)=S^{(1/2)}(AB). These relations may be organized as follows

S(1/2)(BC)=S(1/2)(A),\displaystyle S^{(1/2)}(BC)=S^{(1/2)}(A)\,,\qquad\qquad S(1/2)(AC)=S(1/2)(B),\displaystyle S^{(1/2)}(AC)=S^{(1/2)}(B)\,, (78)
S(1/2)(C)=S(1/2)(AB),\displaystyle S^{(1/2)}(C)=S^{(1/2)}(AB)\,, S(1/2)(ABC)=0.\displaystyle S^{(1/2)}(ABC)=0\,.

Substituting the relations in eq. (78) into eq. (77) we obtain I(1/2)(A:BC)I(1/2)(A:C)=I(1/2)(A:B)I^{(1/2)}(A:BC)-I^{(1/2)}(A:C)=I^{(1/2)}(A:B), which may be re-expressed as

I(1/2)(A:BC)=I(1/2)(A:B)+I(1/2)(A:C).I^{(1/2)}(A:BC)=I^{(1/2)}(A:B)+I^{(1/2)}(A:C)\,. (79)

Having derived the relation in eq. (79), we are now in a position to show that the holographic construction in [27] may be deduced from eq. (76). The corresponding configuration may be briefly described as an interval AA with two auxiliary intervals B1B_{1} and B2B_{2} on either side. In the bipartite limit B1B2AcB_{1}B_{2}\to A^{c}, the system AB1B2AB_{1}B_{2} is in a pure state, and we may apply the result in eq. (79) to this system to obtain the following relation

limB1B2AcI(1/2)(A:B1B2)=limB1B2Ac[I(1/2)(A:B1)+I(1/2)(A:B2)].\lim_{B_{1}B_{2}\to A^{c}}I^{(1/2)}(A:B_{1}B_{2})=\lim_{B_{1}B_{2}\to A^{c}}\left[I^{(1/2)}(A:B_{1})+I^{(1/2)}(A:B_{2})\right]\,. (80)

The holographic entanglement negativity for this configuration may then be computed by substituting eq. (80) into eq. (76) as

=limB1B2Ac12[I(1/2)(A:B1)+I(1/2)(A:B2)].\mathcal{E}=\lim_{B_{1}B_{2}\to A^{c}}\left.\frac{1}{2}\middle[I^{(1/2)}(A:B_{1})+I^{(1/2)}(A:B_{2})\right]\,. (81)

Through the definition of Rényi mutual information, eq. (81) may be expressed as

=limB1B2Ac12[2S(1/2)(A)+S(1/2)(B1)+S(1/2)(B2)S(1/2)(AB1)S(1/2)(AB2)].\mathcal{E}=\lim_{B_{1}B_{2}\to A^{c}}\left.\frac{1}{2}\middle[2S^{(1/2)}(A)+S^{(1/2)}(B_{1})+S^{(1/2)}(B_{2})-S^{(1/2)}(A\cup B_{1})-S^{(1/2)}(A\cup B_{2})\right]. (82)

Since each of the intervals involved on the right hand side of eq. (82) is a single connected subsystem with spherical geometry, eq. (62) applies to each of them, and we have S(1/2)(X)=𝒳2S(X)=32S(X)S^{(1/2)}(X)=\mathcal{X}_{2}\,S(X)=\frac{3}{2}S(X), where XX represents each of the intervals. We may now re-express eq. (82) as

=limB1B2Ac34[2S(A)+S(B1)+S(B2)S(AB1)S(AB2)].\mathcal{E}=\lim_{B_{1}B_{2}\to A^{c}}\left.\frac{3}{4}\middle[2S(A)+S(B_{1})+S(B_{2})-S(A\cup B_{1})-S(A\cup B_{2})\right]. (83)

Utilizing the RT formula S(X)=X/4GNS(X)=\mathcal{L}_{X}/4G_{N} and eq. (83), we may finally arrive at the following expression for the holographic entanglement negativity for a single interval in the AdS3/CFT2AdS_{3}/CFT_{2} framework

=limBAc316GN[2A+B1+B2AB1AB2],\mathcal{E}=\lim_{B\to A^{c}}\left.\frac{3}{16G_{N}}\middle[2\mathcal{L}_{A}+\mathcal{L}_{B_{1}}+\mathcal{L}_{B_{2}}-\mathcal{L}_{A\cup B_{1}}-\mathcal{L}_{A\cup B_{2}}\right], (84)

where BB1B2B\equiv B_{1}\cup B_{2}. Eq. (84) may be readily identified as the holographic conjecture for the entanglement negativity of a single interval in the AdS3/CFT2AdS_{3}/CFT_{2} scenario as proposed in [27]. This concludes the review of the derivation of the holographic proposal for the entanglement negativity of a single interval.

Appendix B Non extremal and extremal AdSd+1AdS_{d+1}-RN

The numerical constants b0b_{0} and b1b_{1} appearing in eqs. (17) and (19) are given as follows

b0=2(d1)d2,b_{0}=\frac{2(d-1)}{d-2}\;, (85)

b1=dd2.b_{1}=\frac{d}{d-2}\;. (86)

The constants 𝒮0\mathcal{S}_{0} and 𝒮1\mathcal{S}_{1} appearing in eqs. (30), (37) and (49) are given as follows

𝒮0=2d2πd12Γ[d22(d1)](d1)Γ[12(d1)](Γ[d2(d1)]Γ[12(d1)])d2,\mathcal{S}_{0}=\frac{2^{d-2}\pi^{\frac{d-1}{2}}\Gamma\left[-\frac{d-2}{2(d-1)}\right]}{(d-1)\Gamma\left[\frac{1}{2(d-1)}\right]}\left(\frac{\Gamma\left[\frac{d}{2(d-1)}\right]}{\Gamma\left[\frac{1}{2(d-1)}\right]}\right)^{d-2}, (87)

𝒮1=Γ[12(d1)]d+12d1πd2Γ[d2(d1)]dΓ[12+1d1](Γ[1d1]Γ[d22(d1)]+21d1(d2)Γ[1+12(d1)]π(d+1)).\mathcal{S}_{1}=\frac{\Gamma\left[\frac{1}{2(d-1)}\right]^{d+1}2^{-d-1}\pi^{-\frac{d}{2}}}{\Gamma\left[\frac{d}{2(d-1)}\right]^{d}\Gamma\left[\frac{1}{2}+\frac{1}{d-1}\right]}\left(\frac{\Gamma\left[\frac{1}{d-1}\right]}{\Gamma\left[-\frac{d-2}{2(d-1)}\right]}+\frac{2^{\frac{1}{d-1}}(d-2)\Gamma\left[1+\frac{1}{2(d-1)}\right]}{\sqrt{\pi}(d+1)}\right). (88)

The numerical constants N0N_{0} and N1N_{1} appearing in eq. (31) are given as follows

N0\displaystyle N_{0} =2(πΓ[d22(d1)]2(d1)Γ[12(d1)])\displaystyle=2\left(\frac{\sqrt{\pi}~{}\Gamma\left[-\frac{d-2}{2(d-1)}\right]}{2(d-1)\Gamma\left[\frac{1}{2(d-1)}\right]}\right) (89)
+201𝑑x(1x2(d1)xd11b0xd+b1x2(d1)1xd11x2(d1)),\displaystyle+2\int_{0}^{1}dx\left(\frac{\sqrt{1-x^{2(d-1)}}}{x^{d-1}\sqrt{1-b_{0}x^{d}+b_{1}x^{2(d-1)}}}-\frac{1}{x^{d-1}\sqrt{1-x^{2(d-1)}}}\right),

N1=01𝑑x(x1x2(d1)1b0xd+b1x2(d1))(1xd21b0xd+b1x2(d1)).N_{1}=\int_{0}^{1}dx\left(\frac{x\sqrt{1-x^{2(d-1)}}}{\sqrt{1-b_{0}x^{d}+b_{1}x^{2(d-1)}}}\right)\left(\frac{1-x^{d-2}}{1-b_{0}x^{d}+b_{1}x^{2(d-1)}}\right). (90)

The function γd(μT)\gamma_{d}\left(\frac{\mu}{T}\right) appearing in eq. (38) is given as follows

γd(μT)=N(1)+d2(d2)16π2(d1)(μT)201𝑑x(x1x2(d1)1xd)(1xd21xd)+𝒪[(μT)4],\gamma_{d}\left(\frac{\mu}{T}\right)=N(1)+\frac{d^{2}(d-2)}{16\pi^{2}(d-1)}\left(\frac{\mu}{T}\right)^{2}\int_{0}^{1}dx\left(\frac{x\sqrt{1-x^{2(d-1)}}}{\sqrt{1-x^{d}}}\right)\left(\frac{1-x^{d-2}}{1-x^{d}}\right)+\mathcal{O}\left[\left(\frac{\mu}{T}\right)^{4}\right], (91)

where the expression for the numerical constant N(ε)N(\varepsilon) is given in eq. (92).

The numerical constant N(ε)N(\varepsilon) appearing in eqs. (50) and (91) is given as follows

N(ε)=2(πΓ[d22(d1)]2(d1)Γ[12(d1)])+201𝑑x(1x2(d1)xd1f(zHx)1xd11x2(d1)).N(\varepsilon)=2\left(\frac{\sqrt{\pi}\Gamma\left[-\frac{d-2}{2(d-1)}\right]}{2(d-1)\Gamma\left[\frac{1}{2(d-1)}\right]}\right)+2\int_{0}^{1}dx\left(\frac{\sqrt{1-x^{2(d-1)}}}{x^{d-1}\sqrt{f(z_{H}x)}}-\frac{1}{x^{d-1}\sqrt{1-x^{2(d-1)}}}\right). (92)

References