## Quantum repeater protocol for deterministic distribution of macroscopic entanglement

Alexey N. Pyrkov, 1, ∗ Ilia D. Lazarev, 1 and Tim Byrnes 2, 3, 4, 5, †

1 Federal Research Center of Problems of Chemical Physics and Medicinal Chemistry RAS, Acad. Semenov av. 1, Chernogolovka, Moscow region, Russia, 142432 2 New York University Shanghai, 567 West Yangsi Road, Shanghai, 200126, China; NYU-ECNU Institute of Physics at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China; Shanghai Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200126, China. 3 State Key Laboratory of Precision Spectroscopy, School of Physical and Material Sciences, East China Normal University, Shanghai 200062, China 4 Center for Quantum and Topological Systems (CQTS), NYUAD Research Institute, New York University Abu Dhabi, UAE. 5 Department of Physics, New York University, New York, NY 10003, USA

(Dated: August 2, 2024)

Distributing long-distance entanglement is a fundamental goal that is necessary for a variety of tasks such as quantum communication, distributed quantum computing, and quantum metrology. Currently quantum repeater schemes typically aim to distribute one ebit at a time, the equivalent of one Bell pair's worth of entanglement. Here we present a method to distribute a macroscopic amount of entanglement across long-distances using a number of operations that scales only linearly with the chain length. The scheme involves ensembles of qubits and entangling them with an S S z z interaction, which can be realized using atomic gas ensembles coupled by a shared optical mode. Using only local measurements on the intermediate ensembles, this leaves the ensembles at the ends of the chain entangled. We show that there are particular 'magic' interaction times that allow for distribution of entanglement with perfect fidelity, with no degradation with chain length. The scheme is deterministic, such that with suitable local conditional unitary corrections, the same entangled state can always be prepared with good approximation.

## I. INTRODUCTION

distant nodes is achieved through entanglement swapping [10-12].

Entanglement distribution in quantum networks is an important task in the field of quantum information, with potential applications such as secure communication, quantum metrology, distributed quantum computation [1-3]. Photons can carry quantum information over long distances with almost negligible decoherence and are compatible with existing telecommunication fiber technology. However, while photons can be transferred relatively efficiently between two locations in the same laboratory or city, distances exceeding this begin to be problematic due to photon loss. For example, in a typical optic fiber, the photon loss rate become significant beyond distances of ∼ 200 km and the maximum distance is estimated as ∼ 600 km [4, 5]. In classical signal transmission line this problem can be solved with optical amplifiers but it is not possible to amplify quantum signal due to the famous no-cloning theorem [6-8]. To overcome this difficulty, Briegel, Zoller and colleagues proposed the quantum repeater protocol [9], aiming to establish long-distance quantum entanglement. In the quantum repeater scheme, the transmission channel is divided into several short-distance links, where entanglement is generated between the nearest neighbor nodes using photon transmission. Then entanglement between

∗ Electronic address: pyrkov@icp.ac.ru

† Electronic address: tim.byrnes@nyu.edu

In the entanglement swapping protocol, the probabilistic nature of generating photonic entanglement (e.g. by parametric down conversion) across the elementary links means that the entanglement may not be present synchronously. One potential solution to this challenge involves storing the generated entanglement in a quantum memory until it has been successfully prepared in both links [13-15]. It has been shown that many physical systems, including single neutral atoms [16, 17], cold atoms [18], trapped ions [19], and the Nitrogen-vacancy centers in diamonds [20-22] can absorb and store incoming photons, for instance using techniques such as electromagnetically induced transparency [23]. This allows them to serve as the nodes of a quantum network with an embedded quantum memory [19, 24, 25].

As a platform to realize a quantum memory, atomic gas ensembles are an attractive candidate due to their ability to perform strong and controllable coupling between the atoms and photons. The primary advantanges of using an ensemble is that there is a collective enhancement due to the large number of atoms with long coherence times. Duan, Lukin, Cirac, and Zoller (DLCZ) proposed a scheme based on single-photon detection to implement a quantum repeater and realize scalable long-distance quantum communication with atomic ensembles [24]. In the scheme, probabilistic write-out photons are generated from two remote atomic ensembles, and the correlated photons are sent to a beam splitter generating a heralded Fock-state type entanglement between the atomic

ensembles. Kuzmich, van der Wal, and co-workers experimentally realized a quantum memory using the DLCZ protocol with trapped cold Cs atoms [26] and Rb atoms [27]. This was followed by several other works studying entanglement generation between atomic ensembles [28-35]. In the continuous variable (CV) framework, Polzik produced entanglement between atomic ensembles using a quantum nondemolition measurement approach [36, 37]. Another prominent design for a CV repeater involves creating entangled two-mode quadrature squeezed states of light using spontaneous parametric downconversion (SPDC) sources. This is followed by CV nonGaussian quantum error correction through the injection and detection of single photons [38]. The process also includes storing heralded link-level entanglement in quantum memories, as well as entanglement swapping through coherent detection [39-41] or other non-Gaussian methods [42].

Such approaches for generating entanglement between ensembles work in either the single excitation (e.g. DLCZ) or Holstein-Primakoff regime (for CV), such that the amount of entanglement is of the order of one ebit, corresponding to the entanglement in a single Bell pair. However, the very large number of atoms in an atomic ensemble provides an opportunity to create entanglement in a different regime, creating a type of macroscopic entanglement. Such an approach was discussed in works such as Refs. [43-47], where a scheme for producing macroscopic entanglement using an S S z z interaction was analyzed and proposed to be used for quantum computation. Here, we define macroscopic entanglement as von Neumann entropy or logarithmic negativity that is of order of log 2 D , where D /greatermuch 1 is the dimensionality of the subsystems that participate in the entanglement. In the case of atomic ensembles, the effective dimensionality is equal to the number of atoms [48], which is in the range of 10 3 to 10 12 . These type of states have been proposed to be used as a resource for several quantum information processing protocols, such as quantum teleportation, remote state preparation, and others [49-52]. Experimentally, the creation of many-particle entanglement localized in a single spatial location [53] and spatially separate regions [54] within one Bose-Einstein condensate (BEC) has been observed. Bell correlations have been experimentally confirmed in a BEC [53] and a thermal atomic ensemble [55]. Recently, an experiment successfully demonstrated the achievement of entanglement between a spatially split BECs [56].

In this paper, we introduce a new protocol for creating long-distance macroscopic entanglement between atomic ensembles. The basic idea is to use entanglement swapping, in a similar way to that proposed in existing quantum repeaters, but perform this at the macroscopic scale. As illustrated in Fig. 1, we consider a chain of qubit ensembles within a quantum network, where nearest neighbors can be entangled with a S S z z interaction. Only local projective measurements are required on the intermediate ensembles leaving the ensembles at the ends of

FIG. 1: Macroscopic entanglement distribution with qubit ensembles. (a) Ensembles of qubits are placed at the nodes of a quantum network. The nearest neighbors of the network are capable of being entangled with a S S z z Hamiltonian, which can generate macroscopic entanglement. The aim is to create entanglement between distant nodes of the network. A possible physical implementation is using atomic ensembles placed in cavities connected by an optical fiber, such as discussed in Ref. [44]. (b) The entanglement distribution protocol. First the ensembles at each node are polarized in the S x direction, thus preparing in the initial state (6). Here we write | + 〉 ⊗ N = | 1 √ 2 , 1 √ 2 〉〉 . Nearest neighbor nodes are entangled with the Hamiltonian (7), producing the state (12). Finally the intermediate nodes are measured, giving the state (18). Dashed lines connecting the nodes denote the entanglement.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00141v1-figure-e4ae9952-6ba0-4f15-9f8b-f0674a0e62c2.jpg)

the chain in an entangled state (Fig. 1(b)). Operations at the ensemble level are only performed for the protocol (i.e. no microscopic operations of single qubits), and the resources required for the protocol only depends upon the chain length and not the number of qubits in an ensemble. The type of state is considerably more complex than the qubit version due to the many-body entangled state created by the S S z z interaction between the atomic ensembles. Despite this, we are able to find an analytic solution at particular entanglement times. We show that there are 'magic' interaction times where entanglement distributed to the ends of the chain do not degrade with the length of the chain.

## II. ENTANGLEMENT DISTRIBUTION PROTOCOL

