Next Article in Journal
Sleep Information Gathering Protocol Using CoAP for Sleep Care
Next Article in Special Issue
Performance Improvement of Plug-and-Play Dual-Phase-Modulated Quantum Key Distribution by Using a Noiseless Amplifier
Previous Article in Journal
Assessing the Role of Strategic Choice on Organizational Performance by Jacquemin–Berry Entropy Index
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Harvesting Large Scale Entanglement in de Sitter Space with Multiple Detectors

Department of Physics, Graduate School of Science, Nagoya University, Chikusa, Nagoya 464-8602, Japan
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(9), 449; https://doi.org/10.3390/e19090449
Submission received: 4 August 2017 / Revised: 23 August 2017 / Accepted: 23 August 2017 / Published: 28 August 2017
(This article belongs to the Special Issue Entropy and Information in the Foundation of Quantum Physics)

Abstract

:
We consider entanglement harvesting in de Sitter space using a model of multiple qubit detectors. We obtain the formula of the entanglement negativity for this system. Applying the obtained formula, we find that it is possible to access to the entanglement on the super horizon scale if a sufficiently large number of detectors are prepared. This result indicates the effect of the multipartite entanglement is crucial for detection of large scale entanglement in de Sitter space.

1. Introduction

All structure in the Universe can be traced back to primordial quantum fluctuations generated during an inflationary phase of the very early universe. To apprehend the history and origin of our universe, it is essential to understand the mechanism and the nature of these fluctuations with quantum origins. To investigate the quantum property of primordial quantum fluctuation, entanglement is a key concept to distinguish quantum nature from the classical one. Thus, it is an important task to analyze the detail of the entanglement of quantum fluctuations generated by inflation. In this direction, detection of the entanglement of the quantum scalar field using a pair of particle detectors were considered [1,2,3,4,5]. The entanglement of the scalar field can be probed by evaluating the entanglement between these two detectors interacting with the field; an initially non-entangled pair of detectors can evolve to be an entangled state through the interaction with the quantum field. As the entanglement cannot be created by local operations, this implies that the entanglement of the quantum field is transferred to the pair of detectors.
In de Sitter spacetime, although detectors can probe the entanglement of the scalar field on a scale smaller than the Hubble horizon, they cannot catch entanglement beyond the Hubble horizon scale [1,2,4,5]. Qualitatively, a similar result is shown for entanglement between two spatial regions defined via averaging (coarse graining of the scalar field) [6,7], and connection to the quantum–classical transition in the early universe was discussed. On the other hand, the result of a recent lattice calculation [8] that simulates the quantum scalar field near the continuous limit shows entanglement between two spatial regions persists even beyond the Hubble horizon scale and the entanglement negativity does not vanish. This discrepancy may come from the efficiency of entanglement detection using a pair of detectors; we expect that the efficiency of detection increases if the degrees of freedom of detectors grows. Hence, we consider a multiple detectors system and investigate how the maximum possible distance of entanglement detection depends on the number of detectors.
In this paper, we consider m + n qubit detectors and investigate detectability of bipartite entanglement on the super horizon scale in de Sitter space. We obtain negativity of this system analytically in the lowest non-trivial order of perturbation with respect to the coupling constant. Using this result, we discuss the possibility of entanglement harvesting beyond the Hubble horizon scale. The structure of the paper is as follows. In Section 2, we introduce our model of the multiple detectors system and the master equation for detectors state. In Section 3, we review the two qubits detectors case and show the maximum possible distance of entanglement detection cannot exceeds the Hubble horizon scale. In Section 4, we evaluate the negativity of m + n detectors case. In Section 5, we discuss the relation to the monogamy inequality. Section 6 is devoted to a summary. We use the unit in which = c = 1 throughout the paper.

2. Model and Strategy

