Quantum fluctuations of spacetime generate quantum entanglement between   gravitationally polarizable subsystems

Shijing Cheng; Hongwei Yu; Jiawei Hu

arXiv:1812.03419·gr-qc·December 11, 2018

Quantum fluctuations of spacetime generate quantum entanglement between gravitationally polarizable subsystems

Shijing Cheng, Hongwei Yu, Jiawei Hu

PDF

TL;DR

This paper investigates how quantum vacuum fluctuations of spacetime in linearized quantum gravity can induce entanglement between two gravitationally polarizable subsystems, highlighting the roles of polarization and boundaries.

Contribution

It demonstrates the conditions under which quantum vacuum gravitational fields generate entanglement between subsystems, considering polarization effects and boundary influences.

Findings

01

Entanglement depends on polarization alignment.

02

Boundaries can enable entanglement for parallel subsystems.

03

Vertical alignment entanglement is unaffected by boundaries.

Abstract

There should be quantum vacuum fluctuations of spacetime itself, if we accept that the basic quantum principles we are already familiar with apply as well to a quantum theory of gravity. In this paper, we study, in linearized quantum gravity, the quantum entanglement generation at the neighborhood of the initial time between two independent gravitationally polarizable two-level subsystems caused by fluctuating quantum vacuum gravitational fields in the framework of open quantum systems. A bath of fluctuating quantum vacuum gravitational fields serves as an environment that provides indirect interactions between the two gravitationally polarizable subsystems, which may lead to entanglement generation. We find that the entanglement generation is crucially dependent on the polarizations, i.e, they cannot get entangled in certain circumstances when the polarizations of the subsystems are…

Equations111

H = H_{S} + H_{F} + H_{I} .

H = H_{S} + H_{F} + H_{I} .

H_{S} = \frac{ω}{2} σ_{3}^{(1)} + \frac{ω}{2} σ_{3}^{(2)},

H_{S} = \frac{ω}{2} σ_{3}^{(1)} + \frac{ω}{2} σ_{3}^{(2)},

H_{I} = - \frac{1}{2} [Q_{ij}^{(1)} (t) E_{ij} (x^{(1)} (t)) + Q_{ij}^{(2)} (t) E_{ij} (x^{(2)} (t))],

H_{I} = - \frac{1}{2} [Q_{ij}^{(1)} (t) E_{ij} (x^{(1)} (t)) + Q_{ij}^{(2)} (t) E_{ij} (x^{(2)} (t))],

V = \int ρ_{m} (x) Φ (x) d^{3} x .

V = \int ρ_{m} (x) Φ (x) d^{3} x .

Φ (x) = Φ (x_{0}) + x_{i} \frac{\partial Φ ( x _{0} )}{\partial x _{i}} + \frac{1}{2} x_{i} x_{j} \frac{\partial ^{2} Φ ( x _{0} )}{\partial x _{i} \partial x _{j}} + \dots,

Φ (x) = Φ (x_{0}) + x_{i} \frac{\partial Φ ( x _{0} )}{\partial x _{i}} + \frac{1}{2} x_{i} x_{j} \frac{\partial ^{2} Φ ( x _{0} )}{\partial x _{i} \partial x _{j}} + \dots,

H_{I} = \frac{1}{2} \int d^{3} x ρ_{m} (x) x_{i} x_{j} \frac{\partial ^{2} Φ}{\partial x _{i} \partial x _{j}} .

H_{I} = \frac{1}{2} \int d^{3} x ρ_{m} (x) x_{i} x_{j} \frac{\partial ^{2} Φ}{\partial x _{i} \partial x _{j}} .

H_{I} = - \frac{1}{2} Q_{ij} E_{ij},

H_{I} = - \frac{1}{2} Q_{ij} E_{ij},

Q_{ij} = \int d^{3} x ρ_{m} (x) (x_{i} x_{j} - \frac{1}{3} δ_{ij} r^{2}),

Q_{ij} = \int d^{3} x ρ_{m} (x) (x_{i} x_{j} - \frac{1}{3} δ_{ij} r^{2}),

E_{ij} = - \frac{\partial ^{2} Φ}{\partial x _{i} \partial x _{j}} + \frac{1}{3} δ_{ij} \nabla^{2} Φ .

E_{ij} = - \frac{\partial ^{2} Φ}{\partial x _{i} \partial x _{j}} + \frac{1}{3} δ_{ij} \nabla^{2} Φ .

Q_{ij}^{(α)} (t) = q_{ij}^{(α)} σ_{-} e^{- iω t} + q_{ij}^{(α) *} σ_{+} e^{iω t},

Q_{ij}^{(α)} (t) = q_{ij}^{(α)} σ_{-} e^{- iω t} + q_{ij}^{(α) *} σ_{+} e^{iω t},

