This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Harvesting correlations from vacuum quantum fields in the presence of a reflecting boundary

Zhihong Liu1 111[email protected], Jialin Zhang1,2 222Corresponding author at [email protected]and Hongwei Yu1,2 333Corresponding author at [email protected] 1 Department of Physics and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, 36 Lushan Rd., Changsha, Hunan 410081, China
2 Institute of Interdisciplinary Studies, Hunan Normal University, 36 Lushan Rd., Changsha, Hunan 410081, China
Abstract

We explore correlations harvesting by two static detectors locally interacting with vacuum massless scalar fields in the presence of an infinite perfectly reflecting boundary. We study the phenomena of mutual information harvesting and entanglement harvesting for two detector-boundary alignments, i.e., parallel-to-boundary and orthogonal-to-boundary alignments. Our results show that the presence of the boundary generally inhibits mutual information harvesting relative to that in flat spacetime without any boundaries. In contrast, the boundary may play a doubled-edged role in entanglement harvesting, i.e., inhibiting entanglement harvesting in the near zone of the boundary while assisting it in the far zone of the boundary. Moreover, there exists an optimal detector energy gap difference between two nonidentical detectors that makes such detectors advantageous in correlations harvesting as long as the interdetector separation is large enough. The value of the optimal detector energy gap difference depends on both the interdetector separation and the detector-to-boundary distance. A comparison of the correlations harvesting in two different alignments shows that although correlations harvesting share qualitatively the same properties, they also display quantitative differences in that the detectors in orthogonal-to-boundary alignment always harvest comparatively more mutual information than the parallel-to-boundary ones, while they harvest comparatively more entanglement only near the boundary.

I Introduction

It has been recognized for a long time that the vacuum state of free quantum fields possesses nonlocal properties and contains correlations between timelike and spacelike separated regions Summers:1985 ; Summers:1987 ; Pozas-Kerstjens:2015 ; Ng:2018 . Such correlations, which may be either quantum entanglement, quantum mutual information or quantum discord M.Naeem:2022 , can be extracted by a pair of initially uncorrelated Unruh-DeWitt (UDW) particle detectors interacting with the vacuum quantum fields. This phenomenon has been known as correlation harvesting, and such harvesting process is referred to as the correlation harvesting protocol Pozas-Kerstjens:2015 .

The correlation harvesting in terms of quantum entanglement has been extensively studied in different circumstances in recent years and has been demonstrated to be sensitive to the properties of spacetime including dimensionality Pozas-Kerstjens:2015 , topology EDU:2016-1 , curvature Steeg:2009 ; Nambu:2013 ; Kukita:2017 ; Ng:2018 ; Ng:2018-2 ; Zhjl:2018 ; Zhjl:2019 ; Robbins:2020jca ; Gallock-Yoshimura:2021 and the presence of boundaries CW:2020 ; Zhjl:2021 ; CW:2019 , the intrinsic motion and energy gap of detectors Salton-Man:2015 ; Zhjl:2020 ; Zhjl:2021 ; Zhjl:2022 ; R.B.Mann:2022 ; Maeso:2022 ; Zhjl:2022.08 , and superpositions of temporal order for detectors Henderson:2020 .

Recently, the study of entanglement harvesting has been extended to nonidentical detectors with different energy gaps in flat spacetime Zhjl:2022.4 , in contrast to previous studies where detectors are usually assumed to be identical. It was argued that the presence of energy gap difference generally enlarges the rangefinding of entanglement harvesting (the harvesting-achievable range for the interdetector separation), and two nonidentical detectors can harvest more entanglement than the identical ones if the interdetector separation is not too small with respect to the interaction duration time. However, the validity of such results in other circumstances such as in curved spacetime or in the presence of boundaries merits further exploration. In fact, it has been recognized that the presence of reflecting boundaries in flat spacetime modifies the fluctuations of quantum fields, resulting in some intriguing quantum effects, such as the Casimir effect Casimir:1948 , the light-cone fluctuations Yu:1999 , the geometric phase Zhjl:2016 and the modification for the radiative properties of accelerated atoms Yu:2005 ; Yu:2006 ; Rizzuto:2009 . Since the dynamic evolution in time of the detectors system is strongly dependent on the fluctuations of quantum fields, modifications of the fields fluctuations caused by the presence of a boundary have been shown to play a significant role in controlling the entanglement creation in the detectors system Zhang:2007 ; Zhang:2007.1 ; Cheng:2018 . And the phenomena of entanglement harvesting by two UDW detectors locally interacting with the fluctuating quantum fields have also been examined in the presence of a reflecting boundary but in a simple scenario: two identical detectors with not too large energy gap in parallel-to-boundary alignment Zhjl:2021 . It was found, through the numerical evaluation, that the reflecting boundary plays a double-edged role in entanglement harvesting, i.e., degrading the harvested entanglement amount in general while enlarging the entanglement harvesting-achievable interdetector separation range. However, what happens to the entanglement harvesting phenomenon in the presence of a boundary when detectors are nonidentical remains to be investigated.

On the other hand, the correlation harvesting in terms of mutual information, a useful measure on information which quantifies the total amount of classical and quantum correlations including entanglement, has also been recently studied. It has been found that the mutual information harvesting also depends upon the intrinsic properties of spacetime, the detectors’ energy gap and noninertial motion Pozas-Kerstjens:2015 ; Simidzija:2018 ; Kendra:2022 ; M.Naeem:2022 . Remarkably, unlike the entanglement harvesting that has a finite harvesting-achievable separation range and cannot occur near a black hole, mutual information harvesting can occur everywhere with an arbitrarily large interdetector separation Pozas-Kerstjens:2015 and even does not vanish near the event horizon of a black hole at an extremely high local Hawking temperature Kendra:2022 . However, in comparison to the entanglement harvesting which is relatively well understood, much remains to be done on understanding the phenomenon of mutual information harvesting, for example, the mutual information harvesting by two identical/nonidentical detectors near a boundary.

In this paper, we will perform a general study of the correlation harvesting phenomenon for two detectors placed near a reflecting boundary, focusing upon entanglement harvesting and mutual information harvesting phenomenon by nonidentical detectors. Our particular interest lies in the influence of energy gaps on mutual information harvesting and the role played by a perfectly reflecting plane boundary, including correlation harvesting by nonidentical detectors with different energy gaps in the scenarios of parallel-to-boundary and orthogonal-to-boundary alignments.

The paper is organized as follows. We begin in section II by briefly reviewing the UDW detector model, the correlation harvesting protocol, and introduce the conventional measures for entanglement and mutual information. In section III, we respectively explore the phenomena of entanglement harvesting and mutual information harvesting for two static detectors aligned parallel to the boundary, where the influence of the boundary and detectors’ energy gaps on correlation harvesting will be studied in detail, and make a qualitative comparison between the phenomena of mutual information harvesting and entanglement harvesting. In section IV, the correlation harvesting phenomenon for detectors orthogonally aligned to the boundary will be studied, and some comparisons of the results between the parallel-to-boundary and the orthogonal-to-boundary alignments are made. Finally, we end with conclusions in section V. Throughout the paper, the natural units =c=kB=1\hbar=c=k_{B}=1 are adopted for convenience.

II The basic formalism

We consider two detectors AA and BB locally interacting with a massless quantum scalar field ϕ[xD(τ)]\phi[x_{D}(\tau)] (D{A,B}D\in\{A,B\}) along their worldlines. The spacetime trajectory of the detector, xD(τ)x_{D}(\tau), is parametrized by its proper time τ\tau. The detector-field coupling is given by the following interaction Hamiltonian

HD(τ)=λχ(τ)[eiΩDτσ++eiΩDτσ]ϕ[xD(τ)],D{A,B},H_{D}(\tau)=\lambda\chi(\tau)\left[e^{i\Omega_{D}\tau}\sigma^{+}+e^{-i\Omega_{D}\tau}\sigma^{-}\right]\phi\left[x_{D}(\tau)\right],~{}~{}D\in\{A,B\}\;, (1)

where λ\lambda is the coupling strength, and χ(τ)=exp[τ2/(2σ2)]\chi(\tau)=\exp[-{\tau^{2}}/(2\sigma^{2})] is the Gaussian switching function with parameter σ\sigma controlling the duration time of the interaction. In principle, all relevant physical quantities can be rescaled with the duration time parameter σ\sigma to be unitless. Here, we use the standard UDW detector model to describe a two-level system with an energy gap ΩD\Omega_{D} between its ground state |0D|0\rangle_{D} and excited state |1D|1\rangle_{D}. The operators σ+=|1D0|D\sigma^{+}=|1\rangle_{D}\langle 0|_{D} and σ=|0D1|D\sigma^{-}=|0\rangle_{D}\langle 1|_{D} are just the SU(2) ladder operators.

Suppose the two detectors are prepared in their ground state and the field is in the Minkowski vacuum state |0M|0\rangle_{M}, then the initial state of the two detectors and field system is given by |Ψi=|0A|0B|0M|\Psi\rangle_{i}=|0\rangle_{A}|0\rangle_{B}|0\rangle_{M}. The time-evolution of the quantum system can be obtained by using the Hamiltonian (1),

|Ψf:=𝒯exp[i𝑑t(dτAdtHA(τA)+dτBdtHB(τB))]|Ψi,|\Psi\rangle_{f}:={\cal{T}}\exp\Big{[}-i\int{dt}\Big{(}\frac{d\tau_{A}}{dt}H_{A}(\tau_{A})+\frac{d\tau_{B}}{dt}{H_{B}}(\tau_{B})\Big{)}\Big{]}|\Psi\rangle_{i}\;, (2)

where 𝒯{\cal{T}} is the time ordering operator and tt is the coordinate time with respect to which the vacuum state of the field is defined. By tracing out the field degrees of freedom in Eq. (2), the density matrix for the final state of the two detectors turns, to leading order in the interaction strength and in the basis {|0A|0B,|0A|1B,|1A|0B,|1A|1B}\{|0\rangle_{A}|0\rangle_{B},|0\rangle_{A}|1\rangle_{B},|1\rangle_{A}|0\rangle_{B},|1\rangle_{A}|1\rangle_{B}\}, out to be Pozas-Kerstjens:2015 ; Zhjl:2019

