Genuine tripartite nonlocality and entanglement in curved spacetime

We study the genuine tripartite nonlocality (GTN) and the genuine tripartite entanglement (GTE) of Dirac fields in the background of a Schwarzschild black hole. We find that the Hawking radiation degrades both the physically accessible GTN and the physically accessible GTE. The former suffers from “sudden death” at some critical Hawking temperature, and the latter approaches to the nonzero asymptotic value in the limit of infinite Hawking temperature. We also find that the Hawking effect cannot generate the physically inaccessible GTN, but can generate the physically inaccessible GTE for fermion fields in curved spacetime. These results show that on the one hand the GTN cannot pass through the event horizon of black hole, but the GTE do can, and on the other hand the surviving physically accessible GTE and the generated physically inaccessible GTE for fermions in curved spacetime are all not nonlocal. Some monogamy relations between the physically accessible GTE and the physically inaccessible GTE are found.


Introduction
The concept of quantum nonlocality was firstly proposed by Einstein, Podolsky and Rosen in 1935 in their famous EPR paradox [1]. Afterwards, Bell established the so-called Bell inequality for conveniently judging the existence of quantum nonlocality [2]. Two quantum systems that admit Bell inequality cannot be regarded as nonlocal, even if they are far apart in space. Contrarily, when the Bell inequality is violated, we say that the two quantum systems are inseparable or have the quantum nonlocality. Quantum nonlocality is a kind of quantum behaviors, which denies the local hidden variable (LHV) model [3][4][5][6][7][8][9]. Quantum nonlocality is a very important quantum resource and has been applied in different fields, such as device-independent quantum computation, a e-mail: hszeng@hunnu.edu.cn (corresponding author) communication complexity, quantum cryptography and randomness generation [10][11][12][13][14][15]. In the case of tripartite quantum systems, a more valuable kind of resource is involved which is called GTN. Svetlichny introduced a kind of GTN and found the so-called Svetlichny inequality to detect its existence [16]. Quantum nonlocality and entanglement are inextricably linked. Quantum nonlocality originates from quantum entanglement, but quantum entanglement does not imply quantum nonlocality.
In this work, we study the properties of GTN and GTE of Dirac fields in the background of a Schwarzschild black hole. Assume that Alice, Bob and Charlie initially share a Greenberger-Horne-Zeilinger-like state. Alice is a Kruskal observer who stays stationarily at an asymptotically flat region, while Bob and Charlie are Schwarzschild observers who hover near the event horizon of the black hole. In addition, there are two imagined observers, anti-Bob and anti-Charlie, in the interior of the event horizon. The Hawking effect would make the information tunnelling from the exterior to the interior of the event horizon, such that a correlated state that involves the above five observers is established. Since the interior of the black hole is causally disconnected from the exterior and observers cannot access to the interior of the event horizon, we thus call the information (including GTN and GTE) that distributes completely in the outside of the event horizon the physically accessible. Otherwise, it is called the physically inaccessible. The main end of this paper is to study the influence of Hawking effect on the physically accessible GTN and GTE, the production of the physically inaccessible GTN and GTE, and the monogamy relationship between the physically accessible and inaccessible information.
The paper is organized as follows. In Sect. 2, we briefly recall the measures of GTN and GTE for tripartite quantum systems. In Sect. 3, we introduce the quantization of Dirac fields in the background of Schwarzschild black hole. Section 4 is the main contribution of our work, where the decay of the physically accessible information and the production of the physically inaccessible information, as well as the monogamy relationship between them, under the influence of Hawking effect, are studied. Finally, the last Sect. 5 is devoted to the conclusion.