We consider the following Hamiltonian for m + n qubit detectors interacting with a scalar field ϕ (see Figure 1):
H tot = H S + H int + H ϕ = α = 1 m + n ω 2 σ 3 ( α ) + g α = 1 m + n ( σ + ( α ) + σ ( α ) ) ϕ ( x α ) + H ϕ ,
where ω represents energy difference between two internal levels | 0 , | 1 and g is a coupling constant between detectors and the scalar field. The tensor products of the operators are defined as
σ j ( 1 ) = σ j 1 1 , σ j ( 2 ) = 1 σ j 1 1 , , σ j ( m + n ) = 1 1 1 σ j ,
with σ 3 = | 1 1 | | 0 0 | , σ + = | 1 0 | , and σ = | 0 1 | .
As a tool of our analysis, we introduce a master equation for the detectors’ state. Regarding the scalar field as an environment, we obtain the reduced density matrix for detectors by tracing out the scalar field degrees of freedom. Provided that the time scale of the environment is shorter than the detector’s time scale and the coupling is weak g 1 , the state of detectors ρ can be shown to obey the following Gorni-Kossakowski-Lindblad-Sudarshan (GKLS) type master equation [9,10,11,12,13]:
ρ t + i [ H eff , ρ ] = L [ ρ ] ,
H eff = H S i 2 α 1 , α 2 = 1 m + n j 1 , j 2 = ± H j 1 j 2 ( α 1 α 2 ) σ j 1 ( α 1 ) σ j 2 ( α 2 ) ,
L [ ρ ] = 1 2 α 1 , α 2 = 1 m + n j 1 , j 2 = ± C j 1 j 2 ( α 1 α 2 ) 2 σ j 2 ( α 2 ) ρ σ j 1 ( α 1 ) σ j 1 ( α 1 ) σ j 2 ( α 2 ) ρ ρ σ j 1 ( α 1 ) σ j 2 ( α 2 ) ,
where H eff is the effective Hamiltonian of the detectors system with quantum corrections. The coefficients H j 1 j 2 ( α 1 α 2 ) , C j 1 j 2 ( α 1 α 2 ) are expressed using the Wightman function of the quantum field as
C j 1 j 2 ( α 1 α 2 ) = 2 g 2 π σ e ω 2 σ 2 e i ω σ j + + d x d y e 1 σ 2 x σ + i 2 ω σ 2 j + 2 1 σ 2 y i 2 ω σ 2 j 2 D ( r c , t + x , y ) ,
H j 1 j 2 ( α 1 α 2 ) = 2 g 2 π σ e ω 2 σ 2 e i ω σ j + + d x d y sgn ( y ) e 1 σ 2 x σ + i 2 ω σ 2 j + 2 1 σ 2 y i 2 ω σ 2 j 2 D ( r c , t + x , y ) ,
with j ± = j 1 ± j 2 and where D ( r c , x , y ) is the Wightman function of the scalar field
D ( r c , x , y ) = ϕ ( t 1 , x α 1 ) ϕ ( t 2 , x α 2 ) , x = ( t 1 + t 2 ) / 2 , y = ( t 1 t 2 ) / 2 ,
where r c = | x α 1 x α 2 | denotes comoving distance between two detectors. The parameter σ , in C j 1 j 2 ( α 1 α 2 ) , H j 1 j 2 ( α 1 α 2 ) , specifies the time scale of coarse graining which is necessary to derive the GKLS master Equation (3). In the limit of σ →∞, this master equation reduces to that with the rotating wave approximation that neglects transition via energy non-conserving processes. The GKLS master equation preserves the trace and complete positivity.
The master Equation (3) was applied to the two detectors system in de Sitter space for the purpose of investigating long time evolution of negativity beyond the Hubble time scale [5]. In the analysis of the present paper, we concentrate on short time evolution from the initial state and do not solve this equation exactly. In such a restricted situation, as we will show in Section 3, prediction by the master Equation (3) coincides with that of the detectors with finite interaction time that is usually imposed by introducing an appropriate switching function of detector.
To examine detection of entanglement of the quantum field, we consider a solution of the master Equation (3) with a separable initial condition and judge the separability of the detectors state after evolution. For Δ t = t t 0 1 / ω , the solution with the initial state ρ 0 = ρ ( t 0 ) is
ρ ( t ) = ρ 0 + Δ t i [ H eff , ρ 0 ] + L [ ρ 0 ] .
For the initial separable state of detectors ρ 0 = | 0 0 0 0 | , H eff and ρ 0 commutes each other and the state after evolution can be written as
ρ ( t ) = ρ 0 + Δ t L [ ρ 0 ] .
Thus the entanglement of the state ρ is completely determined only by the operator L [ ρ 0 ] . From now on, we examine the state (10).
We divide m + n detectors to two groups and assign labels of detectors as
A : α 1 , , m , B : α m + 1 , m + n .
For simplicity of analysis, we assume the distance between two groups is r, and the distance between two detectors belonging to the same group is d (see Figure 1). We denote possible states of m + n detectors after evolution as follows:
  • | 0 : 0 = | 0 , 0 : 0 , 0 : ground state.
  • | i : 0 : i-th detector in group A is excited.
  • | 0 : i : i-th detector in group B is excited.
  • | i 1 : i 2 : i 1 -th detector in group A and i 2 -th detector in group B are excited.
  • | i 1 i 2 : 0 : i 1 -th and i 2 -th detectors ( i 1 i 2 ) in group A are excited.
  • | 0 : i 1 i 2 : i 1 -th and i 2 -th detectors ( i 1 i 2 ) in group B are excited.
  • C ( 0 ) C ( i 1 i 2 ) i 1 = i 2 , C ( r ) C ( i 1 i 2 ) i 1 A , i 2 B .
  • C ( d ) C ( i 1 i 2 ) i 1 i 2 A = C ( i 1 i 2 ) i 1 i 2 B .