h_{ij} = k, λ \sum [a_{k, λ} e_{ij} (k, λ) f_{k} + H . c .],

h_{ij} = k, λ \sum [a_{k, λ} e_{ij} (k, λ) f_{k} + H . c .],

E_{ij} = \frac{1}{2} \ddot{h}_{ij},

E_{ij} = \frac{1}{2} \ddot{h}_{ij},

⟨ E_{ij} (x) E_{k l} (x^{'}) ⟩ = \frac{1}{8 ( 2 π ) ^{3}} \int d^{3} k λ \sum e_{ij} (k, λ) e_{k l} (k, λ) ω^{3} e^{i k \cdot (x - x^{'})} e^{- iω (t - t^{'})},

⟨ E_{ij} (x) E_{k l} (x^{'}) ⟩ = \frac{1}{8 ( 2 π ) ^{3}} \int d^{3} k λ \sum e_{ij} (k, λ) e_{k l} (k, λ) ω^{3} e^{i k \cdot (x - x^{'})} e^{- iω (t - t^{'})},

λ \sum e_{ij} (k, λ) e_{k l} (k, λ) = δ_{ik} δ_{j l} + δ_{i l} δ_{j k} - δ_{ij} δ_{k l} + \hat{k_{i}} \hat{k_{j}} \hat{k_{k}} \hat{k_{l}} + \hat{k_{i}} \hat{k_{j}} δ_{k l} + \hat{k_{k}} \hat{k_{l}} δ_{ij}

λ \sum e_{ij} (k, λ) e_{k l} (k, λ) = δ_{ik} δ_{j l} + δ_{i l} δ_{j k} - δ_{ij} δ_{k l} + \hat{k_{i}} \hat{k_{j}} \hat{k_{k}} \hat{k_{l}} + \hat{k_{i}} \hat{k_{j}} δ_{k l} + \hat{k_{k}} \hat{k_{l}} δ_{ij}

- \hat{k_{i}} \hat{k_{l}} δ_{j k} - \hat{k_{i}} \hat{k_{k}} δ_{j l} - \hat{k_{j}} \hat{k_{l}} δ_{ik} - \hat{k_{j}} \hat{k_{k}} δ_{i l},

\frac{\partial ρ _{tot} ( t )}{\partial t} = - i [H, ρ_{tot} (t)] .

\frac{\partial ρ _{tot} ( t )}{\partial t} = - i [H, ρ_{tot} (t)] .

\frac{\partial ρ ( t )}{\partial t} = - i [H_{eff}, ρ (t)] + D [ρ (t)],

\frac{\partial ρ ( t )}{\partial t} = - i [H_{eff}, ρ (t)] + D [ρ (t)],

H_{eff} = H_{S} - \frac{i}{2} α, β = 1 \sum 2 i, j = 1 \sum 3 H_{ij}^{(α β)} σ_{i}^{(α)} σ_{j}^{(β)},

H_{eff} = H_{S} - \frac{i}{2} α, β = 1 \sum 2 i, j = 1 \sum 3 H_{ij}^{(α β)} σ_{i}^{(α)} σ_{j}^{(β)},

D [ρ (t)] = \frac{1}{2} α, β = 1 \sum 2 i, j = 1 \sum 3 C_{ij}^{(α β)} [2 σ_{j}^{(β)} ρ σ_{i}^{(α)} - σ_{i}^{(α)} σ_{j}^{(β)} ρ - ρ σ_{i}^{(α)} σ_{j}^{(β)}] .

D [ρ (t)] = \frac{1}{2} α, β = 1 \sum 2 i, j = 1 \sum 3 C_{ij}^{(α β)} [2 σ_{j}^{(β)} ρ σ_{i}^{(α)} - σ_{i}^{(α)} σ_{j}^{(β)} ρ - ρ σ_{i}^{(α)} σ_{j}^{(β)}] .

G_{ij k l}^{(α β)} (ω) = \int_{- \infty}^{\infty} d t e^{iω t} G_{ij k l}^{(α β)} (t),

G_{ij k l}^{(α β)} (ω) = \int_{- \infty}^{\infty} d t e^{iω t} G_{ij k l}^{(α β)} (t),

K_{ij k l}^{(α β)} (ω) = \int_{- \infty}^{\infty} d t sgn (t) e^{iω t} G_{ij k l}^{(α β)} (t) = \frac{P}{π i} \int_{- \infty}^{\infty} d λ \frac{G _{ij k l}^{(α β)} ( λ )}{λ - ω},

K_{ij k l}^{(α β)} (ω) = \int_{- \infty}^{\infty} d t sgn (t) e^{iω t} G_{ij k l}^{(α β)} (t) = \frac{P}{π i} \int_{- \infty}^{\infty} d λ \frac{G _{ij k l}^{(α β)} ( λ )}{λ - ω},

C_{ij}^{(α β)}

C_{ij}^{(α β)}

H_{ij}^{(α β)}

A^{(α β)}

A^{(α β)}

A^{(α β)}

G^{(α β)} (ω)

G^{(α β)} (ω)

K^{(α β)} (ω)

⟨ μ ∣ C^{(11)} ∣ μ ⟩ ⟨ ν ∣ (C^{(22)})^{T} ∣ ν ⟩ < ∣ ⟨ μ ∣ Re (C^{(12)}) ∣ ν ⟩ ∣^{2},

⟨ μ ∣ C^{(11)} ∣ μ ⟩ ⟨ ν ∣ (C^{(22)})^{T} ∣ ν ⟩ < ∣ ⟨ μ ∣ Re (C^{(12)}) ∣ ν ⟩ ∣^{2},

A_{3}^{2} + B_{1} B_{2} > A_{1} A_{2} .

A_{3}^{2} + B_{1} B_{2} > A_{1} A_{2} .

q_{33}^{(1)} = 0, q_{13}^{(1)} q_{13}^{(2)} + q_{23}^{(1)} q_{23}^{(2)} = 0, 2 q_{12}^{(1)} q_{12}^{(2)} + q_{11}^{(1)} (q_{11}^{(2)} - q_{22}^{(2)}) = 0.

q_{33}^{(1)} = 0, q_{13}^{(1)} q_{13}^{(2)} + q_{23}^{(1)} q_{23}^{(2)} = 0, 2 q_{12}^{(1)} q_{12}^{(2)} + q_{11}^{(1)} (q_{11}^{(2)} - q_{22}^{(2)}) = 0.

A_{3} = \frac{1}{128 π L ^{5}} [F_{1} sin (ω L) + K_{1} cos (ω L) ω L],

A_{3} = \frac{1}{128 π L ^{5}} [F_{1} sin (ω L) + K_{1} cos (ω L) ω L],

F_{1} = 4 (3 - 3 ω^{2} L^{2} + ω^{4} L^{4}) q_{12}^{(1)} q_{12}^{(2)} + (9 - 5 ω^{2} L^{2} + ω^{4} L^{4}) (q_{11}^{(1)} q_{11}^{(2)} + q_{22}^{(1)} q_{22}^{(2)})

F_{1} = 4 (3 - 3 ω^{2} L^{2} + ω^{4} L^{4}) q_{12}^{(1)} q_{12}^{(2)} + (9 - 5 ω^{2} L^{2} + ω^{4} L^{4}) (q_{11}^{(1)} q_{11}^{(2)} + q_{22}^{(1)} q_{22}^{(2)})

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

Quantum fluctuations of spacetime generate quantum entanglement between gravitationally polarizable subsystems

Shijing Cheng1, Hongwei Yu1111Corresponding author: [email protected] and Jiawei Hu1222Corresponding author: [email protected]

1Department of Physics and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, Changsha, Hunan 410081, China

Abstract

There should be quantum vacuum fluctuations of spacetime itself, if we accept that the basic quantum principles we are already familiar with apply as well to a quantum theory of gravity. In this paper, we study, in linearized quantum gravity, the quantum entanglement generation at the neighborhood of the initial time between two independent gravitationally polarizable two-level subsystems caused by fluctuating quantum vacuum gravitational fields in the framework of open quantum systems. A bath of fluctuating quantum vacuum gravitational fields serves as an environment that provides indirect interactions between the two gravitationally polarizable subsystems, which may lead to entanglement generation. We find that the entanglement generation is crucially dependent on the polarizations, i.e, they cannot get entangled in certain circumstances when the polarizations of the subsystems are different while they always can when the polarizations are the same. We also show that the presence of a boundary may render parallel aligned subsystems entangled which are otherwise unentangled in a free space. However, the presence of the boundary does not help in terms of entanglement generation if the two subsystems are vertically aligned.

I Introduction

Recently, gravitational wave signals from black hole merging systems have been directly detected Abbotta ; Abbotta1 ; Abbotta2 ; Abbotta3 , and this confirms the prediction based on Einstein’s general relativity over a hundred years ago gw . Naturally, one may wonder what happens if gravitational waves are quantized. One of the consequences, if gravity is quantized, is the quantum fluctuations of spacetime itself. A direct result of spacetime fluctuations is the flight time fluctuations of a probe light signal from its source to a detector Ford95 ; Yua ; Yuc . Another expected effect is the Casimir-like force which arises from the quadrupole moments induced by quantum gravitational vacuum fluctuations Quach15 ; Ford15 ; Wu ; Hu ; Wu2 , in close analogy to the Casimir and the Casimir-Polder forces Casimir ; cp .

In the present paper, we are concerned with yet another effect associated with quantum fluctuations of spacetime, that is, quantum entanglement generation by quantum fluctuations of spacetime. Quantum entanglement is crucial to our understanding of quantum theory, and it has many interesting applications in various novel technologies. However, a significant challenge to the realization of these quantum technologies with quantum entanglement as a key resource is the environmental noises that lead to the quantum to classical transition. In the quantum sense, one environment that no physical system can be isolated from is the vacuum that fluctuates all the time. On one hand, there are obviously vacuum fluctuations of matter fields, and on the other hand, there should also be quantum vacuum fluctuations of spacetime itself if we accept that the basic quantum principles we are already familiar with apply as well to a quantum theory of gravity. Currently, there are discussions with vacuum fluctuations of quantized matter fields as inevitable environmental noises that cause quantum decoherence decoherence .

Like matter fields, the fluctuations of gravitational fields, i.e. the fluctuations of spacetime itself, may also cause quantum decoherence. Since gravitation can not be screened unlike electromagntism, gravitational decoherence is universal. Different models for gravitational decoherence have been proposed Anastopoulos ; Kay ; Power ; Kok ; Wang ; Breuer ; blhu ; Blencowe ; Lorenci ; Oniga16 , see Ref. Bassi17 for a recent review. However, decoherence is not the only role the environment plays. In certain circumstances, the indirect interactions provided by the common bath may also create rather than destroy entanglement between the subsystems via spontaneous emission and photon exchange Braun ; Plenio99 ; Kim ; Schneider ; Basharov ; Jakobczyk ; Reznik ; Piani ; Floreanini ; Zhanga ; Zhangb ; Hu15 ; Yang ; Martnez ; Valentini ; Reznik1 ; Floreanini1 . Here the bath can be either in the vacuum state Jakobczyk ; Reznik ; Martnez ; Valentini ; Reznik1 ; Plenio99 ; Floreanini1 , or in a thermal state Braun ; Basharov ; Schneider ; Kim . In Ref. Piani , a general discussion showed that two independent atoms in a common bath can be entangled during a Markovian, completely positive reduced dynamics. In Ref. Floreanini , this approach has been applied in the study of the entanglement generation for two atoms immersed in a common bath of massless scalar fields, and found that entanglement generation can be manipulated by varying the bath temperature and the distance between the two atoms. This work was further generalized to the case in the presence of a reflecting boundary Zhanga ; Zhangb , which showed that the presence of a boundary may offer more freedom in controlling entanglement generation.

In this paper, we study whether a common bath of fluctuating gravitational fields may generate entanglement between subsystems just as matter fields do. We consider a system composed of two gravitationally polarizable two-level subsystems in interaction with a bath of quantum gravitational fields in vacuum. We will show that in certain conditions, the two subsystems may get entangled due to the fluctuating gravitational fields. Like matter fields, the fluctuations of gravitational fields are also expected to be modified if a boundary is present. Although it is generally believed that gravitational waves can hardly be absorbed or reflected, there have been proposals that the interaction between gravitational waves and quantum fluids might be significantly enhanced compared with normal matter (see Ref. Kiefer for a review), and recently there have been interesting conjectures that superconducting films may act as mirrors for gravitational waves due to the so-called Heisenberg-Coulomb effect hc1 ; hc2 . Therefore, we are interested in how the conditions for entanglement generation due to quantum fluctuations of spacetime is modified if such gravitational boundaries exist.

II the basic formalism

The system we consider consists of two gravitationally polarizable two-level subsystems, which are weakly coupled with a bath of fluctuating quantum gravitational fields in vacuum. In principle, the gravitationally polarizable system can be any quantum system that carries nonzero quadrupole polarizability. For example, a Bose-Einstein condensate (BEC) is such a system of which the quadrupole polarizability can be estimated Hu . The Hamiltonian of the whole system takes the following form

[TABLE]

Here $H_{S}$ represents the Hamiltonian of the two gravitationally polarizable two-level subsystems, which can be described as

[TABLE]

where $\sigma_{i}^{(1)}=\sigma_{i}\otimes\sigma_{0}$ , $\sigma_{i}^{(2)}=\sigma_{0}\otimes\sigma_{i}$ , with $\sigma_{i}\ (i=1,2,3)$ being the Pauli matrices, $\sigma_{0}$ the $2\times 2$ unit matrix, and $\omega$ the energy level spacing. $H_{F}$ is the Hamiltonian of the external gravitational fields, the details of which are not needed here. $H_{I}$ describes the quadrupolar interaction between the gravitationally polarizable subsystems and the fluctuating gravitational fields, which can be expressed as

[TABLE]

where $Q_{ij}^{(\alpha)}$ $(\alpha=1,2)$ is the induced quadrupole moment of the $\alpha$ -th subsystem, and $E_{ij}=-\nabla_{i}\nabla_{j}\phi$ with $\phi$ being the gravitational potential. The quadrupolar interaction Hamiltonian (3) can be obtained as follows. The energy of a localized mass distribution $\rho_{m}(x)$ in the presence of an external gravitational potential $\Phi(x)$ is

[TABLE]

When $\Phi(x)$ varies slowly over the region where the mass is located, it can be expanded as

[TABLE]

so the quadrupolar interaction term reads

[TABLE]

Since $\nabla^{2}\Phi=0$ in an empty space, the equation above can be rewritten as

[TABLE]

where

[TABLE]

and

[TABLE]

In general relativity, $E_{ij}$ is defined as the Weyl tensor $C_{i0j0}$ , which can be shown to coincide with the expression given in Eq. (9) in the Newtonian limit. Here $E_{ij}=C_{i0j0}$ and its dual tensor $B_{ij}=-\frac{1}{2}\epsilon_{ikl}{C^{kl}}_{j0}$ are the gravito-electric and gravito-magnetic tensors which satisfy the linearized Einstein field equations organized in a form similar to the Maxwell equations Campbell76 ; matte53 ; campbell71 ; Szekeres ; maartens98 ; Ruggiero02 ; ramos10 ; Ingraham .

In the present paper, we are considering two quantum subsystems in interaction with quantum gravitational vacuum fluctuations, so the quadrupole moment $Q_{ij}^{(\alpha)}$ ought to be a quantum operator. In the interaction picture, the quadrupole operator can be written as

[TABLE]

where $\sigma_{+}=|+\rangle\langle-|$ , $\sigma_{-}=|-\rangle\langle+|$ , are the raising and lowering operators of the subsystem respectively, $q_{ij}^{(\alpha)}=\langle{0|Q_{ij}^{(\alpha)}|1}\rangle$ is the transition matrix element of the quadrupole operator of the $\alpha$ -th subsystem, which is symmetric and traceless, i.e., $q_{ii}^{(\alpha)}=0$ and $q_{ij}^{(\alpha)}=q_{ji}^{(\alpha)}$ , thus leaving only five independent components $q^{(\alpha)}_{11}$ , $q^{(\alpha)}_{22}$ , $q^{(\alpha)}_{12}$ , $q^{(\alpha)}_{13}$ and $q^{(\alpha)}_{23}$ . The quadrupole moments are induced by the fluctuating gravitational fields. In this paper, we assume that the gravitational quadrupole polarizabilities are real, and so are the quadrupole moments. The gravito-electric field $E_{ij}$ is also supposed to be quantized. The metric for a flat spacetime metric with a perturbation can be expressed as $g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}$ , where $\eta_{\mu\nu}$ denotes the flat spacetime metric, and $h_{\mu\nu}$ a linearized perturbation. In the transverse traceless (TT) gauge, the spacetime perturbation can be quantized as Yua

[TABLE]

where $f_{\mathbf{k}}=(2\omega(2\pi)^{3})^{-\frac{1}{2}}e^{i(\mathbf{k}\cdot\mathbf{x}-\omega t)}$ is the field mode and $e_{\mu\nu}(\mathbf{k},\lambda)$ is the polarization tensor with $\omega=|\mathbf{k}|=(k_{x}^{2}+k_{y}^{2}+k_{z}^{2})^{\frac{1}{2}}$ . Here H.c. denotes the Hermitian conjugate, $\lambda$ labels the polarization state, and units in which $\hbar=c=32\pi G=1$ have been adopted. From the definition of $E_{ij}$ , we have

[TABLE]

where a dot means $\frac{\partial}{\partial t}$ . The correlation function can then be obtained as Yua

[TABLE]

where

[TABLE]

with $\hat{k_{i}}={k_{i}}/{k}$ .

We assume that initially the two subsystems are decoupled with the quantum gravitational fields, i.e. $\rho_{\rm tot}(0)=\rho(0)\otimes|0\rangle\langle 0|$ , where $\rho(0)$ denotes the initial state of the two subsystems, and $|0\rangle$ the vacuum state of the gravitational fields. The density matrix of the total system satisfies the Liouville equation

[TABLE]

Under the Born-Markov approximation, the reduced density matrix of the two gravitationally polarizable subsystems $\rho(t)=\mathrm{Tr}_{F}[\rho_{\rm tot}(t)]$ satisfies the Kossakowskl-Lindblad master equation Kossakowski ; Lindblad ,

[TABLE]

where

[TABLE]

and

[TABLE]

Introducing the Fourier and Hilbert transforms of the gravitational field correlation function $G_{ijkl}^{(\alpha\beta)}(t-t^{\prime})=\langle{E_{ij}(t,x_{\alpha})E_{kl}(t^{\prime},x_{\beta})}\rangle$ ,

[TABLE]

where $P$ denotes the principal value. The coefficient matrix $H_{ij}^{(\alpha\beta)}$ and $C_{ij}^{(\alpha\beta)}$ can then be expressed as

[TABLE]

where

[TABLE]

with

[TABLE]

Here the effective Hamiltonian $H_{\rm eff}$ corresponds to the Lamb-like shift and an environment induced direct coupling between the two subsystems. In the following, we neglect this term and concentrate on the effects produced by the dissipative part $\mathcal{D}[\rho(t)]$ .

III The Condition for Entanglement Generation

In this section, we investigate whether the two subsystems become entangled or not at a neighborhood of the initial time with the help of the partial transposition criterion Peres ; Horodecki , i.e., a two-atom state $\rho$ is entangled if and only if the operation of the partial transposition of $\rho$ does not preserve its positivity. We assume that the initial state is $\rho(0)=|+\rangle\langle+|\otimes|-\rangle\langle-|$ , which is pure and separable, where $|-\rangle$ and $|+\rangle$ are the ground and excited state respectively. In this case, the partial transposition criterion can be found to be equivalent to the following condition, i.e., entanglement can be generated at the neighbourhood of $t=0$ if and only if Floreanini ; Piani

[TABLE]

where the subscript $T$ denotes the matrix transposition, and the three-dimensional vectors $\mu_{i}$ , $\nu_{i}$ are $\mu_{i}=\nu_{i}=\{1,-i,0\}$ . For brevity, we make the following substitution in Eq. (21), $A^{(11)}\rightarrow A_{1}$ , $A^{(22)}\rightarrow A_{2}$ , and $A^{(12)},A^{(21)}\rightarrow A_{3}$ , and do the same to the coefficients $B^{(\alpha\beta)}$ . Plugging the new expressions of $C_{ij}^{(\alpha\beta)}$ into Eq. (27), leads to

[TABLE]

Now we investigate the condition for entanglement generation between a pair of gravitationally polarizable subsystems both in a free space and in a space with a Dirichlet boundary. In the presence of a boundary, we will consider cases of subsystems aligned parallel to and perpendicular to the boundary respectively.

III.1 Entanglement generation in the free Minkowski spacetime

We assume the two gravitationally polarizable subsystems are placed at the points ${\bf x}_{1}=(0,0,0)$ and ${\bf x}_{2}=(0,0,L)$ respectively, where $L$ denotes the separation between the two subsystems, as shown in Fig. 1.

The explicit expressions of the coefficients $A_{i}$ and $B_{i}$ in Eq. (28) are given in appendix A, see Eqs. (40)-(A.1). Since $A_{1}=B_{1}$ and $A_{2}=B_{2}$ , it is clear that the criterion for entanglement generation (28) becomes $A_{3}^{2}>0$ . That is, entanglement can be generated as long as $A_{3}$ , which is dependent on $q_{ij}^{(\alpha)}$ and $\omega L$ , is nonzero.

When the polarizations of the two subsystems are the same, for example, when $q^{(\alpha)}_{11}=-q^{(\alpha)}_{22}=|q|$ or $q^{(\alpha)}_{11}=-q^{(\alpha)}_{33}=|q|$ and other components being zero, we plot, in Fig. 2, the coefficient $A_{3}$ as a function of $\omega L$ in the unit of $\Gamma_{0}$ , with $\Gamma_{0}=|q|\omega^{5}/\pi$ . As $\omega L$ grows, $A_{3}$ oscillates around zero with a decreasing amplitude, and it approaches zero as $\omega L\rightarrow\infty$ . That is, when the separation of the two subsystems is infinitely large, they are always separable. Even for finite separations, there are some special values of $\omega L$ that give $A_{3}=0$ , so entanglement can not be generated at those separations.

When the polarizations of the two subsystems are the same, i.e., $q^{(1)}_{ij}=q^{(2)}_{ij}=q_{ij}$ , there are no solutions to $F=K=0$ (see Eqs. (A.1)-(A.1)) for any given $L$ , i.e., irrespective of the polarization of the two subsystems, the entanglement between the two subsystems with the same polarization can always be generated for a given separation. Therefore, when the polarizations of the two subsystems are the same, they can always get entangled if the separation is finite, except for a series of special values of $\omega L$ which are the zero points of $A_{3}$ . Similar conclusions have been drawn in the case of two independent atoms coupled with massless scalar (matter) fields Floreanini . This suggests that there is no difference between spacetime fluctuations and matter fields fluctuations when the entanglement generation is concerned with two subsystems of the same gravitational polarization.

When the polarizations of the two subsystems are different, it can be shown by solving the equations $F=K=0$ that the two subsystems may not get entangled for any given separations when the components of quadrupole moments of the two subsystems satisfy the following conditions simultaneously,

[TABLE]

Note that we do not distinguish the two subsystems here, i.e. the superscripts $(1)$ and $(2)$ in the equations above can be exchanged. Recall that $q_{ij}$ is traceless, so $q_{33}^{(\alpha)}=0$ indicates that $q_{11}^{(\alpha)}=-q_{22}^{(\alpha)}$ . The components of the quadrupole moment can be interpreted as follows. The diagonal components $q_{11}$ , $q_{22}$ and $q_{33}$ respectively represent the contributions to the total mass quadrupole moment from the mass distributed along the $x$ axis, $y$ axis and $z$ axis, and the off-diagonal components $q_{12}$ , $q_{13}$ and $q_{23}$ respectively represent the contributions to the total mass quadrupole moment from the mass distributed in the $xoy$ , $xoz$ and $yoz$ plane. Therefore, when the mass distributions of the two subsystems induced by quantum gravitational vacuum fluctuations satisfy the condition (29), the two subsystems remain disentangled.

As an example, we consider the case when the two subsystems are only polarizable in the $xoy$ plane, i.e. $q^{(1)}_{3i}=q^{(2)}_{3i}=0$ . In this case, we have

[TABLE]

where

[TABLE]

It is obvious that when

[TABLE]

$A_{3}$ is always zero, i.e., the entanglement generation cannot happen. Otherwise, the value of ${A_{3}}$ oscillates around zero as $\omega L$ varies with a changing amplitude. This reveals a clear difference between entanglement generation by vacuum fluctuations of massless scalar (matter) fields and that of the spacetime, that is, entanglement generation may never happen for certain polarizations if only spacetime fluctuations are considered.

III.2 Entanglement generation for two gravitationally polarizable subsystems aligned parallel to the boundary

Now we concentrate on the effect of a Dirichlet boundary (recall our comments in the Introduction for gravitational boundaries ) on entanglement generation between two gravitationally polarizable subsystems in a bath of fluctuating gravitational fields. We assume that the Dirichlet boundary is located at $y$ =0 in the $xoz$ plane, and the two subsystems separated from each other by a distance $L$ is placed along the $z$ direction at a distance of $y$ , see Fig. 3 .

The coefficients $A_{1},~{}A_{2},~{}A_{3}$ and $B_{1},~{}B_{2},~{}B_{3}$ can be calculated from Eq. (23) and Eq. (25) with the Fourier transform of Eq. (43), and the results are lengthy, so we do not give them explicitly here. As before, we will examine when $A^{\prime}_{3}=0$ to determine whether entanglement can be generated or not, since $A^{\prime}_{i}=B^{\prime}_{i}~{}(i=1,2,3)$ still holds in the present circumstances. In this case, the conditions that entanglement cannot be generated for any $L$ are solved as follows

[TABLE]

Now we consider the same example as before, i.e. when the two subsystems are only polarizable in the $xoy$ plane, $q^{(1)}_{3i}=q^{(2)}_{3i}=0$ . In this case, we have

[TABLE]

where the explicit expressions of $M_{1}$ and $N_{1}$ are given in appendix (see Eqs. (A.2)-(A.2)). Obviously, if entanglement generation cannot happen for any given $L$ and $y$ , the corresponding coefficients satisfy $F_{1}=K_{1}=M_{1}=N_{1}=0$ , and then we have

[TABLE]

Comparing the result above with that obtained in the corresponding case without a boundary, we find the condition Eq. (36) is a special case of Eq. (33). That is, if the components of the quadrupole moments satisfy $q_{11}^{(1)}q_{11}^{(2)}=q_{22}^{(1)}q_{22}^{(2)}=-q_{12}^{(1)}q_{12}^{(2)}\neq 0$ , the two subsystems can be entangled in the presence of a boundary, while they remain separable in the free space.

III.3 Entanglement creation for two gravitationally polarizable subsystems aligned vertical to the boundary

Now we investigate the case when the alignment of the two gravitationally polarizable subsystems is vertical to the boundary. We assume that the Dirichlet boundary is placed at $z$ =0 in the $xoy$ plane and the two subsystems are placed on the $z$ -axis with a separation $L$ , the distance from the boundary to the nearer subsystem being $z$ , see Fig. 4.

Following the same procedures, we find that the conditions that entanglement generation cannot happen for any given $L$ and $z$ are the same with that obtained in the case without a boundary, see Eq. (29). That is, when the quadrupole moments of the two subsystems satisfy certain conditions such that entanglement cannot be generated in the free space, it cannot be generated in the presence of a boundary placed vertically to the alignment of the subsystems either.

IV Summary

In this paper, we have investigated the entanglement generation at the neighborhood of the initial time between two independent gravitationally polarizable two-level subsystems in interaction with a bath of fluctuating quantum gravitational fields in vacuum both with and without a boundary. The partial transposition criterion has been applied to determine whether entanglement can be generated or not at the beginning of evolution. In the free space case, when the polarizations of the two subsystems are the same, the two subsystems can always get entangled as long as the separation is finite. This is similar to what happens when the fluctuations of scalar (matter) fields are considered. When the polarizations of the two subsystems are different, they cannot get entangled in certain circumstances, whatever the separation is. This is in sharp contrast with the case of massless scalar (matter) fields. In the presence of a boundary, we find that in some of the cases where entanglement cannot be generated in a free space, the presence of a boundary placed parallel to the alignment of the subsystems may render the subsystems entangled, but this does not happen for a boundary placed vertically.

Acknowledgements.

This work was supported in part by the NSFC under Grants No. 11435006, No. 11690034, and No. 11805063.

Appendix A

The appendix mainly shows some calculations of the coefficients $A_{i}$ and $B_{i}$ in Eq. (28) for the condition of entanglement generation.

A.1 The case in the free Minkowski spacetime

In the free Minkowski spacetime, we assume the two gravitationally polarizable subsystems are placed at the points ${\bf x}_{1}=(0,0,0)$ and ${\bf x}_{2}=(0,0,L)$ respectively, where $L$ denotes the separation between the two subsystems. Straightforward calculations of Eqs. (LABEL:g1) and (19) show that $\mathcal{G}_{ijkl}^{(\alpha\beta)}(-\omega)=0$ , and the nonzero components of $\mathcal{G}_{ijkl}^{(\alpha\beta)}(\omega)$ can be written as

[TABLE]

When $\alpha=\beta$ , we have $f_{ijkl}^{(11)}=f_{ijkl}^{(22)}$ , where

[TABLE]

and when $\alpha\neq\beta$ , we have $f_{ijkl}^{(12)}=f_{ijkl}^{(21)}$ , where

[TABLE]

The coefficients can be obtained from Eq. (23) and Eq. (25) as

[TABLE]

with

[TABLE]

and

[TABLE]

A.2 The case for two subsystems aligned parallel to the boundary

As for two gravitationally polarizable subsystems aligned parallel to the boundary, we assume that the boundary is located at $y$ =0 in the $xoz$ plane, and the two subsystems separated from each other by a distance $L$ is placed along the $z$ direction at a distance of $y$ . With the help of the method of images, the correlation functions can be written in the following form

[TABLE]

and the coefficients $A^{\prime}_{1},~{}A^{\prime}_{2},~{}A^{\prime}_{3}$ and $B^{\prime}_{1},~{}B^{\prime}_{2},~{}B^{\prime}_{3}$ can be calculated from Eqs. (23)-(25) with the Fourier transform of Eq. (43). The results are lengthy, so we do not give them explicitly here, but $A^{\prime}_{i}=B^{\prime}_{i}~{}(i=1,2,3)$ still holds as it does in the free space.

When the two subsystems are only polarizable in the $xoy$ plane, i.e., $q^{(1)}_{3i}=q^{(2)}_{3i}=0$ , we have

[TABLE]

where

[TABLE]

and

[TABLE]

Bibliography59

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1(1) B. P. Abbott et al (LIGO and Virgo Collaborations), Phys. Rev. Lett. 116 , 061102 (2016).
2(2) B. P. Abbott et al (LIGO and Virgo Collaborations), Phys. Rev. Lett. 116 , 241103 (2016).
3(3) B. P. Abbott et al (LIGO and Virgo Collaborations), Phys. Rev. Lett. 118 , 221101 (2017).
4(4) B. P. Abbott et al (LIGO and Virgo Collaborations), Phys. Rev. Lett. 119 , 141101 (2017).
5(5) A. Einstein, Sitzungsber. K. Preuss. Akad. Wiss. 1 , 688 (1916).
6(6) L. H. Ford, Phys. Rev. D 51 , 1692 (1995).
7(7) H. Yu and L. H. Ford, Phys. Rev. D 60 , 084023 (1999).
8(8) H. Yu, N. F. Svaiter, and L. H. Ford, Phys. Rev. D 80 , 124019 (2009).