ρAB:\displaystyle\rho_{AB}: =trϕ(U|ΨfΨf|U)\displaystyle=\operatorname{tr}_{\phi}\big{(}U\ket{\Psi_{f}}\bra{\Psi_{f}}U^{{\dagger}}\big{)}
=(1PAPB00X0PBC00CPA0X000)+𝒪(λ4),\displaystyle=\begin{pmatrix}1-P_{A}-P_{B}&0&0&X\\ 0&P_{B}&C&0\\ 0&C^{*}&P_{A}&0\\ X^{*}&0&0&0\\ \end{pmatrix}+{\mathcal{O}}(\lambda^{4})\;, (3)

where the detector transition probability reads

PD:=λ2𝑑τ𝑑τχ(τ)χ(τ)eiΩD(ττ)W(xD(t),xD(t))D{A,B},P_{D}:=\lambda^{2}\iint d\tau d\tau^{\prime}\chi(\tau)\chi(\tau^{\prime})e^{-i\Omega_{D}(\tau-\tau^{\prime})}W\left(x_{D}(t),x_{D}(t^{\prime})\right)\quad\quad D\in\{A,B\}\;, (4)

and the quantities CC and XX which characterize correlations of the two detectors, are given by

C\displaystyle C :=λ2𝑑τ𝑑τχ(τ)χ(τ)ei(ΩAτΩBτ)W(xA(t),xB(t)),\displaystyle:=\lambda^{2}\iint d\tau d\tau^{\prime}\chi(\tau)\chi(\tau^{\prime})e^{-i(\Omega_{A}\tau-\Omega_{B}\tau^{\prime})}W\left(x_{A}(t),x_{B}(t^{\prime})\right)\;, (5)
X:=\displaystyle X:= λ2dτdτχ(τ)χ(τ)ei(ΩAτ+ΩBτ)[θ(tt)W(xA(t),xB(t))\displaystyle-\lambda^{2}\iint d\tau d\tau^{\prime}\chi(\tau)\chi(\tau^{\prime})e^{-i(\Omega_{A}\tau+\Omega_{B}\tau^{\prime})}\Big{[}\theta(t^{\prime}-t)W\left(x_{A}(t),x_{B}(t^{\prime})\right)
+θ(tt)W(xB(t),xA(t))].\displaystyle+\theta(t-t^{\prime})W\left(x_{B}(t^{\prime}),x_{A}(t)\right)\Big{]}\;. (6)

Here W(x,x):=0|Mϕ(x)ϕ(x)|0MW(x,x^{\prime}):=\langle 0|_{M}\phi(x)\phi(x^{\prime})|0\rangle_{M} is the Wightman function of the quantum field in the Minkowski vacuum state, and θ(t)\theta(t) is the Heaviside theta function. In particular, if detectors are at rest, one has t=τt=\tau, i.e., the coordinate time of the detector is equal to its proper time. The amount of quantum entanglement can be measured by concurrence WW , which, with the density matrix (II), is given by EDU:2016-1 ; Zhjl:2018 ; Zhjl:2019

𝒞(ρAB)=2max[0,|X|PAPB]+𝒪(λ4).\mathcal{C}(\rho_{AB})=2\max\Big{[}0,|X|-\sqrt{P_{A}P_{B}}\Big{]}+\mathcal{O}(\lambda^{4})\;. (7)

Obviously, the concurrence is a competition between the correlation term XX and geometric mean of the detectors’ transition probabilities PAP_{A} and PBP_{B}. The total amount of correlations is characterized by mutual information, which is defined as Nielsen:2000

(ρAB)=S(ρA)+S(ρB)S(ρAB),\mathcal{I}(\rho_{AB})=S(\rho_{A})+S(\rho_{B})-S(\rho_{AB})\;, (8)

where ρA=trB(ρAB)\rho_{A}=\operatorname{tr}_{B}(\rho_{AB}) and ρB=trA(ρAB)\rho_{B}=\operatorname{tr}_{A}(\rho_{AB}) denote the partial traces of detectors’ state ρAB\rho_{AB}, and S(ρ)=tr(ρlnρ)S(\rho)=-\operatorname{tr}(\rho\ln\rho) is the von Neumann entropy. With the above definition (8), the mutual information for the final detector state (II) takes the following form Pozas-Kerstjens:2015

(ρAB)=\displaystyle\mathcal{I}(\rho_{AB})= +ln(+)+ln()PAln(PA)PBln(PB)+𝒪(λ4)\displaystyle\mathcal{L}_{+}\ln(\mathcal{L}_{+})+\mathcal{L}_{-}\ln(\mathcal{L}_{-})-P_{A}\ln(P_{A})-P_{B}\ln(P_{B})+\mathcal{O}(\lambda^{4}) (9)

with

±\displaystyle\mathcal{L}_{\pm} =12(PA+PB±(PAPB)2+4|C|2).\displaystyle=\frac{1}{2}\Big{(}P_{A}+P_{B}\pm\sqrt{(P_{A}-P_{B})^{2}+4|C|^{2}}\Big{)}\;. (10)

In contrast to the concurrence (7), the mutual information is determined by transition probabilities PAP_{A}, PBP_{B} and the correlation term CC rather than the correlation term XX. It can be verified that mutual information is an increasing function of parameter |C||C|, and monotonically decreases as the transition probabilities increase. According to Eq. (9), we can further obtain (ρAB)=0{\cal{I}}(\rho_{AB})=0 if C=0C=0. Especially, when a transition probability is zero, the mutual information must vanish from the positivity condition of the density matrix, (i.e., PAPB|C|2P_{A}P_{B}\geq|C|^{2}).

In what follows, we will explore, with the harvesting protocol, the correlations harvested by two static detectors near a perfectly reflecting boundary, figuring out the role played by the presence of the boundary in both entanglement harvesting and mutual information harvesting. It is supposed that a perfectly reflecting plane boundary is located at z=0z=0, and two nonidentical UDW detectors with different energy gaps are aligned in two different ways: parallel-to-boundary and orthogonal-to-boundary.

III Correlation Harvesting for the detectors aligned parallel to the boundary

In section, we consider that two static detectors separated by a distance LL are aligned parallel to the boundary with a distance Δz\Delta{z} from the boundary (see Fig. (1)). The spacetime trajectories of the detectors can then be written as

xA:={tA,x=0,y=0,z=Δz},xB:={tB,x=L,y=0,z=Δz}.\displaystyle x_{A}:=\{t_{A}\;,x=0\;,y=0\;,z=\Delta{z}\}\;,~{}~{}x_{B}:=\{t_{B}\;,x=L\;,y=0\;,z=\Delta{z}\}\;. (11)
Refer to caption
Figure 1: Two static detectors separated from each other by a distance LL are aligned parallel to the boundary at a distance Δz\Delta{z} away from the boundary.

Let us now calculate the transition probabilities PDP_{D} given by Eq. (4) and the correlation terms CC and XX respectively given by Eq. (5) and Eq. (II). To do this, we need the Wightman function for vacuum massless scalar fields, which is, according to the method of images, given by Birrell:1984

W(xD,xD)=\displaystyle W\left(x_{D},x^{\prime}_{D}\right)= 14π2[1(ttiϵ)2(xx)2(yy)2(zz)2\displaystyle-\frac{1}{4\pi^{2}}\Big{[}\frac{1}{(t-t^{\prime}-i\epsilon)^{2}-(x-x^{\prime})^{2}-(y-y^{\prime})^{2}-(z-z^{\prime})^{2}}
1(ttiϵ)2(xx)2(yy)2(z+z)2].\displaystyle-\frac{1}{(t-t^{\prime}-i\epsilon)^{2}-(x-x^{\prime})^{2}-(y-y^{\prime})^{2}-(z+z^{\prime})^{2}}\Big{]}\;. (12)

Substituting trajectory (11) and Eq. (III) into Eq. (4), one may straightforwardly obtain the transition probabilities Zhjl:2021

PD=\displaystyle P_{D}= λ24π[eΩD2σ2πΩDσErfc(ΩDσ)]λ2σeΔz2/σ28πΔz\displaystyle\frac{\lambda^{2}}{4\pi}\left[e^{-\Omega_{D}^{2}\sigma^{2}}-\sqrt{\pi}\Omega_{D}\sigma{\rm{Erfc}}(\Omega_{D}\sigma)\right]-\frac{\lambda^{2}\sigma{e}^{-\Delta{z}^{2}/\sigma^{2}}}{8\sqrt{\pi}\Delta{z}}
×{Im[e2iΩDΔzErf(iΔzσ+ΩDσ)]sin(2ΩDΔz)},D{A,B},\displaystyle\times\Big{\{}\mathop{\rm Im}\Big{[}e^{2i\Omega_{D}\Delta{z}}{\rm{Erf}}\Big{(}i\frac{\Delta{z}}{\sigma}+\Omega_{D}\sigma\Big{)}\Big{]}-\sin(2\Omega_{D}\Delta{z})\Big{\}},\qquad D\in\{A,B\}\;, (13)

where Erf(x){\rm{Erf}}(x) is the error function and Erfc(x):=1Erf(x)\operatorname{Erfc}(x):=1-{\rm{Erf}}(x). Similarly, the correlation terms CC and XX in this case can also be worked out

C=λ24πeΔΩ2σ24[f(L)f(L2+4Δz2)],C=\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{\Delta\Omega^{2}\sigma^{2}}{4}}\Big{[}f(L)-f(\sqrt{L^{2}+4\Delta{z}^{2}})\Big{]}\;, (14)
X=λ24πe(2ΩA+ΔΩ)2σ24[g(L)g(L2+4Δz2)],X=-\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{(2\Omega_{A}+\Delta\Omega)^{2}\sigma^{2}}{4}}\Big{[}g(L)-g(\sqrt{L^{2}+4\Delta{z}^{2}})\Big{]}\;, (15)

where the auxiliary functions, ff and gg, are defined by

f(L):=σeL2/(4σ2)L{Im[ei(2ΩA+ΔΩ)L/2Erf(iL+2ΩAσ2+ΔΩσ22σ)]sin[(2ΩA+ΔΩ)L2]},f(L):=\frac{{\sigma}e^{-L^{2}/(4\sigma^{2})}}{L}\Big{\{}\mathop{\rm Im}\Big{[}e^{i(2\Omega_{A}+\Delta\Omega){L}/2}{\rm{Erf}}\Big{(}\frac{iL+2\Omega_{A}\sigma^{2}+\Delta\Omega\sigma^{2}}{2\sigma}\Big{)}\Big{]}-\sin\Big{[}\frac{(2\Omega_{A}+\Delta\Omega)L}{2}\Big{]}\Big{\}}\;, (16)

