## Transmission of quantum information through quantum fields in curved spacetimes

Michael Kasprzak 1, ∗ and Erickson Tjoa 2, †

1 Department of Physics, Princeton University, Princeton, New Jersey 08544, USA 2 Max-Planck-Institut für Quantenoptik, Hans-Kopfermann-Straße 1, D-85748 Garching, Germany

We construct a relativistic quantum communication channel between two localized qubit systems, mediated by a relativistic quantum field, that can achieve the theoretical maximum for the quantum capacity in arbitrary curved spacetimes using the Unruh-DeWitt detector formalism. Using techniques from algebraic quantum field theory, we express the quantum capacity of the quantum communication channel purely in terms of the correlation functions of the field and the causal propagator for the wave equation. Consequently, the resulting quantum channel, and hence the quantum capacity, are by construction manifestly covariant, respect the causal structure of spacetime, and are independent of the details of the background geometry, topology, and the choice of Hilbert space (quasifree) representations of the field.

## I. INTRODUCTION

In the spirit of it-from-(qu)bit, quantum information theory has been extensively used to understand generic physical phenomena from quantum many-body physics to high-energy physics, from applied to foundational physics. Relativistic quantum information (RQI) is one of the many attempts to realize this spirit, which thus far focuses on the role of relativity in quantum information-processing tasks, as well as a deeper understanding of fundamental theories such as relativistic quantum field theories (QFT). The role of quantum information theory in relativistic QFTs cannot be understated: a non-exhaustive list includes rigorous study of the entanglement structure, complexity, and measurement theory for QFTs [1-16]; operational approach to particle production due to non-inertial motion and gravitational fields [17-26]; connections between relativistic causality, information-theoretic causality and reference frames [27-31]; quantum information protocols mediated by the relativistic quantum field [32-54]; etc.

[19, 20, 55-59] (see, e.g., [60] for one different approach). Indeed, transmission of classical information has been analyzed by many authors both perturbatively and nonperturbatively and in flat and curved spacetimes using the UDW model [41-43, 49-51, 54, 61]. Transmission of quantum information in the non-perturbative regime has been studied only in flat spacetimes in very specific examples [52, 53]. Note that in all these cases, the communication is encoded through the field amplitude and not, say, through the polarization degrees of freedom of the electromagnetic field, since on its own the polarization is blind to the relativistic features and constraints.

This work is motivated by a somewhat vague, ambitious, and fundamental question:

Problem: What is the information-carrying capacity of quantum fields?

One can think of several ways to answer this question by making the problem statement more precise. For our purposes, we will adopt an operational approach, namely, whether two parties Alice and Bob, each in possession of some quantum-mechanical probes (atoms, qudits, oscillators, etc.) are able to transmit (quantum) information using the relativistic quantum field as a medium - that is, as a relativistic quantum channel mediating the communication between Alice and Bob. A simple model for of such a communication channel is through the Unruh-DeWitt (UDW) particle detector model and its generalizations or variants

∗ mk7592@princeton.edu

† erickson.tjoa@mpq.mpg.de

Here we present the most general formulation for the transmission of quantum information between Alice and Bob, each in possession of a two-level system (a 'qubit detector'), by suitable coupling to a relativistic scalar field in arbitrary curved spacetime. Clearly, our work constitutes a natural extension of the flat spacetime result in [49-53]. However, we will do so by exploiting techniques from algebraic quantum field theory (AQFT) [62-65]. Crucially, what we gain from this will not be a simple rewriting of all the formulae in the literature: we will see that the quality of the relativistic quantum channel mediated by the quantum scalar field can be formulated purely in terms of smeared correlation functions of the field and the detector parameters. Our results make clear what features of the field-mediated communication channel are necessary for perfect transmission of quantum information, and they are by construction manifestly covariant, respect the causal structure of spacetime, and are independent of the details of the background geometry, topology, and (to some extent) the Hilbert space representations of the field. In particular, it automatically works for all quasifree representations that give rise to Fock representations [6265]).

This paper is organized as follows. In Section II we briefly review the scalar field theory in curved spacetimes using the AQFT framework. In Section III we describe our UDW detector setup and review the quantum capacity of a quantum channel. In Section IV we give the protocol that allows the transmission of quantum information through the quantum field. In Section V we

discuss the validity of the channel and study the regimes in which it fails. In Section VI we conclude with some discussions of our results and outlook. We set c = ℏ = 1 and adopt mostly-plus signature for the metric.

## II. SCALAR FIELD THEORY IN CURVED SPACETIME

In this section we review the scalar field theory in curved spacetimes within the AQFT framework to establish notation and convention, adapting the convention given in [58, 66] (for more details and reviews, see also [58, 62-64]). Readers who are familiar with the framework can skip to Section III.

## A. AQFT for scalar field

Consider a Klein-Gordon (KG) scalar field φ in an ( n + 1)-dimensional globally hyperbolic Lorentzian spacetime M with metric tensor g ab , whose equation of motion is given by the KG equation:

$$( - \nabla _ { a } \nabla ^ { a } + m ^ { 2 } + \xi R ) \phi = 0 \,, \quad \quad ( 1 ) \quad \quad \quad$$

where ∇ a is the covariant derivative with respect to the Levi-Civita connection, m is the mass parameter, R is the Ricci scalar curvature, and ξ ≥ 0 controls the nonminimal coupling of φ to R . Global hyperbolicity implies that we have a foliation of spacetime R × Σ and the KG equation is well-posed with initial data given on Cauchy surfaces.

The KG operator ˆ P := -∇ ∇ a a + m 2 + ξR admits a unique advanced and retarded Green's functions G A ( x , y ) , G R ( x , y ) = G A ( y x , ) which satisfies

$$\hat { P } _ { x } G _ { R / A } ( x, y ) = \delta ^ { n + 1 } ( x - y ) \,, \quad \quad ( 2 )$$

where the subscript on ˆ P x means the KG operator is w.r.t. the variable x . We can use these Green's functions to construct the antisymmetric causal propagator (also known as the Pauli-Jordan distribution)

$$E ( x, y ) \coloneqq G _ { R } ( x, y ) - G _ { _ { \Lambda } } ( x, y ) \,. \quad \quad ( 3 )$$

Canonical quantization of the field φ gives rise to a field operator ˆ obeying the canonical commutation relations φ (CCR)

$$[ \hat { \phi } ( x ), \hat { \phi } ( y ) ] = i E ( x, y ) \mathbb { 1 } \,. \quad \quad ( 4 ) \quad \text{fun$$

As ˆ φ is an operator-valued distribution, in the algebraic approach to QFT we consider smeared field operators , i.e., by viewing the field operator as an R -linear map from the test space of smooth compactly supported functions C ∞ 0 ( M ) to a ∗ -algebra A M ( ), i.e., as a map ˆ : φ C ∞ 0 ( M →A M ) ( ) with

$$\hat { \phi } ( f ) \coloneqq \int \mathrm d V \, f ( \mathrm x ) \hat { \phi } ( \mathrm x ) \,. \quad \ \ ( 5 )$$

The algebra of observables A M ( ) is defined to be the collection of sums of products of ˆ ( φ f ) over all functions in C ∞ 0 ( M ) obeying the smeared CCR relation (4), i.e., a CCR algebra [67]. The KG equation is encoded in the fact that we identify ˆ ( φ f ) ≡ 0 if f = ˆ Ph for some h ∈ C ∞ 0 ( M ).

The dynamical content of the field theory is reflected in terms of the solution space of the KG equation. The solution space Sol R ( M ) can be equipped with a symplectic form σ : Sol R ( M × ) Sol R ( M → ) R , defined as

$$\sigma ( \phi _ { 1 }, \phi _ { 2 } ) \coloneqq \int _ { \Sigma _ { t } } \mathrm d \Sigma ^ { a } \left [ \phi _ { 1 } \nabla _ { a } \phi _ { 2 } - \phi _ { 2 } \nabla _ { a } \phi _ { 1 } \right ], \quad ( 6 )$$

where dΣ a = -t a dΣ, -t a is the inward-directed unit normal to the Cauchy surface Σ , t and dΣ = √ h d n x is the induced volume form on Σ t [68, 69]. As is wellknown, this definition is independent of the choice of Cauchy surface. The field operator ˆ ( φ f ) can be expressed as symplectically smeared field operator [62]

$$\hat { \phi } ( f ) \equiv \sigma ( E f, \hat { \phi } ) \,,$$

and the CCR algebra can be written as

$$\Phi _ { \Phi } ^ { ( 1 ) } \Phi _ { \Phi } ^ { ( \Phi ) }, \sigma ( E g, \hat { \phi } ) ] = i \sigma ( E f, E g ) \mathbf 1 = i E ( f, g ) \mathbf 1 \,, \ ( 8 )$$

where σ Ef,Eg ( ) = E f, g ( ) in the second equality follows from Eq. (5) and (7).

Since the smeared field operators are unbounded operators, for free fields it is more convenient technically to work the Weyl algebra W M ( ), whose elements are bounded operators constructed by 'exponentiating' the field operators. The Weyl algebra W M ( ) is a unital C ∗ -algebra generated by elements that formally take the form

$$W ( E f ) \equiv e ^ { i \hat { \phi } ( f ) } \,, \ \ f \in C _ { 0 } ^ { \infty } ( \mathcal { M } ) \,. \quad \ \ ( 9 )$$

These elements satisfy Weyl relations :

$$\dot { \ } _ { \ } W ( E f ) ^ { \dagger } & = W ( - E f ) \,, \\ W ( E ( \hat { P } f ) ) & = \mathbf 1 \,, \\ W ( E f ) W ( E g ) & = e ^ { - \frac { \dagger } { 2 } E ( f, g ) } W ( E ( f + g ) )$$

where f, g ∈ C ∞ 0 ( M ). The third relation enforces relativistic causality (or microcausality ).

In AQFT, a quantum state is defined as a C -linear functional ω : W M → ( ) C (similarly for A M ( )) such that

$$\omega ( \mathbf 1 ) = 1 \,, \ \omega ( A ^ { \dagger } A ) \geq 0 \ \forall A \in \mathcal { W } ( \mathcal { M } ) \,. \quad ( 1 0 )$$

The state ω is pure if it cannot be written as ω = αω 1 + (1 -α ω ) 2 for any α ∈ (0 , 1) and any two algebraic states ω , ω 1 2 ; otherwise we say that the state is mixed.

The way to pass from AQFT to the Hilbert space approach is through Gelfand-Naimark-Segal (GNS) reconstruction theorem [62-64]. The GNS theorem gives us, for a given algebra of observables W M ( ) and a state ω , a

GNS triple ( H ω , π ω , | Ω ω 〉 ), where π ω : W M →B H ( ) ( ω ) is a Hilbert space representation with respect to state ω . In its GNS representation, any state ω is realized as a vector state | Ω ω 〉 ∈ H ω and A ∈ W M ( ) are represented as bounded operators ˆ A := π ω ( A ) ∈ B H ( ω ), and we write ω A ( ) = 〈 Ω ω | ˆ Ω A | ω 〉 . Since there exists infinitely many unitarily inequivalent representations of the CCR algebra, the algebraic approach allows us work with all representations at once and only pick the physically relevant representation at the very end.

For a fixed state ω , we can construct the Wightman n -point functions , defined by

$$\text{W} ( f _ { 1 }, \dots, f _ { n } ) \coloneqq \omega ( \hat { \phi } ( f _ { 1 } ) \dots \hat { \phi } ( f _ { n } ) ) \quad \text{$\text{$\text{$\text{$i$}$} }$} ) \quad \text{$\text{$\text{$\text{$i$} } } ) \text{$\text{$\text{$\text{$i$} } } )$$

where f j ∈ C ∞ 0 ( M ) and the RHS is computed within some GNS representation of A M ( ). The GNS representation of the Weyl algebra W M ( ) allows us to calculate Eq. (11) by differentiation 1 : for example, the smeared Wightman two-point function reads

$$W ( f, g ) \equiv - \frac { \partial ^ { 2 } } { \partial s \partial t } \Big | _ { s, t = 0 } ^ { \omega ( e ^ { i \hat { \phi } ( s f ) } e ^ { i \hat { \phi } ( t g ) } ) } \quad ( 1 2 ) \quad \text{In a} \quad \text{the}$$

where the RHS is calculated in the GNS representation of W M ( ) [63].

The relevance of quasifree states lies in the fact that they are completely specified once we know the Wightman two-point functions. More precisely, any given quasifree state ω is associated with a real-bilinear inner product µ : Sol R ( M × ) Sol R ( M → ) R that satisfies [65]

The consensus within the AQFT community is that physically reasonable states should be Hadamard states [65, 71]. Very roughly speaking, these states respect local flatness and the expectation values of all observables (in particular, the renormalized stress-energy tensor) are finite [65]. A particularly nice subclass of Hadamard states is the family of quasifree states : for these states, all odd-point functions in the sense of (11) vanish and all higher even-point functions can be written as in terms of just two-point functions. In modern terminology, the term quasifree state is synonymous with Gaussian state , where the one-point functions need not vanish and higher-point functions only depend on one- and twopoint functions [70].

$$| \sigma ( E f, E g ) | ^ { 2 } \leq \mu ( E f, E f ) \mu ( E g, E g ) \,, \quad ( 1 3 ) \quad \text{where}.$$

where we recall that Ef,Eg ∈ Sol R ( M ) for all f, g ∈ C ∞ 0 ( M ). The inequality is saturated if ω µ is a pure quasifree state. Any quasifree state can then defined by those that satisfy

$$\omega ( W ( E f ) ) \coloneqq e ^ { - \mu ( E f, E f ) / 2 } \,. \quad \ \ ( 1 4 ) \quad \text{flat} \\ \circ w + \colon$$

1 Strictly speaking, the state on A M ( ) and on W M ( ) are different, but they are in one-to-one correspondence whenever they are related by Eq. (12) (see [70] for good exposition on this).

However, this definition is not useful unless we can compute the induced norm || Ef || := √ µ Ef,Ef ( ).

It turns out that || Ef || is related to the Wightman two-point function as follows [65]: we complexify Sol R ( M ) into the space of complex solutions Sol C ( M ) and define the KG bilinear product by

$$\langle \varphi _ { 1 }, \varphi _ { 2 } \rangle _ { \text{KG} } \coloneqq i \sigma ( \varphi _ { 1 } ^ { * }, \varphi _ { 2 } ) \text{ \quad \ \ } ( 1 5 )$$

for ϕ , ϕ 1 2 ∈ Sol C ( M ). It can be shown that for any quasifree state, there exists a subspace H ⊂ Sol C ( M ) such that ( H 〈· ·〉 , , KG ) is a Hilbert space and an R -linear map K : Sol R ( M → H ) such that for all ϕ , ϕ 1 2 ∈ Sol R ( M )

- (a) K Sol R ( M ) + iK Sol R ( M ) is dense in H ;
- (b) µ ϕ , ϕ ( 1 2 ) = Re 〈 Kϕ ,Kϕ 1 2 〉 KG ;
- (c) σ ϕ , ϕ ( 1 2 ) = 2Im 〈 Kϕ ,Kϕ 1 2 〉 KG ;
- (d) Sol C ( M ∼ H ⊕ H ) = , where H is the complexconjugate Hilbert space of H and 〈 u, v 〉 KG = 0 for all u ∈ H and v ∈ H .

In a more familiar language of canonical quantization, the map K projects to the 'positive-frequency part' of a real solution to the KG equation. The Wightman twopoint function is then given by [65]

$$\mathbb { m } \text{sounds in} \, \sigma \, \sigma \, \sigma \, \sigma \, \{ \sigma \} \\ \mathbb { W } ( f, g ) & = \langle K E f, K E g \rangle _ { \mathbb { K } G } \\ & = \mu ( E f, E g ) + \frac { i } { 2 } E ( f, g ) \,, \quad \ \ ( 1 6 ) \\ \mathbb { m } \text{sounds in} \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sigma \, \sig \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m n } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m } \, \mathcal { m m }$$

where we have used the fact that σ Ef,Eg ( ) = E f, g ( ). Since E f, g ( ) is antisymmetric, it follows that || Ef || 2 = µ Ef,Ef ( ) ≡ W ( f, f ) and hence

$$\omega ( W ( E f ) ) = e ^ { - \frac { 1 } { 2 } W ( f, f ) } \,. \quad \quad ( 1 7 )$$

Indeed, we may very well take Eq. (17) as the definition of quasifree states.

The most important example of a quasifree state is the 'vacuum state' ω , where we can write the (unsmeared) vacuum Wightman function as

$$\mathfrak { r } _ { t } \quad \omega ( \hat { \phi } ( x ) \hat { \phi } ( y ) ) \equiv \mathsf W ( x, y ) = \, \int d ^ { n } k \, u _ { k } ( x ) u _ { k } ^ { * } ( y ) \,, \quad ( 1 8 )$$

where u k ( x ) are called the positive-frequency modes of KG operator with respect to the KG inner product

$$( \phi _ { 1 }, \phi _ { 2 } ) _ { \text{KG} } \coloneqq i \sigma ( \phi _ { 1 } ^ { * }, \phi _ { 2 } ) \,, \quad \quad ( 1 9 )$$

where φ j ∈ Sol C ( M ) are complexified solutions to Eq. (1). The vacuum state ω is not unique and there is a sense in which every GNS vector is a 'vacuum state' in that particular representation. However, once a particular basis of modes { u k } are chosen, it fixes the vacuum state as one whose Wightman function satisfies (18). In flat spacetime, the standard plane-wave basis in the inertial coordinates defines the Minkowski vacuum, which is the only state that respects the full Poincaré symmetry and minimizes the expectation value of the (renormalized) stress-energy tensor [63]. In curved spacetimes, there are multiple inequivalent states that qualify as a vacuum state [62].

## B. Relationship with canonical quantization

The usual Fock space in canonical quantization arises from the GNS reconstruction applied to A M ( ) with respect to some vacuum state ω . Technically speaking, any pure quasifree (Hadamard) state qualifies as a vacuum state, as their GNS representation gives rise to the Fock vacuum in the standard sense [65] (see also [63] for a different definition in Minkowski spacetime). This is the best one can do in generic curved spacetimes as there is no preferred vacua without further assumptions such as time translation symmetry.

In more detail, the GNS Hilbert space H ω for the pure quasifree representation is the Fock space over the one-particle Hilbert space H

$$\mathcal { H } _ { \omega } \equiv \mathfrak { F } ( \mathcal { H } ) = \bigoplus _ { n = 0 } ^ { \infty } \mathcal { H } ^ { \odot n } \,, \quad \quad ( 2 0 ) \quad \quad \,,$$

where H /circledot 0 ∼ = C and H /circledot n means symmetrized direct sum (the n -particle sector of the Fock space). In this representation, we can write

$$\pi _ { \omega } ( \hat { \phi } ( f ) ) = \hat { a } ( ( K E f ) ^ { * } ) + \hat { a } ^ { \dagger } ( K E f ) \,. \quad ( 2 1 ) \quad \text{where}$$

In what follows we drop π ω if it is clear from the context that we are using a Fock representation. Note that if we consider complex smearing functions f : M→ C , we can write

$$\hat { \phi } ( f ) \equiv \hat { \phi } ( \text{Re} f ) + i \hat { \phi } ( \text{Im} f ) \,. \quad \quad ( 2 2 )$$

The operators ˆ( a u ∗ ) , a ˆ ( † v ) are smeared ladder operators 2 obeying the CCR

$$[ \hat { a } ( u ^ { * } ), \hat { a } ^ { \dagger } ( v ) ] = \langle u, v \rangle _ { K G } \, \mathbb { 1 } \, \quad \ \ ( 2 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

on a suitable dense domain of the Fock space.

The standard canonical approach is recovered by working with the unsmeared field operator ˆ ( φ x ) and considering the Fourier mode decomposition

$$\hat { \phi } ( x ) = \int \mathrm d ^ { n } k \, \hat { a } _ { k } u _ { k } ( x ) + \hat { a } _ { k } ^ { \dagger } u _ { k } ^ { * } ( x ) \,, \quad ( 2 4 ) \quad \text{Ess} \,.$$

where u k ( x ) and u ∗ k ( x ) are positiveand negativefrequency modes normalized to Dirac delta functions via the KG inner product:

$$\langle u _ { k }, u _ { k ^ { \prime } } \rangle _ { K G } = \delta ^ { n } ( k - k ^ { \prime } ) \,, \quad \quad ( 2 5 a )$$

$$\stackrel { \theta ^ { \prime } } { \langle u _ { k }, u _ { k ^ { \prime } } ^ { * } \rangle _ { K G } } = 0 \,.$$

$$\stackrel { \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \} } { \langle u _ { k } ^ { * }, u _ { k ^ { \prime } } ^ { * } \rangle _ { K G } = - \delta ^ { n } ( k - k ^ { \prime } ) \,, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

2 There are various conventions on how to label the smeared ladder operators (see [42]). We follow [62] in that it views ˆ a, a ˆ † as being linear in their arguments. The convention in [63] does not include complex conjugation in the argument for ˆ( ) and a · views ˆ as an anti-linear map while ˆ a a † is linear.

These modes are not proper elements of the one-particle Hilbert space but they are convenient to work with. In particular, we see that ˆ a k ≡ ˆ( a u ∗ k ) and ˆ a † k ≡ ˆ ( a † u k ). If we were to pick any other quasifree state, the resulting one-particle structure will differ and this can give rise to Hilbert space representations that are unitarily inequivalent, e.g., by considering Kubo-Martin-Schwinger (KMS) thermal states [63, 65, 72, 73].

Using these modes, we can decompose the solution Ef in the 'eigenmode basis' { u , u k ∗ k }

$$\underset { \text{head} } { \text{ns} } \underset { \text{E} } { E f } \equiv \int \mathrm d ^ { n } k \, \langle u _ { k }, E f \rangle _ { \text{KG} } u _ { k } + \langle u _ { k } ^ { * }, E f \rangle _ { \text{KG} } u _ { k } ^ { * } \,. \quad ( 2 6 )$$

From the CCR and the definition of the Wightman function, we see that

$$i E ( x, x ^ { \prime } ) = \int \mathrm d ^ { n } k \, u _ { k } ( x ) u _ { k } ^ { * } ( x ^ { \prime } ) - u _ { k } ^ { * } ( x ) u _ { k } ( x ^ { \prime } ) \,, \quad ( 2 7 )$$

the positive-frequency part KEf of Ef reads

$$\stackrel { \lambda } { \text{is} } & & K E f = \int \mathrm d ^ { n } k \, \langle u _ { k }, E f \rangle _ { \text{KG} } u _ { k } \equiv \int \mathrm d ^ { n } k \, f _ { k } u _ { k } \quad ( 2 8 )$$

where

$$f _ { k } \coloneqq - i \int \mathrm d V \, f ( \times ) u _ { k } ^ { * } ( \times ) \,. \quad \quad ( 2 9 )$$

Furthermore, since KEf lies in H , it is convenient to use the notation | KEf 〉 ∈ H and | u k 〉 for the (improper) basis u k . In this notation, we write

$$\left | K E f \right \rangle = \int \mathrm d ^ { n } k \, f _ { k } \left | u _ { k } \right \rangle \,, \quad \quad ( 3 0 )$$

This suggests that we can label the ladder operators using the 'momentum space smearing function' f k rather than the 'real space positive-frequency solution' KEf , i.e., we can write

$$\overset { \mathcal { Y } } { \iota } _ { \iota } \quad \hat { \phi } ( f ) & = \hat { a } ( f _ { k } ^ { * } ) + \hat { a } ^ { \dagger } ( f _ { k } ) \equiv \int \mathrm d ^ { n } k \, \left ( \hat { a } _ { k } f _ { k } ^ { * } + \hat { a } _ { k } ^ { \dagger } f _ { k } \right ) \\ & \equiv \hat { \phi } ( f _ { k } ) \,.$$

Essentially, we are using the elements of the one-particle Hilbert space H to label the field observables rather than using the spacetime smearing functions C ∞ 0 ( M ). This notation is particularly convenient if the focus is more about which modes of the field that the detectors are coupled to [59, 74-76].

## III. SETUP

In this section, we introduce the UDW detector model that we will use to couple Alice's and Bob's qubits to the field and briefly review the notion of quantum capacity for a quantum channel and coherent information. These are relevant measures of how much quantum information can be transmitted (on average) through the quantum channel between Alice and Bob.

## A. UDW detector model

In the UDW framework, Alice and Bob carry their own qubit ('detector') with a free Hamiltonian ˆ H ν = 1 2 Ω ˆ , where ν σ z ν ∈ { A,B } , Ω ν the energy gap and ˆ σ z the PauliZ operator. We denote the eigenstates of ˆ σ s by |± 〉 s with s ∈ { x, y, z } , where we do not label the Pauli matrices with ν to reduce notational clutter and the context will make it clear. The excited and ground states of the free Hamiltonian ˆ H ν are given by |± 〉 z with eigenvalues ± Ω ν / 2 respectively.

The UDW interaction between Alice's and Bob's detectors and the field is assumed to be linear in the field as it represents a simplified model of dipole interaction ˆ d E · ˆ in quantum optics, i.e., as a spin-boson type interaction [59]. The general form of such an interaction is given by a linear combination of the form 3

$$\hat { h } _ { I, \nu } ( x ) \coloneqq \sum _ { i } \lambda _ { \nu, i } f _ { \nu, i } ( x ) \hat { m } _ { \nu, i } ( \tau ( x ) ) \otimes \hat { \phi } ( x ) \,, \quad ( 3 2 ) \quad \text{ \quad \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \text{ \ \ } \end{cases}$$

where λ ν,i is the coupling strength, ˆ m ν,i ( τ ) is a Hermitian operator acting on the qubit, and f ν,i ( x ) is a spacetime smearing function defining the duration of the interaction and the spatial profile of the detector. The most commonly used UDW model involves only one term in ˆ h I,ν ( x ) = λ f ν ν ( x ) ˆ m ν ( τ ( x )) ⊗ ˆ ( φ x ) [55], but we will need the more general version in (32) for our purposes. We drop the tensor product symbol if the context is clear.

In the interaction picture, the time evolution operator generated by the interaction Hamiltonian (32) reads [55]

$$\hat { U } = \mathcal { T } \exp \left [ - i \int _ { \mathcal { M } } \mathrm d V \sum _ { \nu } \hat { h } _ { I, \nu } \right ] \quad \text{(3)} \quad \text{for} \,,$$

where d V := d n x √ -g is the invariant volume element and T is the time-ordering prescription 4 . In general, it is not possible to solve the time evolution analytically due to the time-dependence of the interaction, thus many analyses involving the UDW detectors consider perturbative methods in the weak coupling regimes. However, our goal in this paper is to establish a perfect quantum channel between Alice and Bob and as was explained in [53], this requires our model to be in the non-perturbative regime.

In order to construct the perfect quantum channel that can transmit quantum information perfectly, we need to be able to calculate the interaction unitary ˆ U

3 This is a slight modification of the usual spin-boson model since typically spin-boson models are assumed to have timeindependent interactions [59] and to account for the covariant formulation of the UDW model [55].

4 The time-ordering is assumed to be given by some global time function using some spacetime foliation, see [77] for discussions about some subtleties of time-ordering for the spatially smeared UDW model.

non-perturbatively (in the sense of not truncating the Dyson series in the weak coupling regime). For this purpose, we will need to use the so-called delta-coupling regime [58], where the interaction timescale is assumed to be much faster than all the relevant timescales of the problem so that the interaction can be taken to occur at a single instant in time. Delta-coupled detectors have been used in many contexts such as entanglement extraction, relativistic communication, work extraction, and many others (see, e.g., [41-43, 53, 78-82]). As each detector has its own rest frame, we need to consider the Fermi normal coordinates (FNC) associated with each observer [83], whose center-of-mass (COM) trajectory is parametrized by their respective proper times [55].

The physically reasonable assumption we make is that in the FNC, x ≡ ( τ, x ), the spacetime smearing functions can be decomposed as

$$f ( x ( \bar { x } ) ) = \chi ( \tau ) F ( \bar { x } ) \,, \, \text{ \quad \ \ } ( 3 4 )$$

where ( τ, x = 0 ) labels the COM trajectory. In other words, Alice (and Bob) can in their rest frames distinguish the spatial profile (which is assumed to satisfy 'Born rigidity', i.e., time-independent with respect to τ ) from the switching function that controls the duration and strength of the interactions. The Hamiltonian density in this case reads

$$\text{ur} \quad \hat { h } _ { I, \nu } ( \overline { x } ) \coloneqq \sum _ { i } \lambda _ { \nu, i } \chi _ { \nu, i } ( \tau ) F _ { \nu, i } ( \overline { \mathfrak { x } } ) \hat { m } _ { \nu, i } ( \tau ) \hat { \phi } ( x ( \overline { x } ) ) \quad ( 3 5 )$$

and the delta coupling interaction corresponds to the case when

$$\chi _ { \nu, i } ( \tau ) = \delta ( \tau - \tau _ { \nu, i } ) \text{ \quad \ \ } ( 3 6 )$$

for some constants τ ν,i ∈ R .

For our purposes, we consider the scenario where Alice and Bob each interact with the field via delta coupling twice, so that we need to consider multiple interaction times τ ν,i with τ ν, 1 &lt; τ ν, 2 . More generally, if Alice (or Bob) performs a sequence of N delta interactions with τ ν,i &lt; τ ν,i +1 , the corresponding unitary time evolution ˆ U ν can be decomposed into a sequence of simple-generated unitaries [58, 78] (see also [82])

$$\underset { \text{$neq} } { \text{nel} } \quad \hat { U } _ { \nu } = \prod _ { j = 1 } ^ { N } \hat { U } _ { \nu, j } \,, \, \hat { U } _ { j } = \exp \left [ - i \lambda _ { \nu } \hat { m } _ { \nu } ( \tau _ { \nu, j } ) \hat { \phi } ( F _ { \nu, j } ) \right ] \ \ ( 3 7 )$$

The time-ordering T is taken care of since the unitaries ˆ U i are written in time-increasing order. When two or more observers are involved, we can use the global time coordinates to order ˆ U a and ˆ U b : for example, if Alice and Bob are spacelike separated then the ordering does not matter, while Alice unitary acts before Bob ˆ U b if supp( f b ) is in the causal future of supp( f a ). The unitary ˆ U ν can be computed non-perturbatively since ˆ U j is a controlled unitary : for any two Hermitian operators ˆ A,B ˆ such that ˆ = A ∑ r i =1 a j | a j 〉 〈 a j | is finite-dimensional

operator in its spectral decomposition, we can write

$$\exp \left [ i \hat { A } \otimes \hat { B } \right ] = \sum _ { j = 1 } ^ { r } | a _ { j } \rangle \langle a _ { j } | \otimes \exp \left [ i a _ { j } \hat { B } \right ] \,. \quad ( 3 8 ) \quad \text{$\text{$\text{$th$}$ really}$}$$

As we will see, for our setup the full unitary time evolution ˆ U will consist of products of Alice's two simplegenerated unitaries and also Bob's two simple-generated unitaries.

Last but not least, let us recall the unitaries found in [53] for Minkowski spacetime. There Alice performs two simple-generated unitaries, but one of them involves coupling to a conjugate momentum

$$\hat { U } _ { \Lambda } = e ^ { i \hat { \sigma } _ { x } \hat { \pi } ( F _ { \Lambda } ) } e ^ { i \hat { \sigma } _ { z } \hat { \phi } ( F _ { \Lambda } ) } \,, \quad \quad ( 3 9 ) \quad \text{op}$$

where the spatially smeared operators are

$$\hat { \phi } ( F _ { \lambda } ) = \lambda _ { 1 } \, \int \mathrm d ^ { 3 } x \, F _ { \lambda } ( x ) \hat { \phi } ( t _ { \lambda }, x ) \,, \quad \ \ ( 4 0 \mathrm a ) \quad \ \text{ we a}$$

$$\hat { \pi } ( F _ { \lambda } ) = \lambda _ { 2 } \, \int \mathrm d ^ { 3 } x \, F _ { \lambda } ( x ) \partial _ { t } \hat { \phi } ( t _ { \lambda }, x ) \,. \quad \ ( 4 0 b )$$

Note that here they take t a := t A, 1 and t A, 2 = t a + /epsilon1 and consider the limit as /epsilon1 → 0 + so that the unitary generated by the conjugate momentum acts later. Bob's unitary is then constructed as a 'suitable inverse' of Alice's unitary to faithfully extract Alice's qubit.

Our construction requires us to make two modifications to the above unitaries. First, we can generalize the conjugate momentum operator in curved spacetime, i.e.,

$$\hat { \pi } ( x ) \coloneqq \partial _ { t } \hat { \phi } ( x ) \equiv \sqrt { h } t ^ { a } \nabla _ { a } \hat { \phi } ( x ) \,, \quad \ \ ( 4 1 ) \quad \ A$$

where t a is the forward pointing unit normal, orthogonal to Σ t with an induced metric h . Second, it is straightforward to see from integration by parts that

$$\hat { \pi } ( f ) = \int \text{d} V \, f ( x ) \partial _ { t } \hat { \phi } ( x ) = \hat { \phi } ( - \partial _ { t } f ) \,. \quad ( 4 2 ) \quad \text{a p}$$

Thus we can always view the smeared conjugate momentum operator as a smeared field operator with a different smearing function (namely, its time derivative). In Minkowski spacetime, using the spacetime smearing notation we can thus write

$$\hat { \pi } ( F _ { \Lambda } ) \equiv \hat { \pi } ( \delta \cdot F _ { \Lambda } ) = \hat { \phi } ( - \delta ^ { \prime } \cdot F _ { \Lambda } ) \,, \quad ( 4 3 ) \quad \text{$\text{$\text{$\text{thrc}$} }$}$$

where δ ′ ( ) is the distributional derivative of the Dirac t delta function δ t ( ). This second step might appear unconventional since ˆ π is not usually regarded as a spacetime-smeared operator, but this provides a simple way to recast all basic observables in terms of a single field operator ˆ . φ The connection between the modern AQFT approach based on spacetime-smeared operators ( cf. Section II) and the 'old' (non-covariant) AQFT formulations based on two spatially smeared operators { ˆ φ | Σ , π ˆ | Σ } at a fixed Cauchy slice Σ with equal-time CCR [ ˆ ( φ t, x ) , π ˆ( t, y )] = iδ ( n ) ( x -y ) can be made more precise and rigorous (see, e.g., [84, 85], and also [62, 66] for details).

/negationslash

The two modifications above suggest that what we really need for the generalization to curved spacetime is that Alice couples to two smeared field operators ˆ ( φ f 1 ) , φ ˆ ( f 2 ) that do not commute, i.e., E f , f ( 1 2 ) = 0, and not necessarily to the field and its conjugate momentum. Indeed, in the modern AQFT framework the coupling to momentum is somewhat awkward because the perfect quantum channel in [53] requires that Alice couples to ˆ ( φ F a ) more strongly than ˆ( π F a ), but the coupling strengths λ , λ 1 2 would need to have different units, so one needs an additional length scale to compare their strengths. By considering different smeared field operators ˆ ( φ f i ), we can have λ i with the same units and we could use other length scales (e.g., from the effective size of the spatial smearing) if we insist on interpreting one of them as conjugate momentum. In what follows we assume that f i (and hence λ i ) have the same units for simplicity.

## B. Quantum capacity and coherent information

Here we briefly review the notion of quantum capacity of a quantum channel and quantum coherent information, following closely the exposition in [41, 53] (see [8689] for more in-depth details). We denote by D ( H ) the space of density operators associated with some Hilbert space H . Let Φ : D ( H a ) → D ( H b ) be a quantum channel - a completely-positive trace-preserving (CPTP) map - between Alice's and Bob's density operators.

An ideal quantum communication channel is one that is able to send any quantum state from Alice to Bob, or equivalently, one that preserves entanglement reliably since Alice's state can always be viewed as an entangled state with some environment that purifies it. Schematically the idea is the following. First, Alice prepares a purification of the state ρ a ∈ D ( H a ), namely, ρ ea := | ψ ea 〉 〈 ψ ea | ∈ H ⊗H e a with H ∼ H a = e . Then Alice applies an encoding channel E : D ( H ⊗H e a ) → D ( H ⊗H e ⊗ N a ) to map her share of pure state into N quantum systems, i.e.,

$$\rho _ { \mathrm E A, N } \coloneqq ( \mathbf 1 _ { \mathrm E } \otimes \mathcal { E } ) ( \rho _ { \mathrm E A } ) \,. \quad \quad ( 4 4 )$$

After encoding, Alice sends each quantum system through N independent uses of the communication channel Φ, i.e.,

$$\rho _ { _ { E B }, N } \coloneqq \Phi ^ { \otimes N } ( \rho _ { _ { E A }, N } ) \in \mathcal { D } ( \mathcal { H } _ { _ { E } } \otimes \mathcal { H } _ { _ { B } } ^ { \otimes N } ) \,. \quad ( 4 5 )$$

Bob's task is to use a decoding channel D : D ( H ⊗ e H ⊗ N b ) → D ( H ⊗H e b ), such that

$$\rho _ { \text{EB} } \coloneqq ( \mathbf 1 _ { \text{E} } \otimes \mathcal { D } ) ( \rho _ { \text{EB}, N } ) \,. \quad \quad ( 4 6 )$$

This is a good encoding-decoding procedure if for a given /epsilon1 we can have

$$| | \widetilde { \rho } _ { \text{EB} } - \rho _ { \text{EB} } | | _ { 1 } \leq \epsilon \text{ \quad \quad } ( 4 7 )$$

↦

where || · || 1 is the trace norm, ˜ ρ eb is the state obtained if Φ is replaced with a noiseless channel 1 1 a → b : ( ρ ∈ D ( H a )) → ( ρ ∈ D ( H b )).

The rate of the communication channel is the number of qubits transmitted per use of the channel, i.e.,

$$R \coloneqq \frac { 1 } { n } \log _ { 2 } ( \dim \mathcal { H } _ { \Lambda } ) \text{ \quad \ \ } ( 4 8 ) \text{ \quad \ } \text{can},$$

and this rate is said to be achievable if for any δ, /epsilon1 &gt; 0, there exists an ( n, R -δ, /epsilon1 ) for sufficiently large n . The quantum capacity Q [Φ] of the channel Φ is defined by the supremum of all achievable rates (see [86, 87] for more details).

For practical calculations that we need, we consider a useful form of the quantum capacity in terms of an information-theoretic measure called coherent information . Given a quantum channel Φ : D ( H a ) → D ( H b ), let | ψ ea 〉 be a purification of Alice's input state ρ a and define

$$\rho _ { \text{EB} } \coloneqq ( \mathbf 1 \otimes \Phi ) ( | \psi _ { \text{EA} } \rangle \langle \psi _ { \text{EA} } | ) \,. \quad \quad ( 4 9 ) \quad \overset { \dots } { \Gamma } ^ { \quad t }.$$

Using ρ b := Tr ( e ρ eb ) = Φ( ρ a ), the coherent information I c ( ρ , a Φ) is defined to be [90, 91]

$$I _ { c } ( \rho _ { \Lambda }, \Phi ) \coloneqq S ( \Phi ( \rho _ { \Lambda } ) ) - S ( \rho _ { \text{EB} } ) \,. \quad \ \ ( 5 0 ) \quad \text{ for } \mathbb { m } { \mathbb { m } { f o r } }$$

We may sometimes write I c ( ψ ea , Φ) to make clear the joint state on EA . The coherent information depends on the input state ρ a and the channel Φ but is independent of its purification. From this, we can then define the coherent information of a quantum channel I c (Φ) by maximizing over all input states:

$$I _ { c } ( \Phi ) \coloneqq \max _ { \rho _ { \Lambda } } I _ { c } ( \rho _ { \Lambda }, \Phi ) \,. \quad \quad ( 5 1 ) \quad \begin{matrix} I _ { c } \stackrel { \prime } { ( } \\ \ln } \end{matrix}$$

Note that for qubit channels we have that -1 ≤ I c ( ρ, Φ) ≤ 1 and 0 ≤ I c (Φ) ≤ 1 [87]. We can then rewrite the quantum capacity using the Lloyd-ShorDevetak (LSD) formula [91-93], given by [86-89]

$$\mathcal { Q } [ \Phi ] \coloneqq \lim _ { n \to \infty } \frac { 1 } { n } I _ { c } ( \Phi ^ { \otimes n } ) \,, \quad \ \ ( 5 2 ) \quad t u r _ { \text{(su)} }$$

where the RHS is also known as regularized coherent information.

Despite the somewhat physically intuitive expression for Q [Φ], the asymptotic nature of Eq. (52) makes quantum capacity very difficult to compute in most cases unless the quantum channel has very specific properties (e.g., when the quantum channel is additive , namely I c [Φ ⊗ n ] = nI c [Φ], which is false in general [86, 87]). However, for us we do not actually need this: we follow the strategy adopted in [53], which is to construct a quantum channel for qubits that has I c (Φ) → 1 by appropriately adjusting the setup parameters. This works because the quantum coherent information is superadditive [94, 95]:

$$I _ { c } ( \Phi ^ { \otimes N } ) \geq N I _ { c } ( \Phi ) \text{ \quad \ \ } ( 5 3 ) \text{ \quad \ } _ { c } ^ { 1 }$$

which immediately implies that

$$\mathcal { Q } [ \Phi ] \geq I _ { c } ( \Phi ) \geq I _ { c } ( \rho _ { \Lambda }, \Phi ) \,. \quad \quad ( 5 4 )$$

For qubit channels, we have 0 ≤ Q [Φ] ≤ 1, thus we can construct an essentially perfect quantum channel if we can show that I c [Φ] can be made arbitrarily close to 1.

As a final remark, we would like to point out that the existence of many superadditive quantum channels are in some sense a 'curse' to finding a good characterization of quantum capacity - we refer the readers to [96] for one of the more recent summary and progress on this front. From Eq. (54) we see that we need 0 &lt; I c (Φ) ≤ 1 for the channel to be able to transmit some quantum information and our construction capitalizes on this. However, the case for I c (Φ) = 0 is much more subtle, as many quantum channels exhibit strict superadditivity in both weak and strong sense [96]. Weak superadditivity means that there exists N such that I c (Φ ⊗ N ) &gt; NI c (Φ), and strong superadditivity means I c (Φ 1 ⊗ Φ ) 2 &gt; I c (Φ ) + 1 I c (Φ ) for two channels Φ 2 1 , Φ . 2 In both cases, the 'single-letter' formula for coherent information I c (Φ) is not sufficient to calculate the quantum capacity, although some exceptions are known. For instance, all degradable channels (such as some regimes of the dephasing channels) are known to be weakly additive, in which case Q [Φ] = I c (Φ) and hence non-positive channel coherent information would imply zero quantum capacity [86, 87, 97]. All anti-degradable channels (such as all entanglement-breaking channels and certain regimes of quantum erasure channels) are known to have zero quantum capacity, though the argument does not rely on additivity and more towards no-cloning type arguments [97, 98]. Our construction (similarly for [53]) exploits the structure of the detector-field coupling to make the channel coherent information tractable.

The interpretation of the above analysis is that in the limit where Q [Φ] → 1, the quantum channel can effectively transmit quantum information (qubits) reliably in one use of the channel . It is perhaps not surprising that this is possible if and only if the whole procedure is essentially a 'swapping operation', i.e., Alice swaps her qubit 'into the field', and then Bob swaps out the (approximate) embedded qubit in the field state to recover Alice's state. In effect, this suggests that the core of the protocol comprises of two components: (1) Bob being positioned correctly in spacetime in order to recover Alice's qubit, and (2) Alice and Bob being able to implement as high quality 'SWAP gates' as possible . 5

5 In [53], Alice's unitary is called ENCODING gate, and Bob's unitary is called DECODING gate. While the terminology has physical justification, it does not coincide with the encoding and decoding channels E D , even in the case of a single use of the communication channel Φ ( N = 1), as they are not part of Φ whose capacity is being computed. The ENCODING and DECODING gates were part of the definition of the communication channel.

## IV. PROTOCOL FOR THE PERFECT RELATIVISTIC QUANTUM CHANNEL

In this section, we present the protocol for Alice and Bob to implement a perfect quantum channel through the field. While our protocol closely follows [53], our approach naturally extends to arbitrary background spacetimes and choices of vacuum states of the field by recasting the protocol in the AQFT framework.

A brief visual of the quantum channel is depicted in Fig. 1. Mathematically, the quantum channel from Alice to Bob generically takes the form

$$\Phi ( \rho _ { \Lambda } ) = \text{Tr} _ { \Lambda, \phi } \left [ \hat { U } _ { \text{B} } \hat { U } _ { \Lambda } ( \rho _ { \Lambda } \otimes | 0 \rangle \langle 0 | \otimes \rho _ { \text{B} } ) \hat { U } _ { \text{A} } ^ { \dagger } \hat { U } _ { \text{B} } ^ { \dagger } \right ], \ \ ( 5 5 )$$

where | 0 〉 is the field vacuum , 6 ρ b is Bob's initial state, and ˆ U a and ˆ U b are the unitaries implemented by Alice and Bob through the UDW interaction with the field. The task is to construct a perfect quantum channel that can faithfully transmit quantum information to arbitrary accuracy.

More precisely, Alice starts with an arbitrary state | ψ a 〉 := c 1 | + z 〉 + c 2 |- 〉 z , equivalently ρ a = | ψ a 〉 〈 ψ a | , and encodes that state into the field using a 'SWAP' gate. Bob should then be able to decode the state from the quantum field into his own Hilbert space. If the 'SWAP' gates ˆ U a and ˆ U b are implemented correctly, then in the idealized limit we should have a perfect quantum channel Φ( ρ a ) = 1 1 a → b ( ρ a ) ≡ ρ a ∈ H b , noting that this channel is not trivial because it corresponds to sending the same state ρ a from Alice's Hilbert space to Bob's

Essentially, the basic idea is for Alice to (approximately) perform a 'SWAP' gate 7 on her qubit state to the field, effectively encoding it into some superposition of approximately orthogonal coherent states, so that Bob can extract the qubit state by performing another swap operation. The entire problem thus reduces to analyzing under what circumstances the UDW interaction allows for nearly perfect 'SWAP' operations on both parties which will depend on both the coupling strengths and the supports of the interactions. In particular, we will see that Bob can only perform a perfect 'SWAP' gate if he is causally connected to Alice and the details of the interactions determine the quality of the swap operation.

6 The state | 0 〉 can be replaced with any GNS vector state | Ω ω 〉 associated with any of the quasifree representations, which accounts for the existence of many unitarily inequivalent ground states in QFT [62-65, 70].

7 Technically speaking, the qubit and the quantum field have different Hilbert space dimensions but we will call them 'SWAP' gates anyway since it would be the usual SWAP gate if we were to replace the quantum field with a third qubit. The fact that we are 'embedding' the qubit state into the (much higherdimensional) field can be viewed as a form of encoding operation, which justifies calling Alice's unitary ENCODING gate in [53], and similarly Bob's unitary can be viewed as decoding operation from the field to his qubit system.

FIG. 1. A perfect quantum channel from Alice to Bob through a quantum field ˆ consists φ of two unitary operators: Alice encoding her state into the field and then Bob decoding the state. The time in ∆ t between the two gates is a general t that foliates the spacetime.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00518v3-figure-2c87d316-1d83-4ebb-976f-4a985ff8a89f.jpg)

Hilbert space. We start by discussing how Alice encodes her state into the field:

Proposition 1. Let | ψ a 〉 = c 1 | + z 〉 + c 2 |- 〉 z be Alice's initial state and let | 0 〉 be the initial vacuum state of the field. Consider Alice's unitary operation at time t = t a given by 8

$$\hat { U } _ { A } = e ^ { i \hat { \sigma } _ { x } \hat { \phi } ( f _ { 2 } ) } e ^ { i \hat { \sigma } _ { z } \hat { \phi } ( f _ { 1 } ) } \,. \quad \quad \quad ( 5 6 )$$

Then Alice can approximately swap the state | ψ a 〉 into the field provided the following conditions are satisfied

$$E ( f _ { 1 }, f _ { 2 } ) = \frac { \pi } { 4 } \bmod 2 \pi \,, \quad \quad ( 5 7 \tt a ) \\ \tau \, \tau \, \tau \, \tau \, \circ \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \tau \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \end{array}$$

$$E ( f _ { 1 }, f _ { 2 } ) ^ { 2 } \gg \overset { \cdot } { W } ( f _ { 2 }, f _ { 2 } ) \,. \quad \quad ( 5 7 b )$$

The resulting state after Alice's operation reads

$$\lambda & \quad \cdot \quad \cdot \\ \hat { U } _ { A } ( c _ { 1 } | + _ { z } \rangle & + c _ { 2 } | - _ { z } \rangle ) \, | 0 \rangle \\ & = | + _ { y } \rangle \, ( c _ { 1 } | + \alpha \rangle - i c _ { 2 } | - \alpha \rangle ) \,, \quad ( 5 8 )$$

where |± α 〉 are coherent states of the quantum field defined by

$$\left | \pm \alpha \right \rangle \coloneqq e ^ { \pm i \hat { \phi } \left ( f _ { 1 } \right ) } \left | 0 \right \rangle. \text{ \quad \ \ } \text{ (5 9)}$$

Proof. In order to prove this statement, we will explicitly apply Eq. (56) to Alice's qubit and the field and show that the qubit state is encoded under the conditions in Eqs. (57a) and (57b).

The first unitary in Eq. (56) gives

$$e ^ { i \hat { \sigma } _ { z } \hat { \phi } ( f _ { 1 } ) } \left | \psi _ { _ { A } } \right \rangle \left | 0 \right \rangle = c _ { 1 } \left | + _ { z } \right \rangle \left | + \alpha \right \rangle + c _ { 2 } \left | - _ { z } \right \rangle \left | - \alpha \right \rangle. \ \ ( 6 0 )$$

Next, we need to apply the second unitary e iσ ˆ x ˆ ( φ f 2 ) but we first note that

$$\text{ound} _ { \infty \, a i f } \hat { \phi } ( f _ { 2 } ) \left | \pm \alpha \right \rangle = \pm E ( f _ { 1 }, f _ { 2 } ) \left | \pm \alpha \right \rangle + e ^ { \pm i \hat { \phi } ( f _ { 1 } ) } \hat { \phi } ( f _ { 2 } ) \left | 0 \right \rangle \,, \, ( 6 1 )$$

8 The usual convention would be ˆ U a = e -iσ ˆ x ˆ ( φ f 2 ) e -iσ ˆ z ˆ ( φ f 1 ) , but we opt for the positive sign for the exponent for convenience as many expressions appear more symmetric (this was also implicitly done in [53]). This amounts to the replacement ˆ h I →-ˆ h I , or equivalently f i →-f i .

which follows directly from the Baker-CampbellHaussdorff (BCH) formula and the CCR. Now we see that if

$$E ( f _ { 1 }, f _ { 2 } ) ^ { 2 } \gg \langle 0 | \hat { \phi } ( f _ { 2 } ) \hat { \phi } ( f _ { 2 } ) | 0 \rangle = W ( f _ { 2 }, f _ { 2 } ) \quad ( 6 2 ) \quad \text{Pro} _ { \xi }$$

then |± α 〉 become approximate eigenvectors of ˆ ( φ f 2 )

$$\hat { \phi } ( f _ { 2 } ) \left | \pm \alpha \right \rangle \approx \pm E ( f _ { 1 }, f _ { 2 } ) \left | \pm \alpha \right \rangle. \quad \quad ( 6 3 )$$

If we pick the eigenvalues to be E f , f ( 1 2 ) = π 4 mod 2 π , then the action of Alice's full unitary ˆ U a gives

$$\ t n e n \ t h e \, \text{action or Alice} \, & s \, \text{will unitary } \, \upsilon _ { A } \, \text{ gives} \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \ \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \
 \hat { U } _ { A } & ( c _ { 1 } \, | + _ { z } \, \rangle + c _ { 2 } \, | - _ { z } \, \rangle ) \, | 0 \rangle \\ & = e ^ { i \hat { \sigma } _ { z } } \hat { \phi } ( f _ { 2 } ) ( c _ { 1 } \, | + _ { z } \, \rangle \, | + \alpha \rangle + c _ { 2 } \, | - _ { z } \, \rangle \, | - \alpha \rangle ) \\ & \approx c _ { 1 } e ^ { i \frac { \pi } { 4 } \hat { \sigma } _ { z } } \, | + _ { z } \, \rangle \, | + \alpha \rangle + c _ { 2 } e ^ { - i \frac { \pi } { 4 } \hat { \sigma } _ { z } } \, | - _ { z } \, \rangle \, | - \alpha \rangle ) \\ & = | + _ { y } \rangle \, ( c _ { 1 } \, | + \alpha \rangle - i c _ { 2 } \, | - \alpha \rangle ) \,, \quad \quad \quad$$

where we used e i π 4 ˆ σ x | + z 〉 = | + y 〉 and e -iσ ˆ x π 4 |- 〉 z = - | i + y 〉 in the last equality.

The second of Alice's conditions in Eq. (57b) can be interpreted as a 'strong coupling condition' that requires Alice's coupling through the field via the interaction profile f 1 be much stronger than through f 2 . In [53], this corresponds to choosing f 2 such that ˆ φ 2 corresponds to conjugate momentum ˆ( π F a ) and the condition translates to Alice coupling to the field ˆ ( φ F a ) much more strongly than to ˆ( π F a ) ( cf. Eq. (39)). We can interpret Eq. (57a) as a 'fine-tuning condition' that performs near exact π 4 rotations in the Bloch sphere so that Alice's second unitary decouples her detector from the quantum field.

Now we move on to discussing how at a later time, Bob can decode Alice's state from the quantum field into his own qubit detector. Even in the idealized limit where Alice effectively implements a 'SWAP' gate through unitary ˆ U a , Bob cannot simply perform the inverse unitary ˆ U -1 a for two reasons. First, Bob is in the causal future of Alice, so the smearing functions of Bob's interactions cannot be at the same time slice as f i of Alice. More importantly, ˆ U -1 a will simply undo what Alice does resulting in the trivial identity channel 1 1 ( a ρ a ) = ρ a ∈ H a which corresponds no interactions with the field at all, instead of a noiseless communication channel 1 1 a → b ( ρ a ) = ρ a ∈ H b through the field. In other words, we want that Alice's operations defined through ˆ ( φ f i ) be defined at time t = t a , while we want Bob to decode from the field using only field operators ˆ ( φ g i ) localized at time t = t b . Bob's unitary will thus take the form

$$\hat { U } _ { _ { B } } \coloneqq e ^ { - i \hat { \sigma } _ { z } ^ { B } \hat { \phi } ( g _ { 1 } ) } e ^ { - i \hat { \sigma } _ { x } ^ { B } \hat { \phi } ( g _ { 2 } ) } \,, \quad \quad ( 6 5 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where the superscript 'B' emphasizes that Bob acts on H b and g i are chosen such that ˆ U b implements a 'SWAP' gate to decode the qubit from the field at time t b . For this, we need a way to write the operators ˆ ( φ f i ) using different spacetime smearing functions g i defined at a different Cauchy slice. To achieve this we use the following theorem:

Proposition 2. Suppose we have two smeared field operators ˆ ( φ f ) and ˆ ( φ g ) , then the two operators are equal ˆ ( φ f ) = ˆ ( φ g ) if the smearing functions satisfy

$$E f = E g.$$

Proof. The proof is almost immediate from the symplectic definition of smeared operators in Eq. (7). If we assume Ef = Eg 9 , then

$$\hat { \phi } ( f ) = \sigma ( E f, \hat { \phi } ) = \sigma ( E g, \hat { \phi } ) = \hat { \phi } ( g ) \,. \quad ( 6 7 )$$

Given that Alice and Bob's smearing functions f i and g i are related by Eq. (66) and given that Bob's qubit starts in the | + y 〉 state (the state that Alice's qubits is left as after ˆ U A ), we have achieved our goal of creating a perfect quantum channel from Alice to Bob.

Second, despite its simplicity and simple interpretation, Eq. (66) is largely impractical to use for actually calculating Bob's spatial smearing functions G i ( x ) for a given time t = t b . Using the eigenmode basis provided in Eq. (26) and the momentum space smearing functions in Eq. (29), the condition in Eq. (66) can be equivalently written in 'momentum space' as

Some remarks are in order in view of the above Proposition 2. First, intuitively Alice's state is encoded into the field amplitude that propagates into the causal future of Alice's interaction region J + (supp f i ), where J + (supp( f i )) is the causal future of supp( f i ); for Bob to decode that state, his qubit has to be coupled everywhere where the initial data has propagated at a given time t = t b , which in turn forces Bob's detector to interact with the field with the correct support and coupling strengths. In particular, supp( g i ) should essentially be Σ t b ∩ supp( Ef i ), where supp( Ef i ) ⊆ J + (supp f i ). For simplicity, we assume that supp( f 1 ) = supp( f 2 ) and hence supp( g 1 ) = supp( g 2 ), which in the delta-coupling scenario is easily enforced as we are fixing the spatial profile of the interaction on Alice's side - this, in turn, fixes the needed spatial profile for Bob's coupling.

$$f _ { k } = g _ { k }.$$

These equations in momentum space are often easier to work with since in many simple cases (such as all FRW spacetimes) the spatial part of the mode function u k ( x ) ∝ e i k x · and we can use Fourier transforms, which is similar to what was originally done in [53] where

9 The converse statement is false. From the axioms of AQFT, we also have that ˆ ( φ Pf ˆ ) = 0 which is the KG equation. Thus ˆ ( φ g ) is only unique up to homogeneous solution of the KG equation, i.e., g ∼ g + ˆ Pf .

Bob's spatial smearing functions were related to Alice's in momentum space. More generally, in generic curved spacetimes, we do not have a simple plane-wave basis, so it is essential that all the conditions we demand for the perfect quantum channel are given only in terms of the correlation functions of the field without reference to Fourier transforms.

/negationslash

In summary, we have constructed a perfect relativistic quantum channel that allows Alice to send her qubit state to Bob through the quantum field. The quantum channel consists of two steps:

Last but not least, note that there is an asymmetry in how we quantify the quality of Alice's and Bob's respective 'SWAP' gates, i.e., how close their unitaries are in achieving the ideal 'SWAP' operation. On the one hand, Alice needs two smearing functions f 1 and f 2 that satisfy the conditions in Proposition 1 to approximate the ideal 'SWAP'. On the other hand, Bob needs two smearing functions g 1 and g 2 that are related to Alice's via Proposition 2 and hence automatically satisfy the conditions in Proposition 1. However, Proposition 1 alone is not enough for the perfect quantum channel: if Ef i = Eg i , Bob will not be able to decode Alice's state even though by construction Bob's smearing functions g 1 and g 2 satisfy Eqs. (57a) and (57b). In a sense, to decode the state using the 'SWAP' operation with the field Bob needs an extra condition, namely that Bob needs to interact 'at the correct spacetime regions'. The reason for this asymmetry is that the Hilbert space dimension of the field ˆ is φ (much) larger than that of the detectors: it is much simpler for Alice to encode her state in a larger system but Bob will need additional information to decode it.

- (1) At time τ = τ a , Alice encodes her qubit into the field using the unitary ˆ U a such that the conditions in Proposition 1 are satisfied.
- (2) At a later time τ = τ b , Bob decodes the state from the field using the unitary ˆ U b defined in Eq. (65). However, to implement this protocol, Bob's detector must interact across a spacetime region that satisfies the conditions specified in Proposition 2.

As the proper times are not the same in both trajectories, in practice one can use the ambient global coordinate systems and foliation to help with the timeordering of these operations. Once times τ a and τ b are specified, we demand that Alice's spatial smearing functions F i ( x ) fulfill conditions (57a) and (57b), and Bob's functions G i ( x ) fulfill the condition (66) 10 .

10 In principle, Alice or Bob could each need up to four spatial smearing functions since f i or g i could include both δ t ( ) and δ ′ ( ) terms. t For example, in flat spacetime, we can have f 1 = δ t F ( ) 11 ( x ) + δ ′ ( ) t F 12 ( x ) (similarly for f 2 ), but for simplicity these can be chosen to be equal. Both g 1 and g 2 are fixed by f 1 and f 2 in the ideal case.

In order to simplify our subsequent discussions, we make a physically reasonable assumption that Alice's interaction profiles f i are completely fixed except their coupling constants, i.e., in FNC we have

$$f _ { i } = \lambda _ { i } \chi _ { i } ( \tau ) F ( \overline { x } )$$

where χ i are either the Dirac delta function δ or its distributional derivative δ ′ and the spatial profile is the same for both f , f 1 2 (hence they have the same supports supp( f i ) = supp( F ) ⊆ M ). In other words, Alice only varies her interaction profile f i by tuning the coupling strengths λ i without changing the 'shape' of the interaction profile, and this is sufficient to satisfy the perfect 'SWAP' conditions (57b)-(57a) in Proposition 1.

## V. ANALYTIC TESTS OF THE QUANTUM CHANNEL

In this section, we will prove that the protocol given in Section IV produces a perfect quantum channel in the correct limit without any approximations like Eq. (63). This proves that entanglement can always be perfectly transmitted between two localized qubits through a quantum field regardless of the background spacetime geometry. We also discuss the case when Alice and Bob can in principle perfectly implement the SWAP gates with the field but are spacelike separated from each other. For this channel, the quantum and classical channel capacity are both zero, as we expect with relativistic causality.

## A. Ideal case: perfect quantum channel

Recall that our task is to show that the channel coherent information I c (Φ) can be made arbitrarily close to 1, thus proving that the (relativistic) quantum channel can transmit quantum information perfectly. This means computing the coherent information I c ( ρ a , 0 , Φ) and maximizing over all possible input states on Alice's side. Using the definition of coherent information (51), it turns out that it is sufficient to consider the case when Alice's initial state is the maximally mixed state ρ A, 0 = 1 2 1 1 and Bob initial state is | + y 〉 state in order to show that we can attain maximum channel coherent information I c (Φ) → 1 and hence maximal quantum capacity. We stress that this does not mean Alice should send a maximally mixed state to Bob, which is 'garbage' as a message and is even classical, since in practice Alice wants to be able to send pure states properly. Also, in general, it is not true that I c (Φ) is maximized by the maximally mixed state: for example, the completely depolarizing channel Φ( ) ρ = 1 2 has 1 / I c ( 1 1 / 2 Φ) = , -1 but the channel coherent information I c (Φ) = 0 is attained by any pure state.

In order to show that I c (Φ) → 1 is attainable, we need to introduce a qubit of the environment E that purifies

Alice's maximally mixed state 1 1 / 2, such that the joint state of Alice and the environment is the maximally entangled state | ψ ea 〉 = 1 √ 2 ( |- 〉 | z + z 〉 + + | z 〉 |-z 〉 ) ea . As discussed in Section III, a perfect quantum channel will preserve all the entanglement between Alice and the environment EA to Bob and the same environment EB . To compute I c ( 1 1 / 2 Φ), we must evaluate the state , ρ eb in Eq. (49).

Writing out the expression for ρ eb in terms of the unitaries ˆ U a and ˆ U b , we obtain:

$$\rho _ { \text{EB} } = & \text{Tr} _ { \text{A} \phi } [ \hat { U } _ { \text{B} } \hat { U } _ { \text{A} } ( | \psi _ { \text{EA} } \rangle \langle \psi _ { \text{EA} } | \otimes | 0 \rangle \langle 0 | \otimes | + _ { y } \rangle \langle + _ { y } | ) \hat { U } _ { \text{A} } ^ { \dagger } \hat { U } _ { \text{B} } ^ { \dagger } ]. \quad \text{ obtain} \\ \text{Tr} \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \ iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota\iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iia \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \iota \jmath \,,$$

The simplest way of continuing is to decompose the unitaries using projections onto the eigenvectors of the qubit operators. That is, for ˆ U a we write

$$\omega \cdots \, \tilde { x } \cdots \, \tilde { \cdot } \cdots \, \tilde { \cdot } \cdots \, \tilde { \cdot } \cdots \\ \hat { U } _ { \Lambda } & = \sum _ { x, z \in \{ \pm \} } \hat { P } _ { x } \hat { P } _ { z } \otimes e ^ { i x \hat { \phi } ( f _ { 2 } ) } e ^ { i z \hat { \phi } ( f _ { 1 } ) } \\ & = \sum _ { x, z \in \{ \pm \} } \hat { P } _ { x } \hat { P } _ { z } \otimes W ( x E f _ { 2 } ) W ( z E f _ { 1 } ) \quad ( 7 1 ) \quad \text{$(71)$} \quad \text{$begin{matrix}$\theta$} \\ \text{Be main}$} \\ \text{$\text{$\text{$\text{$dow}$} } } }$$

where ˆ P i are the projectors onto the eigenstates of ˆ σ i and in the last line we used the Weyl algebra notation introduced in Eq. (9). For Bob's unitary, we write

$$\hat { U } _ { \text{B} } = \sum _ { x, z \in \{ \pm \} } \hat { P } _ { z } \hat { P } _ { x } \otimes W ( - z E f _ { 1 } ) W ( - x E f _ { 2 } ) \quad ( 7 2 )$$

where the projection operators are understood to act on Bob's system. We have written Bob's unitary as in Eq. (65), the imperfect case is discussed in the following sections. We also expand the Bell state as | ψ ea 〉 = 1 √ 2 ∑ j |-j z 〉 | j z 〉 with j ∈ {±} to get the expression

$$r ^ { \dots \dots \dots } _ { \rho _ { \text{EB} } } & = \frac { 1 } { 2 } \sum _ { j, k, x _ { i }, z _ { i } } \omega ( \hat { O } ) \left \langle k _ { z } | \hat { P } _ { z _ { 1 } } \hat { P } _ { x _ { 1 } } \hat { P } _ { x _ { 4 } } \hat { P } _ { z _ { 4 } } | j _ { z } \right \rangle } \quad \text{while} \quad \text{$full as} \\ & \times | - j _ { z } \right \rangle _ { \text{E} } \left \langle - k _ { z } | \otimes \hat { P } _ { z _ { 3 } } \hat { P } _ { x _ { 3 } } | + _ { y } \right \rangle _ { \text{B} } \left \langle + _ { y } | \hat { P } _ { x _ { 2 } } \hat { P } _ { z _ { 2 } } \quad ( 7 3 ) \quad \text{ity, $b$} \\ & \quad \text{tests}$$

where we write ω O ( ˆ ) ≡ 〈 0 | ˆ 0 O | 〉 and

$$\hat { O } & = W ( - z _ { 1 } E f _ { 1 } ) W ( - x _ { 1 } E f _ { 2 } ) W ( x _ { 2 } E f _ { 2 } ) W ( z _ { 2 } E f _ { 1 } ) \quad \text{ \quad states} \\ & \times W ( - z _ { 3 } E f _ { 1 } ) W ( - x _ { 3 } E f _ { 2 } ) W ( x _ { 4 } E f _ { 2 } ) W ( z _ { 4 } E f _ { 1 } ), \quad \text{ \quad \ } \text{We} \\ & \quad \text{Eq. (}$$

↦

↦

x i stands for x , x 1 2 , x 3 , x 4 , similarly for z i , and all sums are taken over the set {±} . Provided that we can obtain an expression for ω O ( ˆ ), a computer algebra system can straightforwardly sum over the 2 10 terms in the sum. We start by redefining summation indices x 1 → -x 1 , z 1 →-z 1 , x 3 →-x 3 , and z 3 →-z 3 such that

↦

↦

$$\mu _ { 1 } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ \rho _ { \text{EB} } & = \frac { 1 } { 2 } \sum _ { j, k, x _ { i }, z _ { i } } \omega ( \hat { O } ) \left \langle k _ { z } | \hat { P } _ { - z _ { 1 } } \hat { P } _ { - x _ { 1 } } \hat { P } _ { x _ { 4 } } \hat { P } _ { z _ { 4 } } | j _ { z } \right \rangle \quad \text{Note for a} \\ & \times | - j _ { z } \right \rangle _ { \text{E} } \left \langle - k _ { z } | \otimes \hat { P } _ { - z _ { 3 } } \hat { P } _ { - x _ { 3 } } | + _ { y } \right \rangle _ { \text{B} } \left \langle + _ { y } | \hat { P } _ { x _ { 2 } } \hat { P } _ { z _ { 2 } } \ ( 7 5 ) \quad \text{for a} \quad \text{state}$$

and ˆ only contains terms without minus signs O

$$\text{en-} \\ \text{As} \quad & \hat { O } = W$$

In [53], 〈 0 | ˆ 0 O | 〉 is evaluated using the Baker-CampbellHausdorff (BCH) formula and Wick's theorem. In this paper, we adopt a much simpler approach using the AQFT framework. It is straightforward (albeit somewhat tedious) to apply the third Weyl condition W Ef W Eg ( ) ( ) = e -i 2 E f,g ( ) W E f ( ( + )) seven times to g obtain

$$\L \quad \text{s} \quad \text{s} \quad \text{s$$

where only one Weyl operator remains multiplied by a phase factor.

Before we evaluate the expectation value of the remaining Weyl operator, it will be convenient to consider the decomposition of a Wightman two-point function into its real (symmetric) and imaginary (antisymmetric) parts, i.e.,

$$\lambda _ { \nu } \, \text{d} \, \text{d} \$$

where we used the CCR to obtain the second term and H ( f, g ) := ω ( { ˆ ( φ f ) , φ ˆ ( g ) } ) ≡ 1 2 µ Ef,Eg ( ) is the smeared Hadamard two-point distribution of the field ( cf. Eq. (13)). This decomposition is important because only the symmetric part H ( f, g ) is state-dependent , while E f, g ( ) is independent of the state. This is useful as it tells us which parts of our analysis depend on the field state. Furthermore, while we have assumed that the field starts in some vacuum state for simplicity, but the protocol in Section IV automatically extends to all quasi-free states such as thermal states [65], squeezed vacuum state [58], and with slight modifications to more general Gaussian states such as coherent states (see, e.g., [70]).

We can now use the definition of quasifree state ω in Eq. (17) for the vacuum state to evaluate the expectation value ω O ( ˆ ) in terms of Wightman and Hadamard functions giving us

$$\i s u n s \quad \cdots \quad \circ \quad \cdots \\ \$$

Note that this expression is valid automatically for any quasifree states such as squeezed vacuum state and (squeezed) thermal state [58]. The

density matrix in the ordered basis ρ eb in the {| + + z z 〉 be , | + z - 〉 z be , |-z + z 〉 be , |- - 〉 z z be } now reads

$$\begin{array} {$$

where

$$\begin{array} { c c c c c } \mathcal { U } ^ { * } & B ^ { * } & 0 & P _ { - } / & \text{where} \\ \text{Pro} \\ \text{where} \\ P _ { * } = \frac { 1 } { 4 } \left ( 1 \pm e ^ { - 2 W ( f _ { 1 }, f _ { 1 } ) } \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) \right ) \\ X = \frac { 1 } { 4 } \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) \left ( e ^ { - 2 W ( f _ { 2 }, f _ { 2 } ) } + \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) \right ) & \text{for} \, e \\ A = - \frac { i } { 4 } e ^ { - 2 W ( f _ { 1 }, f _ { 1 } ) } \cos ( 2 E ( f _ { 1 }, f _ { 2 } ) ) \\ \times \left ( \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) - e ^ { - 2 W ( f _ { 2 }, f _ { 2 } ) } \sinh ( 4 W ( f _ { 2 }, f _ { 1 } ) ) \right ) & \text{and} \, v \\ B = - \frac { i } { 4 } e ^ { - 2 W ( f _ { 1 }, f _ { 1 } ) } \cos ( 2 E ( f _ { 1 }, f _ { 2 } ) ) \\ \times \left ( \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) + e ^ { - 2 W ( f _ { 2 }, f _ { 2 } ) } \cosh ( 4 W ( f _ { 1 }, f _ { 2 } ) ) \right ) & \text{in or} \\ C = - \frac { 1 } { 4 } e ^ { - 8 W ( f _ { 1 }, f _ { 1 } ) } \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) & \text{in} \\ \times \left ( \sin ( 2 E ( f _ { 1 }, f _ { 2 } ) ) - e ^ { - 2 W ( f _ { 2 }, f _ { 2 } ) } \cosh ( 4 H ( f _ { 1 }, f _ { 2 } ) ) \right ). & \text{c>} \, \end{array} \right ) \\ \text{Now that we have an analytic expression for $\rho_{EB}$, we} \, \text{can use our two conditions in Proposition 1 and ver-}$$

Now that we have an analytic expression for ρ eb , we can use our two conditions in Proposition 1 and verify that the coherent information can be made arbitrarily close to 1. First we use the fine-tuning condition E f , f ( 1 2 ) = π 4 which simplifies ρ eb considerably to

$$\rho _ { \text{EB} } = \begin{pmatrix} \widetilde { P } _ { - } & 0 & 0 & \widetilde { C } \\ 0 & \widetilde { P } _ { + } & \widetilde { P } _ { + } & 0 \\ 0 & \widetilde { P } _ { + } & \widetilde { P } _ { + } & 0 \\ \widetilde { C } & 0 & 0 & \widetilde { P } _ { - } \end{pmatrix} \quad \text{(81)}$$

where

$$\widetilde { P } _ { \pm } & = \frac { 1 } { 4 } \left ( 1 \pm e ^ { - 2 W ( f _ { 2 }, f _ { 2 } ) } \right ) & \quad \text{stat} \quad \text{able} \\ \widetilde { C } & = \frac { 1 } { 4 } e ^ { - 8 W ( f _ { 1 }, f _ { 1 } ) } \left ( e ^ { - 2 W ( f _ { 2 }, f _ { 2 } ) } \cosh ( 4 H ( f _ { 1 }, f _ { 2 } ) ) + 1 \right ). \quad \text{is fi} \\ \cdots \cdots \cdots$$

Next, we use the strong coupling condition E f , f ( 1 2 ) 2 /greatermuch W ( f , f 2 2 ). Recall that we are assuming for simplicity that on Alice's side, her interaction profile will be kept fixed except for the coupling strengths, thus only λ , λ 1 2 are adjustable parameters. In this case, the strong-coupling condition essentially means that λ 1 /greatermuch λ 2 and hence we have the hierarchy of scales

$$W ( f _ { 1 }, f _ { 1 } ) \sim \mathcal { O } ( \lambda _ { 1 } ^ { 2 } ) \gg E ( f _ { 1 }, f _ { 2 } ) ^ { 2 } \sim \mathcal { O } ( \lambda _ { 1 } ^ { 2 } \lambda _ { 2 } ^ { 2 } ) \ ( 8 2 \mathfrak { a } ) \quad \text{perfe}$$

$$E ( f _ { 1 }, f _ { 2 } ) ^ { 2 } \sim \mathcal { O } ( \lambda _ { 1 } ^ { 2 } \lambda _ { 2 } ^ { 2 } ) \gg \mathsf W ( f _ { 2 }, f _ { 2 } ) \sim \mathcal { O } ( \lambda _ { 2 } ^ { 2 } ) \quad ( 8 2 b ) \quad \Big.$$

$$W ( f _ { 1 }, f _ { 1 } ) \sim \mathcal { O } ( \lambda _ { 1 } ^ { 2 } ) \gg \mathsf H ( f _ { 1 }, f _ { 2 } ) \sim \mathcal { O } ( \lambda _ { 1 } \lambda _ { 2 } ) \quad ( 8 2 c )$$

With this hierarchy, we can express λ 1 in units of λ 2 so that λ 1 := cλ 2 for some c ≥ 0, which is also convenient since in general λ i has physical units 11 . We can thus write

$$W ( f _ { 1 }, f _ { 1 } ) = c ^ { 2 } \lambda _ { 2 } ^ { 2 } w _ { 1 } \,, \ \ W ( f _ { 2 }, f _ { 2 } ) = \lambda _ { 2 } ^ { 2 } w _ { 2 } \,, \quad ( 8 3 \text{a} )$$

$$E ( f _ { 1 }, f _ { 2 } ) = c \lambda _ { 2 } ^ { 2 } \text{e} \,, \quad \text{H} ( f _ { 1 }, f _ { 2 } ) = c \lambda _ { 2 } ^ { 2 } \text{h} \,, \quad ( 8 3 \text{b} )$$

where w 1 , w 2 , e h , are fixed by the spacetime smearings. Proposition 1 then translates to

$$\lambda _ { 2 } ^ { 2 } ( c ^ { 2 } \mathrm e ^ { 2 } ) \gg w _ { 2 } \,, \quad \lambda _ { 2 } ^ { 2 } ( \mathrm e ) = \frac { \pi } { 4 } \mod 2 \pi \,. \quad ( 8 4 )$$

Note that the fine-tuning condition can still be satisfied: for example, for a fixed coupling λ 2 we find c such that

$$\lambda _ { 2 } ^ { 2 } ( \text{ce} ) = \frac { \pi } { 4 } \mod 2 \pi \,,$$

and we can always ensure that c is sufficiently large so that

$$\lambda _ { 2 } ^ { 2 } ( c ^ { 2 } e ^ { 2 } ) = c \cdot e \left ( \frac { \pi } { 4 } + 2 \pi n \right ) \gg w _ { 2 } \quad \ \ ( 8 6 )$$

in order to satisfy both conditions in (84).

In terms of the rescaled coupling constants, the strong-coupling condition implies that for fixed λ 2 and c /greatermuch 1, we get

$$\rho _ { \text{EB} } \xrightarrow { c \gg 1 } \begin{pmatrix} \widetilde { P } _ { - } & 0 & 0 & 0 \\ 0 & \widetilde { P } _ { + } & \widetilde { P } _ { + } & 0 \\ 0 & \widetilde { P } _ { + } & \widetilde { P } _ { + } & 0 \\ 0 & 0 & 0 & \widetilde { P } _ { - } \end{pmatrix}.$$

Using negativity N as an entanglement measure [99],

$$\mathcal { N } [ \rho _ { \text{EB} } ] = \frac { | | \rho _ { \text{EB} } ^ { \Gamma _ { \text{E} } } | | _ { 1 } - 1 } { 2 } = \frac { 1 } { 2 } e ^ { - 2 \lambda _ { 2 } ^ { 2 } w _ { 2 } } \quad \quad ( 8 8 )$$

where ρ Γ e eb is the partial transpose of ρ eb and || · || 1 is the trace norm. This means that ρ eb is always an entangled state for any finite coupling λ 2 /lessmuch λ 1 and the channel is able to preserve some quantum entanglement.

Finally, we obtain the perfect quantum channel in the idealized limit where c → ∞ and λ 2 → 0 such that λ 1 is finite and still satisfies (84), i.e.,

$$\rho _ { \text{EB} } = \frac { 1 } { 2 } \begin{pmatrix} 0 & 0 & 0 & 0 \\ 0 & 1 & 1 & 0 \\ 0 & 1 & 1 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} \,,$$

i.e., a Bell state with maximum negativity N [ ρ eb ] = 1 / 2. Therefore, in this limit the coherent information I c (Φ) → 1, which corresponds to a perfect quantum channel that can transmit a maximally entangled state perfectly.

11 The unit of λ i is [Length] n -3 2 , and we have assumed for simplicity that λ i have the same units since f i are chosen to be ( cf. the end of Section III A)

## B. Alice and Bob are spacelike separated

The preceding analysis shows that the ability to implement the approximate 'SWAP' gate depends on how closely one is able to satisfy the fine-tuning and strongcoupling conditions in Proposition 1. While this is straightforward for Alice, Proposition 2 tells us that for Bob to emulate ˆ ( φ f i ) defined at some (global coordinate) time t ( τ a ), Bob's spacetime smearing functions g i at time t ( τ b ) has to be somewhat delocalized over some region of spacetime to catch all the signals from Alice. Since the quantum channel Φ is built out of a relativistic quantum field, we ought to be able to show that there can be no transmission of both classical and quantum information through Φ.

Formally, we start from the CCR which tells us that

$$E ( f _ { i }, g _ { j } ) = 0 \quad \forall i, j = 1, 2 \quad \quad ( 9 0 )$$

when Alice and Bob are spacelike-separated, i.e., supp( g i ) ∩ J ± (supp( f i )) = ∅ . Note that the Hadamard function H ( f , g i j ) is generically nonzero for spacelikeseparated regions due to entanglement contained in the field state [5, 6, 38]. The calculations proceeds similarly as in the previous section, except that Bob will have his own arbitrary smearing functions and their respective UDW interaction unitaries are

$$\hat { U } _ { \text{B} } = \sum _ { x, z \in \{ \pm \} } \hat { P } _ { z } \hat { P } _ { x } \otimes W ( - z E g _ { 1 } ) W ( - x E g _ { 2 } ). \quad ( 9 1 b )$$

The condition for spacelike separation in Eq. (90) will be enforced later in the calculation once everything is written in terms of two-point functions.

We can analyse this scenario very generally by allowing Alice and the environment to start in any arbitrary state ρ ea instead of a maximally entangled state or purifications of any initial state ρ a , 0 . The expression for the density matrix ρ eb now reads

$$\rho _ { \text{EB} } & = \sum _ { x _ { i }, z _ { i } } \omega ( \hat { O } ) \, \text{Tr} _ { \text{A} } \left [ \hat { P } _ { x _ { 4 } } \hat { P } _ { z _ { 4 } } \rho _ { \text{AE} } \hat { P } _ { - z _ { 1 } } \hat { P } _ { - x _ { 1 } } \right ] \\ & \otimes \hat { P } _ { - z _ { 3 } } \hat { P } _ { - x _ { 3 } } | + _ { y } \rangle \langle + _ { y } | _ { \text{B} } \hat { P } _ { x _ { 2 } } \hat { P } _ { z _ { 2 } }$$

where the projection operators act on Alice or Bob and ˆ is now defined to be O

$$\text{the} \\ \text{ly} \\ \text{is} \\ \text{ive}$$

$$\hat { U } _ { \text{A} } = \sum _ { x, z \in \{ \pm \} } \hat { P } _ { x } \hat { P } _ { z } \otimes W ( x E f _ { 2 } ) W ( z E f _ { 1 } ) \quad ( 9 1 \text{a} ) \quad \text{$\text{$\text{$this$e$}$}$}$$

Applying the third Weyl relation and using Eq. (17), this evaluates to

$$\omega ( \hat { O } ) & = e ^ { - \frac { i } { } ( z _ { 1 } - z _ { 4 } ) ( z _ { 1 } + x _ { 4 } ) E ( f _ { 1 }, f _ { 2 } ) } e ^ { - \frac { i } { } ( z _ { 2 } + z _ { 3 } ) ( z _ { 3 } - z _ { 2 } ) E ( g _ { 1 }, g _ { 2 } ) } e ^ { - \frac { i } { } ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } - z _ { 4 } ) E ( f _ { 1 }, g _ { 1 } ) } \\ & \times e ^ { - \frac { i } { } ( z _ { 1 } - z _ { 4 } ) ( z _ { 2 } + z _ { 3 } ) E ( f _ { 1 }, g _ { 2 } ) } e ^ { - \frac { i } { } ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } - z _ { 4 } ) E ( f _ { 2 }, g _ { 1 } ) } e ^ { - \frac { i } { } ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } - z _ { 4 } ) E ( f _ { 2 }, g _ { 2 } ) } \\ & \times e ^ { - \frac { i } { } ( ( z _ { 1 } + z _ { 4 } ) ( z _ { 1 } + z _ { 4 } ) W ( f _ { 1 }, f _ { 1 } ) + ( z _ { 1 } + z _ { 4 } ) ( z _ { 1 } + z _ { 4 } ) W ( f _ { 2 }, f _ { 2 } ) + ( z _ { 2 } + z _ { 3 } ) W ( g _ { 1 }, g _ { 1 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 2 } + z _ { 3 } ) W ( g _ { 2 }, g _ { 2 } ) } \\ & \times e ^ { - \frac { i } { } ( ( x _ { 1 } + z _ { 4 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 1 }, f _ { 2 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 1 }, g _ { 1 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 1 }, g _ { 2 } ) ) } \\ & \times e ^ { - \frac { i } { } ( ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 2 }, g _ { 1 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 2 }, g _ { 2 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 2 } + z _ { 3 } ) H ( g _ { 1 }, g _ { 2 } ) ) } \,.$$

This is the most general form that ω O ( ˆ ) can take in terms of Alice and Bob's two-point functions. We now enforce spacelike-separation using Eq. (90), giving us

$$x \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \\ \omega ( \hat { O } ) & = e ^ { - \frac { 1 } { 2 } ( x _ { 1 } - z _ { 4 } ) ( x _ { 1 } + x _ { 4 } ) E ( f _ { 1 }, f _ { 2 } ) } e ^ { - \frac { 1 } { 2 } ( z _ { 4 } + z _ { 3 } ) ( x _ { 3 } - z _ { 2 } ) E ( g _ { 1 }, g _ { 2 } ) } \\ & \quad \times e ^ { - \frac { 1 } { 2 } ( ( x _ { 1 } + z _ { 4 } ) ( z _ { 1 } + z _ { 4 } ) W ( f _ { 1 }, f _ { 1 } ) + ( x _ { 1 } + z _ { 4 } ) ( x _ { 1 } + z _ { 4 } ) W ( f _ { 2 }, f _ { 2 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 2 } + z _ { 3 } ) W ( g _ { 1 }, g _ { 1 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 2 } + z _ { 4 } ) H ( f _ { 1 }, g _ { 1 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 1 }, g _ { 2 } ) ) } \\ & \quad \times e ^ { - \frac { 1 } { 2 } ( ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 2 }, g _ { 1 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 1 } + z _ { 4 } ) H ( f _ { 2 }, g _ { 2 } ) + ( z _ { 2 } + z _ { 3 } ) ( z _ { 2 } + z _ { 3 } ) H ( g _ { 1 }, g _ { 2 } ) ) } \,,$$

The density matrix ρ eb can now be calculated

$$\rho _ { _ { E B } } = \rho _ { _ { E } } \otimes \begin{pmatrix} \frac { 1 } { 2 } & X \\ X ^ { * } & \frac { 1 } { 2 } \end{pmatrix} \quad \quad ( 9 6 ) \quad \quad,$$

$$\begin{array} { c } \cdot \\ \text{with } \rho _ { E } = \text{Tr} _ { \Lambda } [ \rho _ { E A } ] \text{ and} \\ X = - \frac { i } { 2 } e ^ { - 2 W ( g _ { 1 }, g _ { 1 } ) } \\ \left [ e ^ { - 2 W ( g _ { 2 }, g _ { 2 } ) } \cosh ( 2 H ( g _ { 1 }, g _ { 2 } ) ) + \sin ( 2 E ( g _ { 1 }, g _ { 2 } ) ) \right ]. \\ \end{array}$$

Observe that the density matrix ρ eb depends only on

the initial condition of the environment and Bob's local two-point functions, i.e. not involving Alice's smearing functions at all, and is clearly separable across B and E . Consequently, the quantum channel preserves no entanglement at all from the joint system EA (viewing it as 1 1 ⊗ Φ). In fact, since we start from any arbitrary state ρ ea , we have effectively shown that Φ is a PPT channel [100, 101], and hence also an entanglement-breaking channel (because PPT is equivalent to separability for two qubits), and these are known to fall under the class of anti-degradable channels with zero quantum capacity Q [Φ] = 0 [98].

From an information-theoretic viewpoint, zero quantum capacity says nothing about the ability of the quantum channel to transmit classical information, since the classical channel capacity C [Φ] is an upper bound for quantum capacity, i.e., C [Φ] ≥ Q [Φ]. Quantum channels can also transmit classical bits by, for instance, considering only diagonal density matrix

$$\rho _ { \Lambda, 0 } = \begin{pmatrix} \alpha & 0 \\ 0 & 1 - \alpha \end{pmatrix}, \quad \alpha \in [ 0, 1 ] \,. \quad \begin{pmatrix} \text{def} \\ ( 9 7 ) & \text{ano} \\ \text{the} \end{pmatrix}$$

and indeed it has been shown that for certain deltacoupled UDW models, the classical capacity C [Φ] can be computed exactly and is zero if Alice and Bob are spacelike [41, 42]. Our quantum channel is more complicated as Alice and Bob each perform a sequence of two simple-generated unitaries, but it is not necessary for us to find the explicit formula for C [Φ] to show that it is zero. This follows directly from the fact that the density matrix in Eq. (96) does not depend on ρ a , 0 at all, therefore the quantum channel carries no information whatsoever about Alice's initial state, hence C [Φ] = 0.

We stress that one of the main features of a relativistic quantum channel is precisely that one should be able to show that it respects the causal structure of spacetime by construction. This requires one to be able to show, for instance, that no information transmission is allowed at spacelike separation, or more generally between causally disconnected observers. In fact we can say more: our formalism does not require us to specify anything about the (violation of) strong Huygens' principle [38, 102-106], since everything is encoded in the causal propagator. For example, in (3 + 1)-dimensional Minkowski spacetime, Bob can be timelike-separated from Alice but for a massless scalar field, the causal propagator has zero support in the interior of the light cone (see, e.g., [38, 105]). In this case, if Bob's spacetime smearing functions to the field g i are strictly in the timelike interior and supp( g i ) ∩ J + (supp f j ) = ∅ for any i, j = 1 2, then , Q [Φ] = C [Φ] = 0.

## C. Intermediate scenarios

In the two preceding subsections, we have analyzed two opposite regimes: the ideal case with Ef i = Eg i and the worst case with Alice and Bob spacelike-separated. Let us make some comments on intermediate regimes.

When the coherent information of a quantum channel is not close to 1, e.g., when 0 ≤ I c (Φ) /lessmuch 1, it is in general very difficult to make quantitative claims about the quantum channel capacity for the reasons like superadditivity described in Section III. In fact, even for quantum channels with I c (Φ) = 0, only small subsets of these are known to have zero quantum capacity, namely the anti-degradable channels [86, 98, 107] and the PPT channels [100] (see also [101]). For intermediate scenarios, such as when Bob is only partially covering J ± (supp( f i )) or implementing faulty SWAP gates we are not able to give generic quantitative statements about Q [Φ], as is the case for generic quantum communication channels.

The above discussion is not necessarily the only (or natural) way to view this. One way to think about this from the perspective of the quantum field is this: if Bob manages to decode Alice's state at time t , the 'logical qubit' state c 1 | + α 〉 -ic 2 |-α 〉 in Eq. (64) would have been taken out by Bob, and the field would now change into a new state that has no logical qubit in it (essentially because Bob swapped out the state into his physical detector). Charlie, who couples to the field at any time t + /epsilon1 for any arbitrary /epsilon1 &gt; 0, would no longer have any logical qubit to extract from the field, so Charlie cannot extract another qubit for free. Yet another way to see this is that Bob's interaction requires one to find g i in J + (supp( f j )) that can 'undo' Alice's operation. Therefore, if we choose to pick g i that does not cover the full region needed for the decoding, then Bob fails to implement his ideal ˆ U B and instead performs a different unitary operation. In this interpretation, we would

However, in some situations we may try to use nocloning type argument to conclude that Q [Φ] = 0 in some scenarios as was done in [52, 53] (see also [100, 108, 109]). Suppose Bob is allowed to couple only to at most half of Alice's signal at some future time t that defines the Cauchy slice for Bob's interaction. Consider another observer Charlie who now chooses to couple to the other half of the region at the same time slice so that jointly the interaction support is the full region needed for the perfect quantum channel. Since Charlie's interaction is by construction spacelike separated from Bob, Charlie's interaction is invisible to Bob (and vice versa). If Bob and Charlie were able to individually decode Alice's qubit, this would violate the no-cloning theorem, so we can conclude that the quantum capacity must be zero if the UDW interactions are set up this way 12 .

12 It is worth emphasizing that the protocol we have (and also in [53]) are generically not ideal, since for any finite interaction strengths the 'SWAP' gates are never perfect and the 'logical qubit' in the field is only approximate - consequently one has to be careful about using exact cloning argument. In practice, it would be better to think of the scenario in terms of approximate cloning. Indeed, in Figure 8 of [53] the coherent information vanishes smoothly as one varies Bob's interaction support when analyzing the cloning-type argument.

say that Bob is unable to perform the required 'SWAP' gate to decode Alice's state. Showing that Bob has zero quantum capacity would mean showing that Bob's operations somehow completely depolarize the qubit state and lose all the quantum information in it - we will not pursue this further in this paper and leave this for future work.

## VI. DISCUSSION AND OUTLOOK

We motivated our work through the question: What is the information-carrying capacity of quantum fields? The UDW framework allows us to formulate a precise statement of this question, in terms of quantum capacity of the relativistic quantum channel between Alice and Bob who try to transmit quantum information between their qubit detectors through the quantum field. In this work, we have successfully generalized the results of [49-53], in particular [53], where we construct a quantum channel that can transmit quantum information arbitrarily well through a relativistic quantum scalar field in arbitrary curved spacetimes.

Crucially, in this generalized formulation it becomes clear that there are two essential features of a nearperfect quantum channel that can transmit quantum information arbitrarily well: (1) Alice and Bob's interactions need to satisfy the so-called fine-tuning and strong coupling conditions in Proposition 1 in order to approximate an ideal 'SWAP' gates; (2) Bob needs to couple to the field at appropriate spacetime regions, fixed by Proposition 2, so that he could execute the approximately ideal 'SWAP' gates. Indeed, where Bob should be is fixed by the causal propagator that determines where Alice's information propagates in spacetime. The asymmetry between Alice and Bob can be viewed as the consequence of the asymmetry between encoding the qubit state into the field, with the logical basis being the coherent states |± α 〉 , and decoding the state while taking into account relativistic causality. In this sense, our construction constitutes the most general relativistic qubit channel within the UDW framework that can transmit quantum information arbitrarily well in appropriate regimes and we also provided analyses when the

The features of our construction is that the resulting quantum channel, and hence the quantum capacity, are manifestly covariant, respect the causal structure of spacetime, and are independent of the details of the background geometry, topology, and the choice of Hilbert space (quasifree) representations of the CCR algebra. Furthermore, in this generalization we no longer need to deal with (violation of) the strong Huygens' principle that may occur in curved spacetimes for nonconformally invariant fields (see, e.g.,[49, 50, 102-104]), since this is fully encoded in the causal propagator. We also showed explicitly that the quantum and classical capacity is zero for spacelike separated observers, as we expect from quantum communication channel built from a relativistic quantum field.

quantum channel is sub-optimal.

There are several future directions that merit investigation, at least from a more fundamental standpoint:

- · Is it possible to construct a perfect relativistic quantum channel using higher-dimensional detector models (using qudit detector models [110, 111] rather than qubits)? A priori it is not obvious how to construct logical qudits in the field using only controlled displacement operators and what kind of interaction Hamiltonians are needed to do this, if it exists at all. It is also natural to consider several emitters and receivers, exploiting the fact that one can interpret displacement operations can be viewed as Weyl generators with respect to some spacetime smearing functions [58].
- · As mentioned in [53], one may be interested in more directional coupling, which can be envisaged in (1+1)D models (see, e.g., [37, 52]). Our protocol works automatically simply because only the twopoint functions change to reflect that the leftand right- movers of the field modes are independent. A more interesting setup would be to study how realistic settings involving waveguides and non-trivial dispersion relations can be adapted for this purpose.
- · Is the information-carrying capacity different between different fields with different spins and statistics? For example, how does one compare the scalar field with a fermionic field? The issue one has to deal with is that we have to separate the effects coming from the model dependence of the coupling (e.g., we cannot couple ˆ σ x ⊗ ˆ where ψ ψ is the fermionic field) from truly field-theoretic physics that contribute to the quantum capacity, if we wish to do it in the same way as we have done in this paper (see, e.g., some proposed fermionic couplings in RQI literature [112, 113]).
- · Our construction relies on the fact that the quantum field is non-interacting (Gaussian). There are very few studies involving the UDW model with interacting QFTs such as the φ 4 theory, even perturbatively and for a single qubit (see, e.g., [24] for one example). The effects of interactions are unclear for (long-distance) quantum communication, both perturbatively and non-perturbatively. A possibly relevant testbed would be to consider known rigorous constructions in (1+1)D interacting QFTs (see, e.g., [114-118]).
- · Is it possible to systematically study quantum capacity in the non-relativistic quantum many-body settings? Somewhat related studies to our setup are, for example, [119-121]. It would be interesting to see the similarities and differences when they are compared with setups involving relativistic quantum fields.

We leave these problems open for the future.

## ACKNOWLEDGMENT

M. K. acknowledges support from the Princeton International Internship Program (IIP). E. T. acknowl-

- [1] R. Verch and R. F. Werner, Distillability and positivity of partial transposes in general quantum field systems, Reviews in Mathematical Physics 17 , 545 (2005).
- [3] S. Hollands and A. Ranallo, Channel divergences and complexity in algebraic qft, Communications in Mathematical Physics 404 , 927 (2023).
- [2] S. Hollands and K. Sanders, Entanglement measures and their properties in quantum field theory , Vol. 34 (Springer, 2018).
- [4] K. Sanders, On separable states in relativistic quantum field theory, Journal of Physics A: Mathematical and Theoretical 56 , 505201 (2023).
- [6] S. J. Summers and R. Werner, Bell's inequalities and quantum field theory. I. General setting, J. Math. Phys. 28 , 2440 (1987).
- [5] S. J. Summers and R. Werner, The vacuum violates Bell's inequalities, Phys. Lett. A 110 , 257 (1985).
- [7] H. Casini, M. Huerta, J. M. Magán, and D. Pontello, Entanglement entropy and superselection sectors. Part I. Global symmetries, Journal of High Energy Physics 2020 , 1 (2020).
- [9] R. Longo and F. Xu, Relative entropy in CFT, Advances in Mathematics 337 , 139 (2018).
- [8] H. Casini, M. Huerta, J. M. Magán, and D. Pontello, Entropic order parameters for the phases of QFT, Journal of High Energy Physics 2021 , 1 (2021).
- [10] C. J. Fewster and R. Verch, Quantum fields and local measurements, Communications in Mathematical physics 378 , 851 (2020).
- [12] J. Polo-Gómez, L. J. Garay, and E. Martín-Martínez, A detector-based measurement theory for quantum field theory, Phys. Rev. D 105 , 065003 (2022).
- [11] H. Bostelmann, C. J. Fewster, and M. H. Ruep, Impossible measurements require impossible apparatus, Phys. Rev. D 103 , 025017 (2021).
- [13] I. Jubb, Causal state updates in real scalar quantum field theory, Phys. Rev. D 105 , 025003 (2022).
- [15] L. van Luijk, R. Schwonnek, A. Stottmeister, and R. F. Werner, The Schmidt rank for the commuting operator framework, Communications in Mathematical Physics 405 , 152 (2024).
- [14] N. Pranzini and E. Keski-Vakkuri, Detector-based measurements for QFT: two issues and an AQFT proposal, arXiv preprint arXiv:2310.06596 (2023).
- [16] L. van Luijk, A. Stottmeister, R. F. Werner, and H. Wilming, Embezzlement of entanglement, quantum fields, and the classification of von Neumann algebras, arXiv preprint arXiv:2401.07299 (2024).
- [18] L. C. B. Crispino, A. Higuchi, and G. E. A. Matsas, The Unruh effect and its applications, Rev. Mod. Phys. 80 , 787 (2008).
- [17] S. W. Hawking, Particle creation by black holes, Comm. Math. Phys. 43 , 199 (1975).
- [19] W. G. Unruh, Notes on black-hole evaporation, Phys. Rev. D 14 , 870 (1976).

edges funding from the Munich Center for Quantum Science and Technology (MCQST), funded by the Deutsche Forschungsgemeinschaft (DFG) under Germany's Excellence Strategy (EXC2111 - 390814868).

- [20] B. S. DeWitt, Quantum gravity: the new synthesis, in General Relativity: An Einstein centenary survey , edited by S. W. Hawking and W. Israel (Cambridge University Press, Cambridge, 1979) pp. 680-745.
- [22] B. A. Juárez-Aubry and J. Louko, Quantum fields during black hole formation: how good an approximation is the Unruh state?, Journal of High Energy Physics 2018 , 1 (2018).
- [21] B. A. Juárez-Aubry and J. Louko, Onset and decay of the 1 + 1 Hawking-Unruh effect: what the derivativecoupling detector saw, Class. and Quantum Gravity 31 , 245007 (2014).
- [23] E. Tjoa and R. B. Mann, Unruh-DeWitt detector in dimensionally-reduced static spherically symmetric spacetimes, Journal of High Energy Physics 2022 , 1 (2022).
- [25] G. Kaplanek and E. Tjoa, Effective master equations for two accelerated qubits, Phys. Rev. A 107 , 012208 (2023).
- [24] G. Kaplanek and C. Burgess, Hot accelerated qubits: decoherence, thermalization, secular growth and reliable late-time predictions, Journal of High Energy Physics 2020 , 1 (2020).
- [26] D. Moustos and C. Anastopoulos, Non-markovian time evolution of an accelerated qubit, Phys. Rev. D 95 , 025020 (2017).
- [28] L. Hardy, Towards quantum gravity: a framework for probabilistic theories with non-fixed causal structure, Journal of Physics A: Mathematical and Theoretical 40 , 3081 (2007).
- [27] V. Vilasini and R. Renner, Embedding cyclic causal structures in acyclic spacetimes: no-go results for process matrices, arXiv preprint arXiv:2203.11245 (2022).
- [29] M. Zych, F. Costa, I. Pikovski, and Č. Brukner, Bell's theorem for temporal order, Nature communications 10 , 1 (2019).
- [31] N. Paunković and M. Vojinović, Causal orders, quantum circuits and spacetime: distinguishing between definite and superposed causal orders, Quantum 4 , 275 (2020).
- [30] F. Costa and S. Shrapnel, Quantum causal modelling, New Journal of Physics 18 , 063032 (2016).
- [32] B. Reznik, Entanglement from the vacuum, Found. Phys. 33 , 167 (2003).
- [34] A. Valentini, Non-local correlations in quantum electrodynamics, Phys. Lett. A 153 , 321 (1991).
- [33] B. Reznik, A. Retzker, and J. Silman, Violating Bell's inequalities in vacuum, Phys. Rev. A 71 , 042104 (2005).
- [35] A. Pozas-Kerstjens and E. Martín-Martínez, Harvesting correlations from the quantum vacuum, Phys. Rev. D 92 , 064042 (2015).
- [36] A. Pozas-Kerstjens and E. Martín-Martínez, Entanglement harvesting from the electromagnetic vacuum with hydrogenlike atoms, Phys. Rev. D 94 , 064074 (2016).

- [37] R. H. Jonsson, E. Martín-Martínez, and A. Kempf, Quantum signaling in cavity QED, Phys. Rev. A 89 , 022330 (2014).
- [39] T. R. Perche, J. Polo-Gómez, B. d. S. L. Torres, and E. Martín-Martínez, Fully relativistic entanglement harvesting, Phys. Rev. D 109 , 045018 (2024).
- [38] E. Tjoa and E. Martín-Martínez, When entanglement harvesting is not really harvesting, Phys. Rev. D 104 , 125005 (2021).
- [40] E. Tjoa, Quantum teleportation with relativistic communication from first principles, Phys. Rev. A 106 , 032432 (2022).
- [42] E. Tjoa and K. Gallock-Yoshimura, Channel capacity of relativistic quantum communication with rapid interaction, Phys. Rev. D 105 , 085011 (2022).
- [41] A. G. S. Landulfo, Nonperturbative approach to relativistic quantum communication channels, Phys. Rev. D 93 , 104019 (2016).
- [43] A. Lapponi, D. Moustos, D. E. Bruschi, and S. Mancini, Relativistic quantum communication between harmonic oscillator detectors, Phys. Rev. D 107 , 125010 (2023).
- [45] H.-K. Lo and H. Chau, Why quantum bit commitment and ideal quantum coin tossing are impossible, Physica D: Nonlinear Phenomena 120 , 177 (1998), proceedings of the Fourth Workshop on Physics and Consumption.
- [44] A. Kent, Unconditionally secure bit commitment, Phys. Rev. Lett. 83 , 1447 (1999).
- [46] E. Adlam and A. Kent, Device-independent relativistic quantum bit commitment, Phys. Rev. A 92 , 022315 (2015).
- [48] V. Vilasini, C. Portmann, and L. del Rio, Composable security in relativistic quantum cryptography, New Journal of Physics 21 , 043057 (2019).
- [47] H. Buhrman, N. Chandran, S. Fehr, R. Gelles, V. Goyal, R. Ostrovsky, and C. Schaffner, Positionbased quantum cryptography: Impossibility and constructions, SIAM Journal on Computing 43 , 150 (2014).
- [49] A. Blasco, L. J. Garay, M. Martín-Benito, and E. Martín-Martínez, Violation of the Strong Huygen's Principle and Timelike Signals from the Early Universe, Phys. Rev. Lett. 114 , 141103 (2015).
- [51] P. Simidzija and E. Martín-Martínez, Information carrying capacity of a cosmological constant, Phys. Rev. D 95 , 025002 (2017).
- [50] A. Blasco, L. J. Garay, M. Martín-Benito, and E. Martín-Martínez, Timelike information broadcasting in cosmology, Phys. Rev. D 93 , 024055 (2016).
- [52] R. H. Jonsson, K. Ried, E. Martín-Martínez, and A. Kempf, Transmitting qubits through relativistic fields, Journal of Physics A: Mathematical and Theoretical 51 , 485301 (2018).
- [54] I. B. Barcellos and A. G. S. Landulfo, Relativistic quantum broadcast channel, Phys. Rev. D 109 , 065020 (2024).
- [53] P. Simidzija, A. Ahmadzadegan, A. Kempf, and E. Martín-Martínez, Transmission of quantum information through quantum fields, Phys. Rev. D 101 , 036014 (2020).
- [55] E. Martín-Martínez, T. R. Perche, and B. de S. L. Torres, General relativistic quantum optics: Finite-size particle detector models in curved spacetimes, Phys. Rev. D 101 , 045017 (2020).
- [56] C. Lima, E. Patterson, E. Tjoa, and R. B. Mann, Unruh phenomena and thermalization for qudit detectors, Phys. Rev. D 108 , 105020 (2023).
- [58] E. Tjoa, Nonperturbative simple-generated interactions with a quantum field for arbitrary Gaussian states, Phys. Rev. D 108 , 045003 (2023).
- [57] R. Lopp and E. Martín-Martínez, Quantum delocalization, gauge, and quantum optics: Light-matter interaction in relativistic quantum information, Phys. Rev. A 103 , 013703 (2021).
- [59] E. Tjoa and F. Gray, The Unruh-DeWitt model and its joint interacting Hilbert space, Journal of Physics A: Mathematical and Theoretical 57 , 325301 (2024).
- [61] A. Lapponi, J. Louko, and S. Mancini, Making two particle detectors in flat spacetime communicate quantumly, Phys. Rev. D 110 , 025018 (2024).
- [60] S. Mancini, R. Pierini, and M. M. Wilde, Preserving information from the beginning to the end of time in a robertsonâ€'walker spacetime, New Journal of Physics 16 , 123049 (2014).
- [62] R. M. Wald, Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics , Chicago Lectures in Physics (University of Chicago Press, 1994).
- [64] I. Khavkine and V. Moretti, Algebraic QFT in curved spacetime and quasifree Hadamard states: An introduction, Mathematical Physics Studies , 191-251 (2015).
- [63] C. J. Fewster and K. Rejzner, Algebraic quantum field theory - an introduction (2019), arXiv:1904.04051 [hep-th].
- [65] B. S. Kay and R. M. Wald, Theorems on the uniqueness and thermal properties of stationary, nonsingular, quasifree states on spacetimes with a bifurcate killing horizon, Physics Reports 207 , 49 (1991).
- [67] O. Bratteli, Inductive Limits of Finite Dimensional C*Algebras, Transactions of the American Mathematical Society 171 , 195 (1972).
- [66] T.-P. Hack, Cosmological applications of algebraic quantum field theory in curved spacetimes , Vol. 6 (Springer, 2015).
- [68] E. Poisson, A Relativist's Toolkit: The Mathematics of Black-Hole Mechanics (Cambridge University Press, 2009).
- [70] M. H. Ruep, Weakly coupled local particle detectors cannot harvest entanglement, Classical and Quantum Gravity 38 , 195029 (2021).
- [69] R. Wald, General Relativity (University of Chicago Press, 2010).
- [71] M. J. Radzikowski, Micro-local approach to the Hadamard condition in quantum field theory on curved space-time, Communications in Mathematical Physics 179 , 529 (1996).
- [73] P. C. Martin and J. Schwinger, Theory of manyparticle systems. i, Phys. Rev. 115 , 1342 (1959).
- [72] R. Kubo, Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems, Journal of the physical society of Japan 12 , 570 (1957).
- [74] H. Spohn, Ground state(s) of the spin-boson Hamiltonian, Communications in Mathematical Physics 123 , 277 (1989).
- [75] D. Hasler, B. Hinrichs, and O. Siebert, On existence of ground states in the spin boson model, Communications in Mathematical Physics 388 , 419 (2021).

- [76] M. Fannes, B. Nachtergaele, and A. Verbeure, The equilibrium states of the spin-boson model, Communications in mathematical physics 114 , 537 (1988).
- [78] P. Simidzija, R. H. Jonsson, and E. Martín-Martínez, General no-go theorem for entanglement extraction, Phys. Rev. D 97 , 125002 (2018).
- [77] E. Martín-Martínez, T. R. Perche, and B. d. S. L. Torres, Broken covariance of particle detector models in relativistic quantum information, Phys. Rev. D 103 , 025007 (2021).
- [79] P. Simidzija and E. Martín-Martínez, All coherent field states entangle equally, Phys. Rev. D 96 , 025020 (2017).
- [81] N. K. Kollas and D. Moustos, An exactly solvable relativistic quantum Otto engine, Phys. Rev. D 109 , 065025 (2024).
- [80] K. Gallock-Yoshimura, Relativistic quantum otto engine: instant work extraction from a quantum field, Journal of High Energy Physics 2024 , 1 (2024).
- [82] J. Polo-Gómez and E. Martín-Martínez, Nonperturbative method for particle detectors with continuous interactions, Phys. Rev. D 109 , 045014 (2024).
- [84] J. Dimock, Algebras of local observables on a manifold, Communications in Mathematical Physics 77 , 219 (1980).
- [83] E. Poisson, A. Pound, and I. Vega, The motion of point particles in curved spacetime, Living Reviews in Relativity 14 , 1 (2011).
- [85] B. S. Kay, Linear spin-zero quantum fields in external gravitational and scalar fields, Commun. Math. Phys 62 , 55 (1978).
- [87] M. M. Wilde, Preface to the second edition, in Quantum Information Theory (Cambridge University Press, 2017) p. xi-xii.
- [86] S. Khatri and M. M. Wilde, Principles of quantum communication theory: A modern approach, arXiv preprint arXiv:2011.04672 (2020).
- [88] L. Gyongyosi, S. Imre, and H. V. Nguyen, A survey on quantum channel capacities, IEEE Communications Surveys &amp; Tutorials 20 , 1149 (2018).
- [90] B. Schumacher and M. D. Westmoreland, Relative entropy in quantum information theory, Contemporary Mathematics 305 , 265 (2002).
- [89] A. Holevo, Quantum channel capacities, Quantum Electronics 50 , 440 (2020).
- [91] S. Lloyd, Capacity of the noisy quantum channel, Phys. Rev. A 55 , 1613 (1997).
- [93] I. Devetak and A. Winter, Classical data compression with quantum side information, Phys. Rev. A 68 , 042301 (2003).
- [92] P. W. Shor, The quantum channel capacity and coherent information, in lecture notes, MSRI Workshop on Quantum Computation (2002).
- [94] M. B. Hastings, Superadditivity of communication capacity using entangled inputs, Nature Physics 5 , 255 (2009).
- [96] F. Leditzky, D. Leung, V. Siddhu, G. Smith, and J. A. Smolin, Generic nonadditivity of quantum capacity in simple channels, Phys. Rev. Lett. 130 , 200801 (2023).
- [95] G. Smith and J. Yard, Quantum communication with zero-capacity channels, Science 321 , 1812 (2008).
- [97] T. S. Cubitt, M. B. Ruskai, and G. Smith, The structure of degradable quantum channels, Journal of Mathematical Physics 49 , 102104 (2008).
- [98] M. Horodecki, P. W. Shor, and M. B. Ruskai, Entanglement breaking channels, Reviews in Mathematical Physics 15 , 629 (2003).
- [100] G. Smith and J. A. Smolin, Detecting incapacity of a quantum channel, Phys. Rev. Lett. 108 , 230507 (2012).
- [99] G. Vidal and R. F. Werner, Computable measure of entanglement, Phys. Rev. A 65 , 032314 (2002).
- [101] S. Singh and N. Datta, Detecting positive quantum capacities of quantum channels, npj Quantum Information 8 , 50 (2022).
- [103] S. Sonego and V. Faraoni, Huygens' principle and characteristic propagation property for waves in curved space-times, Journal of Mathematical Physics 33 , 625 (1992), https://doi.org/10.1063/1.529798.
- [102] R. G. McLenaghan, On the validity of Huygens' principle for second order partial differential equations with four independent variables. Part I : derivation of necessary conditions, Annales de l'I.H.P. Physique théorique 20 , 153 (1974).
- [104] V. Faraoni, Massive spin zero fields in cosmology and the tail-free property, Symmetry 11 , 10.3390/sym11010036 (2019).
- [106] R. H. Jonsson, D. Q. Aruquipa, M. Casals, A. Kempf, and E. Martín-Martínez, Communication through quantum fields near a black hole, Phys. Rev. D 101 , 125005 (2020).
- [105] E. Martín-Martínez, Causality issues of particle detector models in QFT and quantum optics, Phys. Rev. D 92 , 104019 (2015).
- [107] M. M. Wolf and D. Pérez-García, Quantum capacities of channels with small environment, Phys. Rev. A 75 , 012303 (2007).
- [109] C. H. Bennett, D. P. DiVincenzo, and J. A. Smolin, Capacities of quantum erasure channels, Phys. Rev. Lett. 78 , 3217 (1997).
- [108] D. Bruß, D. P. DiVincenzo, A. Ekert, C. A. Fuchs, C. Macchiavello, and J. A. Smolin, Optimal universal and state-dependent quantum cloning, Phys. Rev. A 57 , 2368 (1998).
- [110] C. Lima, E. Patterson, E. Tjoa, and R. B. Mann, Unruh phenomena and thermalization for qudit detectors, Phys. Rev. D 108 , 105020 (2023).
- [112] J. Louko and V. Toussaint, Unruh-DeWitt detector's response to fermions in flat spacetimes, Phys. Rev. D 94 , 064027 (2016).
- [111] G. Verdon-Akzam, E. Martín-Martínez, and A. Kempf, Asymptotically limitless quantum energy teleportation via qudit probes, Phys. Rev. A 93 , 022308 (2016).
- [113] B. d. S. L. Torres, T. R. Perche, A. G. S. Landulfo, and G. E. A. Matsas, Neutrino flavor oscillations without flavor states, Phys. Rev. D 102 , 093003 (2020).
- [115] J. Glimm and A. Jaffe, A λφ 4 Quantum Field Theory without Cutoffs: II. The Field Operators and the Approximate Vacuum, Annals of Mathematics 91 , 362 (1970).
- [114] J. Glimm and A. Jaffe, A λφ 4 Quantum Field Theory without Cutoffs. I, Phys. Rev. 176 , 1945 (1968).
- [116] J. Glimm and A. Jaffe, A λφ 4 Quantum Field Theory without Cutoffs: III. The Physical Vacuum, Acta Mathematica 125 , 203 (1970).
- [117] R. Schrader, Yukawa quantum field theory in two space-time dimensions without cutoffs, Annals of Physics 70 , 412 (1972).

- [118] S. J. Summers, A perspective on constructive quantum field theory, arXiv preprint arXiv:1203.3991 (2012).
- [119] V. Giovannetti and R. Fazio, Information-capacity description of spin-chain correlations, Phys. Rev. A 71 , 032314 (2005).
- [120] S. Bose, Quantum communication through an unmodulated spin chain, Phys. Rev. Lett. 91 , 207901 (2003).
- [121] A. Bayat, D. Burgarth, S. Mancini, and S. Bose, Memory effects in spin-chain channels for information transmission, Phys. Rev. A 77 , 050306 (2008).
Michael Kasprzak Erickson Tjoa
Transmission of quantum information through quantum fields in curved spacetimes
We construct a relativistic quantum communication channel between two
localized qubit systems, mediated by a relativistic quantum field, that can
achieve the theoretical maximum for the quantum capacity in arbitrary curved
spacetimes using the Unruh-DeWitt detector formalism. Using techniques from
algebraic quantum field theory, we express the quantum capacity of the quantum
communication channel purely in terms of the correlation functions of the field
and the causal propagator for the wave equation. Consequently, the resulting
quantum channel, and hence the quantum capacity, are by construction manifestly
covariant, respect the causal structure of spacetime, and are independent of
the details of the background geometry, topology, and the choice of Hilbert
space (quasifree) representations of the field.