As we will see, the following coefficients [5] in the master Equation (3) are necessary to calculate the negativity,
C + ( r ) = 2 g 2 e ω 2 σ 2 π + d y e 1 σ 2 ( y + i ω σ 2 ) 2 D ( r c , t + σ , y ) ,
C + + ( r ) = 2 g 2 e ω 2 σ 2 π 0 + d y e y 2 σ 2 D ( r c , t + σ + i ω σ 2 , y ) ,
where the time coarse graining parameter σ must satisfy H σ < 1 to guarantee the assumption to derive the master Equation (3). The parameter σ corresponds to the width of switching function in analysis of the standard particle detector model. By applying the saddle point approximation, which is correct for parameters with 1 / ω < σ < 1 / H , σ < r , these coefficients can be evaluated as
C + ( r ) = 2 g 2 e ω 2 σ 2 σ D ( r c , t + σ , i ω σ 2 ) , C + + ( r ) = 2 g 2 e ω 2 σ 2 e 2 i ω σ σ D ( r c , t + σ + i ω σ 2 , 0 ) .
For the massless conformal scalar field in de Sitter spacetime with a spatially flat time slices, these coefficients are given by [4,5]
| C + + ( r ) | = g 2 σ H 2 e ω 2 σ 2 2 π 2 1 H 2 r 2 , C + ( r ) = g 2 σ H 2 e ω 2 σ 2 2 π 2 1 H 2 r 2 + 4 sin 2 θ ,
where θ H σ 2 ω < π , and r = e N r c e N H 1 denotes the physical separation between detectors at e-folding time of inflation N = H × ( t 0 + σ ) . For the massless minimal scalar field,
| C + + ( r ) | = g 2 σ H 2 e ω 2 σ 2 2 π 2 e 2 i θ H 2 r 2 + 1 1 2 Ei ( i k 0 e N r ) + Ei ( i k 0 e N r ) , C + ( r ) = g 2 σ H 2 e ω 2 σ 2 2 π 2 [ 1 H 2 r 2 + 4 sin 2 θ 1 2 Ei ( i k 0 e N ( r 2 i H 1 sin θ ) ) + Ei ( i k 0 e N ( r 2 i H 1 sin θ ) ) + 1 ] ,
where k 0 is the infrared cutoff corresponding to the comoving size of the inflating universe k 0 = H and Ei ( x ) = x d y y e y is the exponential integral.
As a warm up, we first review the 1 + 1 detectors case, which is often adopted as a model of entanglement harvesting in numerous situations.

3. Negativity for 1 + 1 Detectors System (Two Qubits Case)

For the initial separable state of detectors (We adopt the basis { | 1 : 1 , | 1 : 0 , | 0 : 1 , | 0 : 0 } , which is descending order of states in binary numbering.)
ρ 0 = | 0 : 0 0 : 0 | = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 ,
the state of detectors (10) becomes
ρ = 0 0 0 C + + ( r ) Δ t 0 C + ( 0 ) Δ t C + ( r ) Δ t 0 0 C + ( r ) Δ t C + ( 0 ) Δ t 0 C ( r ) Δ t 0 0 1 2 C + ( 0 ) Δ t .
To quantify entanglement between detectors, we introduce the entanglement negativity [14,15]
E N = λ i < 0 | λ i | ,
where λ i are eigenvalues of partially transposed state ρ PT , which is defined by transposing components belonging to group B only. For the state (18),
ρ PT = 0 0 0 C + ( r ) Δ t 0 C + ( 0 ) Δ t C + + ( r ) Δ t 0 0 C ( r ) Δ t C + ( 0 ) Δ t 0 C + ( r ) Δ t 0 0 1 2 C + ( 0 ) Δ t .
Eigenvalues of this state are
λ = 1 2 Δ t C + ( 0 ) , Δ t C + ( 0 ) ± C + + ( r ) .
We are considering ω Δ t 1 with the weak coupling limit, which implies the coefficient | C | = O ( g 2 ) 1 . Hence, the only eigenvalue which can become negative is Δ t C + ( 0 ) C + + ( r ) and the negativity is given by E N = Δ t max | C + + ( r ) | C + ( 0 ) , 0 . As the sign of negativity does not depends on Δ t , the initial separable state can become entangled one instantly. From now on in this paper, we designate the following quantity as the negativity
E N = | C + + ( r ) | C + ( 0 ) .
By this definition of the negativity, the state is entangled for E N > 0 . For the two qubits case ( m = n = 1 ), E N < 0 implies the state is separable. However, for m + n 3 , we can say nothing about separability from the condition E N < 0 . Comparing the formula of the negativity (22) with that derived previously [4] using a switching function of detectors, Equation (22) exactly coincides with previous one and we can confirm that the time coarse graining parameter σ has a meaning of a width of switching function of detectors.
Figure 2 shows parameter regions where entanglement detection is possible. For parameters ( r , θ ) belonging to the left side regions of each lines, the negativity is positive and a pair of detectors can harvest the entanglement of the scalar field.
For the massless conformal scalar field, the maximum distance of entanglement detection is 2 H 1 and for the massless minimal scalar field, that distance is H 1 . In this paper, we regard the Hubble length H 1 as the horizon radius and define the super horizon scale as r > 2 H 1 . Following this definition of “super horizon scale”, a pair of detectors cannot access entanglement on the super horizon scale for both type of scalar fields.