Measure of GTN
Firstly, we briefly review the concept about GTN. Nonlocality in tripartite systems has been considered as the manifestation of genuine tripartite correlations. Generally, local tripartite correlations shared by Alice, Bob and Charlie can be written as where {λ; p λ } is the probability distribution of some hidden variable that controls the outputs a, b, c ∈ {0, 1} of the local measurements performed by Alice, Bob and Charlie on the two-valued variables x, y, z ∈ {0, 1}, and P λ (a|x) is the conditional probability for obtaining the output a when the measurement setting is x and λ. If the tripartite correlation cannot be written as the form of Eq. (1), then the tripartite system is said to have GTN. In 1987, Svetlichny proposed a hybrid local-nonlocal form of correlation to measure GTN [16]. A tripartite correlation under the definition of Svetlichny is called to be local if it admits the following local LHV model P(a, b, c|x, y, z) = λ p λ P λ (a|x)P λ (b, c|y, z) where λ p λ + μ p μ + ν p ν = 1. This kind of correlation is regarded as Svetlichny local, otherwise it is Svetlichny nonlocal [44][45][46]. Assuming that Alice, Bob and Charlie share a state ρ of some three-qubit system. Alice performs measurements on the observable A = a · σ and A = a · σ , where a = (a 1 , a 2 , a 3 ), a = (a 1 , a 2 , a 3 ) ∈ R 3 are any threedimensional unit vectors, and σ = (σ 1 , σ 2 , σ 3 ) is the vector of Pauli matrices. Similar measurements are performed by Bob and Charlie, respectively, on the observable B, B and observable C, C . Introducing Svetlichny operator to detect the GTN of a given three-qubit state ρ. For the three-qubit X state with density matrix in the orthogonal basis {|0, 0, 0 , |0, 0, 1 , . . . , |1, 1, 1 }, the Svetlichny value can be simply given by [9] where N = n 1 − n 2 − n 3 + n 4 − m 4 + m 3 + m 2 − m 1 .

