Quantum entanglement for continuous variables sharing in an expanding spacetime

Detecting the structure of spacetime with quantum technologies has always been one of the frontier topics of relativistic quantum information. Here, we analytically study the generation and redistribution of Gaussian entanglement of the scalar fields in an expanding spacetime. We consider a two-mode squeezed state via a Gaussian amplification channel that corresponds to the time-evolution of the state from the asymptotic past to the asymptotic future. Therefore, the dynamical entanglement of the Gaussian state in an expanding universe encodes historical information about the underlying spacetime structure, suggesting a promising application in observational cosmology. We find that quantum entanglement is more sensitive to the expansion rate than the expansion volume. According to the analysis of quantum entanglement, choosing the particles with the smaller momentum and the optimal mass is a better way to extract information about the expanding universe. These results can guide the simulation of the expanding universe in quantum systems.


I. INTRODUCTION
Quantum entanglement, predicted by Schrödinger in 1935, plays an important role in quantum information theory.It is considered to be a resource for quantum information tasks, such as computational task, quantum teleportation, quantum communication, quantum control, and quantum simulations [1][2][3][4].Recently, the influence of gravitational effects on quantum entanglement has received considerable attention.Therefore, many efforts have been expended on the study of quantum entanglement of the field modes in a noninertial frame, in the black hole, and in de Sitter space .Besides, many papers have studied the influences of different spacetime on quantum entanglement and explored the properties of spacetime using quantum entanglement.It is clear that these investigations not only contribute to understanding the key questions about quantum information and exploring the structure of spacetime, but also play a central role in the study of the information paradox of the black hole and the entanglement entropy [38][39][40][41].
In this work, we study the redistribution of bipartite entanglement for continuous variables of the scalar fields in the background of an expanding universe.We initially consider a twomode squeezed Gaussian state shared by Alice and Bob in the asymptotic past.A vacuum state in the asymptotic past evolves into a thermal state through the expansion of the universe in the asymptotic future [42][43][44][45].In a quantum information scenario, such a process of the expansion of the universe can be described as a Gaussian channel acting on a two-mode squeezed Gaussian state.There are two reasons why I choose the two-mode squeezed Gaussian state: firstly, it is a typical continuous variable entangled state, which approximates to an arbitrarily good extent the EPR pair; secondly, the state can be produced in the laboratory and exploited for any current realization of quantum information with continuous variables [46,47].In the standard cosmology model, the dynamical quantum entanglement for continuous variables acts as witnesses in the history of the universe from the era of big-bang nucleosynthesis to the era of large-scale structure formation and can encode historical information about the underlying spacetime structure, which suggests a promising application in observational cosmology.Thus, studying the behavior of quantum entanglement is crucial to understanding the history of our universe and its fate.
A two-mode squeezed Gaussian state in the asymptotic past becomes a four-mode Gaussian state in the asymptotic future.Therefore, we evaluate not only the initial bipartite entanglement as influenced by the expansion of the universe, but also, remarkably, the multipartite entanglement, which arises from two bosonic modes and two antibosonic modes.The result shows that the bipartite entanglement initially prepared in the asymptotic past is exactly redistributed into four-partite correlations in the asymptotic future, as a consequence of the monogamy constraints on quantum entanglement distribution.On the other hand, we understand the properties of the expanding spacetime through quantum entanglement for continuous variables.Our results may guide the cosmological observations and the simulation of the expanding spacetime in quantum systems [48][49][50].
The paper is organized as follows.In Sec.II, we briefly introduce the measures of bipartite entanglement for continuous variables and the Coffman-Kundu-Wootters inequality.In Sec.III, we discuss how the expansion of the universe is described by a Gaussian channel.In Sec.IV, we study the redistribution of the two-mode Gaussian entanglement and the generated 1 → 3 entanglement under the influence of the expansion of the underlying spacetime.The last section is devoted to a brief conclusion.