4. Negativity for m + n Detectors System

Let us consider entanglement harvesting using m + n detectors. The initial state of detectors is assumed to be
ρ 0 = | 0 : 0 0 : 0 | ,
and the state after time evolution is obtained by evaluating L [ ρ 0 ] in Equation (10):
ρ = ρ 0 + Δ t 2 α 1 , α 2 = 1 m + n j 1 , j 2 = ± C j 1 j 2 ( α 1 α 2 ) 2 σ j 2 ( α 2 ) ρ 0 σ j 1 ( α 1 ) σ j 1 ( α 1 ) σ j 2 ( α 2 ) ρ 0 ρ 0 σ j 1 ( α 1 ) σ j 2 ( α 2 ) = ρ 0 + Δ t [ i 1 , i 2 A C + ( i 1 i 2 ) | i 1 : 0 i 2 : 0 | + i 1 A i 2 B C + ( i 1 i 2 ) | i 1 : 0 0 : i 2 | + i 1 B i 2 A C + ( i 1 i 2 ) | 0 : i 1 i 2 : 0 | + i 1 , i 2 B C + ( i 1 i 2 ) | 0 : i 1 0 : i 2 | ] Δ t 2 [ 2 i 1 A i 2 B C + + ( i 1 i 2 ) | i 1 : i 2 0 : 0 | + i 1 i 2 A C + + ( i 1 i 2 ) | i 1 i 2 : 0 0 : 0 | + i 1 i 2 B C + + ( i 1 i 2 ) | 0 : i 1 i 2 0 : 0 | + ( m + n ) C + ( 0 ) | 0 : 0 0 : 0 | ] Δ t 2 [ 2 i 1 A i 2 B C ( i 1 i 2 ) | 0 : 0 i 1 : i 2 | + i 1 i 2 A C ( i 1 i 2 ) | 0 : 0 i 1 i 2 : 0 | + i 1 i 2 B C ( i 1 i 2 ) | 0 : 0 0 : i 1 i 2 | + ( m + n ) C + ( 0 ) | 0 : 0 0 : 0 | ] .
After partial transposition of the state with respect to the group B,
ρ PT = ρ 0 + Δ t [ i 1 , i 2 A C + ( i 1 i 2 ) | i 1 : 0 i 2 : 0 | + i 1 A i 2 B C + ( i 1 i 2 ) | i 1 : i 2 0 : 0 | + i 1 B i 2 A C + ( i 1 i 2 ) | 0 : 0 i 2 : i 1 | + i 1 , i 2 B C + ( i 1 i 2 ) | 0 : i 2 0 : i 1 | ] Δ t 2 [ 2 i 1 A i 2 B C + + ( i 1 i 2 ) | i 1 : 0 0 : i 2 | + i 1 i 2 A C + + ( i 1 i 2 ) | i 1 i 2 : 0 0 : 0 | + i 1 i 2 B C + + ( i 1 i 2 ) | 0 : 0 0 : i 1 i 2 | + ( m + n ) C + ( 0 ) | 0 : 0 0 : 0 | ] Δ t 2 [ 2 i 1 A i 2 B C ( i 1 i 2 ) | 0 : i 2 i 1 : 0 | + i 1 i 2 A C ( i 1 i 2 ) | 0 : 0 i 1 i 2 : 0 | + i 1 i 2 B C ( i 1 i 2 ) | 0 : i 1 i 2 0 : 0 | + ( m + n ) C + ( 0 ) | 0 : 0 0 : 0 | ] = Δ t [ i 1 , i 2 A C + ( i 1 i 2 ) | i 1 : 0 i 2 : 0 | + i 1 , i 2 B C + ( i 1 i 2 ) | 0 : i 2 0 : i 1 | i 1 A i 2 B C + + ( i 1 i 2 ) | i 1 : 0 0 : i 2 | i 1 A i 2 B C ( i 1 i 2 ) | 0 : i 2 i 1 : 0 | ] + 1 Δ t ( m + n ) C + ( 0 ) | 0 : 0 0 : 0 | + Δ t i 1 A i 2 B C + ( i 1 i 2 ) | i 1 : i 2 0 : 0 | + i 1 B i 2 A C + ( i 1 i 2 ) | 0 : 0 i 2 : i 1 | Δ t 2 i 1 i 2 A C + + ( i 1 i 2 ) | i 1 i 2 : 0 0 : 0 | + i 1 i 2 B C + + ( i 1 i 2 ) | 0 : 0 0 : i 1 i 2 | Δ t 2 i 1 i 2 A C ( i 1 i 2 ) | 0 : 0 i 1 i 2 : 0 | + i 1 i 2 B C ( i 1 i 2 ) | 0 : i 1 i 2 0 : 0 | ρ 1 + ρ 2 ,
where ρ 1 denotes a part of ρ PT spanned by the basis { | i 1 : 0 , | 0 : i 1 } and ρ 2 is a part of ρ PT spanned by the basis { | 0 : 0 , | i 1 : i 2 , | i 1 i 2 : 0 , | 0 : i 1 i 2 } . As these two sets of basis are orthogonal to each other, a matrix representation of ρ PT with these basis has a block diagonal structure. Hence, to obtain eigenvalues of ρ PT , the only thing we have to do is to consider eigenvalues of ρ 1 and ρ 2 separately.

