We simulate entanglement sharing between two end-nodes of a linear chain quantum network using SeQUeNCe, an open-source simulation package for quantum networks. Our focus is on the rate of entanglement generation between the end-nodes with many repeaters with a finite quantum memory lifetime. Numerical and analytical simulations show limits of connection performance for a given number of repeaters involved, memory lifetimes, the distance between the end-nodes, and an entanglement management protocol. Our findings demonstrate that the performance of quantum connection depends highly on the entanglement management protocol, which schedules entanglement generation and swapping, resulting in the final end-to-end entanglement.

## I. INTRODUCTION

In recent years, much research effort has focused on developing quantum communication networks and sharing quantum entanglement among spatially separated qubits.^{1–10} Small linear quantum networks (with maximum distance in a pair of nodes up to about 100 km) have already been demonstrated.^{11,12} Theoretical development of network protocols has reached its third generation. The generations are classified by methods adopted to suppress loss and operation errors by either heralded entanglement generation or quantum error correction.^{13} The first two generations are based on heralded entanglement generation, requiring signaling back to communicating nodes about the generation status. The third generation of the quantum network is free from this obligation that significantly limits the network throughput. However, the third generation networks make severe demands on the fidelity of quantum gates that are now incompatible with today's hardware characteristics.^{13,14}

Currently, available quantum hardware components are still far from allowing the realization of a fully functional first-generation network for long-distance end-to-end quantum communication.^{15} Key bottlenecks include limited coherence time of qubits and photon loss in the quantum channel medium between the end-nodes. While the latter issue can, in principle, be overcome by using quantum repeaters,^{2,3,16–20} including multimode solutions with multiplexing;^{21–27} the former has been a persistent impediment to progress in quantum communication in particular and quantum information technologies in general. The multi-path routing,^{28,29} involving multiple paths for routing entanglement between a pair of end-users, can enable a long-distance entanglement generation rate with a higher rate than what is possible with a linear repeater chain.^{30}

Quantum memories are essential ingredients for a quantum repeater.^{31,32} Their coherent lifetimes are crucial for creating and maintaining high fidelity entanglement^{33–35} and are dependent on the materials platform. While typical spin qubits coherence time of a few milliseconds up to a second in silicon^{36,37} and diamond^{38–40} is comparable to an average ping time of about 0.1 ms in classical networks, the trapped ions demonstrate memory lifetimes from several minutes to hours.^{41–43} However, the frequency conversion efficiency to telecommunication wavelength (∼1560 *μ*m) remains low.^{44–46} Currently, there is intense research for quantum memory development based on novel two-dimensional materials,^{47} rare Earth ion-doped optical fibers,^{48,49} and quantum dots,^{50} which would operate at the telecommunication wavelength of the quantum memory developments, thus obviating the need for the frequency conversion step. The variety of materials platforms^{51,52} for the quantum network components make the design of quantum networks a challenging engineering task. Preliminary simulations are necessary for designing optical-fiber classical connections specifically for quantum networks.^{53–59} This kind of simulation allows a designer to construct the most tolerant protocol to account for classical network ping, single-photon traveling time fluctuations, and other parameter imperfections in the quantum network hardware, and evaluate if the existing hardware satisfies the error-robustness requirements.

One of the essential resources for a quantum network is the entanglement between any two network nodes, which would allow the transfer of unknown states between the nodes via quantum teleportation. Many protocols have been suggested to generate and distribute such remote entanglement,^{60–63} with photonic networks being the most mature and technologically realistic approach. For example, one family of protocols^{64,65,77} aims to exclude troubles associated with the imperfection of quantum memories in intermediate nodes. However, this requires that all the intermediate links between the end-nodes succeed simultaneously. Such an approach reduces the impact of qubit decoherence, but the need for synchronicity can be demanding. Other protocols rely on qubits to preserve an entangled state for a long time until all the links are successfully established between the end-nodes. Such an approach reduces the number of attempts but requires the qubits to stay in the entangled states with high fidelity for long times.^{20,66,78,79}

Here, we analyze the entanglement generation capabilities between the end-nodes in a simple linear network using SeQEeNCe quantum network simulator.^{54,57} Although several codes are available for simulating quantum networks with similar capabilities, we decided to use SeQUeNCe because it is easy to modify the code for non-developers. We examine the two protocols discussed above in detail and apply them to the generation of elementary links between adjacent nodes. Our results for both scenarios set the requirements on the number of repeaters and memory lifetime for different distances between the end-nodes.

## II. GENERATION OF THE ELEMENTARY QUANTUM LINK

An elementary link between two quantum network nodes is established when quantum entanglement is created between a pair of qubits belonging to the nodes. The entanglement, once created, is a resource used for quantum information transmission from one node to the other. Once a quantum state from the sending node is transmitted, the entanglement is destroyed. An essential parameter characterizing the performance of a given quantum connection is the maximum rate at which entanglement is generated.

Many protocols have been proposed to generate entanglement between two qubits. Here, we choose to analyze the Barrett–Kok protocol's performance^{67} for entanglement generation between two nodes with and without intermediate nodes (quantum repeaters). The Barrett–Kok protocol is robust against losses, and it has been implemented experimentally in several studies.^{68,69} Our analysis can be generalized to other protocols straightforwardly. An elementary link in the Barrett–Kok protocol is shown in Fig. 1. The nodes are located at a distance *L* from each other, and the Bell-state measuring (BSM) station is located right in the middle (in terms of the length of optical communication fiber). We assume that synchronization and scheduling instructions start with basic operations, such as single-photon emission. The instructions are obtained from a single node, called C-node (controlling node), one of the nodes forming a chain and participating in the end-to-end quantum connection. In Fig. 1(a), there are no intermediate nodes between the end-nodes, so an end-node is chosen as a control node. It is most efficient to assign C-node duties to the middle node in the chain, so the maximum ping to the other nodes is minimal. The presence of C-node makes the system hierarchic, i.e., managed from the single-center that allows us to focus on the analysis of the performance of the quantum network and frees us from developing more complex communication protocols specific to peer-to-peer networks.^{70}

The entanglement generation is organized in the following stages. At time *t* = 0, one of the nodes that acts as the C-node sends a message to another node to start the Barrett–Kok protocol. At timestamp $t=L/v$, where *v* is the photon velocity in optical fiber, the protocol start time is received. Once the message is received, immediately (assuming no delay) the nodes excite their qubits, so that the qubit would then be entangled with the photon in the optical cavity it resides in. Photon emission from the cavity happens with the probability $Em$ called memory efficiency. At timestamp $t=2.5L/v$, the leaked photons should reach the BSM station, where the measurement takes place. For simplicity, we assume that the signal speed in the classical channel is *v*, the same as in the quantum channel, and the lengths of the corresponding quantum and classical optical channels are the same. Typically, dark fiber infrastructure for secure communication is installed parallel with the standard optical fibers. Therefore, it is natural to assume that optical paths for the classical and quantum signals are the same. Thus, the result of measurements in BSM is received at timestamp $t=3L/v$. Additionally, BSM measurements in the Barrett–Kok protocol are done two times successively, and once results of the second round of measurements are obtained, both nodes know the final result of the entanglement generation attempt at timestamp $t=4L/v$. If the entanglement is successfully generated, both nodes immediately use it according to their program (e.g., for quantum teleportation or simultaneous measurement in different bases as required in the BB84 protocol^{71}) We increase the counter of generated entanglements by one. Neglecting delays in the measurements and classical messages formation, we assume that at a timestamp of $t=4L/v$, the described steps above are repeated from the beginning.

After escaping their cavities, photons can be lost in the fiber due to attenuation, which we model by the decay probability $e\u2212\gamma L/2$, $\gamma =\alpha /10\xd7log\u200910$, where *α* is the attenuation constant (given in dB per length) and $L/2$ is the fiber length, as shown in Fig. 1(a). The BSM station contains four single-photon detectors, each detecting a photon with probability $Ed$ called detector efficiency. If both photons are detected, two out of four Bell-states can be measured by this scheme. Ideally, the success rate for the Barrett–Kok protocol for entanglement swapping with BSM is $Eb=$ 50%.^{67} By combining probabilities for the length-dependent losses, two detectors, two quantum memories, and one BSM swapping, we arrive at the success probability for establishing entanglement between two nodes:

where $Em$ and $Ed$ are memories and detectors efficiencies, correspondingly.

Knowing the timestamp and an average number of attempts until successful entanglement generation, one can calculate how much time it takes for one successful entanglement generation. Furthermore, entanglement can exist only for a limited time $\tau mem$, called memory lifetime. If the entanglement degrades too fast before the photon reaches BSM, it would be impossible to establish a quantum link between the two nodes. In SeQUeNCe, memory lifetime is a fixed parameter counting lifetime of an entanglement since its inception when a qubit had successfully entangled with the photon. Once two photons with different scheduled expiration timestamps *t*_{1} and *t*_{2} are successfully measured by BSM, the lifetime for the generated entanglement between the two qubits is equal to $min{t1,\u2009t2}$. Defining the entanglement generation rate $R$ as reversed average time spent on the generation of entanglement, one can write the following expression:

where the top line is given by the effective attempt frequency $v/4L$ times *P*_{1} in Eq. (1). We assume that the repetition rate for the quantum state preparation is higher than *v*/*L*, and it is not a limiting factor in determining entanglement generation rate, which is not the case in the small *L* limit.^{72,73} In Fig. 1(b), we find excellent agreement between the analytical expression in Eq. (2) and numerical simulations.

## III. END-TO-END QUANTUM CONNECTION IN A NETWORK

Once an elementary quantum link between the nodes is established, it is possible to join an elementary link in a more extended end-to-end quantum link using the entanglement swapping protocol (see Fig. 2).^{74} When establishing the end-to-end connection between some pairs of distant nodes, some algorithms would choose the optimal path from one end-node to another via intermediate nodes, considered quantum routers.^{29,75,76} Once the path is determined, the problem is reduced to establishing connections between the successive nodes in the chain. Quantum routers give us dual benefits. The first one allows achieving connectivity by a much smaller number of physical communication links between the nodes than $N(N\u22121)/2$, where *N* is the number of nodes in the network. The second one is overcoming the exponential factor in Eq. (2). In the case of a network, *L* in Eq. (2) is the communication path length between the end-nodes. To recover the entanglement/fidelity losses, we divide the distance between the end-nodes A and B into several segments by setting up a number *r* of intermediate nodes (or quantum repeaters), as shown in Fig. 2 for the case of *r* = 3. For simplicity, we consider only a chain-like topology of the networks with equal distances between the adjacent nodes and generation end-to-end quantum connectivity between the end-nodes.

Suppose all the links among adjacent qubits in Fig. 2 have been established by either synchronous or independent generation (discussed below in detail). The next step is to perform entanglement swapping^{74} within each of the repeater nodes simultaneously, which would allow qubits in the end-nodes to become entangled if all the swappings are successful. Two kinds of procedures contribute to the time duration of this process. The first one is the time the middle nodes need to acknowledge the swapping results to the side nodes with entangled qubits. The second one is the time needed for all the nodes involved in the swapping to acknowledge the C-node about the operation results. The total time spent on the entanglement swapping can be estimated as $n\u0303L/2v$, where $n\u0303=log2(r+1)$ is the number of stages at which the total entanglement swapping procedure is done. We assume that the number of repeaters equals to $r=2n\u0303\u22121$. In Fig. 2, we demonstrate a case of *r* = 3, which would require $n\u0303=2$ stages to establish entanglement between the end-nodes. Our numerical simulations show analytical expressions for the entanglement generate rate in Eqs. (3) and (6), derived under the assumption that $r=2n\u0303\u22121$, works very well, even for the cases when $n\u0303$ is not an integer.

The probability of the success of swapping operation is $Es$. If swapping fails at any elementary link, the whole process fails, and all the links between the adjacent nodes must be discarded and regenerated. When a particular swapping fails, the rest of the links are not as *fresh* as newly generated anymore. If the established links had to wait for the regeneration of the broken links, they would have degraded further and become less reliable. Thus, everything has to be done from scratch in this approach, even if a single link fails to establish. The whole process is repeated until entanglement between the end-nodes is established.

## IV. SYNCHRONOUS GENERATION OF ENTANGLEMENTS BETWEEN ADJACENT NODES

Let us consider the case of entanglement generation with *r* repeaters and *r* + 1 elementary links. We try to generate entanglement between adjacent repeater nodes simultaneously and discard an attempt if one of the links fails. In such a scheme, we obtain the *freshest* entanglement because the entanglement between adjacent nodes is established at the same time, and there is no need for waiting. Due to the shorter lengths between neighboring repeater nodes, communication time between repeaters decreases, such that the total time needed for one generation attempt becomes $3L/(v(r+1))+L/v$. This time now determines the effective attempt frequency, which was $v/4L$ in Eq. (2) for a single link. The entanglement swapping is needed only when all adjacent links are successfully generated, which is a rare event: in the synchronous approach, it usually takes repeated attempts before all links are established simultaneously, while in the independent approach, there is usually a long wait before all the links are established. Therefore, we neglect the time needed for entanglement swapping but take into account the time needed for the C-node (in the middle of the repeaters chain) to inform all the nodes about the start time of synchronous entanglement generation and the time needed for the nodes to return their elementary link generation status to the C-node. The probability of successful generation of all the *r* + 1 adjacent links simultaneously equals to $[P1(L/[r+1])]r+1$. The probability of successful swapping of *r* links is given by $Esr$ (see Fig. 2 for the case of *r* = 3). A product of those two probabilities times the effective attempt frequency leads to the following entanglement generation rate:

For simplicity, when analyzing the cases with an arbitrary number of repeaters, we do not account for the memory coherence time in analytical expression Eq. (3). However, one should keep in mind that despite its efficiency in saving *freshness* of the generated entanglements, the described scheme cannot overcome the limit imposed by the distance between the end-nodes, i.e., $L\u22720.5v\tau mem$. Otherwise, the entangled memories would degrade during the photon travel time of about $2L/v/(r+1)+n\u0303L/v$, as in the no-repeater case.

Equation (3) was derived in the limit when the probability of photon loss in a quantum channel is very low. In this case, the most probable outcome of the generation attempt is a failure, and the outcomes whose duration is longer than $3L/(v(r+1))+L/v$ are negligible. However, plots in Fig. 3 show that this assumption is also applicable in the cases of relatively short distances between the end-nodes, i.e., $L\u22431$ km. The minor discrepancies between Eq. (3) and the numerical simulations are caused by the lack of path optimization notifying signal sent to the C-node from the nodes generating elementary links. In other words, we assume that the notification is sent only by the node that is closer to the C-node. Thus, the time delay between the Barrett–Kok protocol operation while generating the elementary links instead of *L*/*v* should be smaller by at least $L/v/(r+1)$. When we simulate cases when $r+1\u22602l,\u2009l\u2208N$, the C-node is not exactly in the middle of the chain. A more precise version of Eq. (3) is significantly more complicated. Thus, we keep the current form, which gives the single swapping round duration $0.5L/v$ as in the case $r+1=2l,\u2009l\u2208N$. In the current implementation of our simulations, we add *L*/*v* delay between plain SeQUeNCe's swapping protocol operation rounds to account for the delays due to the C-node notifying nodes about all the elementary links generation attempts.

As shown in Eq. (3), one feature of synchronous entanglement generation is that it gives an even lower entanglement generation rate than in the no-repeater case. However, such a scheme allows to utilize the intermediate nodes as routers and organize multiple nodes into a quantum network without the need to lay communication fiber between each pair of the network nodes. In addition, the lifetime of quantum memory needs to be much smaller than in the case of the independent generation scheme discussed below. This can be essential if solid-state-based scalable hardware components are employed for quantum repeaters.

## V. INDEPENDENT GENERATION OF ENTANGLEMENTS BETWEEN ADJACENT NODES

Another way to generate entanglement links between the end-nodes involves generation links among adjacent nodes independently, without discarding already established links. In this scheme, a link between the end-nodes can be generated with a much higher probability than in the synchronous protocol. Consequently, some of the links established early on would have to wait for a longer time, making decoherence a much more stringent requirement than in the case of the synchronous protocol.

In this scenario, the tyranny of fiber losses is overcome since the probability for establishing *r* + 1 links is no longer $P1r+1$, but instead, it is proportional to *P*_{1}. Our numerical simulations aim to establish the coefficient of proportionally. We show that the longer the time difference between the completion of the first and the last links, the more probabilistic events are tested before the success. As a result, a better statistic is obtained, which is the case for any Monte-Carlo type simulation. Since the success probability for the generation of an elementary link *P*_{1} is given by Eq. (1), the success probability after *k* attempts is given by

Note that the probability distribution *P _{k}* is normalized, since $\u2211k=1\u221ePk=1$. The average number of attempts until success and its standard deviation are given by

where $fk\xaf$ means averaging the function *f _{k}* over the probability distribution given by Eq. (4). In most cases, $P1\u226a1$, so that $\delta k2\xaf\u22481/P1$.

We simulate a set of *N* processes corresponding to $N=r+1$ elementary links in the chain with *r* repeaters. Each process is described by the distribution in Eq. (4). Using the pseudo-random number generation, we find for each process a number of attempts *k _{i}* until success, where $i=1,\u2026,N$. We define a value $\mu (N)$ as a maximum number of attempts normalized to the statistical average number $k\xaf$, such that $\mu (N)=min{ki}/k\xaf$. Figure 4 shows the result for $\mu (N)$ using elementary link success probability $P1=0.001$. For smaller values of

*P*

_{1}, the result does not depend on

*P*

_{1}.

For $N\u22648$, the obtained dependence can be approximated well by $\mu (N)\u2248N$ dependence. Therefore, in the case when *r* + 1 entanglements are generated by the pairs of nodes independently, the average number of attempts until the last link is generated is given by $\mu (r+1)/P1$. Analogously to Eqs. (2) and (3), one can obtain an entanglement generation rate as a product of probability to generate *N* independent links, $P1(L/[r+1])/\mu (r+1)$, entanglement swapping probability $Esr$, and effective attempt frequency to arrive at

Note that the effective attempt frequency, in this case, is reduced as compared to the synchronous case in Eq. (3) because the time *L*/*v* needed to inform the control node that all independent links are established is much shorter than the time needed to establish the links.

The average age of the oldest entangled memory, when the end-to-end entanglement is established, equals

where the first term in Eq. (7) corresponds to the time needed to generate all adjacent links and the second term for entanglement swapping.

That time is much longer than in the synchronous generation scenario, and from the previous two equations, one finds that $\Delta t\u226bL/v$. Consequently, quantum memory lifetime should be at least $\tau mem\u2273\Delta t$ to maintain quantum states during the entanglement generation time between the end-nodes. However, the simulations in Fig. 5 show that, indeed, after $Rind\u22121$ reaches approximately $\tau mem\u22121$, the entanglement generation rate starts to degrade faster than the analytical expression Eq. (6), although not as rapidly as in the synchronous scheme considered above.

Figure 5 demonstrates that in the limit of long memory lifetime, Eq. (6) works very well for predicting entanglement generation rate as a function of distance, a number of repeaters, and parameters defining hardware performance. Therefore, it can be used for a reverse problem of finding the required parameters of the quantum local area network hardware components for a given distance and entanglement generation rate.

## VI. CONCLUSION

We have simulated two types of protocols of elementary quantum links in a chain-like network supporting quantum entanglement between the end-nodes. The hardware components and basic protocols for entanglement generation and swapping are adopted from the SeQUeNCe package. We have explored the effects of finite memory lifetime on entanglement generation in a quantum network for two entanglement swapping protocols. For the synchronous generation of entanglement between adjacent nodes, the advantages of quantum repeaters are limited. On the other hand, for the independent generation of entanglement between adjacent nodes, additional repeaters enable communication at much longer distances between the end-nodes. We present analytical solutions for the entanglement generation rate for both scenarios, which are almost exact in the limit of infinitely long quantum memory lifetimes. In both cases, our numerical simulations demonstrate that entanglement generation degrades due to a finite quantum memory lifetime, whereas this degradation is less severe in the independent entanglement generation scenario. Our results demonstrate the ultimate performance of a quantum network as a function of the parameters defined by the network's hardware components, the number of repeaters, the distance between the end-nodes, and the lifetime of the quantum memories.

## ACKNOWLEDGMENTS

We gratefully acknowledge Kathy-Anne Soderberg and David Hucul from Air Force Research Lab, Rome, New York, for insightful discussions and funding from the Vice President for Research and Economic Development (VPRED) at the University at Buffalo and State University of New York Research Seed Grant Program, and computational facilities at the Center for Computational Research at the University at Buffalo (http://hdl.handle.net/10477/79221). X.H. acknowledges support by U.S. Army Research Lab via Grant No. W911NF1710257. E.F. acknowledges support by National Science Foundation via Grant No. NSF PHY-1707919.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

## DATA AVAILABILITY

The data that support the findings of this study are available within the article.