II. QUANTIFYING ENTANGLEMENT FOR CONTINUOUS VARIABLES BY THE R ÉNYI-2 ENTROPY
The set of Gaussian states by the definition is the set of states with quasi-probability distributions and Gaussian characteristic functions in the quantum phase space.The properties of the Gaussian state are entirely determined by the first and second statistical moments of the quadrature operators.As the first moments can be arbitrarily adjusted by the local unitary operations, which keep all informationally relevant properties (such as entropy or any measure of correlations) invariant, we can adjust them to be zero.Therefore, the second moments for describing Gaussian state become the unique elements.Based on our research, we initially consider a twomode Gaussian state ρ AB shared by Alice and Bob.We define a vector of quadrature operators as R = ( X1 , P1 , ..., Xn , Pn ) T , which satisfies the canonical commutation relations [ Rk , Rl ] = iΩ kl , with the symplectic matrix Ω = 0 1 −1 0 ⊕n [51][52][53].The elements of a covariance matrix (CM) σ AB can be defined as σ ij = Tr { Ri , Rj } + ρ AB .The CM σ AB of the Gaussian state can be put into a block form For a physical Gaussian state, its CM σ AB must satisfy the uncertainty relation sure of quantum information.Rényi-α entropies are defined as [54,55] S α (ρ) = 1 1 − α ln tr(ρ α ).
In the limit α → 1, it reduces to the von Neumann entropy.For the case α = 2, the special Rényi-2 entropy can be calculated very easily where σ is the CM of the Gaussian state with density matrix ρ.For arbitrary Gaussian states, Rényi-α entropies fulfill the strong subadditivity inequality.This allows us to define relevant Gaussian measures of information and entanglement quantities under a unified approach.
For a N-mode bipartite Gaussian state with CM σ AB in Eq.( 1), the Rényi-2 quantum entanglement E(σ AB ), quantifying the quantum entanglement between Alice and Bob, can be defined as For a pure Gaussian state, the minimum is saturated by σ AB = γ AB , so that E(σ AB ) = 1 2 ln(det γ A ), where γ A is the reduced CM for Alice.For a mixed state, the minimization is over pure N-mode Gaussian states with CM γ AB which satisfies 0 < γ AB ≤ σ AB and det γ AB = 1.
[For two real symmetric matrices M and N, M ≥ N means that the matrix M − N has all nonnegative eigenvalues.]For a special class of two-mode Gaussian states, the closed formulae E can be obtained by the same procedure of entanglement of formation [57].
In multipartite systems, quantifying entanglement is generally very involved.Unlike classical correlations, quantum entanglement is monogamous, which means that it cannot be freely shared among multiple subsystems of a composite system.So far, this fundamental constraint on quantum entanglement sharing has been mathematically demonstrated not only for arbitrary systems of qubits within the discrete-variable scenario but also for a special case of two qubits and an infinitedimensional system and for all N-mode Gaussian states within the continuous variables scenario.
In the general case of a quantum state distributed among N parties, the monogamy constraint can be presented in the form of the Coffman-Kundu-Wootters inequality [56] where the multipartite system has N subsystems S k (k = 1, ..., N), each owned by a corresponding party, and E is a proper quantification of bipartite entanglement.The left-hand side of inequality ( 6), E S i |(S 1 ...S i−1 S i+1 ...S N ) , can quantify the bipartite entanglement between a probe subsystem S i and the remaining N − 1 subsystems.The right-hand side of inequality ( 6), N j =i E S i |S j , can quantify the total bipartite entanglements between subsystem S i and each one of the other subsystems S j =i in the reduced states.The nonnegative difference between these two quantum entanglements, which is minimized over all choices of the probe subsystem, is known as the residual multipartite entanglement.It quantifies the quantum entanglements not encoded in pairwise form, so it includes all manifestations of genuine N-partite entanglement.

III. THE EXPANSION OF THE UNIVERSE DESCRIBED BY GAUSSIAN CHANNEL
Let us start with a 1+1 dimensional Robertson-Walker expanding universe with the metric , where a(t) is the scale factor.With the conformal time η relating to the cosmological time t by η = t 0 dτ a(τ ) , the metric of the Robertson-Walker expanding universe can be rewritten as [42][43][44][45] Here, the conformal scale factor takes the form where the parameters ǫ and σ characterize the volume and the rapidity of the expansion, respectively.It is obvious that the spacetime is flat in the distant past and the far future corresponding to the metric ds 2 = dη 2 − dx 2 when η → −∞ and the metric ds 2 = (1 + 2ǫ)(dη 2 − dx 2 ) when η → +∞, respectively.Therefore, the timelike Killing vector and the particle content of the field are defined in these two limits.
A real scalar field Φ(x, η) in the Robertson-Walker expanding spacetime obeys the Klein-Gordon equation where Having solved the Klein-Gordon equation at the limits of η → ±∞, we obtain a set of modes u in in the distant past ("in" region) and a set of modes u out in the far future ("out" region).Using the inner product, the Bogoliubov transformations between the modes u in k and u out k take the form where the Bogoliubov coefficients are given by with Γ being the gamma function, For convenience, we define and easily obtain where |β k | 2 equals the average number of particles created at "out" mode k.Hence, θ 2 k → 0 means that the average number of particles of the mode k is vanishing, and θ 2 k → 1 means that the average number of particles of the mode k approaches infinity.