and

g(L):=σeL2/(4σ2)L{Im[eiΔΩL/2Erf(iL+ΔΩσ22σ)]+icos(ΔΩL2)}.g(L):=\frac{{\sigma}e^{-L^{2}/(4\sigma^{2})}}{L}\Big{\{}\mathop{\rm Im}\Big{[}e^{i\Delta\Omega{L}/2}{\rm{Erf}}\Big{(}\frac{iL+\Delta\Omega\sigma^{2}}{2\sigma}\Big{)}\Big{]}+i\cos\Big{(}\frac{\Delta\Omega{L}}{2}\Big{)}\Big{\}}\;. (17)

Without loss of generality, we here assume the energy gap of the detector BB is not less than that of detector AA, i.e., ΔΩ:=ΩBΩA0\Delta\Omega:=\Omega_{B}-\Omega_{A}\geq 0 throughout the paper. From Eq. (14) and Eq. (15), it is easy to find that as long as the energy gap difference is large with respect to duration time parameter (ΔΩ1/σ\Delta\Omega\gg 1/\sigma), both the correlation terms CC and XX become vanishingly small, and so mutual information and entanglement can hardly be harvested.

III.1 Entanglement harvesting

In general, it is very difficult to capture the characteristics of the correlation harvesting phenomena with the very complicated analytical results above. So, the numerical evaluations are usually called for investigating the influence of the boundary and the energy gaps on correlation harvesting. However, in some special cases, one can still perform analytical approximations to the afore-results. For two identical detectors, i.e., ΩA=ΩB=Ω\Omega_{A}=\Omega_{B}=\Omega, the approximate results for harvested entanglement can be derived in two cases: small energy gap (Ωσ1\Omega\sigma\ll 1) and large energy gap (Ωσ1\Omega\sigma\gg 1).

III.1.1 Small energy gap

When the detectors system is very close the boundary (i.e., Δz/σ1\Delta{z}/\sigma\ll 1), the transition probabilities can be approximated as

PA=PBλ2Δz22πσ2(13πΩσ2),\displaystyle P_{A}=P_{B}\approx\frac{\lambda^{2}\Delta{z}^{2}}{2\pi\sigma^{2}}\Big{(}\frac{1}{3}-\frac{\sqrt{\pi}\Omega\sigma}{2}\Big{)}, (18)

and the correlation term |X||X| as

|X|{λ2Δz2σ2L3π(1Ω2σ2),Lσ1;2λ2Δz2σ2L4π(1Ω2σ2),Lσ1.|X|\approx\left\{\begin{aligned} &\frac{\lambda^{2}\Delta{z}^{2}\sigma}{2L^{3}\sqrt{\pi}}\Big{(}1-\Omega^{2}\sigma^{2}\Big{)},&\frac{L}{\sigma}\ll 1;\\ &\frac{2\lambda^{2}\Delta{z}^{2}\sigma^{2}}{L^{4}\pi}\Big{(}1-\Omega^{2}\sigma^{2}\Big{)},&\frac{L}{\sigma}\gg 1\;.\end{aligned}\right. (19)

The concurrence, denoted here by 𝒞P(ρAB)\mathcal{C}_{P}(\rho_{AB}) can be approximately expressed as

𝒞P(ρAB){λ2Δz2π(σL3+14Lσ13πσ2+Ω2σ),ΔzσLσ1;max[0,λ2Δz23πσ2+4λ2Δz2σ2L4π]=0,Δzσ1,Lσ1.\mathcal{C}_{P}(\rho_{AB})\approx\left\{\begin{aligned} &\frac{\lambda^{2}\Delta{z}^{2}}{\sqrt{\pi}}\Big{(}\frac{\sigma}{L^{3}}+\frac{1}{4L\sigma}-\frac{1}{3\sqrt{\pi}\sigma^{2}}+\frac{\Omega}{2\sigma}\Big{)},&\quad\frac{\Delta{z}}{\sigma}\ll\frac{L}{\sigma}\ll 1\;;\\ &\max\Bigg{[}0,-\frac{\lambda^{2}\Delta{z}^{2}}{3\pi\sigma^{2}}+\frac{4\lambda^{2}\Delta{z}^{2}\sigma^{2}}{L^{4}\pi}\Bigg{]}=0,&\quad\frac{\Delta{z}}{\sigma}\ll 1,\frac{L}{\sigma}\gg 1\;.\end{aligned}\right. (20)

Eq. (20) shows that the concurrence vanishes in the limit of Δz0\Delta{z}\rightarrow 0, i.e, the two detectors are located at the boundary, as well as when L/σ1L/\sigma\gg 1, i.e., the interdetector separation is very large. This means that entanglement harvesting cannot occur in these cases. While for small interdetector separation (L/σ1L/\sigma\ll 1), the harvested entanglement is an increasing function of the energy gap (Ω\Omega) and the detector-to-boundary distance (Δz\Delta{z}). In the far zone to the boundary (Δz/σ1\Delta{z}/\sigma\gg 1), one can also approximately estimate the concurrence

𝒞P(ρAB){λ22π(σL1π+σ22πΔz2+Ωσ),Δzσ1,Lσ1;0,ΔzσLσ1;\mathcal{C}_{P}(\rho_{AB})\approx\left\{\begin{aligned} &\frac{\lambda^{2}}{2\sqrt{\pi}}\Big{(}\frac{\sigma}{L}-\frac{1}{\sqrt{\pi}}+\frac{\sigma^{2}}{2\sqrt{\pi}\Delta{z}^{2}}+\Omega\sigma\Big{)},&\quad\frac{\Delta{z}}{\sigma}\gg 1,\frac{L}{\sigma}\ll 1\;;\\ &0,&\quad\frac{\Delta{z}}{\sigma}\gg\frac{L}{\sigma}\gg 1\;;\end{aligned}\right. (21)

In the case of L/σ1L/\sigma\ll 1, the concurrence is now a decreasing function of Δz\Delta{z}, which is in contrast to the case of the near zone to the boundary where the concurrence is a growing function of Δz\Delta{z} (see Eq. (20)). Therefore, one can infer that the amount of harvested entanglement ought to peak at a certain distance between the detectors and the boundary. This is what has been demonstrated in the numerical evaluations in Ref. Zhjl:2021 .

III.1.2 Large energy gap

When the detectors’ energy gap is much larger than the Heisenberg energy (Ω1/σ\Omega\gg 1/\sigma), one can easily deduce from Eq. (15) that the concurrence should be vanishingly small. For such identical detectors near the boundary with not too large separation (Δz/σ1,L/σ1\Delta{z}/\sigma\ll 1,L/\sigma\ll 1) the concurrence quantifying the harvested entanglement reads approximately,

𝒞P(ρAB)\displaystyle\mathcal{C}_{P}(\rho_{AB})\approx λ2eΩ2σ2Δz2σL3π,\displaystyle\frac{\lambda^{2}e^{-\Omega^{2}\sigma^{2}}\Delta{z}^{2}\sigma}{L^{3}\sqrt{\pi}}\;, (22)

As a result, the harvested entanglement degrade to zero as the energy gap increases.

Refer to caption
(a) ΩAσ=0.10\Omega_{A}\sigma=0.10
Refer to caption
(b) ΩAσ=1.10\Omega_{A}\sigma=1.10
Figure 2: The concurrence 𝒞P{\cal{C}}_{P} as a function of interdetector separation L/σL/\sigma is plotted for ΩAσ=0.10\Omega_{A}\sigma=0.10, ΔΩ/ΩA={0,0.80,10.00}\Delta\Omega/\Omega_{A}=\{0,0.80,10.00\} in (a) and ΩAσ=1.10\Omega_{A}\sigma=1.10, ΔΩ/ΩA={0,0.50,1.10}\Delta\Omega/\Omega_{A}=\{0,0.50,1.10\} in (b) with fixed Δz/σ=1.00\Delta z/\sigma=1.00.

III.1.3 Numerical estimation

In order to gain a better understanding of the entanglement harvesting for nonidentical detectors with different energy gaps in the presence of the reflecting boundary, we now resort to numerical calculation since analytical approximations are hard to obtain in generic cases. We first show the results from our numerical evaluation on the role played by the reflecting boundary on entanglement harvesting in Figs. (2) and (3).

In Fig. (2), the amount of acquired entanglement is plotted as a function of interdetector separation for various detectors’ energy gaps. Obviously, the harvested entanglement monotonically degrade with the increasing interdetector separation no matter whether the two detectors are identical or not. For a fixed distance Δz\Delta{z} from the boundary, the identical detectors usually have an advantage in acquiring more entanglement than nonidentical detectors with unequal energy gaps (e.g., see Fig. (2b)). In contrast, when the energy gaps are small and the interdetector separation is large enough with respect to the duration time parameter (ΩAσ<ΩBσ1\Omega_{A}\sigma<\Omega_{B}\sigma\ll 1 and L>σL>\sigma), the nonidentical detectors are instead more likely to harvest more entanglement [see Fig. (2a)]. Notice that there is only a finite harvesting-achievable range for the interdetector separation for entanglement harvesting, which can be enlarged by a large energy gap (or energy gap difference).

In Fig. (3), we show how the presence of the boundary influences entanglement harvesting by plotting the concurrence as a function of Δz/σ\Delta{z}/\sigma for various detectors’ energy gaps. It is easy to see that the amount of harvested entanglement is zero in the limit of Δz0\Delta{z}\rightarrow 0 and approaches to their corresponding values in flat spacetime without any boundaries as Δz\Delta{z}\rightarrow\infty as expected. No matter whether the two detectors are nonidentical with different energy gaps or identical, the harvested entanglement measure by concurrence has a peak value when the distance Δz\Delta{z} is comparable to the duration time σ\sigma [e.g., see Fig. (3a)]. The exact position of detectors where the harvested entanglement peaks is not sensitive to the energy gap difference. Moreover, the reflecting boundary plays a double-edged role in entanglement harvesting (inhibiting the entanglement harvesting near the boundary or enhancing it when Δz/σ>1\Delta{z}/\sigma>1 as compared to the case without any boundaries), regardless of whether the detectors are identical or not. Furthermore, both the amount of harvested entanglement and the influence of the boundary can be suppressed by increasing detectors’ energy gaps or the energy gap difference between the two detectors.

Refer to caption
(a) ΩAσ=0.10\Omega_{A}\sigma=0.10
Refer to caption
(b) ΩAσ=1.10\Omega_{A}\sigma=1.10
Figure 3: The concurrence is plotted as a function of the distance Δz/σ\Delta z/\sigma for fixed L/σ=0.50L/\sigma=0.50. We assume ΩAσ=0.10\Omega_{A}\sigma=0.10, ΔΩ/ΩA={0,0.80,10.00}\Delta\Omega/\Omega_{A}=\{0,0.80,10.00\} in (a), and ΩAσ=1.10\Omega_{A}\sigma=1.10, ΔΩ/ΩA={0,0.50,1.10}\Delta\Omega/\Omega_{A}=\{0,0.50,1.10\} in (b). All the colored dashed lines correspondingly represent the results in flat spacetime without any boundaries, i.e., the results in the limit of Δz\Delta{z}\rightarrow\infty.

In order to further study how the entanglement harvesting phenomenon for two nonidentical detectors depends on the energy gap difference, we plot concurrence as a function of the energy gap difference in Fig. (4). When the interdetector separation LL is small relative to the duration time σ\sigma, the energy gap difference would generally hinder harvesting of entanglement, i.e., the harvested entanglement would rapidly degrade to zero with the increasing energy gap difference between two nonidentical detectors [see Fig. (4a)]. However, when the interdetector separation grows to comparable to or larger than the interaction duration time parameter, i.e., LσL\gtrsim\sigma, the concurrence is no longer a monotonically decreasing function of the energy gap difference but may initially increase with the increasing energy gap difference, and then reach a maximum value at certain nonzero energy gap difference before degrade to zero as the energy gap difference further increases [see Fig. (4b)]. This suggests that there exist an optimal energy gap difference that maximizes the amount of entanglement harvested, which we denote by ΔΩ𝒞\Delta\Omega_{\mathcal{C}}. It seems that the value of ΔΩ𝒞\Delta\Omega_{\mathcal{C}} depends upon both the interdetector separation and the detector-to-boundary distance. In order to analyze this clearly, we demonstrate how the value of the optimal energy gap difference is influenced by the presence of the boundary for various interdetector separations in Fig. (5). Obviously, the optimal energy gap difference is an increasing function of Δz/σ\Delta{z}/\sigma, i.e., the presence of the boundary reduces the optimal energy gap difference for entanglement harvesting. It is also easy to find the larger the interdetector separation, the larger the optimal value of the energy gap difference, which is consistent with the conclusion in Ref. Zhjl:2022.4 .

Refer to caption
(a) L/σ=0.10L/\sigma=0.10
Refer to caption
(b) L/σ=1.50L/\sigma=1.50
Figure 4: The concurrence as a function of ΔΩ/ΩA\Delta\Omega/\Omega_{A} with ΩAσ=0.10\Omega_{A}\sigma=0.10 for various Δz/σ\Delta{z}/\sigma. We assume L/σ=0.10L/\sigma=0.10 in (a) and L/σ=1.50L/\sigma=1.50 in (b). Here, the dashed line in all plots indicates the results in flat spacetime without any boundaries. There is a probability that nonidentical detectors may harvest more entanglement than the identical ones for large interdetector separation.
Refer to caption
Figure 5: The plot of ΔΩ𝒞/ΩA\Delta{\Omega_{\mathcal{C}}}/\Omega_{A} versus Δz/σ\Delta{z}/\sigma for two detectors aligned parallel to the boundary. Here, we set ΩAσ=0.10\Omega_{A}\sigma=0.10 and L/σ={0.10,3.00,4.00,5.00}L/\sigma=\{0.10,3.00,4.00,5.00\}. It is apparent that ΔΩ𝒞/ΩA\Delta{\Omega_{\mathcal{C}}}/\Omega_{A} is an increasing function of Δz/σ\Delta{z}/\sigma as long as the interdetector separation is not too small with respect to the interaction duration time.

III.2 Mutual information harvesting

Now, we begin to study mutual information harvesting. For identical detectors, the harvested mutual information can also be approximated in two cases: small energy gap and large energy gap.

III.2.1 Small energy gap

When the detectors system is placed near the boundary (i.e., Δz/σ1\Delta{z}/\sigma\ll 1), the correlation term CC is approximately given by

C{λ2Δz22πσ2(13L215σ2πΩσ2),Lσ1;λ2Δz2π(2σ2L4πΩ4σeL24σ2),Lσ1,C\approx\left\{\begin{aligned} &\frac{\lambda^{2}\Delta{z}^{2}}{2\pi\sigma^{2}}\Big{(}\frac{1}{3}-\frac{L^{2}}{15\sigma^{2}}-\frac{\sqrt{\pi}\Omega\sigma}{2}\Big{)},&\frac{L}{\sigma}\ll 1;\\ &\frac{\lambda^{2}\Delta{z}^{2}}{\pi}\Big{(}\frac{2\sigma^{2}}{L^{4}}-\frac{\sqrt{\pi}\Omega}{4\sigma}e^{-\frac{L^{2}}{4\sigma^{2}}}\Big{)},&\frac{L}{\sigma}\gg 1\;,\end{aligned}\right. (23)

then the harvested mutual information can be approximated as

P(ρAB){ln23λ2Δz2πσ2+λ2L2Δz215πσ4ln(Lσ)ln22λ2ΩΔz2πσ,ΔzσLσ1;12λ2Δz2σ6L8π(2+3πΩσ),Δzσ1,Lσ1.\mathcal{I}_{P}(\rho_{AB})\approx\left\{\begin{aligned} &\frac{\ln 2}{3}\frac{\lambda^{2}\Delta{z}^{2}}{\pi\sigma^{2}}+\frac{\lambda^{2}L^{2}\Delta{z}^{2}}{15\pi\sigma^{4}}\ln\Big{(}\frac{L}{\sigma}\Big{)}-\frac{\ln 2}{2}\frac{\lambda^{2}\Omega\Delta{z}^{2}}{\sqrt{\pi}\sigma},&\quad\frac{\Delta{z}}{\sigma}\ll\frac{L}{\sigma}\ll 1\;;\\ &\frac{12\lambda^{2}\Delta{z}^{2}\sigma^{6}}{L^{8}\pi}\big{(}2+3\sqrt{\pi}\Omega\sigma\big{)},&\quad\frac{\Delta{z}}{\sigma}\ll 1,\frac{L}{\sigma}\gg 1\;.\end{aligned}\right. (24)

As can be seen from Eq. (24), the mutual information harvested by identical detectors, in the limit of Δz0\Delta{z}\rightarrow 0, i.e., the detectors are located at the boundary, must vanish. While, as the detector-to-boundary distance Δz\Delta{z} increases, the mutual information grows like Δz2\Delta{z}^{2} in the near zone of the boundary. It is also interesting to note that the harvested mutual information is a decreasing function of the detectors’ energy gap Ω\Omega in the case of L/σ1L/\sigma\ll 1, while it is an increasing function of the energy gap in the case of L/σ1L/\sigma\gg 1. This is similar to the behavior for the mutual information harvesting in flat spacetime without any boundaries Pozas-Kerstjens:2015 .

For Ωσ1\Omega\sigma\ll 1 and Δz/σ1\Delta{z}/\sigma\gg 1, the harvested mutual information can be approximated as,

P(ρAB){λ22π[ln2+L26σ2ln(Lσ)σ22Δz2ln2πΩσln2],Δzσ1,Lσ1;λ2σ4L2π(1L212Δz2+πΩσL2),ΔzσLσ1.\mathcal{I}_{P}(\rho_{AB})\approx\left\{\begin{aligned} &\frac{\lambda^{2}}{2\pi}\Bigg{[}\ln{2}+\frac{L^{2}}{6\sigma^{2}}\ln\Big{(}\frac{L}{\sigma}\Big{)}-\frac{\sigma^{2}}{2\Delta{z}^{2}}\ln{2}-\sqrt{\pi}\Omega\sigma\ln{2}\Bigg{]},&\quad\frac{\Delta{z}}{\sigma}\gg 1,\frac{L}{\sigma}\ll 1\;;\\ &\frac{\lambda^{2}\sigma^{4}}{L^{2}\pi}\Big{(}\frac{1}{L^{2}}-\frac{1}{2\Delta{z}^{2}}+\frac{\sqrt{\pi}\Omega\sigma}{L^{2}}\Big{)},&\quad\frac{\Delta{z}}{\sigma}\gg\frac{L}{\sigma}\gg 1\;.\end{aligned}\right. (25)

Eq. (25) shows that the amount of harvested mutual information is an increasing function of Δz\Delta{z}, and its boundary dependence is of minus Δz2\Delta{z}^{-2}. The harvested mutual information is also a decreasing function of Ω\Omega in the case of L/σ1L/\sigma\ll 1, but becomes an increasing function of Ω\Omega in the case of L/σ1L/\sigma\gg 1. In the limit of Δz\Delta{z}\rightarrow\infty, the Wightman function (III) reduces to that in flat spacetime, and so the harvested correlations (mutual information or entanglement) approach to the corresponding results in flat spacetime without any boundaries as expected.

III.2.2 Large energy gap

When the detectors’ energy gap is much larger than the Heisenberg energy (Ω1/σ\Omega\gg 1/\sigma), the mutual information for two identical detectors near the boundary with not too large separation (Δz/σ1,L/σ1\Delta{z}/\sigma\ll 1,L/\sigma\ll 1) reads approximately

P(ρAB)λ2eΩ2σ2Δz28πΩ6σ10[2Ω2σ4ln2L2ln(Ωσ)].\displaystyle\mathcal{I}_{P}(\rho_{AB})\approx\frac{\lambda^{2}e^{-\Omega^{2}\sigma^{2}}\Delta{z}^{2}}{8\pi\Omega^{6}\sigma^{10}}\Big{[}2\Omega^{2}\sigma^{4}\ln 2-L^{2}\ln(\Omega\sigma)\Big{]}\;. (26)

One can see that the harvested mutual information degrades to zero as the energy gap increases.

III.2.3 Numerical estimation

In more general cases, the mutual information harvesting behavior can be captured by resorting to numerical estimations. In Figs. (6) and (7), the amount of mutual information is plotted as a function of interdetector separation and detectors-to-boundary distance for various detectors’ energy gaps, respectively. Obviously, similar to the entanglement, the harvested mutual information is also a monotonically decreasing function of the interdetector separation no matter whether the two detectors are identical or not. Interestingly, when the energy gaps are small and the interdetector separation is large enough with respect to the duration time parameter (ΩAσ<ΩBσ1\Omega_{A}\sigma<\Omega_{B}\sigma\ll 1 and L>σL>\sigma), the nonidentical detectors with unequal energy gaps are likely to harvest more mutual information than the identical detectors [see Fig. (6a)], contrary to the usual wisdom that the identical detectors may be advantageous in mutual information harvesting. It should be pointed out that, unlike entanglement harvesting, mutual information can be harvested by two static detectors with an arbitrarily large interdetector separation, and physically, this means that the mutual information harvested outside the harvesting-achievable range of entanglement is either classical correlation or nondistillable entanglement Pozas-Kerstjens:2015 .

Fig. (7) demonstrates that in the limit of Δz\Delta{z}\rightarrow\infty the mutual information approaches to its corresponding values in flat spacetime without any boundaries, while on the boundary the mutual information vanishes as expected. In sharp contrast to the entanglement harvesting phenomenon which displays a peak-value behaviour when the distance Δz\Delta{z} is comparable to the duration time σ\sigma [see Fig. (3)], the harvested mutual information is always a monotonically increasing function of Δz/σ\Delta{z}/\sigma, approaching asymptotically to the value without any boundaries. So, this means that the reflecting boundary always plays an inhibiting role in mutual information harvesting in contrast to its double-edged role in entanglement harvesting (inhibiting the entanglement harvesting near the boundary and enhancing it when Δz/σ>1\Delta{z}/\sigma>1). Moreover, both the amount of harvested mutual information and the influence of the boundary can be suppressed by increasing detectors’ energy gaps or the energy gap difference between the two detectors.

Refer to caption
(a) ΩAσ=0.10\Omega_{A}\sigma=0.10
Refer to caption
(b) ΩAσ=1.10\Omega_{A}\sigma=1.10
Figure 6: The plots of mutual information P{\cal{I}}_{P} versus the interdetector separation L/σL/\sigma. Here, we have set ΩAσ=0.10\Omega_{A}\sigma=0.10, ΔΩ/ΩA={0,0.80,10.00}\Delta\Omega/\Omega_{A}=\{0,0.80,10.00\} in (a) and ΩAσ=1.10\Omega_{A}\sigma=1.10, ΔΩ/ΩA={0,0.50,1.10}\Delta\Omega/\Omega_{A}=\{0,0.50,1.10\} in (b) with fixed Δz/σ=1.00\Delta z/\sigma=1.00.
Refer to caption
(a) ΩAσ=0.10\Omega_{A}\sigma=0.10
Refer to caption
(b) ΩAσ=1.10\Omega_{A}\sigma=1.10
Figure 7: The mutual information is plotted as the function of Δz/σ\Delta z/\sigma with ΩAσ=0.10\Omega_{A}\sigma=0.10, ΔΩ/ΩA={0,0.80,10.00}\Delta\Omega/\Omega_{A}=\{0,0.80,10.00\} in (a), and ΩAσ=1.10\Omega_{A}\sigma=1.10, ΔΩ/ΩA={0,0.50,1.10}\Delta\Omega/\Omega_{A}=\{0,0.50,1.10\} in (b). Here, we have set L/σ=0.50L/\sigma=0.50. The corresponding results in flat spacetime without any boundaries are shown as dashed lines.

In order to further investigate the influence of the energy gap difference on mutual information harvesting, we plot mutual information as a function of the energy gap difference in Fig. (8). Similar to entanglement harvesting, the harvested mutual information is also a monotonically decreasing function of the energy gap difference when the interdetector separation is timelike and very small relative to the duration time σ\sigma [see Fig. (8a)]. However, when the interdetector separation is spacelike (LσL\gg\sigma), the mutual information is no longer a monotonically decreasing functions of the energy gap difference but may peak at a certain nonzero energy gap difference and then degrade to zero with the increase of the energy gap difference [see Fig. (8b)], i.e., there also exists an optimal energy gap difference that maximizes the mutual information. We demonstrate how the optimal energy gap difference, denoted here by ΔΩ\Delta{\Omega_{\mathcal{I}}}, depends upon the presence of the boundary in Fig. (9). Analogous to that in the case of entanglement harvesting, the optimal energy gap difference ΔΩ\Delta{\Omega_{\mathcal{I}}} generally increases as the interdetector separation increases. However, the presence of the boundary can either increase or decrease the optimum energy gap difference for mutual information harvesting. Specifically, when the interdetector separation is large with respect to the duration time, the optimum energy gap difference for mutual information harvesting is an increasing function of the detector-to-boundary distance; however, when the interdetector separation grows to too large, the optimum energy gap difference for mutual information harvesting becomes a decreasing function rather than an increasing function of the detector-to-boundary distance, which is quite different from that in entanglement harvesting.

Refer to caption
(a) ΩAσ=0.10\Omega_{A}\sigma=0.10 and L/σ=0.10L/\sigma=0.10
Refer to caption
(b) ΩAσ=0.10\Omega_{A}\sigma=0.10 and L/σ=5.00L/\sigma=5.00
Figure 8: Mutual information is plotted as a function of ΔΩ/ΩA\Delta\Omega/\Omega_{A} for various fixed Δz/σ\Delta{z}/\sigma with ΩAσ=0.10\Omega_{A}\sigma=0.10 and L/σ=0.10L/\sigma=0.10 in (a) and a spacelike interdetector separation L/σ=5.00L/\sigma=5.00 in (b). Here, the dashed line denotes the results in flat spacetime without any boundaries. Plot (b) shows that for large interdetector separation the nonidentical detectors may harvest more mutual information than identical detectors.
Refer to caption
Figure 9: The plot of ΔΩ/ΩA\Delta{\Omega_{\mathcal{I}}}/\Omega_{A} versus Δz/σ\Delta{z}/\sigma with ΩAσ=0.10\Omega_{A}\sigma=0.10, L/σ={0.10,3.00,5.00,7.00}L/\sigma=\{0.10,3.00,5.00,7.00\} for two detectors aligned parallel to the boundary.

It is well known that once there is a reflecting plane boundary in flat spacetime, the isotropy of spacetime would be lost. Therefore, the orientation of two detectors system in various angular alignments with respect to the boundary may have a non-negligible impact on correlation harvesting. So, in the next section, we will consider the situation in which the two detectors are orthogonally aligned to the boundary plane.

IV Correlation Harvesting for the detectors orthogonally aligned to the boundary

The spacetime trajectories of two static detectors orthogonally aligned to the boundary can be written as [see Fig. (10)]

xA:={tA,x=0,y=0,z=Δz},\displaystyle x_{A}:=\{t_{A}\;,x=0\;,y=0\;,z=\Delta{z}\}\;,~{}~{} xB:={tB,x=0,y=0,z=Δz+L}.\displaystyle x_{B}:=\{t_{B}\;,x=0\;,y=0\;,z=\Delta{z}+L\}\;. (27)
Refer to caption
Figure 10: Two static detectors are separated by a distance LL orthogonally aligned to the boundary plane and Δz\Delta{z} is the distance between the boundary and the detector which is closer.

Substituting Eq. (III) and Eq. (27) into Eq. (4), we can see that PAP_{A} is just the expression of Eq. (III), and PBP_{B} can be obtained by replacing Δz\Delta{z} with Δz+L\Delta{z}+L in Eq. (III). Similarly, the correlation terms CC and XX can be calculated out

C=\displaystyle C= λ24πeΔΩ2σ24[f(L)f(L+2Δz)],\displaystyle\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{\Delta\Omega^{2}\sigma^{2}}{4}}\Big{[}f(L)-f(L+2\Delta{z})\Big{]}\;, (28)
X=λ24πe(2ΩA+ΔΩ)2σ24[g(L)g(L+2Δz)],\displaystyle X=-\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{(2\Omega_{A}+\Delta\Omega)^{2}\sigma^{2}}{4}}\Big{[}g(L)-g(L+2\Delta{z})\Big{]}\;, (29)

with auxiliary functions f(L)f(L) and g(L)g(L) defined by Eq. (16) and Eq. (17), respectively. According to Eq. (7) and Eq. (9), the concurrence and harvested mutual information can be directly obtained.

IV.1 Entanglement harvesting

The approximate expressions for the concurrence for identical detectors can be also obtained in some special cases.

IV.1.1 Small energy gap

For small energy gap (Ωσ1\Omega\sigma\ll 1), the concurrence, denoted here by 𝒞V(ρAB)\mathcal{C}_{V}(\rho_{AB}), can be shown after some long algebraic manipulations to take the following approximate form

𝒞V(ρAB){λ2Δzπ(σL2+14σL3πσ2+LΩ2σ),ΔzσLσ1;λ22π(σL1π+σ22πΔz2+Ωσ),Δzσ1,Lσ1;0,Lσ.\mathcal{C}_{V}(\rho_{AB})\approx\left\{\begin{aligned} &\frac{\lambda^{2}\Delta{z}}{\sqrt{\pi}}\Big{(}\frac{\sigma}{L^{2}}+\frac{1}{4\sigma}-\frac{L}{3\sqrt{\pi}\sigma^{2}}+\frac{L\Omega}{2\sigma}\Big{)},&\quad\frac{\Delta{z}}{\sigma}\ll\frac{L}{\sigma}\ll 1\;;\\ &\frac{\lambda^{2}}{2\sqrt{\pi}}\Big{(}\frac{\sigma}{L}-\frac{1}{\sqrt{\pi}}+\frac{\sigma^{2}}{2\sqrt{\pi}\Delta{z}^{2}}+\Omega\sigma\Big{)},&\quad\frac{\Delta{z}}{\sigma}\gg 1,\frac{L}{\sigma}\ll 1\;;\\ &0,&\quad{L}\gg{\sigma}\;.\end{aligned}\right. (30)

Similar to the case of parallel-to-boundary alignment, Eq. (30) demonstrates that the entanglement cannot be extracted at extremely large interdetector separation LL or vanishing distance Δz\Delta{z} from the boundary. Comparing Eq. (30) with Eqs. (20) and (21), we can find for Δz/σ1\Delta{z}/\sigma\ll 1 the harvested entanglement in the orthogonal-to-boundary scenario is also an increasing function of Δz\Delta{z}, but the boundary dependence behaves like Δz\Delta{z} rather than Δz2\Delta{z}^{2}, while for Δz/σ1\Delta{z}/\sigma\gg 1 the approximate form of concurrence in both the orthogonal-to-boundary and parallel-to-boundary alignment looks the same and behaves as a decreasing function of Δz\Delta{z}. This non-monotonicity of a function with respect to Δz\Delta{z} implies that the harvested entanglement in the orthogonal-to-boundary case also possesses a peak behavior. Beside these, one may find that the detectors orthogonally aligned to boundary can harvest comparatively more entanglement than these in parallel-to-boundary alignment for the case of ΔzLσ\Delta{z}\ll L\ll\sigma. This is physically in accordance with our finding the boundary inhibits the entanglement harvesting for the detectors close to it.

IV.1.2 Large energy gap

For large energy gap (Ω1/σ\Omega\gg 1/\sigma), the concurrence in the case of Δz/σ1\Delta{z}/\sigma\ll 1 and L/σ1L/\sigma\ll 1 can be approximated as

𝒞V(ρAB)λ2eΩ2σ2ΔzσL2π.\mathcal{C}_{V}(\rho_{AB})\approx\frac{\lambda^{2}e^{-\Omega^{2}\sigma^{2}}\Delta{z}\sigma}{L^{2}\sqrt{\pi}}\;. (31)

Similar to the case of the parallel-to-boundary alignment, the harvested entanglement in the orthogonal-to-boundary alignment also degrade with the increasing detectors’ energy gap.

IV.1.3 Numerical estimations

From above approximations, it is easy to infer that the influence of the boundary on the entanglement harvesting phenomenon for two detectors orthogonally aligned to the boundary should be similar to that for two detectors aligned parallel to the boundary, i.e., qualitatively, the reflecting boundary should play a double-edge role (inhibiting entanglement harvesting for Δzσ\Delta{z}\ll\sigma and assisting entanglement harvesting for Δz>σ\Delta{z}>\sigma), and moreover, entanglement harvesting possesses a finite harvesting-achievable range regardless of the presence of the boundary. However, the quantitative details should be slightly different. To better understand it, we define the difference of concurrence between the scenarios of orthogonal-to-boundary and parallel-to-boundary alignments: Δ𝒞:=𝒞V𝒞P\Delta\mathcal{C}:=\mathcal{C}_{V}-\mathcal{C}_{P}.

The difference in the amount of harvested entanglement is plotted as a function of the distance from the boundary in Fig. (11). As we can see from Fig. (11), the detectors in orthogonal-to-boundary alignment near the boundary (Δz/σ1\Delta{z}/\sigma\ll 1) could harvest comparatively more entanglement than that in parallel-to-boundary alignment. However, when the detectors are located far from the boundary (Δz/σ1\Delta{z}/\sigma\gg 1) with not too large interdetector separation, this would be reversed, i.e., the parallel-to-boundary alignment turns out to be favorable for entanglement harvesting. This is because the boundary would play a strongly inhibiting role in entanglement harvesting for two detectors placed near the boundary, and the detectors system in orthogonal-to-boundary alignment has a relatively longer effective distance from the boundary and thus less inhibiting effect, resulting in more entanglement harvested in comparison with those in parallel-to-boundary alignment. While when two detectors are far away from the boundary, the assisting role played by the boundary leads to more entanglement harvested by detectors in parallel-to-boundary alignment in comparison to orthogonal-to-boundary alignment.

Refer to caption
(a) ΩAσ=0.50\Omega_{A}\sigma=0.50 and L/σ=0.50L/\sigma=0.50
Refer to caption
(b) ΩAσ=0.50\Omega_{A}\sigma=0.50 and L/σ=3.00L/\sigma=3.00
Figure 11: The concurrence difference Δ𝒞\Delta{\cal{C}} between orthogonal-to-boundary and parallel-to-boundary alignments, is plotted as a function of the distance Δz/σ\Delta z/\sigma. Here, we have fixed L/σ=0.50L/\sigma=0.50 in plot (a) and L/σ=3.00L/\sigma=3.00 in plot (b). In all plots, we have set ΩAσ=0.50\Omega_{A}\sigma=0.50 with ΔΩ/ΩA={0,0.80,3.00}\Delta\Omega/\Omega_{A}=\{0,0.80,3.00\}.

In Fig (12), we demonstrate how the optimal energy gap difference that maximizes the amount of harvested entanglement depends upon the detector-to-boundary distance. It is easy to find that the optimal energy gap difference ΔΩ𝒞\Delta{\Omega_{\mathcal{C}}} in orthogonal-to-boundary alignment is also an increasing function of the detector-to-boundary distance, but the quantitative detail is slightly different from that in the scenario of parallel-to-boundary alignment [comparing Fig (12) and Fig. (5)].

Refer to caption
Figure 12: The plot of ΔΩ𝒞/ΩA\Delta{\Omega_{\mathcal{C}}}/\Omega_{A} versus Δz/σ\Delta{z}/\sigma with ΩAσ=0.10\Omega_{A}\sigma=0.10 and L/σ={0.10,3.00,4.00,5.00}L/\sigma=\{0.10,3.00,4.00,5.00\} for two detectors in orthogonal-to-boundary alignment.

IV.2 Mutual information harvesting

Now we are in a position of estimating the mutual information harvested by detectors orthogonally aligned to the boundary. Similarly, the analytical approximation of mutual information can also be obtained for identical detectors in some special cases.

IV.2.1 Small energy gap

For small energy gap (Ωσ1\Omega\sigma\ll 1), the mutual information denoted here by V(ρAB)\mathcal{I}_{V}(\rho_{AB}) in orthogonal-to-boundary alignment can be approximated as

V(ρAB){λ2(23πΩσ)Δz26πσ2ln(LΔz),ΔzσLσ1;32λ2(1+πΩσ)σ4Δz2L6πln(σΔz),Δzσ1,Lσ1;λ22π[ln2+L26σ2ln(Lσ)σ22Δz2ln2πΩσln2],Δzσ1,Lσ1;λ2σ4L2π(1L212Δz2+πΩσL2),ΔzσLσ1.\mathcal{I}_{V}(\rho_{AB})\approx\left\{\begin{aligned} &\frac{\lambda^{2}\big{(}2-3\sqrt{\pi}\Omega\sigma\big{)}\Delta{z}^{2}}{6\pi\sigma^{2}}\ln\Big{(}\frac{L}{\Delta{z}}\Big{)},&\quad\frac{\Delta{z}}{\sigma}\ll\frac{L}{\sigma}\ll 1\;;\\ &\frac{32\lambda^{2}\big{(}1+\sqrt{\pi}\Omega\sigma\big{)}\sigma^{4}\Delta{z}^{2}}{L^{6}\pi}\ln\Big{(}\frac{\sigma}{\Delta{z}}\Big{)},&\quad\frac{\Delta{z}}{\sigma}\ll 1,\frac{L}{\sigma}\gg 1\;;\\ &\frac{\lambda^{2}}{2\pi}\Bigg{[}\ln{2}+\frac{L^{2}}{6\sigma^{2}}\ln\Big{(}\frac{L}{\sigma}\Big{)}-\frac{\sigma^{2}}{2\Delta{z}^{2}}\ln{2}-\sqrt{\pi}\Omega\sigma\ln{2}\Bigg{]},&\quad\frac{\Delta{z}}{\sigma}\gg 1,\frac{L}{\sigma}\ll 1\;;\\ &\frac{\lambda^{2}\sigma^{4}}{L^{2}\pi}\Big{(}\frac{1}{L^{2}}-\frac{1}{2\Delta{z}^{2}}+\frac{\sqrt{\pi}\Omega\sigma}{L^{2}}\Big{)},&\quad\frac{\Delta{z}}{\sigma}\gg\frac{L}{\sigma}\gg 1\;.\end{aligned}\right. (32)

In comparison with Eq. (24), the boundary dependence of mutual information, in the case Δz/σ1\Delta{z}/\sigma\ll 1, is Δz2ln(σ/Δz)\Delta{z}^{2}\ln{(\sigma/\Delta{z})} in orthogonal-to-boundary alignment, rather than Δz2\Delta{z}^{2} in parallel-to-boundary alignment. Hence, the detectors orthogonally aligned to the boundary would harvest more mutual information than those aligned parallel to the boundary in the near zone of the boundary. Physically, we can understand it as follows. Since the reflecting boundary always inhibits the harvesting of mutual information as previously discovered, the orthogonal-to-boundary alignment has a comparatively longer effective detector-to-boundary distance than the parallel-to-boundary alignment, so the detectors orthogonally aligned to the boundary could harvest more mutual information. When compared with Eq. (25) for Δz/σ1\Delta{z}/\sigma\gg 1, the approximated expressions of Eq. (32) have the same form as that of Eq. (25), which in the limit of Δz\Delta{z}\rightarrow\infty reduce to the result in flat Minkowski spacetime, which is isotropic, as expected.

IV.2.2 Large energy gap

For large energy gap (Ω1/σ\Omega\gg 1/\sigma), the mutual information in the case of Δz/σ1\Delta{z}/\sigma\ll 1 and L/σ1L/\sigma\ll 1 can be approximated as

V(ρAB)λ2eΩ2σ2Δz24πΩ4σ6ln(L/Δz).\mathcal{I}_{V}(\rho_{AB})\approx\frac{\lambda^{2}e^{-\Omega^{2}\sigma^{2}}\Delta{z}^{2}}{4\pi\Omega^{4}\sigma^{6}}\ln(L/\Delta{z})\;. (33)

Obviously, the harvested mutual information in orthogonal-to-boundary alignment will degrade to zero as the energy gap grows.

IV.2.3 Numerical estimations

With the above discussions, it is easy to infer that the influence of the boundary on the mutual information harvesting phenomenon for two detectors orthogonally aligned to the boundary should be similar to that for two detectors aligned parallel to the boundary, i.e., qualitatively, the reflecting boundary should always play an inhibiting role in mutual information harvesting, and moreover mutual information can be harvested at arbitrarily large interdetector separation. However, the quantitative details may in general be slightly different. To better understand it, we define the difference of mutual information between the scenarios of orthogonal-to-boundary and parallel-to-boundary alignments: Δ:=VP\Delta\mathcal{I}:=\mathcal{I}_{V}-\mathcal{I}_{P}.

The difference in the amount of harvested mutual information is plotted, in Fig. (13), as a function of the distance between the boundary and the detector which is closer. As we can see from Fig. (13), the case of orthogonal-to-boundary alignment in general is more favorable for mutual information harvesting than that of parallel-to-boundary alignment. This is because the boundary always plays an inhibiting role in mutual information harvesting and the case of orthogonal-to-boundary alignment has relatively longer effective distance from boundary.

Refer to caption
(a) ΩAσ=0.50\Omega_{A}\sigma=0.50 and L/σ=0.50L/\sigma=0.50
Refer to caption
(b) ΩAσ=0.50\Omega_{A}\sigma=0.50 and L/σ=3.00L/\sigma=3.00
Figure 13: The mutual information difference Δ\Delta{\cal{I}} is plotted as a function of the distance Δz/σ\Delta z/\sigma. Here, we have fixed L/σ=0.50L/\sigma=0.50 in plot (a) and L/σ=3.00L/\sigma=3.00 in plot (b). In all plots, we have set ΩAσ=0.50\Omega_{A}\sigma=0.50 with ΔΩ/ΩA={0,0.80,3.00}\Delta\Omega/\Omega_{A}=\{0,0.80,3.00\}.

In addition, the boundary dependence for the optimal energy gap difference between two nonidentical detectors that maximizes the amount of mutual information is displayed in Fig. (14). Similar to the case of parallel-to-boundary alignment, the optimal energy gap difference ΔΩ\Delta{\Omega_{\mathcal{I}}} in the case of orthogonal-to-boundary alignment can be either an increasing or a decreasing function of detector-to-boundary distance Δz\Delta{z}. Specifically, when the interdetector separation is large with respect to the duration time, the optimum energy gap difference ΔΩ\Delta{\Omega_{\mathcal{I}}} is an increasing function of Δz\Delta{z}; however, when the interdetector separation is too large, ΔΩ\Delta{\Omega_{\mathcal{I}}} becomes a decreasing function of Δz\Delta{z}.

Refer to caption
Figure 14: The plot of ΔΩ/ΩA\Delta{\Omega_{\mathcal{I}}}/\Omega_{A} versus Δz/σ\Delta{z}/\sigma with ΩAσ=0.10\Omega_{A}\sigma=0.10 and L/σ={0.10,3.00,5.00,7.00}L/\sigma=\{0.10,3.00,5.00,7.00\} for two detectors in orthogonal-to-boundary alignment.

V conclusion

We have performed a detailed study on the phenomenon of correlations harvesting by a pair of UDW detectors interacting with vacuum scalar fields in flat spacetime with the presence of a perfectly reflecting boundary. The influences of the energy gaps, the energy gap difference and the presence of the boundary on the mutual information harvesting and entanglement harvesting for nonidentical detectors are examined in the scenarios of parallel-to-boundary and orthogonal-to-boundary alignments.

For the case of parallel-to-boundary alignment, both the harvested entanglement and mutual information decrease as the interdetector separation increases, and they are also exponentially suppressed by increasing the detectors’ energy gaps. It is worth pointing out the entanglement harvesting phenomenon possesses only a finite harvesting-achievable range for interdetector separation, which can be enlarged by increasing the detectors’ energy gaps (or energy gap difference), while mutual information, in contrast, can be extracted at an arbitrary interdetector separation with an infinite harvesting-achievable range. Moreover, when detectors are located near the boundary (Δz0\Delta{z}\rightarrow 0), both mutual information harvesting and entanglement harvesting are inhibited. However, as the distance to the boundary grows to infinity, the harvested mutual information always increases, finally approaching to the corresponding value in flat spacetime without any boundaries. In contrast, the harvested entanglement will surpass the corresponding value in flat spacetime without any boundaries at a certain distance and peak when the distance becomes comparable to the duration time parameter before turn to decrease and ultimately approach to the corresponding value in flat spacetime without any boundaries when the distance becomes infinite, no matter whether the two detectors are identical or not. Hence, the reflecting boundary, in the sense of correlations harvesting, always plays an inhibiting role in mutual information harvesting but a double-edged role in entanglement harvesting in contrast.

When the interdetector separation is small with respect to the duration time, the harvested correlations are a monotonically decreasing function of the energy gap difference no matter whether the boundary exists or not, i.e., identical detectors will be advantageous to extracting correlations as compared to nonidentical detectors with different energy gaps. However, for not too small interdetector separation, the amount of both harvested entanglement and mutual information possess a maximum value at a certain nonzero energy gap difference, suggesting that there is an optimal value of energy gap difference between two nonidentical detectors which renders the detectors to harvest most correlations. The optimal energy gap difference generally increases as the interdetector separation increases. Interestingly, the presence of the boundary always decreases the value of the optimal energy gap difference for entanglement harvesting, while it can either increase or decrease the optimum energy gap difference for mutual information harvesting. To be specific, the optimum energy gap difference for mutual information harvesting may be decreased by the presence of the boundary when the interdetector separation is large with respect to the duration time, and be increased when the interdetector separation grows to too large.

As for the orthogonal-to-boundary alignment case, the influence of energy gaps and the boundary on both entanglement harvesting and mutual information harvesting is qualitatively similar to that for parallel-to-boundary alignment case, and only quantitative details are slightly different. Specifically, the detectors system in orthogonal-to-boundary alignment always harvest comparatively more mutual information than the parallel-to-boundary ones, while they harvest comparatively more entanglement only near the boundary. This can be understood physically as a result of that the detectors system in orthogonal-to-boundary alignment has a longer effective distance to the boundary.

Finally, we have explored in detail only two alignments of the detectors with respect to the boundary, i.e., the parallel-to-boundary and orthogonal-to-boundary alignments, one may wonder what happens when the detectors are misaligned, i.e., when the line drawn through the detectors is neither parallel nor perpendicular to the boundary. In this case, it is easily seen that the effective distance of the detectors system to the boundary would be shorter than that in the orthogonal-to-boundary alignment but longer than the distance in the parallel-to-boundary alignment. Therefore, one may expect that, on one hand, qualitatively, the influence of energy gaps and the boundary on both entanglement harvesting and mutual information harvesting will be similar, and on the other hand, quantitatively, the detectors in the misaligned scenario would always harvest comparatively more/less mutual information than the parallel-to-boundary/orthogonal-to-boundary alignment, while they could harvest comparatively more/less entanglement than the parallel-to-boundary/orthogonal-to-boundary alignment only near the boundary.

Acknowledgements.
This work was supported in part by the NSFC under Grants No.12075084 and No.12175062, and Postgraduate Scientific Research Innovation Project of Hunan Province under Grant No.CX20220507.

Appendix A The derivations of CC and XX in parallel-to-boundary alignment

A.1 The correlation term CC

Letting u=τu=\tau and s=ττs=\tau-\tau^{\prime}, we have, after carrying out the integration with respect to uu in Eq. (5),

C=\displaystyle C= λ2πσeσ2(ΩAΩB)2/4𝑑ses2+2isσ2(ΩA+ΩB)4σ2W(s).\displaystyle\lambda^{2}\sqrt{\pi}\sigma e^{-\sigma^{2}(\Omega_{A}-\Omega_{B})^{2}/4}\int_{-\infty}^{\infty}dse^{-\frac{s^{2}+2is\sigma^{2}(\Omega_{A}+\Omega_{B})}{4\sigma^{2}}}W(s)\;. (34)

Considering the detectors’ trajectories (11), one can obtain the corresponding Wightman function

W(s)=14π2[1(siϵ)2L21(siϵ)2L24z2].\displaystyle W(s)=-\frac{1}{4\pi^{2}}\Big{[}\frac{1}{(s-i\epsilon)^{2}-L^{2}}-\frac{1}{(s-i\epsilon)^{2}-L^{2}-4z^{2}}\Big{]}\;. (35)

Then the correlation term CC takes the form

C=C0Cz,C=C_{0}-C_{z}\;, (36)

with

C0\displaystyle C_{0} =λ2σ4π3/2eΔΩ2σ2/4es2/(4σ2)(siϵ)2L2ei(2ΩA+ΔΩ)s/2𝑑s,\displaystyle=-\frac{\lambda^{2}\sigma}{4\pi^{3/2}}e^{-\Delta\Omega^{2}\sigma^{2}/4}\int_{-\infty}^{\infty}\frac{e^{-s^{2}/(4\sigma^{2})}}{(s-i\epsilon)^{2}-L^{2}}e^{-i(2\Omega_{A}+\Delta\Omega)\cdot{s}/2}ds\;, (37)

and

Cz\displaystyle C_{z} =λ2σ4π3/2eΔΩ2σ2/4es2/(4σ2)(siϵ)2L24z2ei(2ΩA+ΔΩ)s/2𝑑s,\displaystyle=-\frac{\lambda^{2}\sigma}{4\pi^{3/2}}e^{-\Delta\Omega^{2}\sigma^{2}/4}\int_{-\infty}^{\infty}\frac{e^{-s^{2}/(4\sigma^{2})}}{(s-i\epsilon)^{2}-L^{2}-4z^{2}}e^{-i(2\Omega_{A}+\Delta\Omega)\cdot{s}/2}ds\;, (38)

where ΔΩ:=ΩBΩA\Delta\Omega:=\Omega_{B}-\Omega_{A}. Using the Sokhotski formula,

1x±iϵ=𝒫1xiπδ(x),\frac{1}{x\pm{i}\epsilon}={\cal{P}}\frac{1}{x}\mp{i}\pi\delta(x)\;, (39)

we have

C0=\displaystyle C_{0}= λ2σ4π3/2eΔΩ2σ2/4𝒫es2/(4σ2)s2L2ei(2ΩA+ΔΩ)s/2𝑑s\displaystyle-\frac{\lambda^{2}\sigma}{4\pi^{3/2}}e^{-\Delta\Omega^{2}\sigma^{2}/4}{\cal{P}}\int_{-\infty}^{\infty}\frac{e^{-s^{2}/(4\sigma^{2})}}{s^{2}-L^{2}}e^{-i(2\Omega_{A}+\Delta\Omega)\cdot{s}/2}ds
iλ2σ4π1/2eΔΩ2σ2/4es2/(4σ2)ei(2ΩA+ΔΩ)s/2sgn(s)δ(s2L2)𝑑s\displaystyle-\frac{i\lambda^{2}\sigma}{4\pi^{1/2}}e^{-\Delta\Omega^{2}\sigma^{2}/4}\int_{-\infty}^{\infty}e^{-s^{2}/(4\sigma^{2})}e^{-i(2\Omega_{A}+\Delta\Omega)\cdot{s}/2}{\rm{sgn}}(s)\delta(s^{2}-L^{2})ds\;
=\displaystyle= λ2σeΔΩ2σ2/44π3/2𝒫ei(2ΩA+ΔΩ)s/2es2/(4σ2)s2L2𝑑sλ2σeΔΩσ4+L24σ24πLsin[(2ΩA+ΔΩ)L2],\displaystyle-\frac{\lambda^{2}\sigma{e}^{-\Delta\Omega^{2}\sigma^{2}/4}}{4\pi^{3/2}}{\cal{P}}\int_{-\infty}^{\infty}\frac{e^{-i(2\Omega_{A}+\Delta\Omega)\cdot{s}/2}e^{-s^{2}/(4\sigma^{2})}}{s^{2}-L^{2}}ds-\frac{\lambda^{2}\sigma{e}^{-\frac{\Delta\Omega\sigma^{4}+L^{2}}{4\sigma^{2}}}}{4\sqrt{\pi}L}\sin\big{[}\frac{(2\Omega_{A}+\Delta\Omega){L}}{2}\big{]}\;, (40)

where sgn(s)\rm{sgn}(s) represents the sign function. The integration in the last line of Eq. (A.1) can be performed straightforwardly by using the convolution of two functions in Fourier transforms, i.e.,

λ2σ4π3/2eΔΩ2σ2/4𝒫es2/(4σ2)s2L2ei(2ΩA+ΔΩ)s/2𝑑s\displaystyle-\frac{\lambda^{2}\sigma}{4\pi^{3/2}}e^{-\Delta\Omega^{2}\sigma^{2}/4}{\cal{P}}\int_{-\infty}^{\infty}\frac{e^{-s^{2}/(4\sigma^{2})}}{s^{2}-L^{2}}e^{-i(2\Omega_{A}+\Delta\Omega)\cdot{s}/2}ds
=\displaystyle= λ2σ4πLe(L2+ΔΩ2σ4)/(4σ2)Im[ei(2ΩA+ΔΩ)L/2Erf(iL+2ΩAσ2+ΔΩσ22σ)].\displaystyle\frac{\lambda^{2}\sigma}{4\sqrt{\pi}L}e^{-(L^{2}+\Delta\Omega^{2}\sigma^{4})/(4\sigma^{2})}\mathop{\rm Im}\Big{[}e^{i(2\Omega_{A}+\Delta\Omega){L}/2}{\rm{Erf}}\Big{(}\frac{iL+2\Omega_{A}\sigma^{2}+\Delta\Omega\sigma^{2}}{2\sigma}\Big{)}\Big{]}. (41)

Then C0C_{0} can be expressed in the form

C0=λ24πeΔΩ2σ24f(L)C_{0}=\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{\Delta\Omega^{2}\sigma^{2}}{4}}f(L)\; (42)

with

f(L):=σeL2/(4σ2)L{Im[ei(2ΩA+ΔΩ)L/2Erf(iL+2ΩAσ2+ΔΩσ22σ)]sin[(2ΩA+ΔΩ)L2]}.f(L):=\frac{{\sigma}e^{-L^{2}/(4\sigma^{2})}}{L}\Big{\{}\mathop{\rm Im}\Big{[}e^{i(2\Omega_{A}+\Delta\Omega){L}/2}{\rm{Erf}}\Big{(}\frac{iL+2\Omega_{A}\sigma^{2}+\Delta\Omega\sigma^{2}}{2\sigma}\Big{)}\Big{]}-\sin\Big{[}\frac{(2\Omega_{A}+\Delta\Omega)L}{2}\Big{]}\Big{\}}\;. (43)

Comparing the integration in Eq. (37) with that in Eq. (38), we easily obtain the expression of CzC_{z} from C0C_{0} with LL replaced by L2+4Δz2\sqrt{L^{2}+4\Delta{z}^{2}}

Cz=λ24πeΔΩ2σ24f(L2+4Δz2).C_{z}=\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{\Delta\Omega^{2}\sigma^{2}}{4}}f(\sqrt{L^{2}+4\Delta{z}^{2}})\;. (44)

Therefore, combining Eqs. (42) and (44) can verify the validity of Eq. (14).

A.2 The correlation term XX

Upon considering W(xA(τ),xB(τ))=W(xB(τ),xA(τ))W(x_{A}(\tau^{\prime}),x_{B}(\tau))=W(x_{B}(\tau^{\prime}),x_{A}(\tau)), the expression of Eq. (II) can be simplified as

X=\displaystyle X= λ2𝑑uχ(u)χ(us)ei(ΩA+ΩB)u0𝑑s[eiΩAsW(s)+eiΩBsW(s)]\displaystyle-\lambda^{2}\int_{-\infty}^{\infty}du\chi(u)\chi(u-s)e^{-i(\Omega_{A}+\Omega_{B})u}\int_{0}^{\infty}ds\Big{[}e^{i\Omega_{A}s}W(-s)+e^{i\Omega_{B}s}W(-s)\Big{]}
=\displaystyle= λ2πσeσ2(ΩA+ΩB)240𝑑ses24σ2[eisΩAΩB2W(s)+eisΩAΩB2W(s)],\displaystyle-\lambda^{2}\sqrt{\pi}\sigma e^{-\frac{\sigma^{2}(\Omega_{A}+\Omega_{B})^{2}}{4}}\int_{0}^{\infty}dse^{-\frac{s^{2}}{4\sigma^{2}}}\Big{[}e^{is\frac{\Omega_{A}-\Omega_{B}}{2}}W(-s)+e^{-is\frac{\Omega_{A}-\Omega_{B}}{2}}W(-s)\Big{]}\;, (45)

where u=τu=\tau and s=ττs=\tau-\tau^{\prime} are introduced.

Inserting the Wightman function (35) into Eq. (A.2), we have

X=X0Xz,X=X_{0}-X_{z}\;, (46)

with

X0:=λ2σ2π3/2eσ2(2ΩA+ΔΩ)2/40es2/(4σ2)cos(sΔΩ/2)(siϵ)2L2𝑑s,\displaystyle X_{0}:=\frac{\lambda^{2}\sigma}{2\pi^{3/2}}e^{-\sigma^{2}(2\Omega_{A}+\Delta\Omega)^{2}/4}\int_{0}^{\infty}\frac{e^{-s^{2}/(4\sigma^{2})}\cos(s\Delta\Omega/2)}{(-s-i\epsilon)^{2}-L^{2}}ds\;, (47)

and

Xz=\displaystyle X_{z}= λ2σ2π3/2eσ2(2ΩA+ΔΩ)2/40es2/(4σ2)cos(sΔΩ/2)(siϵ)2L24Δz2𝑑s.\displaystyle\frac{\lambda^{2}\sigma}{2\pi^{3/2}}e^{-\sigma^{2}(2\Omega_{A}+\Delta\Omega)^{2}/4}\int_{0}^{\infty}\frac{e^{-s^{2}/(4\sigma^{2})}\cos(s\Delta\Omega/2)}{(-s-i\epsilon)^{2}-L^{2}-4\Delta{z}^{2}}ds\;. (48)

Using the Sokhotski formula (39), we find, after carrying out the fourier transformation, Eq. (47) takes the form Zhjl:2022.4

X0=\displaystyle X_{0}= σλ2eL2/(4σ2)4Lπe(2ΩA+ΔΩ)2σ24{Im[eiΔΩL/2Erf(iL+ΔΩσ22σ)]+icos(ΔΩL2)}\displaystyle-\frac{{\sigma}\lambda^{2}e^{-L^{2}/(4\sigma^{2})}}{4L\sqrt{\pi}}e^{-\frac{(2\Omega_{A}+\Delta\Omega)^{2}\sigma^{2}}{4}}\Big{\{}\mathop{\rm Im}\Big{[}e^{i\Delta\Omega{L}/2}{\rm{Erf}}\Big{(}\frac{iL+\Delta\Omega\sigma^{2}}{2\sigma}\Big{)}\Big{]}+i\cos\Big{(}\frac{\Delta\Omega{L}}{2}\Big{)}\Big{\}}
=\displaystyle= λ24πe(2ΩA+ΔΩ)2σ24g(L)\displaystyle-\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{(2\Omega_{A}+\Delta\Omega)^{2}\sigma^{2}}{4}}g(L)\; (49)

with

g(L):=σeL2/(4σ2)L{Im[eiΔΩL/2Erf(iL+ΔΩσ22σ)]+icos(ΔΩL2)}.g(L):=\frac{{\sigma}e^{-L^{2}/(4\sigma^{2})}}{L}\Big{\{}\mathop{\rm Im}\Big{[}e^{i\Delta\Omega{L}/2}{\rm{Erf}}\Big{(}\frac{iL+\Delta\Omega\sigma^{2}}{2\sigma}\Big{)}\Big{]}+i\cos\Big{(}\frac{\Delta\Omega{L}}{2}\Big{)}\Big{\}}\;. (50)

Similarly, XzX_{z} takes the following form

Xz=λ24πe(2ΩA+ΔΩ)2σ24g(L2+4Δz2).X_{z}=-\frac{\lambda^{2}}{4\sqrt{\pi}}e^{-\frac{(2\Omega_{A}+\Delta\Omega)^{2}\sigma^{2}}{4}}g(\sqrt{L^{2}+4\Delta{z}^{2}})\;. (51)

Combining Eq. (A.2) with Eq. (51) gives Eq. (15).

References