4.1. Eigenvalues of ρ 1

The eigenvalue equation is
ρ 1 | λ = λ | λ .
From the configuration of detectors we are considering, we can assume the following form of the eigenvector
| λ = α i A | i : 0 + β i B | 0 : i ,
where α , β are coefficients to be determined. By applying ρ 1 ,
ρ 1 | λ = α i 1 , i 2 A C + ( i 1 i 2 ) | i 1 : 0 α i 1 A i 2 B C ( i 1 i 2 ) | 0 : i 2 + β i 1 , i 2 B C + ( i 1 i 2 ) | 0 : i 2 β i 1 A i 2 B C + + ( i 1 i 2 ) | i 1 : 0 = i 1 A α i 2 A C + ( i 1 i 2 ) β i 2 B C + + ( i 1 i 2 ) | i 1 : 0 + i 1 B α i 2 A C ( i 1 i 2 ) + β i 2 B C + ( i 1 i 2 ) | 0 : i 1 ,
and the eigenvalue equation is reduced to be
λ α = α i 2 A C + ( i 1 i 2 ) β i 2 B C + + ( i 1 i 2 ) with i 1 A ,
λ β = α i 2 A C ( i 1 i 2 ) + β i 2 B C + ( i 1 i 2 ) with i 1 B .
Thus,
λ α = α C + ( 0 ) + α ( m 1 ) C + ( d ) β n C + + ( r ) , λ β = β C + ( 0 ) + β ( n 1 ) C + ( d ) α m C ( r ) .
By eliminating α , β , we obtain
λ C + ( 0 ) ( m 1 ) C + ( d ) λ C + ( 0 ) ( n 1 ) C + ( d ) = m n C + + ( r ) 2 ,
and eigenvalues are
λ = C + ( 0 ) + m + n 2 1 C + ( d ) ± m n C + + ( r ) 2 + m n 2 2 C + ( d ) 2 1 / 2 .
These quantities are O ( g 2 ) .

4.2. Eigenvalues of ρ 2

We will show that eigenvalues of ρ 2 are positive up to O ( g 2 ) and they do not contribute to the negativity in the present order of calculation. We assume the form of the eigenvector as
| λ = | 0 : 0 + i 1 A i 2 B α i 1 i 2 | i 1 : i 2 + i 1 i 2 A β i 1 i 2 | i 1 i 2 : 0 + i 1 i 2 B γ i 1 i 2 | 0 : i 1 i 2 ,
where α , β , γ are coefficients to be determined. By applying ρ 2 ,
ρ 2 | λ = 1 Δ t ( m + n ) C + ( 0 ) | 0 : 0 + Δ t 2 2 i 1 A i 2 B C + ( i 1 i 2 ) | i 1 : i 2 i 1 i 2 A C + + ( i 1 i 2 ) | i 1 i 2 : 0 i 1 i 2 B C ( i 1 i 2 ) | 0 : i 1 i 2 + Δ t i 1 A i 2 B α i 1 i 2 C + ( i 1 i 2 ) | 0 : 0 Δ t 2 i 1 i 2 A β i 1 i 2 C ( i 1 i 2 ) | 0 : 0 + i 1 i 2 B γ i 1 i 2 C + + ( i 1 i 2 ) | 0 : 0 .
From this, we have the following equations
λ = 1 Δ t ( m + n ) C ( 0 ) + Δ t 2 2 i 1 A i 2 B α i 1 i 2 C + ( i 1 i 2 ) i 1 i 2 A β i 1 i 2 C ( i 1 i 2 ) i 1 i 2 B γ i 1 i 2 C + + ( i 1 i 2 ) ,
λ α i 1 i 2 = Δ t C + ( i 1 i 2 ) ( i 1 A , i 2 B ) ,
λ β i 1 i 2 = Δ t 2 C + + ( i 1 i 2 ) ( i 1 i 2 A ) ,
λ γ i 1 i 2 = Δ t 2 C ( i 1 i 2 ) ( i 1 i 2 B ) .
After eliminating coefficients α , β , γ , we obtain
λ 2 1 Δ t ( m + n ) C + ( 0 ) λ ( Δ t ) 2 i 1 A i 2 B C + ( i 1 i 2 ) 2 + 1 4 i 1 i 2 A C + + ( i 1 i 2 ) 2 + 1 4 i 1 i 2 B C + + ( i 1 i 2 ) 2 = 0 ,
and eigenvalues are
λ = 1 2 1 Δ t ( m + n ) C + ( 0 ) ± 1 2 1 Δ t ( m + n ) C + ( 0 ) 2 + terms C 2 1 / 2 = 1 Δ t ( m + n ) C + ( 0 ) + O ( g 4 ) , O ( g 4 ) .
Thus, up to O ( g 2 ) , ρ 2 does not have negative eigenvalues.