The annihilation and creation operators satisfy
where b in,k and b † in,k are the bosonic annihilation and creation operators acting on the states in the asymptotic past, b out,k and b † out,k are the the bosonic annihilation and creation operators acting on the states in the asymptotic future, and b out,−k and b † out,−k are the the antibosonic annihilation and creation operators, respectively.We use b in,k |0 k in = 0 to find the relation between the "in" vacuum state and the "out" vacuum state.If substituting b in,k with Eq.( 14), we obtain According to the the normalization condition, the "in" vacuum state can be expressed in the asymptotic future as where where ] is a two-mode squeezing operator.The squeezing parameter r k is defined as cosh(r k ) = |α k |.It is worth emphasizing that the squeezing operator U k is a Gaussian operation, which preserves the Gaussianity from the input states.Therefore, Eq. (18) shows that the expansion of a Robertson-Walker spacetime can be described by a Gaussian (a bosonic amplification) channel.In the phase space, the two-mode squeezing operator U k corresponds to the symplectic transformation where I 2 denotes the unity matrix in 2 × 2 space, and Z 2 denotes the third Pauli matrix.

IV. THE EFFECT OF THE EXPANDING UNIVERSE ON GAUSSIAN ENTANGLEMENT
In this paper, we initially consider a pure two-mode squeezed Gaussian state shared by Alice and Bob in the distant past.It has covariance matrix [60] where s is the squeezing parameter.If Alice and Bob undergo the expanding universe associating to the symplectic transformation in Eq.( 19), the initial two-mode squeezed state σ in AB becomes a four-mode Gaussian state in the asymptotic future.In other words, the initial two-mode squeezed state in the asymptotic past is transformed, via the expansion of a Robertson-Walker spacetime, to a four-mode Gaussian state in the asymptotic future.Therefore, a complete description of the quantum system involves four modes: the bosonic mode A described by Alice; the bosonic mode B described by Bob; the antibosonic mode Ā described by anti-Alice; the antibosonic mode B described by anti-Bob.The covariance matrix describing the complete system thus becomes [60]  where S A, Ā and S B, B given by Eq.( 19) are the phase space representation of the two-mode squeezing operation, and I Ā B denotes the 4 × 4 identity matrix.

A. Bipartite Gaussian entanglement
Because Alice and Bob cannot detect the antibosonic modes, we should take the trace over the modes Ā and B. By performing this operation on Eq.( 21), we obtain the reduced state between the modes A and B Employing Eq.( 5), we can obtain an analytic expression of Gaussian entanglement From Eq.( 23), we can see that quantum entanglement E(σ out AB ) depends not only on the squeezing parameter s, but also on the volume ǫ and the rapidity σ of the expanding universe, meaning that quantum entanglement encodes historical information about the cosmological parameters.
In Fig. 1 (a), we plot the Gaussian entanglement E(σ out AB ) between Alice and Bob as a function of the expansion rate σ for different momentums k. and a monotonically increasing function of the momentum k.We find that as the expansion rate σ increases, quantum entanglement decreases to an asymptotic value dependent on the momentum k.At the limit of k → ∞, we obtain where ln[cosh(2s)] is the initial entanglement in Eq. (20).Therefore, we cannot extract any information about the spacetime at the limit of k → ∞.Fig. 1 (b) shows that quantum entanglement E(σ out AB ) decreases with the increase of the expansion volume ǫ.We can also see that quantum entanglement between the modes A and B first decreases to the minimum value and then increases to the initial value with the growth of the mass m.Therefore, quantum entanglement of the bosonic fields with the larger mass m is insensitive to the expansion volume ǫ.This means that choosing the bosons with the appropriate mass makes quantum entanglement more sensitive to the cosmological parameters.
We calculate quantum entanglements in all possible bipartite divisions of the four-mode quantum system to explore the distribution of Gaussian entanglement in an expanding spacetime.
Firstly, tracing over the modes in A and B, we obtain the covariance matrix σ out Ā B between the modes Ā and B σ out Secondly, taking the trace over the modes B and Ā, we get the covariance matrix between Alice and anti-Bob σ out ĀB = σ out A B can be directly obtained by a simple calculation.According to Eq.( 5), we find E(σ out Ā B ) = E(σ out A B ) = 0, meaning that the expansion of the universe cannot generate quantum entanglement between anti-Alice and anti-Bob (or Alice and anti-Bob).
Finally, we fix our eyes on quantum entanglement between the modes A and Ā. Tracing over the modes B and B, we obtain the covariance matrix σ out A Ā for Alice and anti-Alice Using Eqs.( 5) and ( 26), we obtain an analytic expression of quantum entanglement between the modes A and Ā Unlike the Gaussian entanglement E(σ out AB ), the Gaussian entanglement E(σ out A Ā) does not depend on the initial squeezing parameter s.What we need to notice is that In Fig. 2 (a), we plot the Gaussian entanglement E(σ out A Ā) as functions of the expansion rate σ and the momentum k.Fig. 2 (b) shows how the mass m and the expansion volume ǫ influence the Gaussian entanglement E(σ out A Ā). Fig. 2 (a) shows that the Gaussian entanglement E(σ out A Ā) increases with the expansion rate σ, which means that the expansion of the universe can generate quantum entanglement between Alice and anti-Alice.We find that the particles with the smaller momentum k help the expansion of the universe to generate greater quantum entanglement.From Fig. 2 (b), we can see that quantum entanglement E(σ out A Ā) increases with the growth of the expansion volume ǫ.When we choose the particles with the suitable mass m, the expansion of the underlying spacetime has a more pronounced effect on quantum entanglement.From Fig. 2, we can also see that quantum entanglement is sensitive to the expansion rate but not to the expansion volume.
Combining Fig. 1 and Fig. 2, we come to three conclusions: (i) with the increase of the expansion rate and the expansion volume, the initial Gaussian entanglement between Alice and Bob decreases, and at the same time, the Gaussian entanglement between Alice and anti-Alice (or Bob and anti-Bob) increases, which means that the expansion of the universe redistributes the initial entanglement; (ii) quantum entanglement is more sensitive to the expansion rate than the expansion volume; (iii) the information about the expanding universe can be better extracted by choosing the particles with the smaller momentum and the appropriate mass.Studying the properties of quantum entanglement can help us to understand the history of the expanding universe.
Therefore, these results can help us better simulate the expanding spacetime and the production of cosmological particles with quantum systems in a laboratory setting [48][49][50].

B. Generation of a four-mode Gaussian entanglement
The four-mode Gaussian state σ out AB Ā B of Eq.( 21) is entirely inseparable, meaning that it contains genuine multipartite Gaussian entanglement distributed among all the four parties involved.Let us now compute the bipartite entanglements in the 1 → 3 partitions of the Gaussian state σ out AB Ā B .Employing Eq.( 5), the 1 → 3 entanglements are found to be For any nonzero value of s, σ and ǫ, each single party is in an entangled Gaussian state with the block of the remaining three parties, with respect to all possible global splitting of the modes.
In Fig. 3    entanglement E(σ out Ā|AB B ) increases from zero to the maximum and then reduces to zero.This implies that, for the bosons with the larger mass m, the effect of the expansion of the universe on the 1 → 3 entanglement is not apparent.
The residual multipartite entanglement from the monogamy inequality in Eq.( 6) is an entanglement monotone under Gaussian local operations and classical communication for a four-mode pure state σ out AB Ā B .The residual entanglement of the Gaussian state σ out AB Ā B is defined as We verify that the second quantity can achieve the minimum.Therefore, we obtain an analytic expression for the residual entanglement E res (σ out AB Ā B ) From Eq.( 31), we find that the residual entanglement E res (σ out AB Ā B ) is always greater than or equal to zero for any parameters, which proves that the Coffman-Kundu-Wootters inequality in Eq.( 6) is still correct in an expanding spacetime.In addition, it precisely quantifies the multipartite quantum correlations that cannot be stored in the bipartite form.
In Fig. 4 (a), we plot the residual entanglement E res (σ out AB Ā B ) as functions of the expansion rate σ and the momentum k.Fig. 4 (b) shows how the mass m and the expansion volume ǫ influence the residual entanglement E res (σ out AB Ā B ). From Fig. 4, we can see that the residual entanglement monotonically increases with the increase of the expansion rate and the expansion volume.This indicates that the expansion of the universe can generate the residual entanglement.We can also see that the residual entanglement is more influenced by the expansion rate than the expansion volume.In addition, the residual entanglement is more sensitive to the cosmological parameters when we choose the particles with the smaller momentum and the suitable mass in an expanding spacetime.