Measure of GTE
GTE can be defined by its opposite of biseparability. We call a tripartite pure state | is biseparable, if it has a bipartition of the form | = | A ⊗ | B , where | A and | B are the monomeric or bipartite pure states. Obviously, a biseparable pure state has at least one pure marginal. If the tripartite state | is not biseparable with respect to any of its bipartition, then it is called GTE. One can then define the so-called genuine tripartite concurrence C(| ) = min χ i ∈χ 2[1 − Tr(ρ 2 for describing the degrees of the GTE for the pure state | , where χ = {A i |B i } denotes the set of all possible bipartitions of the tripartite system, and ρ A χ i is the reduced density operator of system A corresponding to the bipartition χ i [42,49]. The GTE for a mixed state ρ can be obtained by a convex roof construction where the infimum takes over all possible decompositions For three-qubit X states given by Eq. (5), the GTE is given by where ν i = 4 j =i √ n j m j [43].

Quantization of Dirac fields in a Schwarzschild black hole
The Dirac equation under a general background spacetime can be written as [48] [γ a e a μ (∂ μ + μ )] = 0, where γ a are the Dirac matrices, the four-vectors e a μ is the inverse of the tetrad e a μ , and μ = 1 8 [γ a , γ b ]e a ν e bν;μ are the spin connection coefficients. The metric of the Schwarzschild black hole can be written as where M denotes the mass of the black hole. For simplicity, we takeh, G, c and k as unity in this paper. Solving the Dirac equation of Eq. (9) near the event horizon of black hole, a set of positive frequency outgoing solutions for the inside and outside regions of the event horizon can be obtained [20,24,50] where φ(r ) denotes four-component Dirac spinor, ω is a monochromatic frequency, and u = t − r * with the tortoise coordinate r * = r + 2M ln r −2M 2M . The modes + k,in and + k,out are usually called Schwarzschild modes. According to future-directed timelike Killing vector under each region, particles and antiparticles will be classified.
Making an analytic continuation for Eqs. (11) and (12) according to suggestion of Damour and Ruffini, we obtain a complete basis of positive energy modes, i.e., the Kruskal modes [51]. Then, we can use Schwarzschild mode and Kruskal mode to expand the Dirac field, respectively, leading to the Bogoliubov transformations between annihilation operator and creation operator under the Schwarzschild and Kruskal coordinates [52,53]. After properly normalizing the state vector, the vacuum state and excited state of the Kruskal particle in the single-mode approximation are given by where T = 1 8π M is the Hawking temperature, {|n I } and {|n I I } are the number states for the particle outside the region and the antiparticle inside the region of the event horizon, respectively.

Evolution of the GTN and GTE in Schwarzschild black hole
Consider a Greenberger-Horne-Zeilinger-like (GHZ-like) state of the Dirac fields shared by Alice, Bob and Charlie in the asymptotically flat region where α is the state parameter that runs from 0 to 1. Now, we assume that Alice still stays at an asymptotically flat region, while Bob and Charlie hover outside the event horizon of the black hole, then we can rewrite Eq. (14) in terms of Kruskal modes for Alice and Schwarzschild modes for Bob and Charlie as Generally, this is a 5-partite entangled state consisted by subsystems: subsystem A observed by Alice, subsystems B I and C I observed by Bob and Charlie outside the event horizon of black hole, and subsystems B I I and C I I observed by anti-Bob and anti-Charlie inside the event horizon, respectively. Since the interior region of black hole is causally disconnected from the exterior region, and Alice, Bob and Charlie cannot access the modes inside the event horizon, we thus call the modes B I and C I outside the event horizon the accessible modes, and the modes B I I and C I I inside the event horizon the inaccessible modes. Taking trace over the inaccessible modes on state AB I B I I C I C I I , we obtain the reduced density operator ρ AB I C I where the matrix elements are written by According to Eqs. (6) and (8), we obtain the Svetlichny value and GTE for state ρ AB I C I , (17) and respectively. Obviously, the Svetlichny value and GTE depend not only on the state parameter α, but also on the Hawking temperature T , meaning that the Hawking radiation of black hole will affect the physically accessible GTN and GTE between Alice, Bob and Charlie. In Fig. 1, we plot the Svetlichny value S(ρ AB I C I ) and the GTE C(ρ AB I C I ) as functions of the Hawking temperature T for different initial parameter α. We find that S(ρ AB I C I ) is larger than 4 at first and then smaller than 4 with the increase of Hawking temperature T . The critical Hawking temperature for S(ρ AB I . This implies that the thermal noise introduced by Hawking temperature destroys the physically accessible GTN between Alice, Bob and Charlie, and takes place "sudden death" at the critical temperature T c . However, the physically accessible GTE C(ρ AB I C I ) is a monotonic decreasing function of T , and has the asymptotic value α √ 1 − α 2 in the limit of infinite Hawking temperature. The GTE C(ρ AB I C I ) never takes place "sudden death" in the finite Hawking temperature. These results suggest that the GTE in the initial state of Eq. (14) can be distinguished into two different parts: nonlocal and local. The nonlocal GTE is destroyed completely by the Hawking effect after the temperature T > T c , and finally only the local GTE (whole or partial) is preserved. At this stage, no quantum information tasks based on GTN can work, but tasks based on GTE can still work. In other words, GTE is more suitable for relativistic quantum information tasks than GTN.
By comparing Fig. 1a, b, we find that both the critical temperature T c for the sudden death of GTN and the asymptotic value of C(ρ AB I C I ) in the infinite temperature depend on the initial GTE in Eq. (14), i.e., parameter α. The more the initial GTE is, the longer for the death time of GTN is (i.e., the larger the T c is ), and the larger for the asymptotic value of C(ρ AB I C I ) is. In Fig. 1a ), the critical Hawking temperature (about 1.13) and the asymptotic value (0.5) of C(ρ AB I C I ) are the maximal. For α = 1 √ 6 ( Fig. 1b), they reduce to about 0.34 and 0.37 respectively.
To further inspect the behaviors of GTN and GTE of the tripartite subsystem AB I C I in the regions of the initial parameter α ∈ (0, 1/ √ 2) and α ∈ (1/ √ 2, 1), we plot Fig. 1c, d. We can see from Fig. 1c that the GTE C(ρ AB I C I ) is symmetrical with respect to α 2 = 1/2 in the parameter region α 2 ∈ [0, 1]. In fact, this symmetry can also be verified via the fact that Eq. (18) is covariant under the exchange α 2 ←→ 1 − α 2 . From Fig. 1d, we see that the Svetlichny value S(ρ AB I C I ) is obviously asymmetrical with respect to α 2 = 1/2. However, the GTN (i.e., the part of S(ρ AB I C I ) ≥ 4) is symmetrical with respect to α 2 = 1/2. This symmetry also can be verified analytically via Eq. (17).
. Thus the second term in the right hand side of Eq. (17) becomes as 4| cos 2 ζ cos 2 (2η) + sin 2 ζ | ≤ 4. This means that, in the inspection of the GTN of the tripartite subsystem AB I C I , the second term in the right hand side of Eq. (17) can be ignored and only the first term need to be considered. Therefore, the GTN is symmetrical with respect to α 2 = 1/2. Note that both the GTE and GTN in Fig. 1c, d decrease when Hawking temperature increases, which are consistent with the results from Fig. 1a, b. Naturally, in the region α ∈ (0, 1/ √ 2), GTN and GTE change slowly with the increase of the α; in the region α ∈ (1/ √ 2, 1), GTN and GTE change steeply with the increase of the α.
We can also make the similar discussions for other tripartite subsystems. Tracing over the modes B I I and C I on the state AB I B I I C I C I I , we obtain the reduced density operator ρ AB I C I I as  with the matrix elements given by The Svetlichny value and the GTE for this state are given by and respectively. According to the exchange symmetry for Bob and Charlie, we can get S(ρ AB I I C I ) = S(ρ AB I C I I ), C(ρ AB I I C I ) = C(ρ AB I C I I ). (22) Thus the analysis for tripartite system AB I I C I is the same as for tripartite system AB I C I I .
In Fig. 2, we plot the Svetlichny value S(ρ AB I C I I ) and the GTE C(ρ AB I C I I ) as functions of the Hawking temperature T for different initial parameter α. It is shown that C(ρ AB I C I I ) increases from zero and approaches to the asymptotic value α √ 1 − α 2 in the infinite Hawking temperature. This means that the Hawking effect can generate physically inaccessible GTE between modes A, B I and C I I , even though they are separated by the event horizon of black hole. Physically, it can be regarded as a kind of entanglement transfer: Initially, there is GTE between modes A, B I and C I . Lately, the Hawking effect produce entanglement between the modes C I and C I I , which is equivalent to an interaction between modes C I and C I I . This interaction transfers some information from mode C I to mode C I I . Therefore, the GTE between the modes A, B I and C I I is established. Comparing Fig. 2a, b, we find that the produced GTE C(ρ AB I C I I ) depends on the initial accessible GTE in Eq. (14). Under given Hawking temperature, more initially accessible GTE will produce more C(ρ AB I C I I ). In the limit of infinite Hawking temperature, the asymptotic value of C(ρ AB I C I I ) is 0.5 in Fig. 2a and 0.37 in Fig. 2b respectively. The figure shows that S(ρ AB I C I I ) is always smaller than 4 for any T , thus the physically inaccessible GTN between modes A, B I and C I I cannot be produced. It also means that the produced inaccessible GTE C(ρ AB I C I I ) is local. Similar analysis is also valid for can pass through the event horizon of black hole, while the nonlocality cannot. Now, we discuss the physically inaccessible GTN and GTE between the modes A, B I I and C I I . Tracing over the modes B I and C I on the state AB I B I I C I C I I , we obtain the reduced density operator ρ AB I I C I I as with the matrix elements given by

S(ρ AB I I C I I
and

C(ρ AB I I C I I
respectively. In Fig. 3, we plot the Svetlichny value S(ρ AB I I C I I ) and the GTE C(ρ AB I I C I I ) as functions of the Hawking temperature T for different initial parameter α. We find the similar result as observed in the tripartite system of AB I C I I : the Hawking effect can generate physically inaccessible GTE between the modes A, B I I and C I I , but cannot generate the physically inaccessible GTN, i.e., entanglement can pass through the event horizon of black hole, and nonlocality cannot. The more the initial GTE in Eq. (14) is, the more the produced C(ρ AB I I C I I ) is. The produced GTE has the asymptotic value C(ρ AB I I C I I ) = α √ 1 − α 2 in the limit of infinite Hawking temperature, which is 0.5 for Fig. 3a and 0.37 for Fig. 3b respectively. The mechanism for the production of this GTE is also the result of entanglement transfer.
Finally, we investigate the physically inaccessible GTN and GTE for the system ρ AB I B I I . In the bases |000 , |100 , |010 , |001 , |101 , |111 , |110 , and |011 for A, B I and B I I , the density operator ρ AB I B I I has its matrix expression with the matrix elements given by The corresponding Svetlichny value and the GTE read and respectively. According to the exchange symmetry for Bob and Charlie, we can get the Svetlichny value and GTE between modes A, C I and C I I ,

S(ρ AC I C I I ) = S(ρ AB I B I I ), C(ρ AC I C I I ) = C(ρ AB I B I I ).
(29) Thus, we just need to analyze the GTN and GTE for the tripartite system AB I B I I . In Fig. 4a, b, we plot the Svetlichny value S(ρ AB I B I I ) and GTE C(ρ AB I B I I ) between Alice, Bob and anti-Bob as functions of the Hawking temperature T for different initial parameter α. It is shown that Hawking effect can also produce physically inaccessible GTE between modes A, B I and B I I , but cannot produce the physically inaccessible GTN between them. The produced C(ρ AB I B I I ) increases monotonically from zero and reaches the asymptotic value C(ρ AB I B I I ) = α 2 for T → ∞. Comparing Fig. 4a, b, we find that C(ρ AB I B I I ) also depends on the initially accessible GTE in Eq. (14). Under given Hawking temperature, more initially accessible GTE can produce more C(ρ AB I B I I ).
Besides GTN and GTE for tripartite systems, we can also study the bipartite nonlocality and entanglement for the considered system under the influence of Hawking effect. The bipartite nonlocality and entanglement may be described by the CHSH inequality and concurrence respectively for qubit systems. We review these concepts in the Appendix A and present all the pairwise bipartite nonlocality and entanglement for our considered system in the appendix B. In Fig.  4c, d, we plot the maximal Bell signal B(ρ B I B I I ) and bipartite concurrence C(ρ B I B I I ) between the modes B I and B I I as functions of the Hawking temperature T for different initial parameter α. We find that B(ρ B I B I I ) is always less than 2, meaning that Hawking effect cannot generate Bell nonlocality between Bob and anti-Bob. However, Hawking effect can generate entanglement between Bob and anti-Bob, i.e., entanglement C(ρ B I B I I ) can pass through the event horizon of black hole. Actually, this is just the physical nature of the Hawking radiation-produce entangled pairs of particle and antiparticle between the causally disconnected regions. The same analysis is also valid for the reduced state ρ C I C I I . In fact, from the calculation in Appendix B, we can find that all the reduced bipartite subsystems in state AB I B I I C I C I I have no nonlocality, meaning that nonlocality cannot pass through the event horizon of black hole. Except for C(ρ B I B I I ) and C(ρ C I C I I ), there is no pairs of entanglement in any other reduced bipartite subsystems. This result is easy to understand: Initially, there is no entanglement between modes A and B, modes A and C, modes B and C. The Hawking radiation is essentially a local operation on B I B I I or C I C I I , it of course can only produce quantum entanglement between modes B I and B I I or C I and C I I , and cannot produce any other kinds of bipartite entanglement.
From the discussions of Figs. 1, 2, 3 and 4, we find that the physically accessible GTN, i.e., S(ρ AB I C I ) − 4 reduces with Hawking temperature and suffers from a "sudden death" at some critical Hawking temperature, but the physically inaccessible GTN is never generated by the Hawking effect, i.e., S(ρ AB I C I I ), S(ρ AB I I C I ), S(ρ AB I I C I I ), S(ρ AB I B I I ), and S(ρ AC I C I I ) are all less than 4. This result suggests that the GTN can not be redistributed. However, the GTE behaves differently. When the physically accessible GTE reduces with Hawking temperature, at the same time, the physically inaccessible GTE is generated. This result implies that the GTE may be redistributable through the Hawking effect. To manifest this inference, we try to find some monogamy relation for the GTE. Through careful inspection, we find three monogamy relations between the physically accessible GTE and the physically inaccessible GTE, (32) where 2α √ 1 − α 2 is the initial GTE in state of Eq. (14). These monogamy relations reflect the restrictions in the redistribution of entanglement from physically accessible to physically inaccessible patterns. Especially the Eq. (30) shows that the total sum of the physically accessible GTE C(ρ AB I C I ) and the physically inaccessible GTE C(ρ AB I I C I I ) is equal to the initial GTE. The monogamy relations are important for understanding the transfer of quantum information in relativistic spacetime.
Besides the above monogamy relations, we also find that the physically accessible entanglement fulfills the following Coffman-Kundu-Wootters monogamy inequality where (i, j, k) denote all the permutations of the three modes A, B I , C I . These Coffman-Kundu-Wootters monogamy inequalities reflect the distribution of the physically accessible entanglement in the environment of Schwarzschild black hole.

Conclusions
The effect of Hawking radiation on the GTN and GTE for Dirac fields in Schwarzschild spacetime has been investigated. It has been shown that Hawking effect degrades both the physically accessible GTN and the physically accessible GTE, where the former takes place "sudden death" at some critical Hawking temperature, and the latter approaches to the nonzero asymptotic value in the infinite Hawking temperature. This means that on the one hand the surviving physically accessible GTE is not nonlocal, and on the other hand the GTE has more resistance to the Hawking noise than GTN. Further investigation has demonstrated that the GTN is not redistributable, but GTE can be redistributed through Hawking effect. With the growth of the Hawking temperature, the physically accessible GTN decreases, but no physically inaccessible GTN is generated. The GTE however behaves differently: With the loss of the physically accessible GTE, the physically inaccessible GTE is generated continually through Hawking effect. Further, the physically accessible GTE and the physically inaccessible GTE fulfil some monogamy relations. All these phenomena suggest that the GTE is redistributable. We can regard the redistribution of entanglement as a kind of phenomenon of information tunnelling, i.e. the flow of quantum entanglement can pass through the event horizon of black hole, but the flow of quantum nonlocality can not. This result has been demonstrated also by the pairwise bipartite nonlocality and entanglement in the underlying system considered in this paper.
Note that the authors in reference [34] studied the tripartite entanglement in environment of Schwarzschild black hole used the measure of π -tangle. They found that the physically accessible π -tangle decreases with Hawking temperature and approaches to a nonzero asymptotic value for the infinite Hawking temperature. Similar Coffman-Kundu-Wootters monogamy inequalities to Eq. (33) for π -tangle was found. These results agree with ours.
Quantum entanglement and nonlocality are the important manifestations of quantum correlation. They have potential applications in various fields of science. We expect our research can present helps for these applications and enriches the theory of the relativistic quantum information science.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data included in this study are available upon request by contact with the corresponding author Hao-Sheng Zeng.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: CHSH inequality and concurrence for bipartite systems
In this appendix, we review the concepts of CHSH inequality and concurrence for bipartite systems. The CHSH inequality is an important tool for judging quantum nonlocality. It is considered to be both necessary and sufficient conditions of the separability for bipartite pure states. The key element for the CHSH inequality is the Bell operator defined by [54] where a, a , b and b are unit vectors in R 3 , and σ = (σ 1 , σ 2 , σ 3 ) is the vector of Pauli matrices. For any bipartite mixed state ρ, the well-known CHSH inequality can be expressed as CHSH inequality holds for any state ρ that admits local hidden variable model, and the violation of it implies the Bell nonlocality of the underlying state. In practical applications, we need to find the maximal Bell signal B(ρ), which can be equivalently expressed for two-qubit systems as where Z i and Z j are the two largest eigenvalues of U (ρ) = T T ρ T ρ , and the correlation matrix is defined by T = (t i j ) with t i j = Tr[ρσ i ⊗ σ j ]. Bell nonlocality can be witnessed by the maximum violation of CHSH inequality.
For two-qubit X-state, where ρ i j is the element of density matrix ρ X [43].

Appendix B: Pairwise quantum nonlocality and entanglement in state AB I B I I C I C I I
We now study the pairwise bipartite quantum nonlocality and entanglement in the state AB I B I I C I C I I . By tracing over the irrelevant modes on state AB I B I I C I C I I , we obtain all the pairwise density operators as