4.3. Negativity

As the ρ 2 does not have negative eigenvalues, the negativity of the bipartite state ρ is given by the eigenvalue of ρ 1 and we obtain the following key formula of the negativity in this paper
E N = m n C + + ( r ) 2 + m n 2 2 C + ( d ) 2 1 / 2 C + ( 0 ) + m + n 2 1 C + ( d ) .
For a fixed value of the total number of detectors m + n ,
E N m = n m 2 m n C + + ( r ) 2 + m n 2 2 C + ( d ) 2 1 / 2 C + + ( r ) 2 + C + ( d ) 2 ,
and if the total number of detectors is even, m = n provides a maximum value of the negativity
E N = n C + + ( r ) C + ( 0 ) + ( n 1 ) C + ( d ) n .
Figure 3 shows an example of detection of super horizon scale entanglement with ten detectors. In this case, ten detectors catch nonzero negativity on the super horizon scale r = 2.8 H 1 .
Using Equation (14), the negativity for the massless conformal scalar field is
E N n 1 H 2 r 2 1 n 1 4 sin 2 θ + n 1 4 sin 2 θ + H 2 d 2 .
For a given n and separation d of detectors in the same group, we introduce a ratio δ = d / r with 0 δ 1 . Then, the maximum distance of entanglement detection is given by
r max = 2 sin θ H 1 δ 2 δ 2 1 / 2 n 1 / 4 n + 4 δ 2 ( 1 δ 2 ) 2 1 / 2 n 1 / 2 1 / 2 .
For δ = 1 , we have
r max = 2 H 1 n 1 / 4 sin θ ,
and r max can become super horizon scale as n increases. For δ < 1 , in the limit of large n 1 , r max approaches the following asymptotic value
r max 2 H 1 sin θ 1 δ 2 .
For the massless minimal scalar field, although we cannot obtain the analytic expression of r max , it is possible to obtain the asymptotic formula. For n 1 with δ 1 ,
r max ( 2 sin 2 θ ) 1 / 2 n ln n 1 / 4 H 1 for δ = 1 , ( sin 2 θ ) 1 / 2 δ ( ln δ ) 1 / 4 H 1 for δ 1 .
For δ 1 with θ = π / 4 , r max can exceeds 2 H 1 for δ greater than ≈0.96. Figure 4 shows r max as a function of number of detectors n.
For both types of scalar fields, r max approaches constant values as n for δ = d / r < 1 . These constant values become larger than 2 H 1 if δ is sufficiently close to unity. If we take δ = 1 , r max grows as the number of detectors increases and becomes infinity as n . Therefore, it is possible to detect the super horizon scale entanglement if we prepare sufficiently large number of detectors.

5. Monogamy Inequality

