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.


I. 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, 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.
The combination of quantum information science, relativity theory and quantum field theory gives us a deeper understanding of quantum mechanics. It is necessary to understand quantum effects in the relativistic framework, because the world is essentially noninertial or/and curved.
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 Sec. II, we briefly recall the measures of GTN and GTE for tripartite quantum systems. In Sec. III, we introduce the quantization of Dirac fields in the background of Schwarzschild black hole. Sec. IV 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 section is devoted to the conclusion.

II. MEASURES OF GTN AND GTE
A. 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 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 ), are any three-dimensional 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 then the Svetlichny inequality fulfills for any state ρ that admits LHV model, where S is taken over all possible Svetlichny operators, i.e., the measurements are taken over all the directions of three-dimensional space.
Equivalently, if a tripartite state violates the Svetlichny inequality for some Svetlichny operator, then this state is genuinely nonlocal. For convenience, we introduce the Svetlichny value, i.e., the maximum expectation value of Svetlichny operator [8,9,47] S(ρ) = max S tr(Sρ) 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] S(ρ X ) = max{8 where N = n 1 − n 2 − n 3 + n 4 − m 4 + m 3 + m 2 − m 1 . 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 ρ = i p i |Ψ i Ψ i |. For three-qubit X states given by Eq.(5), the GTE is given by where

III. 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 Γ µ = can be written as where M denotes the mass of the black hole. For simplicity, we take , 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 II } are the number states for the particle outside the region and the antiparticle inside the region of the event horizon, respectively.

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 II and C II 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 II and C II inside the event horizon the inaccessible modes. Taking trace over the inaccessible modes on state Ψ AB I B II C I C II , 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 , 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 . 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.1(a) and (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.1(a)(α = 1 √ 2 ), 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.1(b)), 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.1(c) and (d). We can see from Fig.1(c) 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.1(d), 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). Setting cos 2 η = (1 + e − ω T ) −1 , sin 2 η = (1 + e ω T ) −1 and cos 2 ζ = α 2 , . 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.1 We can also make the similar discussions for other tripartite subsystems. Tracing over the modes B II and C I on the state Ψ AB I B II C I C II , we obtain the reduced density operator ρ AB I C II 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 II C I ) = S(ρ AB I C II ), C(ρ AB II C I ) = C(ρ AB I C II ).
Thus the analysis for tripartite system AB II C I is the same as for tripartite system AB I C II .
In Fig.2, we plot the Svetlichny value S(ρ AB I C II ) and the GTE C(ρ AB I C II ) as functions of the Hawking temperature T for different initial parameter α. It is shown that C(ρ AB I C II ) 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 II , 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 II , which is equivalent to an interaction between modes C I and C II . This interaction transfers some information from mode C I to mode C II . Therefore, the GTE between the modes A, B I and C II is established. Comparing Fig.2 (a) and (b), we find that the produced GTE C(ρ AB I C II ) depends on the initial accessible GTE in Eq. (14). Under given Hawking temperature, more initially accessible GTE will produce more C(ρ AB I C II ). In the limit of infinite Hawking temperature, the asymptotic value of C(ρ AB I C II ) is 0.5 in Fig.2(a) and 0.37 in Fig.2(b) respectively. The figure shows that S(ρ AB I C II ) is always smaller than 4 for any T , thus the physically inaccessible GTN between modes A, B I and C II cannot be produced. It also means that the produced inaccessible GTE C(ρ AB I C II ) is local. Similar analysis is also valid for the tripartite system of modes A, B II and C I . The different behaviors between S(ρ AB I C II ) and C(ρ AB I C II ) under the Hawking effect suggest that the information flows of different quantum resources inside and outside the event horizon of a black hole are completely different: The entanglement 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 II and C II . Tracing over the modes B I and C I on the state Ψ AB I B II C I C II , we obtain the reduced density operator ρ AB II C II as with the matrix elements given by

The corresponding Svetlichny value and GTE read
and respectively.
In Fig.3, we plot the Svetlichny value S(ρ AB II C II ) and the GTE C(ρ AB II C II ) 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 II : The Hawking effect can generate physically inaccessible GTE 3.5

4.
0.092 between the modes A, B II and C II , 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 II C II ) is. The produced GTE has the asymptotic value C(ρ AB II C II ) = α √ 1 − α 2 in the limit of infinite Hawking temperature, which is 0.5 for Fig.3(a) and 0.37 for Fig.3(b) 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 II . In the bases |000 , |100 , |010 , |001 , |101 , |111 , |110 , and |011 for A, B I and B II , the density operator ρ AB I B II 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 II , Thus, we just need to analyze the GTN and GTE for the tripartite system AB I B II .
In Fig.4 (a) and (b), we plot the Svetlichny value S(ρ AB I B II ) and GTE C(ρ AB I B II ) 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 II , but cannot produce the physically inaccessible GTN between them. The produced C(ρ AB I B II ) increases monotonically from zero and reaches the asymptotic value C(ρ AB I B II ) = α 2 for T → ∞. Comparing Fig.4 (a) and (b), we find that C(ρ AB I B II ) 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 II ).
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.4 (c)-(d), we plot the maximal Bell signal B(ρ B I B II ) and bipartite concurrence C(ρ B I B II ) between the modes B I and B II as functions of the Hawking temperature T for different initial parameter α. We find that B(ρ B I B II ) 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 II ) 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 II . In fact, from the calculation in Appendix B, we can find that all the reduced bipartite subsystems in state Ψ AB I B II C I C II have no nonlocality, meaning that nonlocality cannot pass through the event horizon of black hole. Except for C(ρ B I B II ) and C(ρ C I C II ), 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 II or C I C II , it of course can only produce quantum entanglement between modes B I and B II or C I and C II , and cannot produce any other kinds of bipartite entanglement.
From the discussions of Figs.1-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 II ), S(ρ AB II C I ), S(ρ AB II C II ), S(ρ AB I B II ), and S(ρ AC I C II ) 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, 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 II C II ) 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.
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 ij ) with t ij = Tr[ρσ i ⊗ σ j ]. Bell nonlocality can be witnessed by the maximum violation of CHSH inequality.
For two-qubit X-state, with all elements ρ ij being real, the three eigenvalues corresponding to the matrix U(ρ) = T T ρ T ρ are, As Z 1 is greater than Z 2 , so the maximal Bell signal reads with B 1 = 2 √ Z 1 + Z 2 and B 2 = 2 √ Z 1 + Z 3 [55,56].
Appendix B: Pairwise quantum nonlocality and entanglement in state Ψ AB I B II C I C II We now study the pairwise bipartite quantum nonlocality and entanglement in the state Ψ AB I B II C I C II . By tracing over the irrelevant modes on state Ψ AB I B II C I C II , we obtain all the pairwise density operators as