We now describe concretely the entanglement distribution protocol using qubit ensembles. The basic overview of the scheme is shown in Fig. 1(b). In order to keep our presentation independent of a particular physical implementation, we consider generically a system consisting of N qubits located at each node of the quantum network. We assume that the qubits within the ensemble cannot be controlled individually, but the ensembles may be controlled by total spin operations as we describe below.

Each of the ensembles is initially prepared in a spin coherent state, which we denote

$$| \alpha, \beta \rangle \rangle \equiv \prod _ { l =$$

where | 0 〉 l and | 1 〉 l denotes the computational states for the l th qubit in the ensemble. Here, α and β are arbitrary complex numbers satisfying | α | 2 + | β | 2 = 1.

Each ensemble is controlled using total spin operators

$$S ^ { \gamma } = \sum _ { l = 1 } ^ { N } \s$$

where γ ∈ { x, y, z } and σ γ l is the Pauli spin operator for the l th qubit within the ensemble. The spin operations conserve the particle number, and we assume that N is the same for each ensemble for simplicity.

Since the spin operators and the initial state are completely symmetric under particle interchange, one can write an equivalent representation using the Schwinger boson formalism [48]. The spin coherent state in the bosonic formulation is

$$| \alpha, \beta \rangle \rangle \equiv \frac { 1 } { \$$

where the bosonic creation and annihilation operators a, b obey commutation relations [ a, a † ] = [ b, b † ] = 1. The spin operators are written as

$$\begin{array} {$$

where the number operators for the two qubit states are defined as n a = a a, n † b = b b † . The eigenstates of S z operator, which we call the number states, are defined as

$$| k \rangle = \frac { ( a ^ { \dagger } ) ^ { k }$$

where |∅〉 is the vacuum state. The eigenvalues are S z | k 〉 = (2 k -N k ) | 〉 . A mapping of the number states in the qubit formulation is given in Ref. [57].

We consider a chain of M ensembles prepared in the initial state

$$| \phi _ { 0 } \rangle = | \frac { 1 } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 1 } \otimes \dots \otimes | \frac { 1 } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M }. \quad ( 6 )$$

We assume a scenario where the nearest neighbors of the chain can be entangled using a S S z z interaction. More distant pairs cannot be directly entangled. The Hamiltonian that we will consider in the chain is

$$H & = \sum _ { j = 1 } ^ { M - 1 } ( - 1 ) ^ { j } n _ { j } ^ { a } n _ { j + 1 } ^ { a }, & ( 7 ) \\ \intertext { a - a - a - a } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i } \lambda _ { i }$$

where j labels the M ensembles. The signs in the Hamiltonian are chosen in this way so that the simplest analytical formulas can be obtained, as shown later. This is the same as a S S z z interaction up to local rotations because

$$n _ { j } ^ { a } n _ { j + 1 } ^ { a } = \frac { 1 } { 4 } ( S _ { j } ^ { z } + \hat { N } _ { j } ) ( S _ { j + 1 } ^ { z } + \hat { N } _ { j + 1 } ). \quad ( 8 )$$

Numerous proposals for producing a S S z z interaction between qubits have been proposed, particular in the context of atomic ensembles and BECs. In Ref. [44] a scheme mediated by photons in an optical fiber was proposed. This was later adapted to a geometric phase gate [47]. Another approach is to use spin-dependent forces on neighboring BECs [58]. Such an interaction is also produced when one-axis squeezed states are split [59]. This type of entanglement was experimentally realized in BECs [54, 56]. In this paper, we assume that such an interaction can be induced using suitable methods, and the interaction time t is a controllable parameter.

Evolving the Hamiltonian for a time t between the ensembles gives the unitary evolution

$$U = \prod _ { j = 1 } ^ { M - 1 } e ^ { - i ( - 1 ) ^ { j } n _ { j } ^ { a } n _ { j + 1 } ^ { a } t }. \quad \quad ( 9 ) \\ - \hat { \alpha } \dots \dots \dots \dots \dots \dots$$

The state after applying the gates is

$$| \psi _ { M } \rangle = U | \phi _ { 0 } \rangle. \text{ \quad \ \ } ( 1 0 )$$

To understand the type of state that is produced by this interaction, consider first the M = 2 case which gives

$$| \psi _ { 2 } \rangle & = e ^ { i n _ { 1 } ^ { a } n _ { 2 } ^ { a } t } | \frac { 1 } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 1 } | \frac { 1 } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 } \\ & = \frac { 1 } { \sqrt { 2 ^ { N } } } \sum _ { k _ { 1 } } \sqrt { C _ { N } ^ { k _ { 1 } } } | k _ { 1 } \rangle | \frac { e ^ { i k _ { 1 } t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 }. \quad ( 1 1 ) \\ As discussed in Ref. \ [ 4 3 ]. \text{ this produces a tvpe of entan-}$$

As discussed in Ref. [43], this produces a type of entangled state where number states in the first ensemble are correlated with spin coherent states around the equator of the Bloch sphere of the second. The unitary operation U can be viewed as a conditional rotation of ensemble 2 with an angle dependent on the number state on ensemble 1.

For a chain of M ensembles, the wavefunction after the entangling gate is

$$\begin{array} { c c c c } \cdots & \circ & \circ \cdots \\ | \psi _ { M } \rangle = \frac { 1 } { \sqrt { 2 ^ { \frac { M N } { 2 } } } } \sum _ { k _ { 1 }, k _ { 3 }, \dots, k _ { M - 1 } = 0 } ^ { N } \left ( \prod _ { i \in \text{odd} } \sqrt { C ^ { k _ { i } } _ { N } } \right ) \\ & \times | k _ { 1 } \rangle _ { 1 } | \frac { e ^ { i ( k _ { 1 } - k _ { 3 } ) t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 } | k _ { 3 } \rangle _ { 3 } | \frac { e ^ { i ( k _ { 3 } - k _ { 5 } ) t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 4 } & \text{which of me} \\ & \cdots \otimes | k _ { M - 1 } \rangle _ { M - 1 } | \frac { e ^ { i k _ { M - 1 } t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M } & ( 1 2 ) \\ \text{for even } M. \text{ We have} \end{array}$$

for even M . We have

$$\text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{$w$} \, \text{\quad} | \psi _ { M } \rangle = \frac { 1 } { \sqrt { 2 \frac { ( M + 1 ) N } { 2 } } } \sum _ { k _ { 1 }, k _ { 3 }, \dots, k _ { M } = 0 } ^ { N } \left ( \prod _ { i \in \text{odd} } \sqrt { C _ { N } ^ { k _ { i } } } \right ) _ { i } \Big | \Psi _ { \bar { q } } \\ \times | k _ { 1 } \rangle _ { 1 } | \frac { e ^ { i ( k _ { 1 } - k _ { 3 } ) t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 } | k _ { 3 } \rangle _ { 3 } | \frac { e ^ { i ( k _ { 3 } - k _ { 5 } ) t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 4 } \\ \cdots \otimes | \frac { e ^ { i ( k _ { M - 2 } - k _ { M } ) t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M - 1 } | k _ { M } \rangle _ { M } \quad ( 1 3 ) \\ \text{for $M$ is odd. In the extensions} \, ( 1 2 ) \, \text{and} \, ( 1 3 ) \, \text{the odd}$$

for M is odd. In the expressions (12) and (13), the odd numbered ensembles are expanded to number states and the remaining ensembles on even sites are conditionally rotated. We note that the wavefunction can be always written in a way such that on adjacent sites, there is a conditional rotation of the ensemble written as a spin coherent state, dependent on the number states on alternating nodes. Hence the parity of the number of ensembles M plays important role. Depending on the parity of M , the last ensemble is written in terms of a spin coherent state or a number state. Henceforth, we will show expressions for even M and defer the expressions for odd M to Appendix A.

Next we perform a measurement of the intermediate ensembles in the x -basis. The measurement operator for this is written

$$M _ { q } ^ { ( j ) } = | q \rangle _ { j } ^ { ( x ) } \langle q | _ { j } ^ { ( x ) }, \quad \quad ( 1 4 ) \quad \text{$\text{wh}$}$$

where the number state in the x -basis is | k 〉 ( x ) = e -iS y π/ 4 | k 〉 and q ∈ [0 , N ] labels the measurement outcome. We note that in conventional entanglement swapping, two qubits from different Bell pairs are measured in the Bell basis to create long-distance entanglement. Such entangled basis measurements are implementable in a photonic context but are less convenient in the context of atomic ensembles. Here, we perform only local measurements on intermediate nodes such as to transfer the entanglement between the first and last ensembles in the chain.

Later, we show that it is beneficial to introduce a slight angular offset from the x -axis of the measurement. For this purpose we perform a equatorial rotation before performing the measurement, given by

$$V ^ { ( j ) } ( \phi ) = e ^ { i n _ { j } ^ { a } \phi }.$$

Hence the final state is the projected state

$$| \Psi _ { \vec { q } } \rangle = \left ( \bigotimes _ { j = 2 } ^ { M - 1 } M _ { q _ { j } } ^ { ( j ) } V ^ { ( j ) } ( \phi ) \right ) | \psi _ { M } \rangle, \quad \ ( 1 6 ) \\ \text{ich is an unnormalized state.  Here we wrote the full set}$$

which is an unnormalized state. Here we wrote the full set of measurement outcomes of the intermediate ensembles as

$$\vec { q } = ( q _ { 2 }, q _ { 3 }, \dots, q _ { M - 1 } ). \quad \quad ( 1 7 )$$

The final unnormalized state after all the measurements is

$$| \Psi _ { \vec { q } } \rangle = \frac { 1 } { \sqrt { 2 ^ { \frac { M N } { 2 } } } } \sum _ { k _ { 1 }, k _ { 3 }, \dots, k _ { M - 1 } = 0 } ^ { N } \sqrt { C _ { N } ^ { k _ { 1 } } } \left ( \prod _ { j = 2 } ^ { M - 1 } \Omega _ { q _ { j } } ^ { ( j ) } \right ) \\ & \times | k _ { 1 } \rangle | \frac { e ^ { i k _ { M - 1 } t } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M } \quad ( 1 8 ) \\, \quad, \, \cdot \, \cdot \,,$$

where we defined

$$\text{odd} \\ \text{all} y \\ \text{ays} \\ \text{is a} \quad \text{The explicit expressions for the matrix elements are}$$

 The explicit expressions for the matrix elements are

$$\underset { \substack { \text{pm} \\ \text{ter-} \\ \text{em-}} } } { \text{pm} } \quad \langle q | ^ { ( x ) } | \frac { e ^ { i \alpha } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle = i ^ { N - q } e ^ { i N \alpha / 2 } \sqrt { C _ { N } ^ { q } } \cos ^ { q } \frac { \alpha } { 2 } \sin ^ { N - q } \frac { \alpha } { 2 } } \\ \underset { \text{th-} } { \text{th-} } }$$

and

$$\text{low} \\ \text{dd} \\ \text{ate} ^ { \left \langle q \right | ^ { ( x ) } | k \right \rangle } = \frac { 1 } { \sqrt { q! ( N - q )! 2 ^ { N } } } \sum _ { l = 0 } ^ { q } \sum _ { m = 0 } ^ { N - q } C _ { q } ^ { l } C _ { N - q } ^ { m } \\ \text{for} \quad ( - 1 ) ^ { N - q - m } \sqrt { ( l + m )! ( N - l - m )! } \delta _ { k, l + m } \quad ( 2 1 ) \\ \text{ where } \delta _ { i, j } \text{ is Kronecker delta } [ 4 8 ]. \text{ The probability of} \\ \text{ obtaining} \sigma \text{ this outcome labeled by} \, \vec { a } \text{ is}$$

where δ i,j is Kronecker delta [48]. The probability of obtaining this outcome labeled by /vector q is

$$p _ { \vec { q } } = \langle \Psi _ { \vec { q } } | \Psi _ { \vec { q } } \rangle.$$

## III. ENTANGLEMENT DISTRIBUTION

The entanglement for the states can be calculated using the von Neumann entropy

$$E = - \sum _ { l = 0 } ^ { N } \lambda _ { l } \log _ { 2 } \lambda _ { l }, \quad \quad ( 2 3 ) \\ \int _ { \mu \nu } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots } \int _ { \dots }$$

where λ l are the eigenvalues of the density matrix

$$\rho _ { M } = \frac { \text{Tr} _ { 1 } | \Psi _ { \vec { q } } \rangle \langle \Psi _ { \vec { q } } | } { p _ { \vec { q } } } = \frac { 1 } { p _ { \vec { q } } } \sum _ { k _ { 1 } = 0 } ^ { N } \langle k _ { 1 } | \Psi _ { \vec { q } } \rangle \langle \Psi _ { \vec { q } } | k _ { 1 } \rangle. \quad ( 2 4 )$$

FIG. 2: Entanglement as measured by von Neumann entropy (23) between the first and last nodes of the chain, after projecting out the intermediate ensembles. The states used are (18) for even M and (A1) for odd M . The measurement outcome q j = N is used for all curves with offset angle φ = 0. The interaction time t is with respect to the Hamiltonian (7) and different lines show various chain lengths M . Each ensemble has N = 20 bosons.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00141v1-figure-ceff1030-5853-432e-aadd-dbb176f2edef.jpg)

The dynamics of entanglement versus time and the number of nodes in chain is shown in Fig. 2. Here we show the case when all measurement results are chosen equal to q j = N and the offset angle φ = 0. We later generalize our analysis to other measurement outcomes, but this case shall be an interesting particular case to examine first. We firstly note that the form of the graph is similar to the characteristic 'devil's crevasse' shape giving sharp dips at fractional multiples of π [43]. The fact that this type of dependence is seen can be understood from the fact that the underlying entanglement between nearest neighbor nodes is generated using S S z z interactions. It naturally follows that the more entanglement that is present between the neighboring nodes, the larger the amount of entanglement that is distributed to the ends of the chain after the projection is applied.

An important question is the performance with scaling the protocol with chain length M , such that long distance entanglement distribution can be achieved. Typically, we observe from Fig. 2 that the entanglement degrades with increasing chain length (see for example short interaction times t &lt; π/ 4). What is interesting is that for some particular times the entanglement does not degrade, and remains constant regardless of the chain length. We call these special points in the graph 'magic times', as these regions are of particular interest in the context of entanglement distribution. It is also interesting to see that as the chain length increases, in the vicinity of the magic times the entanglement curve flattens out and remains constant for a range of times beyond the magic times. This is important as it suggests that there is some tolerance allowable for the interaction times, which is beneficial if some errors exist in controlling t .

To understand the nature of the state that is generated at the magic times, we evaluate the fidelity with respect to the M = 2 state wavefunction given by (11). The fidelity is given by

$$\begin{array} { c } \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \pi \end{array} \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \pi \end{array} \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 } \end{array} \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 \\ \cd.. \\ \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 \\ \cd.. \\ \cdots } \end{array} \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 \\ \cots \\ \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 \\ \cots \\ \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 \\ \cd: \\ \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 \\ \cd. \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cd.. \\ \cdots \\ \cd.. \\ \cdots \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd: \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.: \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd... \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.; \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cdots \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd. \\ \cd.. \\ \cd. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd. \\ \cd.. \\ \cd: \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd. \\ \cd.. \\ \cd.} \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.} \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd. \\ \cd.. \\ \cd... \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd.. \\ \cd..$$