We expect that detectability of large scale entanglement on the super horizon scale is related to multipartite entanglement. To make the connection clear, we check the monogamy inequality of negativity [16] for the present m + n detectors system. For a tripartite system A B 1 B 2 , the negativity between A and B 1 , between A and B 2 and between A and B 1 B 2 should obey the following monogamy inequality
E N 2 ( A : B 1 ) + E N 2 ( A : B 2 ) E N 2 ( A : B 1 B 2 ) .
This inequality implies
j 1 , j 2 ( E N ( A j 1 : B j 2 ) ) 2 ( E N ( A : B ) ) 2 ,
where A = j = 1 , , m A j and B = j = 1 , , n B j . For the present detectors model with m = n , the left hand side of the inequality (51) is
n 2 × | C + + ( r ) | C + ( 0 ) 2 ,
and the right hand side of the inequality (51) is
n 2 | C + + ( r ) | C + ( 0 ) + ( n 1 ) C + ( d ) n 2 = n 2 C + + ( r ) C + ( 0 ) + 1 1 n C + ( 0 ) C + ( d ) 2 .
As C + ( 0 ) C + ( d ) , which can be directly confirmed from Equations (14) and (15), the inequality (51) definitely holds.
The difference between both side of the inequality (51) can be interpreted as the residual entanglement and regarded as quantifying degrees of multipartite entanglement [16]. This quantity is
n 2 1 1 n 2 C + ( 0 ) C + ( d ) 2 + 2 1 1 n C + + ( r ) C + ( 0 ) C + ( 0 ) C + ( d ) .
If we take d = 0 , this difference becomes zero and r max reduces to 2 H 1 for the massless conformal scalar case, which is the same value attained by a pair of detectors. The residual entanglement becomes maximum for d = r , in which case r max can become larger than the Hubble horizon scale provided that sufficiently large number of detectors are prepared. Thus, the effect of multipartite entanglement is crucial for detection of the bipartite entanglement on the super horizon scale.

6. Summary

We investigated detection of entanglement of the scalar field on the super horizon scale in de Sitter space using multiple detectors. For this purpose, we obtained the formula of the negativity for an m + n qubit detectors system. The maximum possible distance of detecting nonzero values of negativity is bounded by n 1 / 4 H 1 for the massless conformal scalar field and ( n / ln n ) 1 / 4 H 1 for the massless minimal scalar field. For both type of scalar fields, these bounds grow as the number of detectors increases. Thus, it is possible to detect entanglement on the super horizon scale if we prepare sufficient large number of detectors.
As a practical method to confirm entanglement of a detectors system, the test of Bell-CHSH [17] inequality for a pair of detectors is usually accepted. In our previous studies [2,4], we have confirmed that there is no violation of Bell-CHSH inequality on the super horizon scale. However, there is a possibility that effect of multipartite entanglement can violate Bell-like inequalities. Bell-Mermin-Klyshoko (BMK) inequalities [18,19,20] is such a candidate which can capture multipartite entanglement. Several authors discuss cosmological implication of this inequalities [21]. It may be interesting task to evaluate degrees of violation of these inequalities for the present detectors model.

Acknowledgments

We would like to thank S. Ishizaka for informing us monogamy relation of negativity for qubit system. This work was supported in part by the JSPS KAKENHI Grant Number 16H01094.

Author Contributions

Both authors contributed equally to this work. Both authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Steeg, G.V.; Menicucci, N.C. Entangling power of an expanding universe. Phys. Rev. D 2009, 79. [Google Scholar] [CrossRef] [Green Version]
  2. Nambu, Y.; Ohsumi, Y. Classical and quantum correlations of scalar field in the inflationary universe. Phys. Rev. D 2011, 84. [Google Scholar] [CrossRef]
  3. Martín-Martínez, E.; Menicucci, N.C. Cosmological quantum entanglement. Class. Quantum Gravity 2012, 29. [Google Scholar] [CrossRef]
  4. Nambu, Y. Entanglement Structure in Expanding Universes. Entropy 2013, 15, 1847–1874. [Google Scholar] [CrossRef]
  5. Kukita, S.; Nambu, Y. Entanglement dynamics in de Sitter spacetime. arXiv, 2017; arXiv:1706.09175. [Google Scholar]
  6. Nambu, Y. Entanglement of quantum fluctuations in the inflationary universe. Phys. Rev. D 2008, 78. [Google Scholar] [CrossRef]
  7. Nambu, Y.; Ohsumi, Y. Entanglement of a coarse grained quantum field in the expanding universe. Phys. Rev. D 2009, 80. [Google Scholar] [CrossRef]
  8. Matsumura, A.; Nambu, Y. Large Scale Quantum Entanglement in de Sitter Spacetime. arXiv, 2017; arXiv:1707.08414. [Google Scholar]
  9. Lidar, D.A.; Bihary, Z.; Whaley, K. From completely positive maps to the quantum Markovian semigroup master Equation. Chem. Phys. 2001, 268, 35–53. [Google Scholar] [CrossRef]
  10. Schaller, G.; Brandes, T. Preservation of positivity by dynamical coarse graining. Phys. Rev. A 2008, 78. [Google Scholar] [CrossRef]
  11. Benatti, F.; Floreanini, R.; Marzolino, U. Entangling two unequal atoms through a common bath. Phys. Rev. A 2010, 81. [Google Scholar] [CrossRef]
  12. Majenz, C.; Albash, T.; Breuer, H.-P.; Lidar, D.A. Coarse graining can beat the rotating-wave approximation in quantum Markovian master equations. Phys. Rev. A 2013, 88. [Google Scholar] [CrossRef]
  13. Nambu, Y.; Kukita, S. Derivation of Markovian Master Equation by Renormalization Group Method. J. Phys. Soc. Jpn. 2016, 85. [Google Scholar] [CrossRef]
  14. Peres, A. Separability Criterion for Density Matrices. Phys. Rev. Lett. 1996, 77, 1413–1415. [Google Scholar] [CrossRef] [PubMed]
  15. Horodecki, M.; Horodecki, P.; Horodecki, R. Separability of mixed states: Necessary and sufficient conditions. Phys. Lett. A 1996, 223, 1–8. [Google Scholar] [CrossRef]
  16. Ou, Y.-C.; Fan, H. Monogamy inequality in terms of negativity for three-qubit states. Phys. Rev. A 2007, 75. [Google Scholar] [CrossRef]
  17. Clauser, J.; Horne, M.; Shimony, A.; Holt, R. Proposed Experiment to Test Local Hidden-Variable Theories. Phys. Rev. Lett. 1969, 23, 880–884. [Google Scholar] [CrossRef]
  18. Mermin, N.D. Extreme quantum entanglement in a superposition of macroscopically distinct states. Phys. Rev. Lett. 1990, 65, 1838–1840. [Google Scholar] [CrossRef] [PubMed]
  19. Belinskii, A.V.; Klyshko, D.N. Interference of light and Bell’s theorem. Physics-Uspekhi 1993, 36, 653–693. [Google Scholar] [CrossRef]
  20. Gisin, N.; Bechmann-Pasquinucci, H. Bell inequality, Bell states and maximally entangled states for n qubits. Phys. Lett. A 1998, 246, 1–6. [Google Scholar] [CrossRef]
  21. Kanno, S.; Soda, J. Infinite violation of Bell inequalities in inflation. arXiv, 2017; arXiv:1705.06199. [Google Scholar]