V. CONCLUTIONS
We have studied the redistribution of quantum entanglement for continuous variables in an expanding spacetime.We consider four modes: the bosonic mode A observed by Alice; the bosonic mode B observed by Bob; the antibosonic mode Ā observed by anti-Alice; the antibosonic mode B observed by anti-Bob.We get a phase space description for a quantum state evolution under the influence of the expansion of a Robertson-Walker spacetime.When the quantum state evolves from the asymptotic past to the asymptotic future, the dynamical entanglement contains historical information concerning the expanding spacetime.We find that quantum entanglement is more sensitive to the expansion rate than the expansion volume.We show the redistribution of the initial entanglement: quantum entanglement between the modes A and B decreases with the growth of the expansion rate and the expansion volume; at the same time, quantum entanglement between the modes A and Ā (or B and B ) can be generated by the expansion of the underlying spacetime.This means that the loss of quantum entanglement can be interpreted as a redistribution of the initial entanglement into multipartite quantum correlations.We also find that the 1 → 3 entanglement increases with the expansion rate and the expansion volume.The residual entanglement is always greater than or equal to zero, which directly proves that the Coffman-Kundu-Wootters inequality is true in an expanding spacetime.According to the analysis of quantum entanglement, choosing the particles with the smaller momentum and the optimal mass is a better way to extract information about the expanding universe.These results can guide the simulation of the expanding universe in different quantum systems [48][49][50].

FIG. 1 :
FIG. 1: (a) Quantum entanglement E(σ out AB ) between the modes A and B as a function of the expansion rate σ for different momentums k, with the fixed values of m = ǫ = s = 1.(b) Quantum entanglement E(σ out AB ) as a function of the mass m for different expansion volumes ǫ, with the fixed values of k = σ = s = 1.
Fig.1 (b) shows how the mass m influences the Gaussian entanglement E(σ out AB ) for different expansion volumes ǫ.From Fig.1 (a), we can see that quantum entanglement E(σ out AB ) is a monotonically decreasing function of the expansion rate σ

FIG. 2 :
FIG. 2: (a) Quantum entanglement E(σ out A Ā) between the mode A and the mode Ā as functions of the expansion rate σ and the momentum k for the fixed values of m = ǫ = s = 1.(b) Quantum entanglement E(σ out A Ā) as functions of the mass m and the expansion volume ǫ for the fixed values of k = σ = s = 1.
(a) and (d), we plot the 1 → 3 bipartite entanglement as a function of the expansion rate σ for different squeezing parameters s.We find that the 1 → 3 bipartite entanglement monotonically increases with the increase of the expansion rate σ.In addition, as the squeezing parameter s increases, the influence of the expansion rate σ on the generated entanglement E(σ out Ā|AB B ) is more obvious.The 1 → 3 bipartite entanglement is plotted in Fig.3 (b) and (e) as a function of the expansion volume ǫ for different momentums k.It is shown that the 1 → 3 bipartite entanglement increases with the increase of the expansion volume ǫ.However, the influence of the expansion volume on the 1 → 3 entanglement is less than that of the expansion rate.From Fig.3 (b) and (e), we can see that, for a fixed expansion volume ǫ, the 1 → 3 entanglement decreases with the increase of the momentum k.

Fig. 3 (
Fig.3 (c) and (f) show how the mass m influences the 1 → 3 entanglement.From Fig.3 (c), we find that quantum entanglement E(σ out A|B Ā B ) first increases from the initial entanglement to the maximum and then reduces to the initial entanglement.From Fig.3(f), we can see that quantum

FIG. 3 : 1 →
FIG. 3: 1 → 3 entanglement (a) and (d) as a function of the expansion rate σ for different squeezing parameters s. 1 → 3 entanglement (b) and (e) as a function of the expansion volume ǫ for different momentums k. 1 → 3 entanglement (c) and (f) as a function of the mass m.The other parameters are set to 1.

FIG. 4 :
FIG. 4: (a) The residual entanglement E res (σ out AB Ā B ) as functions of the expansion rate σ and the momentum k for the fixed values of m = ǫ = s = 1.(b) The residual entanglement E res (σ out AB Ā B ) as functions of the mass m and the expansion volume ǫ for the fixed values of k = σ = s = 1.
n k denotes the boson number, and n −k denotes the antiboson number.This means that an initial vacuum state |0 k in evolves into a two-mode squeezed state in the asymptotic future.By rotating the squeezing angle and giving up the phase angle, we obtain