For an odd number of nodes in Fig. 3, we evaluate the fidelity with respect to the M = 3 state wavefunction

$$\Phi _ { \ell } ^ { y } \Big | \Psi _ { N } ^ { ( M = 3 ) } \rangle = \frac { \sum _ { k _ { 1 }, k _ { 3 } = 0 } ^ { N } \sqrt { C _ { N } ^ { k _ { 1 } } C _ { N } ^ { k _ { 3 } } } \Omega _ { N } ^ { ( 2 ) } | k _ { 1 } \rangle | k _ { 3 } \rangle } { \sqrt { \sum _ { k _ { 1 }, k _ { 3 } = 0 } ^ { N } C _ { N } ^ { k _ { 1 } } C _ { N } ^ { k _ { 3 } } | \Omega _ { N } ^ { ( 2 ) } | ^ { 2 } } }. \quad ( 2 6 ) \\ \Phi _ { \ell }. \quad \text{Explicit expressions for the fidelity are given in Appendix} \Lambda \text{ } \text{the odd} \text{ in notation } \text{time to them in} \text{ } \text{thumb in}$$

Explicit expressions for the fidelity are given in Appendix A. The fidelity versus the interaction time is shown in Fig. 3 for various chain lengths. We see that at the magic times the fidelity of the state is equal to 1, for all chain lengths. This shows that the two types of states corresponding to the M = 2 state Eq. (11) the M = 3 state Eq. (26) can be prepared at these magic times without degradation with the length of the chain.

These results suggest that it is possible to distribute the same type of macroscopic entanglement as that generated using the S S z z interaction, as long as the interaction times are in the vicinity of the magic times. However, for this to be considered a practical protocol, we must also analyze outcomes beyond q j = N , as this is only a single measurement outcome from many. For long chain lengths, there is an exponentially small probability of obtaining q j = N , and would be extremely inefficient if only this was the only successful outcome. If one simply considers other measurement outcomes we do not have the stationary entanglement property as seen in Fig. 2, and the particular state that is generated is dependent on the particular measurement outcome. We now show that by modifying the measurement slightly it is possible to produce the same entanglement between the end nodes of the chain for any measurement outcome.

## IV. SPIN-CAT STATES

To understand what is special about the entanglement at the magic times we perform an analysis of the wavefunction at these times. For the purposes of our analysis we define the magic times to be at evolution times

$$t _ { L } = \frac { 2 \pi } { L }$$

where L takes an integer value. We note that there are additional times that have the invariance property with M , e.g. t = 4 π/ 5 (see Fig. 3). For the simplicity and

FIG. 3: The fidelity between the projected state | Ψ /vector q 〉 and the M = 2 state (11) for even cases and the M = 3 state (26) for odd cases. Expressions for the fidelity are given by (25) for even M and (A2) for odd M . The measurement outcome q j = N is used for all curves with offset angle φ = 0. The interaction time t is with respect to the Hamiltonian (7) and different lines show various chain lengths M . Each ensemble has N = 20 bosons.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00141v1-figure-2c7cfd9c-1c58-42d6-af93-781892b3c79a.jpg)

generality of the analysis we shall only consider the times (27). For all times L ∈ [1 , √ N ] the entanglement exhibits the stationary behavior.

The reason why times with fractional multiplies of 2 π are special is discussed in Ref. [43] for the M = 2 case. At these times we see that the the spin coherent states factorize as

$$| \psi _ { 2 } ( t = \frac { 2 \pi } { L } ) \rangle = \frac { 1 } { \sqrt { 2 ^ { N } } } \sum _ { k = 0 } ^ { N } \sqrt { C _ { N } ^ { k } } | k \rangle _ { 1 } | \frac { e ^ { 2 \pi i k / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 }, \quad \begin{array} { c } \sin _ { c } \\ \text{coher} \\ \text{Bloch} \\ \text{for th} \end{array}$$

$$= \frac { 1 } { \sqrt { 2 ^ { N } } } \sum _ { m = 0 } ^ { L - 1 } \left ( \sum _ { \substack { k = 0 \\ \{ k ( \text{mod} L ) = m \} } } ^ { N } \sqrt { C _ { N } ^ { k } } | k \rangle \right ) | \frac { e ^ { 2 \pi i m / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 }. \\ = \quad \text{W}.$$

Here we made the replacement

$$\sum _ { k = 0 } ^ { N } \to \sum _ { m = 0 } ^ { L - 1 } \sum _ { \substack { k = 0 \\ \{ k ( \bmod L ) = m \} } } ^ { N } \quad \ \ ( 3 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where the second summation only runs over those k satisfying k (mod L ) = m .

It is then natural to define the set of states

$$| \mathcal { C } _ { m } \rangle = \frac { 1 } { \sqrt { Z _ { m } } } \sum _ { \substack { k = 0 \\ \{ k ( \bmod L ) = m \} } } ^ { N } \sqrt { C _ { N } ^ { k } } | k \rangle \quad \quad ( 3 1 ) \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{ } \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{ } \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{$$

where the normalization Z m ≈ 2 N /L is a good approximation for N /greatermuch 1. These states form an orthonormal set 〈C m |C n 〉 = δ mn . Taking a Fourier transform of these states gives a spin coherent state

$$\frac { 1 } { \sqrt { L } } \sum _ { m = 0 } ^ { L - 1 } e ^ { - 2 \pi i n m / L } | \mathcal { C } _ { m } \rangle = | \frac { e ^ { - 2 \pi i n / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle. \quad ( 3 2 ) \\ \intertext { t n l e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e f e e }$$

The inverse Fourier transform is

$$| \mathcal { C } _ { m } \rangle = \frac { 1 } { \sqrt { L } } \sum _ { n = 0 } ^ { L - 1 } e ^ { 2 \pi i n m / L } | \frac { e ^ { - 2 \pi i n / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle. \quad ( 3 3 ) \\ \mathfrak { w } _ { - } \dots \Delta \Delta + \alpha \cdots \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \alpha \cdots \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta$$

We see that |C m 〉 can be written as a sum of spin coherent states. This is a spin-cat state which has been considered in various contexts before to encode quantum information [60-63].

The M = 2 state at the magic times can be written

$$\underset { \text{the} \\ \text{the} \\ \text{and} } \quad | \psi _ { 2 } ( t = \frac { \pi } { 2 L } ) \rangle = \frac { 1 } { \sqrt { L } } \sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \frac { e ^ { 2 \pi i m / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 2 }. \quad ( 3 4 ) \\ \text{the} \\ \text{the} \\ \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{end} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad  \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad  \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{or} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \text{and} \quad \\ \text{like}$$

For well-separated spin coherent states

$$\underset { \text{lines} } { \text{lines} } \quad | \langle \langle \frac { e ^ { 2 \pi i n / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } | \frac { e ^ { 2 \pi i m / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle | & = \cos ^ { N } ( n - m ) \pi / L \\ & \approx e ^ { - N ( n - m ) ^ { 2 } \pi ^ { 2 } / 2 L ^ { 2 } } \\ \text{of } 2 \pi \, \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 \, \ }$$

Hence we can consider the state (34) to be an entangled state of L terms where on the first ensemble there is a spin-cat state and on the second ensemble there is a spin coherent state at various positions on the equator of the Bloch sphere. This fact was already noted in Ref. [43] for the case L = 2.

## V. ENTANGLEMENT DISTRIBUTION AT MAGIC TIMES

We now return to the general case of M ensembles and evaluate the quantum state at the magic times. Using the same methods as in the previous section we obtain the wavefunction before measurement for even M as

$$k \text{sat} ^ { \quad } _ { \quad } \colon & \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \col. \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \col: \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colots \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \color \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \col.: \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \col� \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \colon \text{quad} \col:( 3 6 ) \\ & \text{Written in this form we may deduce why there is the en-}$$

Written in this form we may deduce why there is the entanglement invariance property that was observed in Fig. 2. Let us make the measurement of the intermediate ensembles in two steps, first projecting the even numbered

nodes within j ∈ [2 , M -1], then the odd nodes. Splitting the measurements into two steps does not affect the final result as all the measurements commute. Again set the rotation parameter φ = 0 for now. Considering first for simplicity obtaining the outcome

$$| q = N \rangle ^ { ( x ) } = | \frac { 1 } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle \quad \quad ( 3 7 ) \quad I _ { \text{me} }$$

and assume that L &lt; √ N such that the spin coherent states are approximately orthogonal as in (35). Then this enforces m 1 = m 3 = · · · = m M -1 for even M . The resulting state for even M is approximately

$$\bigotimes _ { j = 2, 4, \dots M - 2 } | q = N \rangle _ { j } ^ { ( x ) } \langle q = N | _ { j } ^ { ( x ) } | \psi _ { M } \rangle \, \infty \quad \begin{array} { c } \underset { \text{com} \\ \text{unil} \\ \text{easil} \end{array} } \\ \sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \mathcal { C } _ { m } \rangle _ { 3 } \dots | \mathcal { C } _ { m } \rangle _ { M - 1 } | \frac { e ^ { 2 i m \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M }, \quad ( 3 8 ) \quad \text{(38)} \quad \text{an a}.$$

which is valid for N /greatermuch 1. This correlates all the remaining ensembles, in particular the first and last ensemble such that they have the same label m .

$$\sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \mathcal { C } _ { m } \rangle _ { 3 } \dots | \mathcal { C } _ { m } \rangle _ { M - 1 } | \frac { e ^ { 2 i m \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M }, \quad ( 3 8 ) \\ \text{which is valid for $N\gg 1$. This correlates all the remain-}$$

Now looking at the remaining intermediate ensembles, the spin-cat states can be written as a superposition of spin coherent states as shown in (33). Measurement of the cat states will collapse the superposition as

$$\langle q = N | ^ { ( x ) } | \mathcal { C } _ { m } \rangle \approx \frac { 1 } { \sqrt { L } } \, \quad \quad ( 3 9 ) \, \quad \underset { \text{jec} } { \text{jec} }$$

where only the n = 0 term in (33) contributes for wellseparated spin coherent states. The correlations between the first and last ensembles are hence preserved and we obtain

$$\underset { j = 2 } { \overset { M - 1 } { \bigotimes } } | q = N \rangle _ { j } ^ { ( x ) } \langle q = N | _ { j } ^ { ( x ) } | \psi _ { M } \rangle \underset { \pi } { \infty } \\ \sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \frac { e ^ { 2 i m \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M }, \quad ( 4 0 ) \quad \text{whe} \\ \mathfrak { h } \text{ is valid for } N \gg 1. \ \text{ This is precisely the same} \\ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots$$

which is valid for N /greatermuch 1. This is precisely the same quantum state as (34). This is obtained for arbitrary chain length M . This explains the invariance property of the entanglement and fidelity observed in Figs. 2 and 3 for the even case.

For odd M , we can use the same manipulations to obtain

$$\stackrel { \sim } { \bigotimes _ { j = 2 } ^ { M - 1 } } | q = N \rangle _ { j } ^ { ( x ) } \langle q = N | _ { j } ^ { ( x ) } | \psi _ { M } \rangle \, \infty \, \sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \mathcal { C } _ { m } \rangle _ { M }. \\ ( 4 1 ) \quad \text{Figur} \\. \quad. \. \. \.$$

This is the same state as (26) which can be written as

$$| \Psi _ { N } ^ { ( M = 3 ) } \rangle & \approx \frac { 1 } { \sqrt { L } } \sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \mathcal { C } _ { m } \rangle _ { 3 }, & \quad ( 4 2 ) \quad \text{respec} \\ 1 & \stackrel { L - 1 } { = } \ e ^ { 2 i n \pi / L } \quad _ { e ^ { 2 i n \pi / L } } 1 \quad _ { e ^ { - 2 i n \pi / L } } 1 \quad \text{oscilla}$$

$$= \frac { 1 } { \sqrt { L } } \sum _ { n = 0 } ^ { \stackrel { \prime \prime - \upsilon } { L - 1 } } | \frac { e ^ { 2 i n \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 1 } | \frac { e ^ { - 2 i n \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { 3 } \quad \text{$\sin} \\ ( 4 3 ) \quad \text{where}$$

## VI. MEASUREMENTS OF SPIN COHERENT STATES

In the previous section we showed that for the measurement outcome q j = N at the magic times, an entangled state consisting of spin-cat states is generated. We now wish to generalize this result for arbitrary measurement outcomes. One can see from the form of (36) that if one were able to measure in the spin coherent state basis, a similar result would be obtained for all measurement outcomes. This is however experimentally challenging, since unlike the optical case, homodyne measurements are not easily realizable with spin ensemble systems. Therefore, an alternative is required to implement the measurements of the last section. In this section we show that standard spin measurements are sufficient in achieving the same collapse as described in the previous section.

To illustrate the measurement scheme, first consider performing a measurement of the spin-cat state (33) in the x basis. Plotting the positions of the spin coherent states involved in this superposition immediately identifies a problem: spin coherent states labeled with n and L -n have the same x position (Fig. 4(a)). Hence projection in the x basis will collapse the superposition, but will catch two spin coherent states together, unless they are located at the extremal positions 〈 S x 〉 = ± N . To avoid this issue, it is beneficial to slightly offset the spin coherent states such that one measures instead the state

$$e ^ { i \phi n ^ { a } } | \mathcal { C } _ { m } \rangle = \frac { 1 } { \sqrt { L } } \sum _ { n = 0 } ^ { L - 1 } e ^ { 2 \pi i n m / L } | \frac { e ^ { - i ( 2 \pi n / L - \phi ) } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle,$$

where a suitable angle choice is

$$\phi = \frac { \pi } { 2 L }.$$

This slight rotation offsets the spin coherent states slightly such that a projection on the x -axis results in a collapse of the state to uniquely one of the spin coherent states (Fig. 4(b)).

The effect of the measurement is also evident by examining the probability of the projection

$$p _ { q } = | \langle q | ^ { ( x ) } e ^ { i \phi n ^ { a } } | \mathcal { C } _ { m } \rangle | ^ { 2 }. \quad \quad ( 4 6 )$$

Figure 4(c)(d) shows the probability of the rotated spincat state for the same states as shown in Figs. 4(a)(b) respectively. We see peaks in the probability distribution at the expectation value of S x each of the spin coherent states involved in the spin-cat state. When the projection involves two spin coherent states, there is a rapid oscillation of the probability distribution (Fig. 4(c)). By applying the rotation φ this ambiguity is removed and smooth Gaussian peaks are seen for each of the spin coherent states (Fig. 4(d)).

FIG. 4: Visualization and measurement of a spin-cat state. (a)(b) Visualization of the spin-cat states (33) and (44) respectively. Figures are plotted by centering a Gaussian e -N [( 〈 S x 〉-cos θ ) 2 +( 〈 S y 〉-sin θ ) 2 ] / 4 , which is the Q-function distribution for a spin coherent state. (c)(d) The probability distribution (46) for φ = 0 , π/ 2 L respectively. For all plots L = 4 , N = 100 is used.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00141v1-figure-58f0bc20-de9b-4221-af29-da8b9d06b7ca.jpg)

For a given measurement outcome q , one can deduce which spin coherent state has been projected using the peak position of the Gaussians:

$$q _ { \text{peak} } = N \cos ^ { 2 } ( \pi ( \bar { n } - 1 / 4 ) / L ), \quad \ \ ( 4 7 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where the spin coherent state is labeled by ¯ n ∈ [0 , L -1]. For a particular outcome q , by finding the closest q peak one can deduce which spin coherent state labeled by ¯ n the collapse occurred at.

## VII. APPROXIMATE WAVEFUNCTION FOR ARBITRARY MEASUREMENT OUTCOMES

Using the technique in the previous section, we may now examine the remaining projection outcomes in addition to the q j = N case. Examining (36), first consider performing a measurement on the even numbered ensembles. Assuming that the collapse occurs such as to pick out one and only one of the spin coherent states we obtain the relation

$$m _ { j - 1 } - m _ { j + 1 } = \bar { m } _ { j } ( \text{mod} L ). \quad \ \ ( 4 8 ) \quad \ \ |,$$

The state is then

$$& \text{The state is then} \\ & \bigotimes _ { j = 2, 4, \dots M - 2 } | q _ { j } \rangle _ { j } ^ { ( x ) } \langle q _ { j } | _ { j } ^ { ( x ) } | \psi _ { M } \rangle \, \infty \quad \sum _ { m = } ^ { L - } \\ & \sum _ { m = 0 } ^ { L - 1 } | \mathcal { C } _ { m } \rangle _ { 1 } | \mathcal { C } _ { m - \bar { m } _ { 2 } } \rangle _ { 3 } | \mathcal { C } _ { m - \bar { m } _ { 2 } - \bar { m } _ { 4 } } \rangle _ { 5 } \dots | \mathcal { C } _ { m - \bar { m } _ { e v e n } } \rangle _ { M - 1 } \quad \text{Fig} \\ & \otimes | \frac { e ^ { 2 i ( m + \bar { m } _ { e v e n } ) \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M }. \quad \text{(49)} \quad \text{sis is}$$

where we defined

$$\bar { m } _ { \text{even} } = \sum _ { j = 2, 4, \dots, M - 2 } \bar { m } _ { j }. \quad \quad ( 5 0 ) \\. \quad \dots \quad \dots \quad \dots \quad \dots$$

This correlates the first and last ensemble up to an offset in m . In (49), the proportionality factor is

$$\begin{array} { c c c } \int m. \, \Pi \, ( \mathcal { O } ), \, \text{the po/policency trace} \, \text{ is} \\ \int _ { 1 } & \frac { 1 } { \sqrt { L ^ { M / 2 } } } \prod _ { j = 2, 4, \dots, M - 2 } \langle q _ { j } | ^ { ( x ) } | \frac { e ^ { 2 i ( m _ { j - 1 } - m _ { j + 1 } ) \pi / L + i \phi } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle \\ \Big \rfloor & = \frac { e ^ { i N \tilde { m } _ { \text{even} } \pi / L + i N \phi / 2 } } { \sqrt { L ^ { M / 2 } } } \\ & \times \prod _ { j = 2 } ^ { M - 2 } \sqrt { C _ { N } ^ { q _ { j } } } \cos ^ { q _ { j } } ( \frac { \tilde { m } _ { j } \pi } { L } + \frac { \phi } { 2 } ) \sin ^ { N - q _ { j } } ( \frac { \tilde { m } _ { j } \pi } { L } + \frac { \phi } { 2 } ) } \, \quad ( 5 1 ) \\ & \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array} \, \Big \rceil \,.$$

$$& \times \prod _ { j = 2 } ^ { M - 2 } \sqrt { C _ { N } ^ { q _ { j } } } \cos ^ { q _ { j } } ( \frac { \bar { m } _ { j } \pi } { L } + \frac { \phi } { 2 } ) \sin ^ { N - q _ { j } } ( \frac { \bar { m } _ { j } \pi } { L } + \frac { \phi } { 2 } ) \ ( 5 1 ) \\ & \dots \. \. \. \. \. \.$$

which has no dependence on m hence does not affect the wavefunction.

Measurement of the odd numbered intermediate ensembles will give factors

$$\stackrel { \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \lor \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \OC } \text{$\dagger$} \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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\text{$\text{$\test$\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{$\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{$\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{$$$$

where in the second line we used the result of the previous section that the measurement will collapse the spin-cat state to one of the spin coherent states with label ¯. n Here the important factor that affects the wavefunction is the phase dependence on m . The total m dependence is given by projecting the odd sites is

$$\text{OR} & \quad \langle q _ { 3 } | ^ { ( x ) } | \mathcal { C } _ { m - \bar { m } _ { 2 } } \rangle \langle q _ { 5 } | ^ { ( x ) } | \mathcal { C } _ { m - \bar { m } _ { 2 } - \bar { m } _ { 4 } } \rangle \dots \langle q _ { M - 1 } | ^ { ( x ) } | \mathcal { C } _ { m - \bar { m } _ { e v e n } } \rangle \\ \text{ES} & \quad \in e ^ { 2 \pi i m \bar { n } _ { o d d } / L } \\ \text{$e$ max}$$

where

$$\bar { n } _ { \text{odd} } = \sum _ { j = 3, 5, \dots, M - 1 } \bar { n } _ { j }. \quad \quad ( 5 4 ) \\ \cdot \. \. \. \. \. \. \. \. \. \. \.$$

We thus finally obtain the final state after projection

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ | \Psi _ { \vec { q } } ^ { \text{approx} } \rangle = \bigotimes _ { j = 2 } ^ { M - 1 } | q _ { j } \rangle _ { j } ^ { ( x ) } \langle q _ { j } | _ { j } ^ { ( x ) } | \psi _ { M } \rangle \infty \\ \sum _ { m = 0 } ^ { L - 1 } e ^ { 2 \pi i m \bar { n } _ { \text{odd} } / L } | \mathcal { C } _ { m } \rangle _ { 1 } | \frac { e ^ { 2 i ( m - \bar { m } _ { \text{even} } ) \pi / L } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle _ { M }. \quad ( 5 5 ) \\. \quad \text{Figure 5 show a comparison of the approximate wave-}$$

Figure 5 show a comparison of the approximate wavefunction (55) with the exact wavefunction (18). In Fig. 5(a) the amplitude of the coefficients in the number basis is shown. Excellent agreement is seen, with the coefficents having nearly perfect agreement. The fidelity of

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00141v1-figure-f3b7f340-438d-4192-bb51-80eedef7cf06.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00141v1-figure-0f13f17e-6770-4ff1-badd-b7bed222accf.jpg)

F

FIG. 5: Comparison of the exact wavefunction (18) with the approximate wavefunction (55). (a) The wavefunction amplitude |〈 Ψ /vector q | k , k 1 M 〉| for N = 10 , M = 3 , L = 3 and q 2 = N/ 2 for K = ( N + 1) k 1 + k M + 1. The approximate wavefunction is shown as the dashed line, the exact wavefunction is shown as a solid line. (b) The fidelity F = |〈 Ψ Ψ /vector q | approx /vector q 〉| 2 for M = 3 , L = 3 for three outcomes q 2 = 0 , N/ 2 , N as marked. All cases implement the measurement scheme where the intermediate ensembles are measured using the offset angle φ = π/ L 2 .

the approximate wavefunction improves with N , as seen in Fig. 5(b). The reason that the fidelities improve with N is due to the measurement scheme that we implement as described in the previous section. In order to collapse the state to a single spin coherent state, the Gaussians for each one should be well-separated, as discussed in Fig. 4. This requires L &lt; √ N typically. With increasing N , the Gaussian become better separated and approaches the ideal collapse to a single spin coherent state.

We note that in Fig. 3 perfect fidelity F = 1 was obtained, whereas in Fig. 5 only near-unit fidelity is obtained. The reason is that in Fig. 3 the offset angle φ = 0 and only the outcome q j = N was considered. For this case, the projection on the spin coherent state is perfect, as can be seen in Fig. 4(a) and 4(c). However, there is no guarantee that the outcome q j = N is obtained, in fact for long chains it has an exponentially small probability. Using the offset measurement of Fig. 4(b)(d) is advantageous in this sense, as all measurement outcomes lead to the generation of a nearly ideal state (55). The small price to pay for this is a non-unit fidelity, although as can be seen in Fig. 5(b) the fidelity quickly approaches 1.

Therefore, using the offset measurement include φ = π/ L 2 , one may obtain consistently obtain (55) for all measurement outcomes. The final step to turn this into a deterministic scheme for preparing entanglement is to peform conditional rotations to eliminate the ¯ n odd and ¯ m even terms. Specifically, knowledge of ¯ n odd and ¯ m even allows for a correction of the state to the state (41) using the rotations

$$e ^ { i 2 \pi n ^ { \alpha } / L } | \mathcal { C } _ { m } \rangle = e ^ { i 2 \pi m / L } | \mathcal { C } _ { m } \rangle \quad \quad ( 5 6 ) \quad \text{ when} \\ _ { i \alpha ^ { \prime } } \quad \tilde { \epsilon } \quad _ { i ( \alpha ) + \alpha ^ { \prime } ) } \quad \tilde { \epsilon } \quad \text{ and}$$

$$e ^ { i \phi n ^ { a } } | \frac { e ^ { i \phi ^ { \prime } } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle = | \frac { e ^ { i ( \phi + \phi ^ { \prime } ) } } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle. \quad \ \ ( 5 7 ) ^ { \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Then by applying the corrections

$$e ^ { - 2 \pi i \bar { n } _ { \text{odd} } n _ { 1 } ^ { a } / L } e ^ { 2 \pi i \bar { m } _ { \text{even} } n _ { M } ^ { a } / L } | \Psi _ { \vec { q } } \rangle \approx | \psi _ { 2 } ( t = \frac { \pi } { 2 L } ) \rangle \quad ( 5 8 ) \quad \text{We}$$

we obtain a deterministic preparation of the state (34) which is the aim of this study.

## VIII. RESOURCE SCALING AND EXPERIMENTAL IMPLEMENTATION

We now comment on the resource costs required to execute the protocol. This question is intimately connected to the experimental implementation, since the resources required depend on the way that the operations are implemented. The key aspect that makes this scheme efficient is that the protocol has no dependence on N , the number of qubits in the ensembles, in any of the operations that need to be performed. Despite this, the amount of entanglement that is distributed is of order log 2 N , as can be seen in Fig. 2. The reason for this is that only collective spin operations are used to control the ensembles in all the steps of the protocol.

To make this point clear, let us review the steps of the protocol and examine them each:

$$\i t _ { \quad s e } \quad 1. \text{ Prepare the initial state } | \phi _ { 0 } \rangle = | \frac { 1 } { \sqrt { 2 } }, \frac { 1 } { \sqrt { 2 } } \rangle \rangle ^ { \otimes M }.$$

- 2. Apply the Hamiltonian (7) for a time t .
- 3. Apply the rotation (15) to the intermediate qubits.
- 4. Measure the intermediate qubits in the x -basis (14).
- 5. Based on the measurement results, apply the correction (58).

Specifically, the initial state preparation in Step 1 can be routinely performed in atomic ensembles using optical pumping to populate one of the atomic levels, followed by a π/ 2 pulse implemented with radio frequency (RF) pulses. Both of these operations are independent of N since the optical and RF pulses applied to the ensembles are common to all atoms in the ensemble. The unitary operations in Steps 2, 3, 5 are all in terms of total spin operators S γ , hence are at the ensemble level. Again these can be implemented using RF pulses on the hyperfine ground states of the atoms. The S S z z interaction of Step 2 can be implemented using techniques discussed in Ref. [44, 47], where the interaction is mediated using an optical channel. Finally, the measurement in Step 4 is performed using atomic counting of the ensemble, which is routinely performed in works such as Refs. [64-66]. All these operations only require operations at the ensemble level, not at the microscopic level of the individual qubits. For this reason there is no N dependence to the scheme. The only scaling is with respect to the chain length M , where the M -1 links of the chain must be entangled and M -2 intermediate nodes must be projected out.

## IX. CONCLUSIONS

We have proposed a protocol to distribute macroscopic entanglement on ensembles of qubits, by first generating

entanglement in a chain configuration then projecting out the intermediate ensembles to leave the ends of the chain entangled. Such a protocol is non-trivial in the ensemble case since the entanglement is of a complex many-body nature, which causes the entanglement tends to diminish with longer chain lengths. We have found 'magic' interaction times where the entanglement can be transferred without degradation. The nature of the states at these magic times was found due to an encoding of the entangled states in terms of generalized spin-cat states. This allowed us to derive the explicit form of the states for any ensemble size. The final wavefunction (18), which has the approximate form (55), takes a similar form to the M = 2 ensemble case, up to simple spin rotations. By correcting for these known rotations, it allows for the deterministic preparation of the equivalent state generated by an S S z z interaction.

The amount of entanglement that is generated using this scheme is of order log 2 N , and hence is macroscopic by our definition. The resources required to distribute this do not involve N , and only depend upon M , the chain length. The reason is that only total spin operations are used in the protocol from start to finish. In this sense this has a considerable scaling advantage over preparing the entanglement microscopically. The parallelization occurs since the resources to control an ensemble are the same as that of a qubit. The use of total spin operators means that the control pulses illuminate the entire ensemble simultaneously, which is the origin of the parallelization. The type of entanglement that is distributed is suitable for various protocols that are based on the S S z z interaction, since this is ultimately the form of entanglement that is distributed. This includes schemes such as quantum teleportation [49, 67], Deutsch-Jozsa algorithm [61]. Such tasks could be implemented on a distributed quantum computing platform. The distribution of entanglement at macroscopic quantities is a task that we anticipate will become more important as the quantum internet gains more prominence. Our protocol offers the possibility of distribute entanglement at scale, which can serve as a fundamental primitive in quantum network infrastructure.

## X. ACKNOWLEDGMENTS

A. N. P and I. D. L. are supported by the Russian Science Foundation (grant no 23-21-00507). T.B. is supported by the National Natural Science Foundation of China (62071301); NYU-ECNU Institute of Physics at NYU Shanghai; Shanghai Frontiers Science Center of Artificial Intelligence and Deep Learning; the Joint Physics Research Institute Challenge Grant; the Science and Technology Commission of Shanghai Municipality (19XD1423000,22ZR1444600); the NYU Shanghai Boost Fund; the China Foreign Experts Program (G2021013002L); the NYU Shanghai Major-Grants Seed Fund; Tamkeen under the NYU Abu Dhabi Research In- stitute grant CG008; and the SMEC Scientific Research Innovation Project (2023ZKZD55).

## Appendix A: Expressions for odd M

Here we show the corresponding expressions for odd M for the even M expressions given in the main text.

The odd M counterpart of (18) is

$$\text{ns} \\ \text{for} \\ \text{ich} \\ \text{to} \\ \text{By} \\ \text{de} \\ \text{ted} \\ \text{~} \\ \text{ip} \\ \text{pic} \\ \text{in order to calculate the fidelity for odd case, we take} \\ \text{$\bullet$} \text{ and} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$} \text{$\bullet$}$$

In order to calculate the fidelity for odd case, we take the overlap with the projected M = 3 state (26). The fidelity in this case equals

$$\intertext { \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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test{$\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{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\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{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$$

when we measure N on all intermediate ensembles.

The odd M counterpart of (36) is

$$\text{seq on} \\ \text{orm of} \\ \text{hemes} \\ \text{-Jozsa} \\ \text{1 on a} \\ \text{strib-a task} \\ \text{at task} \\ \text{-octol} \\ ; \text{scale,} \\ \text{scale,} \\ \text{antum}$$

Following similar steps to that in Sec. VII we get the odd M counterpart of (55)

$$\begin{array} { c c c } \cdot & \cdot & \cdot \\ | \Psi _ { \vec { q } } ^ { \text{approx} } \rangle = \bigotimes _ { j = 2 } ^ { M - 1 } | q _ { j } \rangle _ { j } ^ { ( x ) } \langle q _ { j } | _ { j } ^ { ( x ) } | \psi _ { M } \rangle \in \\ \sum _ { m = 0 } ^ { L - 1 } e ^ { 2 \pi i m \bar { n } _ { \text{odd} } / L } | \mathcal { C } _ { m } \rangle _ { 1 } | \mathcal { C } _ { m - \bar { n } _ { \text{even} } } \rangle _ { M }. \quad & ( A 4 ) \\ \end{array}$$

where for (A4) we use the modified definitions

$$\bar { m } _ { \text{even} } & = \sum _ { j = 2, 4, \dots, M - 1 } \bar { m } _ { j } \\ \bar { n } _ { \text{odd} } & = \sum _ { j = 3, 5, \dots, M - 2 } \bar { n } _ { j }.$$

- [1] S.-H. Wei, B. Jing, X.-Y. Zhang, J.-Y. Liao, C.-Z. Yuan, B.-Y. Fan, C. Lyu, D.-L. Zhou, Y. Wang, G.-W. Deng, et al., Laser &amp; Photonics Reviews 16 , 2100219 (2022).
- [2] S. Wehner, D. Elkouss, and R. Hanson, Science 362 (2018).
- [3] K. Azuma, S. E. Economou, D. Elkouss, P. Hilaire, L. Jiang, H.-K. Lo, and I. Tzitrin, Reviews of Modern Physics 95 , 045006 (2023).
- [4] S. Pirandola, R. Laurenza, C. Ottaviani, and L. Banchi, Nat. Commun. 8 , 15043 (2017).
- [5] Y. Mao, P. Zeng, and T.-Y. Chen, Advanced Quantum Technologies 4 , 2000084 (2021).
- [6] J. L. Park, Found. Phys. 1 , 23 (1970).
- [7] W. K. Wootters and W. H. Zurek, Nature 299 , 802 (1982).
- [8] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University Press, 2010).
- [9] H.-J. Briegel, W. D¨ ur, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 81 , 5932 (1998).
- [10] C. H. Bennett, G. Brassard, C. Cr´ epeau, R. Jozsa, A. Peres, and W. K. Wootters, Physical review letters 70 , 1895 (1993).
- [11] M. Zukowski, A. Zeilinger, M. A. Horne, and A. K. Ekert, Physical review letters 71 , 4287 (1993).
- [12] J.-W. Pan, D. Bouwmeester, H. Weinfurter, and A. Zeilinger, Phys. Rev. Lett. 80 , 3891 (1998).
- [13] A. I. Lvovsky, B. C. Sanders, and W. Tittel, Nat. Photonics 3 , 706 (2009).
- [14] C. Simon, M. Afzelius, J. Appel, A. Boyer de la Giroday, S. J. Dewhurst, N. Gisin, C. Y. Hu, F. Jelezko, S. Kr¨ll, o J. H. M¨ller, et al., Eur. Phys. J. D u 58 , 1 (2010).
- [15] K. Heshami, D. England, P. Humphreys, P. Bustard, V. Acosta, J. Nunn, and B. Sussman, Journal of Modern Optics 63 , 2005 (2016).
- [16] A. Reiserer and G. Rempe, Rev. Mod. Phys. 87 , 1379 (2015).
- [17] W. Rosenfeld, F. Hocke, F. Henkel, M. Krug, J. Volz, M. Weber, and H. Weinfurter, Phys. Rev. Lett. 101 , 260403 (2008).
- [18] N. Sangouard, C. Simon, H. de Riedmatten, and N. Gisin, Rev. Mod. Phys. 83 , 33 (2011).
- [19] L.-M. Duan and C. Monroe, Rev. Mod. Phys. 82 , 1209 (2010).
- [20] L. Childress, M. V. G. Dutt, J. M. Taylor, A. S. Zibrov, F. Jelezko, J. Wrachtrup, P. R. Hemmer, and M. D. Lukin, Science 314 , 281 (2006).
- [21] M. V. G. Dutt, L. Childress, L. Jiang, E. Togan, J. Maze, F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, Science 316 , 1312 (2007).
- [22] S. Hermans, M. Pompili, H. Beukers, S. Baier, J. Borregaard, and R. Hanson, Nature 605 , 663 (2022).
- [23] M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. Mod. Phys. 77 , 633 (2005).
- [24] L. M. Duan, M. D. Lukin, J. I. Cirac, and P. Zoller, Nature 414 , 413 (2001).
- [25] S. Ritter, C. N¨lleke, C. Hahn, A. Reiserer, A. Neuzner, o M. Uphoff, M. M¨cke, E. Figueroa, J. Bochmann, and u G. Rempe, Nature 484 , 195 (2012).
- [26] A. Kuzmich, W. P. Bowen, A. D. Boozer, A. Boca, C. W. Chou, L.-M. Duan, and H. J. Kimble, Nature 423 , 731

(2003).

- [27] C. H. van der Wal, M. D. Eisaman, A. Andr´ e, R. L. Walsworth, D. F. Phillips, A. S. Zibrov, and M. D. Lukin, Science 301 , 196 (2003).
- [28] C. W. Chou, H. de Riedmatten, D. Felinto, S. V. Polyakov, S. J. van Enk, and H. J. Kimble, Nature 438 , 828 (2005).
- [29] Z.-S. Yuan, Y.-A. Chen, B. Zhao, S. Chen, J. Schmiedmayer, and J.-W. Pan, Nature 454 , 1098 (2008).
- [30] Y. Yu, F. Ma, X.-Y. Luo, B. Jing, P.-F. Sun, R.-Z. Fang, C.-W. Yang, H. Liu, M.-Y. Zheng, X.-P. Xie, et al., Nature 578 , 240 (2020).
- [31] Y.-F. Pu, S. Zhang, Y.-K. Wu, N. Jiang, W. Chang, C. Li, and L.-M. Duan, Nat. Photonics 15 , 374 (2021).
- [32] J. Simon, H. Tanji, J. K. Thompson, and V. Vuleti´ c, Phys. Rev. Lett. 98 , 183601 (2007).
- [33] X.-H. Bao, A. Reingruber, P. Dietrich, J. Rui, A. D¨ck, u T. Strassel, L. Li, N.-L. Liu, B. Zhao, and J.-W. Pan, Nat. Phys. 8 , 517 (2012).
- [34] N. Sangouard, C. Simon, B. Zhao, Y.-A. Chen, H. de Riedmatten, J.-W. Pan, and N. Gisin, Phys. Rev. A 77 , 062301 (2008).
- [35] J.-L. Liu, X.-Y. Luo, Y. Yu, C.-Y. Wang, B. Wang, Y. Hu, J. Li, M.-Y. Zheng, B. Yao, Z. Yan, et al., Nature 629 , 579 (2024).
- [36] B. Julsgaard, A. Kozhekin, and E. S. Polzik, Nature 413 , 400 (2001).
- [37] H. Krauter, D. Salart, C. Muschik, J. M. Petersen, H. Shen, T. Fernholz, and E. S. Polzik, Nature Physics 9 , 400 (2013).
- [38] J. Dias and T. C. Ralph, Physical Review A 97 , 032335 (2018).
- [39] J. Dias, M. S. Winnel, N. Hosseinidehaj, and T. C. Ralph, Physical Review A 102 , 052425 (2020).
- [40] M. Ghalaii and S. Pirandola, Physical Review A 102 , 062412 (2020).
- [41] F. Furrer and W. J. Munro, Physical Review A 98 , 032335 (2018).
- [42] K. P. Seshadreesan, H. Krovi, and S. Guha, Physical Review Research 2 , 013310 (2020).
- [43] T. Byrnes, Phys. Rev. A 88 , 023609 (2013).
- [44] A. N. Pyrkov and T. Byrnes, New Journal of Physics 15 , 093019 (2013).
- [45] D. Rosseau, Q. Ha, and T. Byrnes, Phys. Rev. A 90 , 052315 (2014).
- [46] T. Byrnes, K. Wen, and Y. Yamamoto, Phys. Rev. A 85 , 040306(R) (2012).
- [47] M. I. Hussain, E. O. Ilo-Okeke, and T. Byrnes, Physical Review A 89 , 053607 (2014).
- [48] T. Byrnes and E. O. Ilo-Okeke, Quantum atom optics: Theory and applications to quantum technology (Cambridge university press, 2021).
- [49] A. N. Pyrkov and T. Byrnes, Phys. Rev. A 90 , 062336 (2014).
- [50] T. Byrnes, D. Rosseau, M. Khosla, A. Pyrkov, A. Thomasen, T. Mukai, S. Koyama, A. Abdelrahman, and E. Ilo-Okeke, Optics Communications 337 , 102 (2015).
- [51] M. Chaudhary, M. Fadel, E. O. Ilo-Okeke, A. N. Pyrkov, V. Ivannikov, and T. Byrnes, Phys. Rev. A 103 , 062417 (2021).

- [52] T. Byrnes, K. Yan, and Y. Yamamoto, New Journal of Physics 13 , 113025 (2011).
- [53] R. Schmied, J.-D. Bancal, B. Allard, M. Fadel, V. Scarani, P. Treutlein, and N. Sangouard, Science 352 , 441 (2016).
- [54] M. Fadel, T. Zibold, B. D´camps, and P. Treutlein, Scie ence 360 , 409 (2018).
- [55] N. J. Engelsen, R. Krishnakumar, O. Hosten, and M. A. Kasevich, Phys. Rev. Lett. 118 , 140401 (2017).
- [56] P. Colciaghi, Y. Li, P. Treutlein, and T. Zibold, Physical Review X 13 , 021031 (2023).
- [57] T. Byrnes, arXiv preprint arXiv:2307.00875 (2023).
- [58] P. Treutlein, T. W. H¨nsch, J. Reichel, A. Negretti, M. A. a Cirone, and T. Calarco, Phys. Rev. A 74 , 022312 (2006).
- [59] Y. Jing, M. Fadel, V. Ivannikov, and T. Byrnes, New Journal of Physics 21 , 093038 (2019).
- [60] G. S. Agarwal, R. R. Puri, and R. P. Singh, Physical Review A 56 , 2249 (1997).
- [61] H. Semenenko and T. Byrnes, Physical Review A 93 , 052302 (2016).
- [62] W. Qin, A. Miranowicz, H. Jing, and F. Nori, Physical Review Letters 127 , 093602 (2021).
- [63] S. Omanakuttan, V. Buchemmavari, J. A. Gross, I. H. Deutsch, and M. Marvian, arXiv preprint arXiv:2401.04271 (2024).
- [64] D. B. Hume, I. Stroescu, M. Joos, W. Muessel, H. Strobel, and M. K. Oberthaler, Physical Review Letters 111 , 253001 (2013).
- [65] H. Ott, Reports on Progress in Physics 79 , 054401 (2016).
- [66] A. Qu, B. Evrard, J. Dalibard, and F. Gerbier, Physical Review Letters 125 , 033401 (2020).
- [67] A. N. Pyrkov and T. Byrnes, New Journal of Physics 16 , 073038 (2014).
Alexey N. Pyrkov Ilia D. Lazarev Tim Byrnes
Quantum repeater protocol for deterministic distribution of macroscopic entanglement
Distributing long-distance entanglement is a fundamental goal that is
necessary for a variety of tasks such as quantum communication, distributed
quantum computing, and quantum metrology. Currently quantum repeater schemes
typically aim to distribute one ebit at a time, the equivalent of one Bell
pair's worth of entanglement. Here we present a method to distribute a
macroscopic amount of entanglement across long-distances using a number of
operations that scales only linearly with the chain length. The scheme involves
ensembles of qubits and entangling them with an $S^z S^z$ interaction, which
can be realized using atomic gas ensembles coupled by a shared optical mode.
Using only local measurements on the intermediate ensembles, this leaves the
ensembles at the ends of the chain entangled. We show that there are particular
``magic'' interaction times that allow for distribution of entanglement with
perfect fidelity, with no degradation with chain length. The scheme is
deterministic, such that with suitable local conditional unitary corrections,
the same entangled state can always be prepared with good approximation.