Figure 1. Setup of a m + n detectors system. The group A consists of m detectors and the group B consists of n detectors. The separation between two groups is r and size of each group is d. We assume d r .
Figure 1. Setup of a m + n detectors system. The group A consists of m detectors and the group B consists of n detectors. The separation between two groups is r and size of each group is d. We assume d r .
Entropy 19 00449 g001
Figure 2. The entanglement detection is possible for parameters ( r , θ ) in regions left side of each solid lines. Left panel: massless conformal scalar field. Right panel: massless minimal scalar field with e-foldings N = 0 , 2 , 10 .
Figure 2. The entanglement detection is possible for parameters ( r , θ ) in regions left side of each solid lines. Left panel: massless conformal scalar field. Right panel: massless minimal scalar field with e-foldings N = 0 , 2 , 10 .
Entropy 19 00449 g002
Figure 3. Dependence of number m of group A on negativity (the massless conformal scalar field with θ = π / 2 ). The total number of detectors is m + n = 10 with r = 2.8 H 1 (super horizon scale). The maximum value of the negativity is attained for m = n = 5 .
Figure 3. Dependence of number m of group A on negativity (the massless conformal scalar field with θ = π / 2 ). The total number of detectors is m + n = 10 with r = 2.8 H 1 (super horizon scale). The maximum value of the negativity is attained for m = n = 5 .
Entropy 19 00449 g003
Figure 4. r max as a function of number of detectors n. Left panel: the massless conformal scalar field with θ = π / 2 . Right panel: the massless minimal scalar field with θ = π / 4 , N = 10 . For both types of scalar fields, with d / r = 1 , r max exceeds the horizon scale 2 H 1 as the number of detectors increases.
Figure 4. r max as a function of number of detectors n. Left panel: the massless conformal scalar field with θ = π / 2 . Right panel: the massless minimal scalar field with θ = π / 4 , N = 10 . For both types of scalar fields, with d / r = 1 , r max exceeds the horizon scale 2 H 1 as the number of detectors increases.
Entropy 19 00449 g004

Share and Cite

MDPI and ACS Style

Kukita, S.; Nambu, Y. Harvesting Large Scale Entanglement in de Sitter Space with Multiple Detectors. Entropy 2017, 19, 449. https://doi.org/10.3390/e19090449

AMA Style

Kukita S, Nambu Y. Harvesting Large Scale Entanglement in de Sitter Space with Multiple Detectors. Entropy. 2017; 19(9):449. https://doi.org/10.3390/e19090449

Chicago/Turabian Style

Kukita, Shingo, and Yasusada Nambu. 2017. "Harvesting Large Scale Entanglement in de Sitter Space with Multiple Detectors" Entropy 19, no. 9: 449. https://doi.org/10.3390/e19090449

APA Style

Kukita, S., & Nambu, Y. (2017). Harvesting Large Scale Entanglement in de Sitter Space with Multiple Detectors. Entropy, 19(9), 449. https://doi.org/10.3390/e19090449

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop