Entangling detectors in anti-de Sitter space

We examine in $(2+1)$-dimensional anti-de Sitter (AdS) space the phenomena of entanglement harvesting - the process in which a pair of detectors (two-level atoms) extract entanglement from a quantum field through local interactions with the field. We begin by reviewing the Unruh-DeWitt detector and its interaction with a real scalar field in the vacuum state, as well as the entanglement harvesting protocol in general. We then examine how the entanglement harvested by a pair of such detectors depends on their spacetime trajectory, separation, spacetime curvature, and boundary conditions satisfied by the field. The harvested entanglement is interpreted as an indicator of field entanglement between the localized regions where the detectors interact with the field, and thus this investigation allows us to probe indirectly the entanglement structure of the AdS vacuum. We find an island of separability for specific values of the detectors' energy gap and separation at intermediate values of the AdS length for which entanglement harvesting is not possible; an analogous phenomena is observed in AdS$_4$, to which we compare and contrast our results. In the process we examine how the transition probability of a single detector, as a proxy for local fluctuations of the field, depends on spacetime curvature, its location in AdS space, and boundary conditions satisfied by the field.


I. INTRODUCTION
The theory of quantum fields propagating on curved spacetimes investigates situations where both quantum and gravitational effects are important, but the gravitational field can be treated classically and the backreaction of the fields on the spacetime can be ignored [1,2]. This area of research has given us insights into the early Universe and black hole physics, leading to valuable clues of what we should expect from a full-fledged quantum theory of gravity.
Physically motivated detector models, such as the Unruh-DeWitt detector [3,4], provide an operational way in which to probe properties of quantum fields on curved spacetimes. These detectors move along a classical trajectory through a given spacetime while their internal degrees of freedom interact locally with a quantum field living on the spacetime. Most often the interaction coupling the detector and field is inspired by the light-matter interaction widely used in the field of quantum optics [5][6][7][8]. By measuring the internal degrees of freedom of a detector, or collection of detectors, after their interaction with the field has ceased, one effectively performs a measurement on the quantum field itself, allowing one to probe properties of the field such as local fluctuations and their correlations.
Entanglement is a particularly interesting property of a quantum field theory that is well suited to being probed by detectors, and is now finding applications in disparate areas of physics such as the study of critical phenomena in condensed matter systems [9][10][11], in the description of non-classical states of light within the field of quantum optics [5,12], and offering an explanation for the origin of black hole entropy [13][14][15]. It was first discovered within the context of algebraic quantum field theory that the vacuum state of a free quantum field is entangled as seen by the local inertial observers in flat spacetime, even if these observers are localized in spacelike-separated regions [16][17][18]. This result is surprising because it suggests that observers can violate a Bell-like inequality by simply observing local vacuum fluctuations of a quantum field -the vacuum is a resource for entanglement.
It was later pointed out that this property of the electromagnetic vacuum could in principle be exploited to generate entanglement between a pair of atoms [19]. Extensive research has since investigated this process of localized detectors (two-level atoms) extracting entanglement/non-local correlations from the vacuum state of a quantum field [8,[20][21][22][23][24][25][26][27][28][29], which has become known as the entanglement harvesting protocol [30]. The harvested entanglement can depend on the properties of the detectors used, however, it also gives an indication of the entanglement structure of the field itself.
A particular interesting application of the entanglement harvesting protocol is to examine how vacuum entanglement is affected by spacetime structure, such as the expansion rate of the Universe [22] and global spacetime topology [25]. However, apart from these investigations, little is known about entanglement harvesting in curved spacetimes. It has been shown that a weakly curved background can either increase or decrease harvested entanglement depending on the choice of vacuum state [31], and recently the entanglement harvesting protocol was investigated in a black hole spacetime [32], illustrating substantive differences differences from the corresponding situation in Minkowski space. The presence of the black hole was found to inhibit entanglement: as two detectors of fixed proper separation approach the black hole, their concurrence decreases to zero at a finite distance away from the horizon, this distance depends on the energy gaps of the detectors, the mass of the black hole, and the curvature of the spacetime.
Motivated by the above, in this article we present a detailed study of entanglement harvesting in anti-de Sitter (AdS) space. Quantum fields on AdS space are of particular interest because of the rapid development of the anti-de Sitter/conformal field theory correspondence which posits a connection between conformal field theories and the bulk geometry of AdS [33]. We note that entanglement harvesting in the bulk may have implications for entanglement entropy in the context of the AdS/CFT correspondence conjecture. If one considers quantum corrections to the well-known Ryu-Takayanagi prescription [34] then the entanglement of quantum fields in the bulk across the Ryu-Takayanagi surface contributes to the entanglement entropy of the boundary CFT [35]. Understanding the nature of entanglement in the bulk is likely to be important for understanding the quantum corrections to the entanglement in the boundary.
We shall investigate how entanglement harvesting depends on the AdS length and the influence of the boundary conditions satisfied by the field at infinity, comparing our results to the flat spacetime counterpart and to the results obtained for the BTZ black hole [32]. For this reason we shall consider (2 + 1)-dimensional AdS spacetime. This lower-dimensional setting has the advantage of yielding insight into the significant physical processes with relative computational ease and efficiency. We find a number of novel features that are particularly visible in the static case, most prominently an island of separability in parameter space at intermediate values of the AdS length, whose origin remains to be understood. This same phenomenon has been observed in a concurrent study in (3+1)-dimensional AdS spacetime [36], and we shall also, where appropriate, compare our results to those obtained in that setting.
We begin in Sec. II by reviewing the basic formalism of the field-detector coupling and describing in general the entanglement harvesting protocol in curved spacetimes. In Sec. III we investigate the influence of the structure of AdS space on the entanglement harvesting protocol, including the AdS length and boundary conditions satisfied by the field. We consider two different arrangements of detectors: a pair of static non-geodesic detectors and a pair of detectors undergoing circular geodesic orbits. This allows for a comparison between entanglement harvesting along geodesic and non-geodesic detector trajectories. In addition the transition probability of a single detector will be derived in all cases. Finally, we conclude in Sec. V by summarizing the results presented and outlining prospects for future research.
Throughout this article we adopt the mostly positive convention for the metric signature (−1, 1, 1) and employ natural units = c = 8G = 1.

II. UNRUH-DEWITT DETECTORS IN ADS3
In this section we introduce the Unruh-DeWitt detector as a simplified model of a two-level atom interacting locally with a quantum field. We review the derivation of the joint state of two such detectors after their interaction with the field has ceased, stating explicitly the matrix elements appearing in this state to leading order in the interaction strength in terms of the vacuum Wightman function. We then introduce scalar field theory on (2+1)-AdS space, stating explicitly the associated vacuum Wightman function.

A. The Unruh-DeWitt detector
Consider a detector moving along the spacetime trajectory x D (τ ), parametrized in terms of the detectors proper time τ , with its ground state |0 D and excited state |1 D separated by an energy gap Ω D . As a simplified model of the light-matter interaction [6,7], suppose the detector couples to a real scalar field φ(x) through the interaction Hamiltonian (1) where λ is the strength of the interaction, χ D (τ ) is a switching function specifying the duration of the detector's interaction with the field, and σ ± denote SU (2) ladder operators acting on the detector's Hilbert space.
Consider now two such detectors, labeled by A and B, which both begin (τ A , τ B → −∞) in their ground state and then couple to the same scalar field, which begins in an appropriately defined vacuum state |0 . In such a scenario the initial state of the detectors and field together is |Ψ i := |0 A |0 B |0 . In the far future (τ A , τ B → ∞), after the interactions between the field and detectors has ceased, the joint state of the detectors and field together is |Ψ f = U |Ψ i , where the evolution operator U is given by where T denotes the time ordering operator and t is the time coordinate with respect to which the vacuum state of the field is defined. The final state of the detectors alone is obtained by tracing out the field degrees of free- which we have expressed in the computational basis is the Wightman function associated with the field, and P D is the probability that a detector has transitioned from its ground to excited state as a result of its interaction with the field; in some contexts this probability is referred to as the response function of the detector [1]. For spatially fixed detectors A and B we have where the quantities γ D will be functions of their spatial location. In what follows we will assume both detectors have the same proper energy gap Ω and switching parameter σ in their own rest frame, i.e, Ω D = Ω A = Ω B = Ω and σ D = σ A = σ B = σ.
To quantify the entanglement among the detectors resulting from their interaction with the field we employ the concurrence as a measure of entanglement [37]. The concurrence evaluated for the final state of the detectors given in Eq. (3) yields [25,38] From Eq. (7) we see that concurrence, and consequently the entanglement harvested, is a competition between two quantities: the off diagonal matrix element |X|, which leads to non-local correlations between the detectors, and the geometric mean of the detector's transition probabilities P A and P B (the "noise"). This process of detectors extracting entanglement from the vacuum through the local interactions described above is known as the entanglement harvesting protocol. Assuming that any interaction between the detectors mediated by the field is negligible, which is certainly the case if the detectors are spacelike separated, the resulting entanglement of the final state of the detectors in Eq. (3) is a consequence of existing entanglement between the regions of the field in which the detectors interacted, some of which has been transferred to the detectors. In this way, the entanglement harvesting protocol, specifically quantifying the entanglement of the state given in Eq. (3) through an entanglement measure like the concurrence, provides an indication of the entanglement structure of the vacuum state. We will exploit this fact in Sec. III to examine the entanglement structure of the AdS vacuum and its dependence on spacetime curvature and boundary conditions satisfied by the field as seen by a pair of detectors.

B. The AdS3 Wightman function
AdS space is a maximally symmetric spacetime of constant negative curvature that solves Einstein's equations with negative cosmological constant Λ := −1/ 2 . The (2+1)-dimensional AdS (AdS 3 ) geometry can be described by a 3-dimensional hyperboloid embedded in a flat 4-dimensional geometry with coordinates (X 1 , X 2 , T 1 , T 2 ) [39]. Global coordinates (τ , ρ, φ) covering AdS 3 are obtained via the transformation in which the induced metric on AdS 3 becomes whereτ ∈ R, ρ ∈ R + , and φ ∈ [0, 2π) [40]. Under the coordinate transformation t := τ and r := sinh ρ, the AdS 3 metric may be cast in the familiar form It will be useful for what follows to introduce the function d(R 1 , R 2 ) denoting the proper distance between the spacetime points x 1 = (t, R 1 , φ) and x 2 = (t, R 2 , φ), We shall consider a massless conformally coupled real scalar field on AdS 3 , since its vacuum Wightman function takes the particularly simple form [39] W (ζ) where is the square distance between x and x in the embedding space R 2,2 and the parameter ζ ∈ {1, 0, −1} specifies either Dirichlet (ζ = 1), transparent (ζ = 0), or Neumann (ζ = −1) boundary conditions satisfied by the field at spatial infinity. Furthermore, it is the vacuum state described by this Wightman function from which the Hartle-Hawking vacuum on the BTZ black hole may be constructed [41]; the BTZ black hole being asymptotically AdS space. Thus our investigations in this paper will allow for a proper comparison between entanglement harvesting between detectors in AdS 3 space and detectors outside the BTZ black hole [32] to investigate the effects of spacetime horizons. Although the general method for computing Wightman functions is to carry out a mode sum over basis functions, the structure of the Wightman function in AdS 3 is simple enough that a combination of analytic and numerical integration is possible. Our methods are complementary to a concurrent study of entanglement harvesting in AdS 4 [36], in which mode sums were carried out to evaluate the relevant quantities of interest.

III. STATIC DETECTORS
In this section we evaluate the transition probability of a static detector at different locations in AdS 3 , and answer the question of how much entanglement two such detectors can harvest from the AdS 3 vacuum.
Suppose detectors A and B are kept at distinct fixed positions R A and R B , respectively, along a common radial direction. The spacetime trajectories describing such detectors parameterized by their proper time τ D is where γ D := (R D / ) 2 + 1. We take the switching functions of the detectors to be Gaussians described by where t 0 > 0 corresponds to detector A interacting with the field before detector B and vice versa for t 0 < 0, where the time delay t 0 is defined in the field frame. Each detector interacts with the field for an approximate amount of proper time σ.
To compute the transition probability P D and matrix element X we evaluate the Wightman function in Eq. (14) along the static detectors' trajectories, substitute the resulting expression into Eqs. (4) and (6), and express the result in a form that lends itself to being evaluated numerically, all of which is done explicitly in Appendix A. The resulting expression for the transition probability is where we have defined The resulting expression for the matrix element X is where In what follows, both the transition probability and matrix element X given above were evaluated numerically in Mathematica to a precision of 10 −17 using the double exponential and double exponential oscillatory integrations methods.

A. The transition probability
We now evaluate the transition probability of a static detector, given in Eq. (18), for a wide range of scenarios. We begin by noting that the transition probability in Eq. (18) can be evaluated perturbatively in σ/ . We find that the leading contributions take a very simple form The first term appearing above is simply the transition probability in (2+1)-dimensional Minkowski space (flat space). The proper distance d(0, R D ) first enters at order (σ/ ) 4 and its appearance can be understood as the manifestation of redshift (time dilation) due to the detector being on a static trajectory.
The perturbative expression in Eq. (22) is particularly useful for understanding how the transition probability in AdS 3 deviates from the flat space result. We see that the leading order correction to the flat space result is proportional to ζ. This means that the most important factor in the deviation from flat spacetime is due to the boundary conditions satisfied by the field. A small amount of negative curvature will result in the detector clicking more (less) than it would in flat spacetime if the field satisfies Neumann (Dirichlet) boundary conditions. If the field satisfies transparent boundary conditions then there is no correction to the flat space result at order σ/ .
The next order correction is also quite interesting. It is independent of the boundary conditions satisfied by the field and so can be thought of as a "universal" contribution due purely to negative curvature. Since this term in the perturbative expansion is always negative for detectors starting in their ground state, we can interpret this as the statement that negative curvature tends to decrease the transition probability of a detector. Conversely, since this correction will always be positive for detectors that are initially excited (Ω D < 0), this means excited detectors in negatively curved spaces are more likely to relax to their ground state than they would be in flat spacetime.
The perturbative expansion, while insightful, cannot be used when σ/ is large, and there are many interesting phenomena in this regime. To explore this regime we must resort to a full numerical evaluation of the integral in Eq. (18) which is used in producing Figs. 1-3.
As seen upon comparison of the plots in Fig. 1, for all boundary conditions, the transition probability P D is only sensitive to the proper distance d(0, R D ) the detector is from the origin when the switching parameter σ is greater than or comparable to the AdS length . As /σ grows, regardless of the boundary conditions, P D approaches its corresponding value in Minkowski space, indicated by the dashed line. Furthermore, for all boundary conditions, the transition probability of a detector located at the origin vanishes as /σ → 0, whereas for a detector positioned some proper distance away from the origin, the transition probability remains almost the same for any . The most notable property of a detector located at the origin is that for a field satisfying Dirichlet boundary conditions (ζ = 1), its transition probability reaches a maximum near /σ ≈ 0.7.
In Fig. 2 the transition probability is plotted as a function of the detector's energy gap Ωσ for fixed AdS length /σ = 1/2, where Ω < 0 corresponds to the detector initially prepared in its excited state. We observe that for Ω > 0, the transition probability decays to zero exponentially with increasing Ωσ without oscillation, regardless of the boundary condition. Moreover, the decay is much faster for the Dirichlet condition (ζ = 1) than for the other boundary conditions. However for Ω < 0, the transition probability oscillates as a function of Ωσ, this oscillation being suppressed by both the transparent boundary condition and for detectors located at large distances d(0, R D ) from the origin. This feature is evidently dependent on dimensionality, and not on the general structure of AdS spacetime, since the analogous graph in (3 + 1)-dimensions exhibits a transition rate (modulo oscillations) that is roughly proportional to |Ω| for large negative energy gap [36]. It is well-known that the detector response exhibits different dependence on Ω as a function of spacetime dimension, and so this difference is not surprising. The dimension dependence arises due  , and (f) the detector is located a proper distance d(0, RD)/σ = 1 away from the origin. Negative energy gaps Ωσ < 0 correspond to a detector that is prepared in its excited state |1 D prior to interacting with the field.
to the short distance behaviour of the Wightman function [42] or, equivalently, due to the energy scaling of the density of states of the field quanta sensed by the detector which goes as Ω d−3 (see, e.g., [43]).
In Fig. 3, we illustrate the behaviour of the transition probability for a detector located at different spatial positions in AdS 3 . The transition probability, regardless of the boundary conditions, is hardly influenced by the change of position if /σ is not vanishingly small. This is to be expected: for large AdS lengths the detector is highly localized in spacetime and so any change in the AdS curvature negligibly affects the detector. It is only when the AdS length is comparable to the width of the switching function that we see a dependence of transition probability of detector B on its proper distance from the origin.

B. Entanglement harvesting with static detectors
We now consider harvesting entanglement with detector A at the origin and detector B placed at various fixed proper distances from A, and for now we assume no relative time delay in the detectors' switching functions (t 0 = 0). The transition probabilities P A and P B and matrix element X can be obtained numerically using Eqs. (18) and (20), after which the concurrence given in Eq. (7) can be easily evaluated as a measure of the resulting entanglement between the detectors. We find that the effects of curvature, spatial separation, and detector energy gap are notably more dramatic for entanglement harvesting as compared to the transition probability, which we shall now demonstrate.
We first consider how the concurrence depends on the AdS length /σ. From Fig. 4 we see that the concurrence vanishes for all boundary conditions as /σ → 0. The concurrence attains its maximum value in the interval 0.5 /σ 1, and this maximum is largest if the field  satisfies Dirichlet boundary conditions (ζ = 1). We also note that as grows, the concurrence asymptotes to the flat space value for all boundary conditions, as expected.
Similar peaking behaviour occurs as the energy gap is changed for fixed . In Fig. 5 we observe a peak in concurrence for positive Ωσ for each boundary condition, the peak again being largest for the Dirichlet boundary condition for the value of /σ considered. For initially excited detectors (Ω < 0) the decrease in concurrence with increasing |Ω| is much more rapid than for detectors initially in their ground states. Hence, similar to what we found for (2 + 1)-dimensional flat space (see Appendix B), it is easier to perform entanglement harvesting using detectors in their ground-state instead of their excited-state.
Next, we analyze the dependence of the concurrence on the proper distance d(0, R A ) detector A is from the origin for a fixed proper separation d(R A , R B ) of the detectors. Figure 6 shows that the large d(0, R A ) behaviour changes depending on the proper separation of the detectors. When the detectors are close together, the concurrence approaches a constant value, but for larger separations, the concurrence rapidly falls to zero with increasing distance from the origin. We also note that the concurrence is maximum when the detectors are close to the origin. In this region they will experience a smaller acceleration and a lower transition probability, which likely leads to the larger value of the concurrence as seen from Eq. (7).
Finally, we consider the behaviour of the concurrence for different proper separations d(R A , R B ) between the two detectors. From the plots shown in Fig. 6 we see that that the larger the value of d(R A , R B ), the less the concurrence, commensurate with analogous results in flat spacetime and consistent with ones expectations from the fall off of the Wightman function in Eq. (14) as the spatial distance between x and x grows. For Dirichlet boundary conditions (ζ = 1) the decay rate is slowest and the complete elimination of entanglement (sudden death) occurs at the largest proper separation, with the opposite true for Neumann boundary conditions (ζ = −1).
Clearly the concurrence exhibits interesting behaviour as a function of AdS length , and the detectors' energy gap Ω and proper separation d(R A , R B ). To illustrate this more fully, we now provide density plots of the concurrence as a function of the detectors' energy gap Ω and proper separation for different boundary conditions in Fig. 7. The shaded areas in each plot are regions in the parameter space where no entanglement harvesting is possible. We do not plot values of /σ > 20 as we have found no appreciable quantitive changes for these values of /σ. We see from Fig. 7 several common features for all boundary conditions over the range of /σ. Regions of large concurrence are always in the lower-left corner (smaller energy gap and smaller detector separation), and regions of zero concurrence are in the lower-right corner (smaller energy gap and larger separation). We note also that entanglement harvesting is generally possible at any detector separation for sufficiently large energy gap, albeit in minuscule amounts compared to smaller energy gaps and detector separation. We also note that the boundary of entanglement in the {d(R A , R B )/σ, Ωσ} parameter space is approximately a straight line for large enough , a feature observed in (3 + 1)-dimensions [36], and also in flat space in all spatial dimensionalities D ≤ 3, noted in [28], and which we have observed in (2 + 1)dimensions (though we have not displayed the graph).
Distinct features arise when we consider the boundary of the zero-entanglement region (shaded region in Fig. 7).  Fig. 7); we will refer to this region as a separabil-ity island. In order to see this separability island clearly, we provide a close-up view of this region for Dirichlet boundary conditions (ζ = 1) in Fig. 8; other boundary conditions yield qualitatively similar plots. We see that the island is in the region 3.4 Ωσ 3.85. This behaviour is due to the matrix element X. Over this region, as d(R A , R B ) increases, X decreases to a minimum and then begins to increase again, while P B remains nearly constant throughout. Considering the expression for X, and recalling that ∆ T = 0 in Eq. (20), specifically the term appearing in the numerator e −a X y 2 cos (β X u), and the cosine factor cos β X in Eq. (21), we can deduce that when |β X | a X and the corresponding X must oscillate. Otherwise, were |β X | a X , then the Gaussian factor e −a X y 2 would play a dominant role in the integration, suppressing the resulting decrease/increase behaviour.
For larger AdS length, /σ 5/2, γ B approaches γ A = 1 and in general |β X | a X , which results in the oscillatory behaviour of X being suppressed. Of course, γ B can increase for detectors with very large proper detector separation d(R A , R B ) so that |β X | a X ; then the factor K X ∝ γ −1/2 B e −Ω 2 σ 2 /2 will result in |X| being extremely small, and tiny fluctuations of |X| will have no impact on the concurrence.
Conversely, when the AdS length is very small, /σ 1, it is easy to deduce that γ B ∝ e d(R A ,R B )/ γ A = 1 and |β X | a X because |β X | ∝ Ω and a X ∝ 2 /(2σ 2 ). Therefore, the result of the integration of X can oscillate as a function of the position of detector B. However, the factor K X then becomes approximately where we have used Eq. (13) in arriving at the above expression. Despite oscillations in X, the above behaviour of K X highly suppresses their intensity, and so such behaviour is negligible. Similar conclusions can also be obtained for transparent (ζ = 0) and Neumann (ζ = −1) boundary conditions. Summarizing, the shape of the separability island in the {d(R A , R B )/σ, Ωσ} parameter space depends sensitively on the AdS length , and vanishes for values of | /σ − 2.5| 1. Outside this region entanglement harvesting is possible because of the influence of K X and a X on the matrix element X. We also note that a similar island has been observed in the same general region of parameter space in (3 + 1)-dimensions [36]. It is clear that this phenomenon merits further study.

C. Entanglement harvesting with a time delay
We now consider the case when the switching functions of the two static detectors are offset by some t 0 = 0 in the coordinate time t. It is clear from the definition of X in Eq. (20) that the concurrence will not be symmetric under the transformation t 0 → −t 0 ; in other words, the amount of entanglement harvested will depend on whether detector A or B interacts with the field first. Indeed, this can be seen in Fig. 9 where we consider the dependence of concurrence on detector separation and switching delay when the detectors' gaps are equal to Ωσ = 2. The concurrence is non-zero for greater detector separation when detector A switches first. This effect is most noticeable in the case of Dirichlet (ζ = 1) boundary conditions. At moderate AdS length ( /σ = 5) we note the formation of two "peninsulas" in the {d(R A , R B )/σ, t 0 /σ)} parameter space, which is largest around the region of small time delay and detector separation. When the AdS length is large ( /σ = 20), the peninsulas are nearly symmetric about t 0 = 0, look very similar for all three boundary conditions, and approach the flat space behaviour. In the limit of large , β X vanishes as −1 so that (∆ T + iβ X ) → ∆ T . Thus in the limit → ∞, X will be even in t 0 , recovering the Minkowski result.
For small AdS length and small detector separation, the concurrence is also largest around the region of small time delay, but for moderate detector separation, the concurrence is much larger when detector A interacts with the field first. Again, this effect is seen for all three boundary conditions, but is most exaggerated for Dirichlet (ζ = 1) boundary conditions. Additionally, in this regime, and when t 0 > 0, we find oscillations dependent on time delay in the concurrence that occur for all three boundary conditions, although they are much more pronounced for Neumann (ζ = −1) and transparent (ζ = 0) boundary conditions where the concurrence goes to zero in the troughs for some parameters. In the case of Neumann boundary conditions (ζ = −1), these oscillations can be seen more clearly in Fig. 10.
It can be seen from the definition of X in Eq. (20) that it will be invariant under the transformation t 0 → −t 0 and Ω → −Ω, however, as depicted in figure 10, the concurrence is not. As illustrated in Fig. 2, the transition probability of a detector initialized in its excited state (Ω < 0) is much higher than when the detector is initialized in its ground state, meaning in the former case, the term √ P A P B dominates over |X| in Eq. (7) and the concurrence is zero.
Finally, in Fig. 11 we consider small energy gaps, and find that the concurrence is nearly symmetric about t 0 = 0 for all three boundary conditions. Here, β X is small so X is approximately even in t 0 (the large y be- haviour of the integrand in Eq. (20) is exponentially suppressed by the Gaussian factor). Unlike in Fig. 9, there are no time delay dependent oscillations, nor is entanglement harvesting possible for moderate or large detector separation regardless of which detector interacts with the field first. However, the concurrence is maximum around the region of small time delay for all boundary conditions. Outside of this region, there are notable differences in behaviour for each boundary condition. For Neumann boundary conditions (ζ = −1) there is little variation in the concurrence as |t 0 | increases. However, for Dirichlet boundary conditions (ζ = 1) there are two regions around |t 0 |/σ = 3 where no entanglement harvesting is possible when the detectors are very close together, but possible again when their separation increases. In the case of transparent boundary conditions (ζ = 0), entanglement harvesting is only possible for three small regions of the parameter space.
Let us emphasize that the time delay asymmetries observed in the harvested entanglement do not, strictly speaking, reflect something specific to (2 + 1)-AdS spacetime. Instead these asymmetries provide another example of the effects of a relative red shift (or relative time dilation) on the ability of atoms to harvest entanglement from the quantum vacuum. Note that the asymmetries can all be traced back to the non-vanishing of β X , which is proportional to γ A − γ B . If one traces the origin of the β X term, it arises from the term exp (−iΩ A τ A − iΩ B τ B ) in the expression for X (6) which does not explicitly depend on any properties of the Wightman function. As a result, we expect that similar asymmetries should arise whenever two detectors have a relative red-shift factor when compared to the frame in which the time delay is measured. Characterizing the differences in these asymmetries between equivalent detector setups in different spacetimes could then provide a probe of the underlying spacetime geometry, topology, etc.
When a time delay is present, the switching functions are modified, yielding a non-zero ∆ T . This is symmetric in t 0 , but the presence of a differential redshift breaks this symmetry. If γ A < γ B and t 0 > 0, then entanglement can be harvested. However if t 0 < 0, a simple transformation of the integrand in Eq. (20) converts the problem into one with Ω < 0 -the problem effectively becomes that of harvesting entanglement from initially excited detectors, for which the √ P A P B term appearing in Eq. (7) is large and the harvested entanglement suppressed. For the case at hand γ A = 1 < γ B and so entanglement harvesting is suppressed for t 0 < 0, the effect diminishing as the AdS length increases, causing γ B → 1. It also diminishes for small gap (Fig. 11), which likewise suppresses the asymmetry.
This asymmetry has also been observed in AdS 4 [36], though the effect is considerably less pronounced. We attribute this to Hugyens' principle being operative in this case, which has the effect of suppressing harvested entanglement within the light cones of the detectors.

IV. DETECTORS ON CIRCULAR GEODESICS
From the analysis presented so far we have seen that redshift effects, manifested through γ A and γ B , can significantly affect entanglement harvesting when considering two static detectors in AdS 3 . While these effects have led to a rich structure, it is certainly worth understanding what happens in the case where there is no relative redshift between the detectors. In this section we consider such a detector configuration, with detector A and detector B in circular geodesic motion about the origin. The trajectories of such detectors are where the angular velocity of the detectors in the coordinate frame (t, r, φ) is 1/ . As can be seen from the detector trajectories, the proper times of each detector are equal and coincide with the coordinate time t. As a consequence the energy gaps of both detectors are equal in the coordinate frame. Additionally, we take the two detectors to have identical switching functions given by Eq. (17).
Since the detectors have the same proper time, their transition probabilities are equal, P A = P B =P D , and whereã := 2 /4σ 2 ; we will used a tilde over quantities to remind us that detectors are moving along a circular geodesic centred at the origin. Note that transition probability of the detectors in this case is identical to that of a static detector located at the origin.
The matrix elementX appearing in the final joint state of these two detectors is where we have defined A detailed derivation of Eqs. (25) and (26) can be found in Appendix A.
Both the transition probabilityP D and matrix element X given above were evaluated numerically in Mathematica to a precision of at least 25 significant digits using the method "DoubleExponential" or "DoubleExpo-nentialOscillatory" depending on the oscillatory nature of the integrand. Mathematica's "PrincipalValue" procedure was implemented to handle the principal value terms appearing in the matrix elements.
We first consider the scenario when detector A is static and located at the origin, detector B orbits around it, and there is no relative time delay in the switching functions (t 0 = 0). In Figs. 12 and 13 we show how the concurrence depends on the energy gap of the detectors Ωσ and AdS length /σ, respectively. In Fig. 12 we see that the concurrence asymptotes to a constant value for large /σ, corresponding to the flat space result; similar behaviour was observed for static detectors in Fig. 4. This is expected because in the limit /σ → ∞, AdS 3 approaches Minkowski space. On the other hand, when the detector energy gap Ωσ is small, we observe a maximum in the concurrence for /σ 1, the precise location and magnitude of which depends on the boundary conditions satisfied by the field. In particular, the maximum is largest for Neumann boundary conditions (ζ = −1). Increasing the energy gap Ωσ results in the broadening of the peak and, for Dirichlet boundary conditions (ζ = 1), its eventual disappearance.
Comparing Figs. 13 and 5, we also see similar behaviour. The concurrence peaks for small positive detector energy gap, Ωσ > 0, and rapidly decays for Ωσ < 0. This behaviour is present regardless of the AdS length; however we find that for very small AdS length, /σ = 1/2, the largest peak corresponds to Neumann boundary conditions (ζ = −1), but for larger AdS lengths /σ ≥ 1, the largest peak is for Dirichlet boundary conditions (ζ = 1). We also find that as the AdS length becomes large, the harvested entanglement for all three boundary conditions approach the flat space results.
In addition, we find that the case where both detectors A and B move in circular geodesic orbits around the origin is mathematically equivalent to the case when detector A is fixed at the origin and detector B is in a circular geodesic orbit around the origin. The coordinates of the detectors only appear in the termα X ; however when the proper distance between the two detectors is kept constant, thenα X = sech d(R A , R B )/ and the concurrence is only dependent on the proper distance between the two detectors. This means that the calculation is insensitive to the proper distance of either detector from the origin, and without loss of generality we can simply consider the scenario where detector B orbits detector A located at the origin.
Having presented results for both static detectors and detectors moving on circular geodesics, it is worth discussing the similarities and differences between the two cases. We plot the difference in the concurrence for these two cases in Fig. 14. First, we note that in the flat space (large ) limit that the static and circular geodesic trajectories both reduce to the same flat space detector configuration. That is, taking the large limit in each case yields two static detectors with proper separation d(R A , R B ) = R B − R A in flat spacetime. Thus, we expect, and indeed observe, that any differences between the static and circular geodesic cases disappear at large AdS length /σ.
At smaller AdS length, there are differences that emerge between the two cases, and we expect that these are mostly due to the lack of redshift effects for detectors moving along circular geodesics. A particularly notable difference is that, at least for small detector energy gaps, detectors on static trajectories harvest more entanglement than detectors on circular geodesic trajectories.
Another significant difference between these two cases is the lack of a separability island for detectors moving on circular geodesics, which was observed for static detectors. For static detectors and particular choices of AdS length, there was a region in the {d(R A , R B )/σ, Ωσ} parameter space (see Fig. 8) where the entanglement vanished. For detectors on circular geodesics there is no such region, suggesting that the island of no entanglement owes its existence largely to the redshift effects present in the static detector case. Instead, for detectors on circular geodesics we observe a "peninsula" of large entanglement in the {Ωσ, /σ} parameter space for small /σ, most prominent in Fig. 12(b) for Neumann boundary conditions (ζ = −1) (clearly displayed in Fig. 15) but also present for all boundary conditions in Fig. 12. We see that, for detectors with small positive energy gaps, the concurrence will vanish and then reappear as the AdS length increases. Finally, we allow the switching functions of the detectors to be offset by some t 0 = 0. Unlike, in the case of static trajectories, definition ofX ensures that the concurrence is symmetric under the transformation t 0 → −t 0 . This is unsurprising, as the two detectors have the same proper time which is equal to the coordinate time. We depict the dependence of the concurrence on detector separation and switching offset in Fig. 16, taking the detectors' energy gaps to be Ωσ = 2, symmetric about the t 0 = 0 axis for all three boundary conditions. We find the concurrence is highest in the region around t 0 = 0 and low-to-moderate detector separation. This maxima is present for all boundary conditions, but most pronounced for Dirichlet boundary conditions (ζ = 1) When the AdS length is small, /σ = 1, we note that there are two secondary maxima present in the case of transparent boundary conditions (ζ = 0) around t 0 /σ ≈ 6.5. Additionally in this regime, entanglement harvesting is possible at large detector separation when there is no relative time delay in the detectors' switching functions, most notably for Dirichlet conditions (ζ = 1). At larger AdS length, /σ = 5 and /σ = 20, entangle- ment harvesting is possible for larger detector separations when the time delay is non-zero, noted by the two "peninsulas" of concurrence in the two rightmost columns in Fig. 16. Again, this effect is present for all three boundary conditions, but most exaggerated for the Dirichlet case (ζ = 1) .
Comparing Fig. 16 to Fig. 9, we note that when the AdS length is large /σ, /σ = 20, the plots look very similar for all boundary conditions as a consequence of both the static and circular trajectories having the same large (flat space) limit. When the AdS length is smaller, the effect of the relative redshift between the detectors is striking. In addition to the asymmetry in the static case, we also find that for small AdS length, when t 0 is tuned correctly, entanglement harvesting is possible for much larger detector separations as compared to the circular geodesic case. Additionally, the time dependent oscillations found when /σ = 1 for static trajectories are no longer present for the circular geodesic trajectories.

V. CONCLUSION
To investigate and quantify field entanglement in an operational manner, one must make measurements of the field by interacting an appropriate measuring apparatus with the field, and then analyze the measurement outcomes for indications of field entanglement. The entanglement harvesting protocol constructs such a measurement: two detectors interact locally with the field and after their interaction become entangled. As the detectors do not interact directly, and assuming any interaction mediated by the field is negligible, the amount of entanglement present in the final state of the detectors is entanglement that has been extracted from the localized regions in which the detectors interacted with the field. In this way the entanglement harvesting protocol can be used to probe the entanglement structure of the vacuum state of a quantum field theory.
In this article we performed a detailed study of both the transition probability of a single detector and the entanglement harvesting protocol for detectors in AdS 3 space interacting with to a real, massless, conformally coupled scalar field beginning in the vacuum state. We examined how this transition probability and protocol depend on the detectors' locations and trajectories in AdS space, spacetime curvature (in AdS 3 the scalar curvature is R = −6/ 2 ), and the boundary conditions satisfied by the field at spatial infinity. The parameter space is a rather large one to explore, and so we considered identical detectors with equal energy gaps and switching widths in their own rest frame. Our results complement the earlier work done in flat spacetime [8,28], and provide a useful benchmark for further investigations of quantum information and detector physics in other asymptotically AdS spacetimes.
Beginning with a study of the transition probability P D of a single detector, we found that it is most sensitive to changes in parameters at small values of the AdS length /σ. For Dirichlet boundary conditions (ζ = 1), the transition probability of a detector at any position approaches zero as /σ → 0, whereas for Neumann boundary conditions (ζ = −1), the transition probability likewise decreases to zero for a detector fixed at the origin, but increases a little with decreasing /σ if the detector is located elsewhere. For transparent boundary conditions (ζ = 0), the transition probability of a detector positioned far away from the origin remains constant as a function of /σ.
For entanglement harvesting, we found that for static detectors there is an optimal AdS length /σ and detector energy gap Ωσ at which the concurrence reaches a maximum value. One unexpected phenomenon is the appearance of "separability islands" for a range of small but finite /σ ≈ 2.5. In this region the detector transition probability P D remains approximately constant as the proper separation of the detectors changes, but the non-local correlation |X| attains a local minimum whose origin remains to be understood.
We also observed a strong effect on the efficacy of entanglement harvesting if the detectors switch on at different times (in other words, the peak of their switching functions is offset by a time t 0 ). For t 0 < 0 entanglement harvesting is significantly suppressed compared to t 0 > 0, an effect due to locally different detector gaps induced by differing location-dependent redshifts. A similar effect is seen in AdS 4 [36] and we expect it to be present whenever the proper times of the detectors differ and a time delay is present in the switching functions.
Finally, we considered a scenario in which both detectors move along circular geodesic orbits about the origin. This has the effect of removing any relative redshift between the two detectors and so constitutes an important comparative setting. We found a number of similarities with entanglement harvesting by static detectors: there is a maximum in the concurrence as a function of both the detector energy gap and the AdS length in both cases, though the quantitative details differ slightly.
The entanglement harvesting protocol provides an operational way in which to probe the entanglement structure of a quantum field. We hope that our investigations here will inspire further studies of entanglement harvesting in other curved spacetimes to better understand how field entanglement depends on spacetime structure. In addition, we hope that connections between entanglement harvesting and other methods used to study entanglement in field theory, such as algebraic and path integral methods, will be made in the near future.

Static detectors
Beginning with the definition of the transition probability P D in Eq. (4), we can express P D in terms of the integration variables u := τ D and s := τ D − τ D and evaluate the integral over u where in the last equality we have expressed the remaining integral in terms of the integration variable ∆t := t − t = s/γ D . Upon substituting the AdS 3 Wightman function given in Eq. (14) into Eq. (A1), it is seen that the transition probability can be expressed as a difference of two terms and Using Eqs. (15) and (10), and the detector's trajectory given in Eq. (16), we may express the denominators appearing in the integrands defining P ± D as where we have made explicit the i dependence indicating the appropriate branch cut [41] and defined α ± D := −(R D / ) 2 ± 1 /γ 2 D ; note that from the definition of γ D below Eq. (16) it is seen that α − D = −1. Let us first express P − D in terms of the dimensionless integration variable y := ∆t/ where β D := γ D Ω and a D := 2 γ 2 D /4σ 2 . Before we conitnue, we note that if we only consider the principal value of the square root, then integrals of the form are infinite, and the Cauchy principle value of integration cannot be applied to correct this. However, P D is calculated using the Wightman function, which is a tempered distribution (i.e. P − D must be finite). In order to correct this, we require that Direct application of Sokhotsky's formula yields the identity which when combined with Eq. (A7) allows for the simplification of P − D to dy e −a D y 2 sin(β D y) sin(y/2) + π n∈Z (−1) n cos(2nπβ D )e −4n 2 π 2 a D .
Turning our attention to P + D , we note that it may also be rewritten in the form Note that |α + D | = 1 for finite , in which case the singularities appearing in the above integrand are integrable and may take → 0. Again, by requiring that the Wightman function is a tempered distribution, we use Eq. (A9), which implies . . . We now evaluate the matrix element X defined in Eq. (6). Taking the switching function to be the Gaussian functions given in Eq. (17), X may be simplified to where in arriving at the last equality we have introduced the integration variables u := t and s := t − t and carried out the integration over u. Upon substituting the spacetime trajectories of the static detectors given in Eq. (16) into the Wightman function in Eq. (14), the denominators become where α ± X := −R A R B / 2 ± 1 /γ A γ B . Using Eqs. (A17) and (14), X simplifies to where we have introduced the integration variable y := s/ and defined Through the methods similar to those used to obtain P + D , the matrix element X can be brought into the form given in Eq. (20).

Detectors on circular geodesics
As discussed in Sec. IV, we evaluate the transition probability of detector B orbiting around the origin on a circular geodesic. Similar to Eq. (A1), the transition probability can be expressed as where in the last equality we have used the integration variable y := ∆t/ and definedã = 2 /4σ 2 . Again, the appropriate branch of the square root function appearing in Eq. (A20) is dictated by the requirement that the Wightman function is a tempered distribution, and so just as in Eq. (A10) we have −1 + cos(y − i ) = i √ 2 sin y/2 − i , 1 + cos(y − i ) = √ 2 cos y/2 − i .
Finally, using Eq. (A11) the transition probability of a detector on a circular geodesic can be brought to the form stated in Eq. (25). We now compute the matrix element X defined in Eq. (6) for the case where detector A and detector B orbit around the origin on a circular geodesic given by Eq. (24). Through the same manipulations leading Eqs. (A16) and (A18),X may be written as where in the last equality we have used the integration variable y := s/ and defined