Abstract
Crucial problems of the quantum Internet are the derivation of stability properties of quantum repeaters and theory of entanglement rate maximization in an entangled network structure. The stability property of a quantum repeater entails that all incoming density matrices can be swapped with a target density matrix. The strong stability of a quantum repeater implies stable entanglement swapping with the boundness of stored density matrices in the quantum memory and the boundness of delays. Here, a theoretical framework of noise-scaled stability analysis and entanglement rate maximization is conceived for the quantum Internet. We define the term of entanglement swapping set that models the status of quantum memory of a quantum repeater with the stored density matrices. We determine the optimal entanglement swapping method that maximizes the entanglement rate of the quantum repeaters at the different entanglement swapping sets as function of the noise of the local memory and local operations. We prove the stability properties for non-complete entanglement swapping sets, complete entanglement swapping sets and perfect entanglement swapping sets. We prove the entanglement rates for the different entanglement swapping sets and noise levels. The results can be applied to the experimental quantum Internet.
Similar content being viewed by others
Introduction
The quantum Internet allows legal parties to perform networking based on the fundamentals of quantum mechanics1,2,3,4,5,6,7,8,9,10,11,12. The connections in the quantum Internet are formulated by a set of quantum repeaters and the legal parties have access to large-scale quantum devices13,14,15,16,17 such as quantum computers18,19,20,21,22,23,24,25,26,27,28,29. Quantum repeaters are physical devices with quantum memory and internal procedures3,4,5,6,9,14,15,16,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48. An aim of the quantum repeaters is to generate the entangled network structure of the quantum Internet via entanglement distribution49,50,51,52,53,54,55,56,57,58,59,60,61. The entangled network structure can then serve as the core network of a global-scale quantum communication network with unlimited distances (due to the attributes of the entanglement distribution procedure). Quantum repeaters share entangled states over shorter distances; the distance can be extended by the entanglement swapping operation in the quantum repeaters3,5,6,10,14,15,16. The swapping operation takes an incoming density matrix and an outgoing density matrix; both density matrices are stored in the local quantum memory of the quantum repeater62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79. The incoming density matrix is half of an entangled state such that the other half is stored in the distant source node, while the outgoing density matrix is half of an entangled state such that the other half is stored in the distant target node. The entanglement swapping operation, applied on the incoming and outgoing density matrices in a particular quantum repeater, entangles the distant source and target quantum nodes. Crucial problems here are the size and delay bounds connected to the local quantum memory of a quantum repeater and the optimization of the swapping procedure such that the entanglement rate of the quantum repeater (outgoing entanglement throughput measured in entangled density matrices per a preset time unit) is maximal. These questions lead us to the necessity of strictly defining the fundamental stability and performance criterions80,81,82,83,84,85,86,87,88 of quantum repeaters in the quantum Internet.
Here, a theoretical framework of noise-scaled stability analysis and entanglement rate maximization is defined for the quantum Internet. By definition, the stability of a quantum repeater can be weak or strong. The strong stability implies weak stability, by some fundamentals of queueing theory89,90,91,92,93. Weak stability of a quantum repeater entails that all incoming density matrices can be swapped with a target density matrix. Strong stability of a quantum repeater further guarantees the boundness of the number of stored density matrices in the local quantum memory. The defined system model of a quantum repeater assumes that the incoming density matrices are stored in the local quantum memory of the quantum repeater. The stored density matrices formulate the set of incoming density matrices (input set). The quantum memory also consists of a separate set for the outgoing density matrices (output set). Without loss of generality, the cardinality of the input set (number of stored density matrices) is higher than the cardinality of the output set. Specifically, the cardinality of the input set is determined by the entanglement throughput of the input connections, while the cardinality of the output set equals the number of output connections. Therefore, if in a given swapping period, the number of incoming density matrices exceeds the cardinality of the output set, then several incoming density matrices must be stored in the input set (Note: The logical model of the storage mechanisms of entanglement swapping in a quantum repeater is therefore analogous to the logical model of an input-queued switch architecture89,90,91). The aim of entanglement swapping is to select the density matrices from the input and output sets, such that the outgoing entanglement rate of the quantum repeater is maximized; this also entails the boundness of delays. The maximization procedure characterizes the problem of optimal entanglement swapping in the quantum repeaters.
Finding the optimal entanglement swapping means determining the entanglement swapping between the incoming and outgoing density matrices that maximizes the outgoing entanglement rate of the quantum repeaters. The problem of entanglement rate maximization must be solved for a particular noise level in the quantum repeater and with the presence of various entanglement swapping sets. The noise level in the proposed model is analogous to the lost density matrices in the quantum repeater due to imperfections in the local operations and errors in the quantum memory units. The entanglement swapping sets are logical sets that represent the actual state of the quantum memory in the quantum repeater. The entanglement swapping sets are formulated by the set of received density matrices stored in the local quantum memory and the set of outgoing density matrices, which are also stored in the local quantum memory. Each incoming and outgoing density matrix represent half of an entangled system, such that the other half of an incoming density matrix is stored in the distant source quantum repeater, while the other half of an outgoing density matrix is stored in the distant target quantum repeater. The aim of determining the optimal entanglement swapping method is to apply the local entanglement swapping operation on the set of incoming and outgoing density matrices such that the outgoing entanglement rate of the quantum repeater is maximized at a particular noise level. As we prove, the entanglement rate maximization procedure depends on the type of entanglement swapping sets formulated by the stored density matrices in the quantum memory. We define the logical types of the entanglement swapping sets and characterize the main attributes of the swapping sets. We present the efficiency of the entanglement swapping procedure as a function of the local noise and its impacts on the entanglement rate. We prove that the entanglement swapping sets can be defined as a function of the noise, which allows us to define noise-scaled entanglement swapping and noise-scaled entanglement rate maximization. The proposed theoretical framework utilizes the fundamentals of queueing theory, such as the Lyapunov methodology89, which is an analytical tool used to assess the performance of queueing systems89,90,91,92,93,94, and defines a fusion of queueing theory with quantum Shannon theory95,96,97,98,99,100,101,102,103 and the theory of quantum Internet.
The novel contributions of our manuscript are as follows:
-
1.
We define a theoretical framework of noise-scaled entanglement rate maximization for the quantum Internet.
-
2.
We determine the optimal entanglement swapping method that maximizes the entanglement rate of a quantum repeater at the different entanglement swapping sets as a function of the noise level of the local memory and local operations.
-
3.
We prove the stability properties for non-complete entanglement swapping sets, complete entanglement swapping sets and perfect entanglement swapping sets.
-
4.
We prove the entanglement rate of a quantum repeater as a function of the entanglement swapping sets and the noise level.
This paper is organized as follows. In Section 2, the preliminary definitions are discussed. Section 3 proposes the noise-scaled stability analysis. In Section 4, the noise-scaled entanglement rate maximization is defined. Section 5 provides a performance evaluation. Finally, Section 6 concludes the results. Supplemental Information is included in the Appendix.
System Model and Problem Statement
System model
Let \(V\) refer to the nodes of an entangled quantum network \(N\), which consists of a transmitter node \(A\in V\), a receiver node \(B\in V\), and quantum repeater nodes \({R}_{i}\in V\), \(i\,=\,\mathrm{1,}\ldots ,q\). Let \(E=\{{E}_{j}\}\), \(j\,=\,\mathrm{1,}\ldots ,m\) refer to a set of edges (an edge refers to an entangled connection in a graph representation) between the nodes of \(V\), where each \({E}_{j}\) identifies an \({{\rm{L}}}_{l}\)-level entanglement, \(l\,=\,\mathrm{1,}\ldots ,r\), between quantum nodes \({x}_{j}\) and \({y}_{j}\) of edge \({E}_{j}\), respectively. Let \(N\,=\,(V,{\mathscr{S}})\) be an actual quantum network with \(|V|\) nodes and a set \({\mathscr{S}}\) of entangled connections. An \({{\rm{L}}}_{l}\)-level, \(l=1,\ldots ,r\), entangled connection \({E}_{{{\rm{L}}}_{l}}(x,y)\), refers to the shared entanglement between a source node \(x\) and a target node \(y\), with hop-distance
since the entanglement swapping (extension) procedure doubles the span of the entangled pair in each step. This architecture is also referred to as the doubling architecture10,14,15,16.
For a particular \({{\rm{L}}}_{l}\)-level entangled connection \({E}_{{{\rm{L}}}_{l}}(x,y)\) with hop-distance (1), there are \(d{(x,y)}_{{{\rm{L}}}_{l}}-1\) intermediate nodes between the quantum nodes \(x\) and \(y\).
Figure 1 depicts a quantum Internet scenario with an intermediate quantum repeater \({R}_{j}\). The aim of the quantum repeater is to generate long-distance entangled connections between the distant quantum repeaters. The long-distance entangled connections are generated by the \({U}_{S}\) entanglement swapping operation applied in \({R}_{j}\). The quantum repeater must manage several different connections with heterogeneous entanglement rates. The density matrices are stored in the local quantum memory of the quantum repeater. The aim is to find an entanglement swapping in \({R}_{j}\) that maximizes the entanglement rate of the quantum repeater.
Entanglement fidelity
The aim of the entanglement distribution procedure is to establish a \(d\)-dimensional entangled system between the distant points \(A\) and \(B\), through the intermediate quantum repeater nodes. Let \(d\,=\,2\), and let \(|{\beta }_{00}\rangle \) be the target entangled system \(A\) and \(B\), \(|{\beta }_{00}\rangle =\frac{1}{\sqrt{2}}(|00\rangle +|11\rangle ),\) subject to be generated. At a particular density \(\sigma \) generated between \(A\) and \(B\), the fidelity of \(\sigma \) is evaluated as
Without loss of generality, an aim of a practical entanglement distribution is to reach \(F\ge 0.98\) in (2) for a given \(\sigma \)10,11,12,14,15,16,17,30.
Entanglement purification and entanglement throughput
Entanglement purification69,104,105 is a probabilistic procedure that creates a higher fidelity entangled system from two low-fidelity Bell states. The entanglement purification procedure yields a Bell state with an increased entanglement fidelity \(F{\prime} \),
where \({F}_{in}\) is the fidelity of the imperfect input Bell pairs. The purification requires the use of two-way classical communications10,11,12,14,15,16,17,30.
Let \({B}_{F}({E}_{{{\rm{L}}}_{l}}^{i})\) refer to the entanglement throughput of a given \({{\rm{L}}}_{l}\) entangled connection \({E}_{{{\rm{L}}}_{l}}^{i}\) measured in the number of \(d\)-dimensional entangled states established over \({E}_{{{\rm{L}}}_{l}}^{i}\) per sec at a particular fidelity \(F\) (dimension of a qubit system is \(d\,=\,2\))10,11,12,14,15,16,17,30.
For any entangled connection \({E}_{{{\rm{L}}}_{l}}^{i}\), a condition \(C\) should be satisfied, as
where \({B}_{F}^{\ast }({E}_{{{\rm{L}}}_{l}}^{i})\) is a critical lower bound on the entanglement throughput at a particular fidelity \(F\) of a given \({E}_{{{\rm{L}}}_{l}}^{i}\), i.e., \({B}_{F}({E}_{{{\rm{L}}}_{l}}^{i})\) of a particular \({E}_{{{\rm{L}}}_{l}}^{i}\) has to be at least \({B}_{F}^{\ast }({E}_{{{\rm{L}}}_{l}}^{i})\).
Definitions
Some preliminary definitions for the proposed model are as follows.
Definition 1
(Incoming and outgoing density matrix). In a \(j\)-th quantum repeater \({R}_{j}\), an \(\rho \) incoming density matrix is half of an entangled state \(|{\beta }_{00}\rangle =\frac{1}{\sqrt{2}}(|00\rangle +|11\rangle )\) received from a previous neighbor node \({R}_{j-1}\). The \(\sigma \) outgoing density matrix in \({R}_{j}\) is half of an entangled state \(|{\beta }_{00}\rangle \) shared with a next neighbor node \({R}_{j+1}\).
Definition 2
(Entanglement Swapping Operation). The \({U}_{S}\)entanglement swapping operation is a local transformation in a \(j\)-th quantum repeater \({R}_{j}\) that swaps an incoming density matrix \(\rho \) with an outgoing density matrix \(\sigma \) and measures the density matrices to entangle the distant source and target nodes \({R}_{j-1}\) and \({R}_{j+1}\).
Definition 3
(Entanglement Swapping Period). Let \(C\) be a cycle with time \({t}_{C}\,=\,\mathrm{1/}{f}_{C}\) determined by the \({o}_{C}\) oscillator in node \({R}_{j}\), where \({f}_{C}\) is the frequency of \({o}_{C}\). Then, let \({\pi }_{S}\) be an entanglement swapping period in which the set \({{\mathscr{S}}}_{I}({R}_{j})={\cup }_{i}{\rho }_{i}\) of incoming density matrices is swapped via \({U}_{S}\) with the set \({{\mathscr{S}}}_{O}({R}_{j})={\cup }_{i}{\sigma }_{i}\) of outgoing density matrices, defined as \({\pi }_{S}=x{t}_{C}\), where \(x\) is the number of C.
Definition 4
(Complete and Non-Complete Swapping Sets). Set \({{\mathscr{S}}}_{I}({R}_{j})\) formulates a complete set \({{\mathscr{S}}}_{I}^{\ast }({R}_{j})\) if set \({{\mathscr{S}}}_{I}({R}_{j})\) contains all the \(Q={\sum }_{i\mathrm{=1}}^{N}\,|{B}_{i}|\) incoming density matrices per \({\pi }_{S}\) that is received by \({R}_{j}\) during a swapping period, where \(N\) is the number of input entangled connections of \({R}_{j}\) and \(|{B}_{i}|\) is the number of incoming densities of the \(i\)-th input connection per \({\pi }_{S}\); thus, \({{\mathscr{S}}}_{I}({R}_{j})={\cup }_{i\mathrm{=1}}^{Q}{\rho }_{i}\) and \(|{{\mathscr{S}}}_{I}({R}_{j})|=Q\). Set \({{\mathscr{S}}}_{O}({R}_{j})\) formulates a complete set \({{\mathscr{S}}}_{O}^{\ast }({R}_{j})\) if \({{\mathscr{S}}}_{O}({R}_{j})\) contains all the \(N\) outgoing density matrices that are shared by \({R}_{j}\) during a swapping period \({\pi }_{S}\); thus, \({{\mathscr{S}}}_{O}({R}_{j})={\cup }_{i\mathrm{=1}}^{N}{\sigma }_{i}\) and \(|{{\mathscr{S}}}_{O}({R}_{j})|=N\).
Let \({\mathscr{S}}({R}_{j})\) be an entanglement swapping set of \({R}_{j}\), defined as
Then, \({\mathscr{S}}({R}_{j})\) is a complete swapping \({{\mathscr{S}}}^{\ast }({R}_{j})\) set, if
with cardinality
Otherwise, \({\mathscr{S}}({R}_{j})\) formulates a non-complete swapping set \({\mathscr{S}}({R}_{j})\), with cardinality
Definition 5
(Perfect Swapping Sets). A complete swapping set \({{\mathscr{S}}}^{\ast }({R}_{j})\) is a perfect swapping set
at a given \({\pi }_{S}\), if
holds for the cardinality.
Definition 6
(Coincidence set). In a given \({\pi }_{S}\), the coincidence set \({{\mathscr{S}}}_{{R}_{j}}^{({\pi }_{S})}\,(({R}_{i},{\sigma }_{k}))\) is a subset of incoming density matrices in \({{\mathscr{S}}}_{I}({R}_{j})\) of \({R}_{j}\) received from \({R}_{i}\) that requires the outgoing density matrix \({\sigma }_{k}\) from \({{\mathscr{S}}}_{O}({R}_{j})\) for the entanglement swapping. The cardinality of the coincidence set is
Definition 7
(Coincidence set increment in an entanglement swapping period) Let \(|B({R}_{i}({\pi }_{S}),{\sigma }_{k})|\) refer to the number of density matrices arriving from \({R}_{i}\) for swapping with \({\sigma }_{k}\) at \({\pi }_{S}\). This means the increment of the \({Z}_{{R}_{j}}^{({\pi }_{S}^{{\prime} })}\,(({R}_{i},{\sigma }_{k}))\) cardinality of the coincidence set is
where \({\pi }_{S}^{{\prime} }\) is the next entanglement swapping period. The derivations assume that an incoming density matrix \(\rho \)chooses a particular output density matrix \(\sigma \) for the entanglement swapping with probability \({\rm{\Pr }}(\rho ,\sigma )=x\ge 0\)(Bernoulli i.i.d.).
Definition 8
(Incoming and outgoing entanglement rate) Let \(|{B}_{{R}_{i}}({\pi }_{S})|\) be the incoming entanglement rate of \({R}_{j}\) per a given \({\pi }_{S}\), defined as
where \(|B({R}_{i}({\pi }_{S}),{\sigma }_{k})|\) refers to the number of density matrices arriving from \({R}_{i}\) for swapping with \({\sigma }_{k}\) per \({\pi }_{S}\).
Then, at a given \(|{B}_{{R}_{i}}({\pi }_{S})|\), the \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|\), the outgoing entanglement rate of \({R}_{j}\) is defined as
where \(L\) is the loss, \(0 < L\le N\), and \(D({\pi }_{S})\) is a delay measured in entanglement swapping periods caused by the optimal entanglement swapping at a particular entanglement swapping set.
Definition 9
(Swapping constraint). In a given \({\pi }_{S}\), each incoming density in \({{\mathscr{S}}}_{I}({R}_{j})\) can be swapped with at most one outgoing density, and only one outgoing density is available in \({{\mathscr{S}}}_{O}({R}_{j})\) for each outgoing entangled connection.
Definition 10
(Weak stable (stable) and strongly stable entanglement swapping). Weak stability (stability) of a quantum repeater \({R}_{j}\) entails that all incoming density matrices can be swapped with a target density matrix in \({R}_{j}\). A \(\zeta ({\pi }_{S})\) entanglement swapping in \({R}_{j}\) is weak stable (stable) if, for every \(\varepsilon > 0\), there exists a \(B > 0\), such that
where \({{\mathscr{S}}}_{I}^{({\pi }_{S})}({R}_{j})\) is the set of incoming densities of \({R}_{j}\) at \({\pi }_{S}\), while \(|{{\mathscr{S}}}_{I}^{({\pi }_{S})}({R}_{j})|\) is the cardinality of the set.
For a strongly stable entanglement swapping in \({R}_{j}\), the weak stability is satisifed and the cardinality of \({{\mathscr{S}}}_{I}^{({\pi }_{S})}({R}_{j})\) is bounded. A \(\zeta ({\pi }_{S})\) entanglement swapping in \({R}_{j}\) is strongly stable if
Noise-scaled entanglement swapping sets
Proposition 1 (Noise Scaled Swapping Sets). Let \(\gamma \) be a noise coefficient that models the noise of the local quantum memory and the local operations, \(0\le \gamma \le 1\). For \(\gamma \,=\,0\), the swapping set at a given \({\pi }_{S}\) is complete swapping set \({{\mathscr{S}}}^{\ast }({R}_{j})\), while for any \(\gamma \, > \,0\), the swapping set is non-complete swapping set \({\mathscr{S}}({R}_{j})\) at a given \({\pi }_{S}\).
In realistic situations, \(\gamma \) corresponds to the noises and imperfections of the physical devices and physical-layer operations (quantum operations, realization of quantum gates, storage errors, losses from local physical devices, optical losses, etc) in the quantum repeater that lead to the loss of density matrices. For further details on the physical-layer aspects of repeater-assisted quantum communications in an experimental quantum Internet setting, we suggest6.
Figure 2 illustrates the perfect swapping set, complete swapping set and non-complete swapping set. For both sets, the incoming densities are stored in incoming set \({{\mathscr{S}}}_{I}({R}_{j})\). Its cardinality depends on the incoming entanglement throughputs of the incoming connections. The outgoing set \({{\mathscr{S}}}_{O}({R}_{j})\) is a collection of outgoing density matrices. The outgoing matrix is half of an entangled state and the other half is shared with a distant target node.
The input set \({{\mathscr{S}}}_{I}({R}_{j})\) and output se \({{\mathscr{S}}}_{O}({R}_{j})\) of \({R}_{j}\) consist of the incoming and outgoing density matrices. For a non-complete entanglement swapping set, the noise is non-zero; therefore, loss is present in the quantum memory. As a convention of our model (see the swapping constraint in Definition 9), any density matrix loss is modeled as a “double loss” that affects both sets \({{\mathscr{S}}}_{I}({R}_{j})\) and \({{\mathscr{S}}}_{O}({R}_{j})\). Because of a loss, the \({U}_{S}\) swapping operation cannot be performed on the incoming and outgoing density matrices.
Problem statement
The problem formulation for the noise-scaled entanglement rate maximization is given in Problems 1–3.
Problem 1 Determine the entanglement swapping method that maximizes the entanglement rate of a quantum repeater at the different entanglement swapping sets as a function of the noise level of the local memory and local operations.
Problem 2 Prove the stability for non-complete entanglement swapping sets, complete entanglement swapping sets and perfect entanglement swapping sets.
Problem 3 Determine the outgoing entanglement rate of a quantum repeater as a function of the entanglement swapping sets and the noise level.
Problem 4 Define the optimal entanglement swapping period length as a function of the noise level at the different entanglement swapping sets.
The resolutions of Problems 1–4 are proposed in Proposition 1 and Theorems 1–4.
Entanglement Swapping Stability at Swapping Sets
This section presents the stability analysis of the quantum repeaters for the different entanglement swapping sets.
Proposition 2 (Noise-scaled weight coefficient). Let \(\omega (\gamma ({\pi }_{S}))\) be a weight at a non-complete swapping set \({\mathscr{S}}({R}_{j})\) at \(\gamma ({\pi }_{S}) > 0\), where \(\gamma ({\pi }_{S})\) is the noise at a given \({\pi }_{S}\). For a \({{\mathscr{S}}}^{\ast }({R}_{j})\) complete swapping set, \(\gamma ({\pi }_{S})\,=\,0\), and \(\omega (\gamma ({\pi }_{S})\,=\,0)={\omega }^{\ast }({\pi }_{S})\) is a maximized weight at \({\pi }_{S}\). For any non-complete set \({\mathscr{S}}({R}_{j})\), \(\omega (\gamma ({\pi }_{S}))\ge {\omega }^{\ast }({\pi }_{S})-f(\gamma ({\pi }_{S}))\), where f(·) is a sub-linear function, and \(\omega (\gamma ({\pi }_{S})) < {\omega }^{\ast }({\pi }_{S})\). At a perfect swapping set \(\hat{{\mathscr{S}}}({R}_{j})\), the weight is \(\hat{\omega }({\pi }_{S})\le {\omega }^{\ast }({\pi }_{S})\).
Proof. At a given \({\pi }_{S}\), let \({\zeta }_{ik}({\rho }_{A},{\sigma }_{k})\) be constant, with respect to the swapping constraint, defined as
thus \({\zeta }_{ik}({\rho }_{A},{\sigma }_{k})=1\), if incoming density \({\rho }_{A}\) is selected from set \({{\mathscr{S}}}_{{R}_{j}}^{({\pi }_{S})}\,(({R}_{i},{\sigma }_{k}))\) for the swapping with \({\sigma }_{k}\), and 0 otherwise. The aim is therefore to construct a
feasible entanglement swapping method for all input and output neighbors of \({R}_{j}\), for all \({\pi }_{S}\) entanglement swapping periods.
Then, from \({Z}_{{R}_{j}}^{({\pi }_{S})}\,(({R}_{i},{\sigma }_{k}))\) (see Definition 7) and (17), a \(\omega ({\pi }_{S})\) weight coefficient can be defined for a given entanglement swapping \(\zeta ({\pi }_{S})\) at a given \({\pi }_{S}\), as
where \(\langle \,\cdot \,\rangle \) is the inner product, \({Z}_{{R}_{j}}({\pi }_{S})\) is a matrix of all coincidence set cardinalities for all input and output connections at \({\pi }_{S}\), defined as
while \(\zeta ({\pi }_{S})\) is as given in (18).
For a perfect and complete entanglement swapping set, at \(\gamma ({\pi }_{S})\,=\,0\), let \({\zeta }^{\ast }({\pi }_{S})\) refer to the entanglement swapping, with \(\omega \,(\gamma ({\pi }_{S})=0)={\omega }^{\ast }({\pi }_{S})\), where \({\omega }^{\ast }({\pi }_{S})\) is the maximized weight coefficient defined as
where \({\zeta }^{\ast }({\pi }_{S})\) is an optimal entanglement swapping method at \(\gamma ({\pi }_{S})=0\) (in general not unique, by theory) with the maximized weight. By some fundamental theory89,90,91,93, it can be verified that for a non-complete set with entanglement swapping \(\chi ({\pi }_{S})\) at \(\gamma ({\pi }_{S}) > 0\), defined as
with norm a \(|\chi ({\pi }_{S})|\)89,91, defined as
with \(|\chi ({\pi }_{S})|\le 1\), the relation
holds for the weights.
Then, let \({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}))\) be a Lyapunov function89,90,91,93 of \({Z}_{{R}_{j}}({\pi }_{S})\) (see (20)) as
Then it can be verified89,90,91,93 that
holds, as \({Z}_{{R}_{j}}({\pi }_{S})\) is sufficiently large89,90,91,93, where \(\varepsilon > 0\). Since (26) is analogous to the condition on strong stability given in (16), it follows that as (21) holds for all \({\pi }_{S}\) entanglement swapping periods, the entanglement swapping at \(\gamma ({\pi }_{S})=0\) in \({R}_{j}\) is a strongly stable entanglement swapping with maximized weight coefficients for all periods.
Since for any complete and perfect entanglement swapping set, the noise is zero, the \(\hat{\omega }({\pi }_{S})\) weight coefficient of at a perfect swapping set \(\hat{{\mathscr{S}}}({R}_{j})\) is also a maximized weight with \(f(\gamma ({\pi }_{S}))=0\) as
As the noise is nonzero, \(\gamma ({\pi }_{S}) > 0\), an \(\zeta ({\pi }_{S})\ne {\zeta }^{\ast }({\pi }_{S})\) entanglement swapping cannot reach the \({\omega }^{\ast }({\pi }_{S})\) maximized weight coefficient (21), thus
where the non-zero noise \(\gamma ({\pi }_{S})\) decreases \({\omega }^{\ast }({\pi }_{S})\) to \(\omega (\gamma ({\pi }_{S})) < {\omega }^{\ast }({\pi }_{S})\) as
where \(f(\,\cdot \,)\) is a sub-linear function, as
and
where \(\gamma ({\pi }_{S})\ge \gamma ({\pi }_{S}\,=\,0)\), and \(c > 0\) is a constant.
For a non-complete set, (26) can be rewritten as
where \(\varepsilon > 0\), which represents the stability condition given in (15). As follows, an entanglement swapping at a non-complete entanglement swapping set is stable. \(\blacksquare \)
These statements are further improved in Theorems 1 and 2, respectively.
Non-complete swapping sets
Theorem 1
(Noise-scaled stability at non-complete swapping sets) An \(\zeta ({\pi }_{S})\) entanglement swapping at \(\gamma ({\pi }_{S})\mathrm{ > 0}\) is stable for any non-complete entanglement swapping set \({\mathscr{S}}({R}_{j})\).
Proof . Let \(\gamma ({\pi }_{S}) > 0\) be the noise at a given \({\pi }_{S}\), and let \(\zeta ({\pi }_{S})\) be the actual entanglement swapping at any non-complete entanglement swapping set \({\mathscr{S}}({R}_{j})\). Using the formalism of93, let \({\mathscr{L}}(X)\) refer to a Lyapunov function of a \(M\times M\) size matrix \(X\), defined as
where \({x}_{ik}\) is the \((i,k)\)-th element of \(X\).
Let \({C}_{1}\) and \({C}_{2}\) constants, \({C}_{1} > 0\), \({C}_{2} > 0\). Then, an \(\zeta ({\pi }_{S})\) entanglement swapping with \(\gamma ({\pi }_{S}) > 0\) is stable if only
where \({\Delta }_{{\mathscr{L}}}\) is a difference of the Lyapunov functions \({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}^{{\prime} }))\) and \({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}))\), where \({\pi }_{s}^{{\prime} }\) is a next entanglement swapping period, defined as
and
To verify (34), first \({\Delta }_{{\mathscr{L}}}\) is rewritten via (25) as
where \({Z}_{{R}_{j}}({\pi }_{s}^{{\prime} })\) can be evaluated as
where \(|\bar{B}({R}_{i}({\pi }_{s}^{{\prime} }),{\sigma }_{k})|\le 1\) is the normalized number of arrival density matrices from \({R}_{i}\) for swapping with \({\sigma }_{k}\) at a next entanglement swapping period \({\pi }_{s}^{{\prime} }\), defined as
where \(|{B}_{{R}_{j}}({\pi }_{S})|={\sum }_{i,k}\,|B({R}_{i}({\pi }_{S}),{\sigma }_{k})|\) is a total number of incoming density matrices of \({R}_{j}\) from the \(N\) quantum repeaters.
Using (38), the result in (37) can be rewritten
thus (34) can be rewritten as
where \({\mathbb{E}}(|\bar{B}({R}_{i}({\pi }_{S}),{\sigma }_{k})|)\) is the expected normalized number of density matrices arrive from \({R}_{i}\) for swapping with \({\sigma }_{k}\) at \({\pi }_{S}\).
Let \(\omega ({\pi }_{S}(\gamma ({\pi }_{S}) > 0))\) be the weight coefficient of the \(\zeta ({\pi }_{S}(\gamma ({\pi }_{S}) > 0))\) actual entanglement swapping at the noise \(\gamma ({\pi }_{S}) > 0\) at a given \({\pi }_{S}\), as
and let \({\alpha }_{ik}\) be defined as
Then, by some fundamental theory89,90,91,93,
where \({\nu }_{z}\ge 0\) is a constant, and \({\zeta }_{z}({\pi }_{S}(\gamma ({\pi }_{S}) > 0))\) is a \(z\)-th entanglement swapping at a noise \(\gamma ({\pi }_{S}) > 0\), while \({\omega }_{z}^{{\prime} }({\pi }_{S}(\gamma ({\pi }_{S}) > 0))\) is the weight of \({\zeta }_{z}({\pi }_{S}(\gamma ({\pi }_{S}) > 0))\).
Then, using (44), \({\mathbb{E}}({\Delta }_{{\mathscr{L}}}|{Z}_{{R}_{j}}({\pi }_{S}))\) from (41) can be rewritten as
where C1 is set as93
Therefore, there as \({C}_{2}\) is selected such that
then
holds, by theory89,90,91,93, since (48) can be rewritten as
where \(\varepsilon \) is defined as
therefore the stability condition in (15) is satisfied via
that concludes the proof. \(\blacksquare \)
As a corollary, (34) is satisfied for the \(\zeta ({\pi }_{S})\) entanglement swapping method with any \(\gamma ({\pi }_{S})\, > \,0\) non-zero noise, at a given \({\pi }_{S}\).
Complete and perfect swapping sets
Lemma 1 extends the results for entanglement swapping at complete and perfect swapping sets.
Lemma 1
(Noise-scaled stability at perfect and complete swapping sets) An \({\zeta }^{\ast }({\pi }_{S})\) optimal entanglement swapping at \(\gamma ({\pi }_{S})\,=\,0\) is strongly stable for any complete \({{\mathscr{S}}}^{\ast }({R}_{j})\) and perfect \(\hat{{\mathscr{S}}}({R}_{j})\) entanglement swapping set.
Proof. Let \({\omega }^{\ast }(\gamma ({\pi }_{S})=0)\) be the weight coefficient of \({\zeta }^{\ast }({\pi }_{S})\) at \(\gamma ({\pi }_{S})\,=\,0\) at a given \({{\mathscr{S}}}^{\ast }({R}_{j})\), as given in (21), with \({\zeta }^{\ast }({\pi }_{S})\) as
where \({\mathscr{S}}(\zeta ({\pi }_{S}))\) is the set of all possible \(N!\) entanglement swapping operations at a given \({\pi }_{S}\), and at \(N\) outgoing density matrices, \(|{{\mathscr{S}}}_{O}({R}_{j})|=N\). For a \(\hat{{\mathscr{S}}}({R}_{j})\) perfect entanglement swapping set, \({\mathscr{S}}(\zeta ({\pi }_{S}))\) is the set of all possible \(N!\) entanglement swapping operations, since \(|{{\mathscr{S}}}_{I}({R}_{j})|=|{{\mathscr{S}}}_{O}({R}_{j})|=N\)89,91).
Then, let \({\Delta }_{{\mathscr{L}}}\) be a difference of the Lyapunov functions \({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}))\) and \({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}^{{\prime} }))\),
and
where \({\pi }_{S}^{{\prime} }\) is a next entanglement swapping period; as
Then, for any complete swapping set \({{\mathscr{S}}}^{\ast }({R}_{j})\), from (45), \({\mathbb{E}}({\Delta }_{{\mathscr{L}}}|{Z}_{{R}_{j}}({\pi }_{S}))\) at \({\zeta }^{\ast }({\pi }_{S})\) is as
where \({C}_{1}\) is set as in (46), and by some fundamentals of queueing theory89,90,91,93, the condition in (16) can be rewritten as
where \(\varepsilon \) as given in (50), while \(|{Z}_{{R}_{j}}^{\ast }({\pi }_{S})|\) is the cardinality of the coincidence sets of \({{\mathscr{S}}}^{\ast }({R}_{j})\) at a given \({\pi }_{S}\), as
Thus, (57) can be rewritten as
By similar assumptions, for any \(\hat{{\mathscr{S}}}({R}_{j})\) perfect entanglement swapping with cardinality \(|{\hat{Z}}_{{R}_{j}}({\pi }_{S})|\) of the coincidence sets of \(\hat{{\mathscr{S}}}({R}_{j})\) at a given \({\pi }_{S}\), the condition in (16) can be rewritten as
Thus, from (59) and (60), it follows that \({\zeta }^{\ast }({\pi }_{S})\) (52) is strongly stable for any complete and perfect entanglement swapping set, which concludes the proof. \(\blacksquare \)
Noise-Scaled Entanglement Rate Maximization
This section proposes the entanglement rate maximization procedure for the different entanglement swapping sets.
Since the entanglement swapping is stable for both complete and non-complete entanglement swapping, this allows us to derive further results for the noise-scaled entanglement rate. The proposed derivations utilize the fundamentals of queueing theory (Note: in queueing theory, Little’s law defines a connection between the \(L\) average queue length and the \(W\) average delay as \(L=\lambda W\), where \(\lambda \) is the arrival rate. The stability property is a required preliminary condition for the relation). The derivations of the maximized noise-scaled entanglement rate assume that an incoming density matrix \(\rho \) chooses a particular output density matrix \(\sigma \) for the entanglement swapping with probability \(Pr(\rho ,\sigma )=x\ge 0\).
Preliminaries
Let \(|{Z}_{{R}_{j}}({\pi }_{S})|\) be the cardinality of the coincidence sets at a given \({\pi }_{S}\), as
and let \(|{B}_{{R}_{i}}({\pi }_{S})|\) be the total number of incoming density matrices in \({R}_{j}\) per a given \({\pi }_{S}\), as
where \(|B({R}_{i}({\pi }_{S}),{\sigma }_{k})|\) refers to the number of density matrices arrive from \({R}_{i}\) for swapping with \({\sigma }_{k}\) per \({\pi }_{S}\).
From (61) and (62), let \(D({\pi }_{S})\) be the delay measured in entanglement swapping periods, as
where \({D}_{{R}_{j}}^{({\pi }_{S})}(({R}_{i},{\sigma }_{k}))\) is the delay for a given \({R}_{i}\) at a given \({\pi }_{S}\), as
At delay \(D({\pi }_{S})\ge 0\), the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) number of swapped density matrices per \({\pi }_{S}\) is
where \(0\, < \,L\le N\) is the number of lost density matrices in \({{\mathscr{S}}}_{O}({R}_{j})\) of \({R}_{j}\) per \({\pi }_{S}\) at a non-zero noise \(\gamma ({\pi }_{S})\, > \,0\), \(L\,=\,0\) if \(\gamma ({\pi }_{S})\,=\,0\), and \({\tilde{\pi }}_{S}\) is the extended period, defined as
with \({\pi }_{S}/({\tilde{\pi }}_{S})\le 1\); thus, (65) identifies the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) outgoing entanglement rate per \({\pi }_{S}\) for a particular entanglement swapping set.
Non-complete swapping sets
Theorem 2
(Entanglement rate decrement at non-complete swapping sets). For a non-complete entanglement swapping set \({\mathscr{S}}({R}_{j})\), \(\gamma > 0\), the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) outgoing entanglement rate is \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|=(1-\frac{L}{N})\frac{1}{1+D({\pi }_{S})}(|{B}_{{R}_{j}}({\pi }_{S})|)\) per \({\pi }_{S}\), where \(|{B}_{{R}_{j}}({\pi }_{S})|\) is the total incoming entanglement throughput at a given \({\pi }_{S}\), and \(D({\pi }_{S})\le |{Z}_{{R}_{j}}({\pi }_{S})|/\)\(|{B}_{{R}_{j}}({\pi }_{S})|\), where \(\mathop{\mathrm{lim}}\limits_{{\pi }_{S}\to \infty }{\mathbb{E}}[|{Z}_{{R}_{j}}({\pi }_{S})|]\le \frac{(N-L)\beta }{{C}_{1}}+\xi (\gamma )\), where \(\beta ={\sum }_{i,k}\,(|\bar{B}({R}_{i}({\pi }_{S}),{\sigma }_{k})|-{|\bar{B}({R}_{i}({\pi }_{S}),{\sigma }_{k})|}^{2})\) and \(\xi (\gamma )=\frac{(N-L)}{2{C}_{1}}f(\gamma ({\pi }_{S}))\), where \({C}_{1}\, > \,0\) is a constant.
Proof. The entanglement rate decrement for a non-complete entanglement swapping set is as follows.
After some calculations, the result in (45) can be rewritten as
where
where \(|\bar{B}({R}_{i}({\pi }_{S}),{\sigma }_{k})|\) refers to the normalized number of density matrices arrive from \({R}_{i}\) for swapping with \({\sigma }_{k}\) at \({\pi }_{S}\) as
while \(\beta \) is defined as89,93
From (67), \({\mathbb{E}}({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{s}^{{\prime} })))\) can be evaluated as
thus, after \(P\) entanglement swapping periods \({\pi }_{S}\,=\,\mathrm{0,}\ldots ,P-1\), \({\mathbb{E}}({\mathscr{L}}({Z}_{{R}_{j}}(P{\pi }_{S})))\) is yielded as
where \({\pi }_{S}(P)\) is the \(P\)-th entanglement swapping, while \({\mathbb{E}}({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}(0))))\) identifies the initial system state, thus
by theory93.
Therefore, after \(P\) entanglement swapping periods, the expected value of \(|{Z}_{{R}_{j}}({\pi }_{S})|\) can be evaluated as
where \({\mathscr{L}}({Z}_{{R}_{j}}({\pi }_{S}(P)))\ge 0\), thus, assuming that the arrival of the density matrices can be modeled as an i.i.d. arrival process, the result in (74) can be rewritten as
Then, since for any noise \(\gamma ({\pi }_{S})\) at \({\pi }_{S}\), the relation
holds, thus for any entanglement swapping period \({\pi }_{S}\), from the sub-linear property of \(f(\,\cdot \,)\), the relation
follows for \(f(\gamma ({\pi }_{S}))\).
Therefore, \(\omega (\gamma ({\pi }_{S}))\) from (29) can be rewritten as
such that for \(f(|{Z}_{{R}_{j}}({\pi }_{S})|)\), the relation
holds, which allows us to rewrite (75) in the following manner:
where
where
Therefore, the \(D({\pi }_{S})\) delay for any non-complete swapping set is as
As a corollary, the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) outgoing entanglement rate per \({\pi }_{S}\) at (84), is as
that concludes the proof. \(\blacksquare \)
Complete swapping sets
Theorem 3
(Entanglement rate decrement at complete swapping sets). For a complete entanglement swapping set \({{\mathscr{S}}}^{\ast }({R}_{j})\), \(\gamma \,=\,0\), the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) outgoing entanglement rate is \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|=\frac{1}{1+{D}^{\ast }({\pi }_{S})}|{B}_{{R}_{j}}({\pi }_{S})|\) per \({\pi }_{S}\), where \({D}^{\ast }({\pi }_{S})=|{Z}_{{R}_{j}}^{\ast }({\pi }_{S})|/|{B}_{{R}_{j}}({\pi }_{S})| < D({\pi }_{S})\), with \(\mathop{\mathrm{lim}}\limits_{{\pi }_{S}\to \infty }{\mathbb{E}}\,[|{Z}_{{R}_{j}}^{\ast }({\pi }_{S})|]\le \frac{N\beta }{{C}_{1}}\).
Proof. Since for any complete swapping set,
it follows that
therefore (82) can be rewritten for a complete swapping set as
Therefore, the \({D}^{\ast }({\pi }_{S})\) decrement for any complete swapping set is as
with relation \({D}^{\ast }({\pi }_{S}) < D({\pi }_{S})\), where \(D({\pi }_{S})\) is as in (84).
As a corollary, the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) entanglement rate at (89), is as
which concludes the proof. \(\blacksquare \)
Perfect swapping sets
Lemma 2
(Entanglement rate decrement at perfect swapping sets). For a perfect entanglement swapping set \(\hat{{\mathscr{S}}}({R}_{j})\), \(\gamma =0\), the \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|\) outgoing entanglement rate is \(|{B}_{{R}_{j}}^{{\rm{{\prime} }}}({\pi }_{S})|=\frac{1}{1+\hat{D}({\pi }_{S})}|{B}_{{R}_{j}}({\pi }_{S})|\) per \({\pi }_{S}\), where \(\hat{D}({\pi }_{S})=|{\hat{Z}}_{{R}_{j}}({\pi }_{S})|/|{B}_{{R}_{j}}({\pi }_{S})|\), with \({\mathbb{E}}[|{\hat{Z}}_{{R}_{j}}({\pi }_{S})|]=N\).
Proof. The proof trivially follows from the fact, that for a \(\hat{{\mathscr{S}}}({R}_{j})\) perfect entanglement swapping set, the \({Z}_{{R}_{j}}^{({\pi }_{S})}(({R}_{i},{\sigma }_{k}))\) cardinality of all \(N\) coincidence sets \({{\mathscr{S}}}_{{R}_{j}}^{({\pi }_{S})}(({R}_{i},{\sigma }_{k}))\), \(i\,=\,\mathrm{1,}\ldots ,N\), \(k\,=\,\mathrm{1,}\ldots ,N\), is \({Z}_{{R}_{j}}^{({\pi }_{S})}(({R}_{i},{\sigma }_{k}))\,=\,1\), thus
As follows,
and the \({B}_{{R}_{j}}^{{\prime} }({\pi }_{S})\) entanglement rate at (92), is as
The proof is concluded here. \(\blacksquare \)
Sub-linear function at a swapping period
Theorem 4
For any \(\gamma ({\pi }_{S})\, > \,0\), at a given \({\pi }_{S}^{\ast }=(1+h){\pi }_{S}\) entanglement swapping period, \(h\, > \,0\), \(f(\gamma ({\pi }_{S}))\) can be evaluated as \(f(\gamma ({\pi }_{S}))=2{\pi }_{S}^{\ast }N\).
Proof. Let \({\zeta }^{\ast }({\pi }_{S})\) be the optimal entanglement swapping at \(\gamma ({\pi }_{S})\,=\,0\) with a maximized weight coefficient \({\omega }^{\ast }({\pi }_{S})\), and let \(\chi ({\pi }_{S})\) be an arbitrary entanglement swapping with defined as
where \({\zeta }^{\ast }({\pi }_{S}-{\pi }_{S}^{\ast })\) an optimal entanglement swapping at an \(({\pi }_{S}-{\pi }_{S}^{\ast })\)-th entanglement swapping period, while \({\pi }_{S}^{\ast }\) is an entanglement swapping period, defined as
where \(h > 0\).
Then, the \(\omega ({\pi }_{S}-{\pi }_{S}^{\ast })\) weight coefficient of entanglement swapping \(\chi ({\pi }_{S})\) (94) at an \(({\pi }_{S}-{\pi }_{S}^{\ast })\)-th entanglement swapping period is as
while \(\omega ({\pi }_{S})\) at \({\pi }_{S}\) is as
It can be concluded, that the difference of (96) and (97) is as
thus (96) is at most \({\pi }_{S}^{\ast }N\) more than (97), since the weight coefficient of \(\chi ({\pi }_{S})\) can at most decrease by \(N\) every period (i.e., at a given entanglement swapping period at most \(N\) density matrix pairs can be swapped by \(\chi ({\pi }_{S})\)).
On the other hand, let
be the weight coefficient of \({\zeta }^{\ast }({\pi }_{S})\) at \({\pi }_{S}\). It also can be verified, that the corresponding relation for the difference of (96) and (97) is as
From (98) and (99), for the \(\omega ({\pi }_{S})\) coefficient of \(\chi ({\pi }_{S})\) at \({\pi }_{S}\), it follows that
therefore, function \(f(\gamma ({\pi }_{S}))\) for any non-zero noise, \(\gamma ({\pi }_{S})\, > \,0\), is evaluated as (i.e., \(\chi ({\pi }_{S})\ne {\zeta }^{\ast }({\pi }_{S})\))
while for any \(\gamma ({\pi }_{S})\,=\,0\) (i.e., \(\chi ({\pi }_{S})={\zeta }^{\ast }({\pi }_{S})\))
which concludes the proof. \(\blacksquare \)
Performance Evaluation
In this section, a numerical performance evaluation is proposed to study the delay and the entanglement rate at the different entanglement swapping sets.
Entanglement swapping period
In Fig. 3(a), the values of \({\pi }_{S}^{\ast }\) as a function of \(h\) and \({\pi }_{S}\) are depicted. \({\pi }_{S}\in [\mathrm{1,10}]\), \(h\in [\mathrm{0,1}]\). In Fig. 3(b), the values of \(f(\gamma ({\pi }_{S}))=2{\pi }_{S}^{\ast }N\) are depicted as a function of \(h\) and \(N\).
Delay and entanglement rate ratio
In Fig. 4(a), the values of \(D({\pi }_{S})\) are depicted as a function of the incoming entanglement rate \(|{B}_{{R}_{j}}({\pi }_{S})|\) for the non-complete \({\mathscr{S}}({R}_{j})\), complete \({{\mathscr{S}}}^{\ast }({R}_{j})\) and perfect \(\hat{{\mathscr{S}}}({R}_{j})\) entanglement swapping sets. The \(D({\pi }_{S})\) delay values are evaluated as \(D({\pi }_{S})=|{Z}_{{R}_{j}}({\pi }_{S})|/|{B}_{{R}_{j}}({\pi }_{S})|\) via the maximized values of (84), (92) and (92), i.e., the delay depends on the cardinality of the coincidence set and the incoming entanglement rate (This relation also can be derived from Little’s law; for details, see the fundamentals of queueing theory89,90,91,93).
In Fig. 4(b), the ratio \(r\) of the \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|\) outgoing and \(|{B}_{{R}_{j}}({\pi }_{S})|\) incoming entanglement rates is depicted as a function of the incoming entanglement rate. For the non-complete entanglement swapping set, the loss is set as \(\tilde{L}=0.2N\).
The highest \(D({\pi }_{S})\) delay values can be obtained for non-complete entanglement swapping set \({\mathscr{S}}({R}_{j})\), while the lowest delays can be found for the perfect entanglement swapping set \(\hat{{\mathscr{S}}}({R}_{j})\). For a complete entanglement swapping set \({{\mathscr{S}}}^{\ast }({R}_{j})\), the delay values are between the non-complete and perfect sets. This is because for a non-complete set, the losses due to the \(\gamma ({\pi }_{S})\, > \,0\) non-zero noise allow only an approximation of the delay of a complete set; thus, \({D}_{{R}_{j}}({\pi }_{S}) > {D}_{{R}_{j}}^{\ast }({\pi }_{S})\). For a complete entanglement swapping set, while the noise is zero, \(\gamma ({\pi }_{S})\,=\,0\), the \(|{Z}_{{R}_{j}}^{\ast }({\pi }_{S})|\) cardinality of the coincidence set is high, \(|{Z}_{{R}_{j}}^{\ast }({\pi }_{S})| > |{\hat{Z}}_{{R}_{j}}({\pi }_{S})|\); thus, the \({D}_{{R}_{j}}^{\ast }({\pi }_{S})\) delay is higher than the \({\hat{D}}_{{R}_{j}}({\pi }_{S})\) delay of a perfect entanglement swapping set, \({D}_{{R}_{j}}^{\ast }({\pi }_{S}) > {\hat{D}}_{{R}_{j}}({\pi }_{S})\). For a perfect set, the noise is zero, \(\gamma ({\pi }_{S})\,=\,0\) and the cardinality of the coincidence sets is one for all inputs; thus, \(|{Z}_{{R}_{j}}^{\ast }({\pi }_{S})| > |{\hat{Z}}_{{R}_{j}}({\pi }_{S})|=N\). Therefore, the \({\hat{D}}_{{R}_{j}}({\pi }_{S})\) delay is minimal for a perfect set.
From the relation of \({D}_{{R}_{j}}({\pi }_{S}) > {D}_{{R}_{j}}^{\ast }({\pi }_{S}) > {\hat{D}}_{{R}_{j}}({\pi }_{S})\), the corresponding delays of the entanglement swapping sets, the relation for the decrease in the outgoing entanglement rates is straightforward, as follows. As \(r\to 1\) holds for the ratio \(r\) of the \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|\) outgoing and \(|{B}_{{R}_{j}}({\pi }_{S})|\) incoming entanglement rates, then \(|{B}_{{R}_{j}}^{{\prime} }({\pi }_{S})|\to |{B}_{{R}_{j}}({\pi }_{S})|\), i.e., no significant decrease is caused by the entanglement swapping operation.
The highest outgoing entanglement rates are obtained for a perfect entanglement swapping set, which is followed by the outgoing rates at a complete entanglement swapping set. For a non-complete set, the outgoing rate is significantly lower due to the losses caused by the non-zero noise in comparison with the perfect and complete sets.
Conclusions
The quantum repeaters determine the structure and performance attributes of the quantum Internet. Here, we defined the theory of noise-scaled stability derivation of the quantum repeaters and methods of entanglement rate maximization for the quantum Internet. The framework characterized the stability conditions of entanglement swapping in quantum repeaters and the terms of non-complete, complete and perfect entanglement swapping sets in the quantum repeaters to model the status of the quantum memory of the quantum repeaters. The defined terms are evaluated as a function of the noise level of the quantum repeaters to describe the physical procedures of the quantum repeaters. We derived the conditions for an optimal entanglement swapping at a particular noise level to maximize the entanglement throughput of the quantum repeaters. The results are applicable to the experimental quantum Internet.
Ethics statement
This work did not involve any active collection of human data.
Data availability
This work does not have any experimental data.
References
Pirandola, S. & Braunstein, S. L. Unite to build a quantum internet. Nature 532, 169–171 (2016).
Pirandola, S. End-to-end capacities of a quantum communication network. Commun. Phys. 2, 51 (2019).
Wehner, S., Elkouss, D. & Hanson, R. Quantum internet: A vision for the road ahead. Science 362, 6412 (2018).
Pirandola, S., Laurenza, R., Ottaviani, C. & Banchi, L. Fundamental limits of repeaterless quantum communications. Nature Communications 8, 15043, https://doi.org/10.1038/ncomms15043 (2017).
Pirandola, S. et al. Theory of channel simulation and bounds for private communication. Quantum Sci. Technol. 3, 035009 (2018).
Pirandola, S. Capacities of repeater-assisted quantum communications. arXiv:1601.00966 (2016).
Pirandola, S. Bounds for multi-end communication over quantum networks. Quantum Sci. Technol. 4, 045006 (2019).
Pirandola, S. et al. Advances in Quantum Cryptography. arXiv:1906.01645 (2019).
Laurenza, R. & Pirandola, S. General bounds for sender-receiver capacities in multipoint quantum communications. Phys. Rev. A 96, 032318 (2017).
Van Meter, R. Quantum Networking. ISBN 1118648927, 9781118648926, John Wiley and Sons Ltd (2014).
Lloyd, S. et al. Infrastructure for the quantum Internet. ACM SIGCOMM Computer Communication Review 34, 9–20 (2004).
Kimble, H. J. The quantum Internet. Nature 453, 1023–1030 (2008).
Gyongyosi, L. & Imre, S. Optimizing High-Efficiency Quantum Memory with Quantum Machine Learning for Near-Term Quantum Devices, Scientific Reports, Nature, https://doi.org/10.1038/s41598-019-56689-0 (2019).
Van Meter, R., Ladd, T. D., Munro, W. J. & Nemoto, K. System Design for a Long-Line Quantum Repeater. IEEE/ACM Transactions on Networking 17(3), 1002–1013 (2009).
Van Meter, R., Satoh, T., Ladd, T. D., Munro, W. J. & Nemoto, K. Path selection for quantum repeater networks, Networking Science, Volume 3, Issue 1–4, pp 82–95, (2013).
Van Meter, R. & Devitt, S. J. Local and Distributed Quantum Computation. IEEE Computer 49(9), 31–42 (2016).
Gyongyosi, L., Imre, S. & Nguyen, H. V. A Survey on Quantum Channel Capacities, IEEE Communications Surveys and Tutorials, https://doi.org/10.1109/COMST.2017.2786748 (2018).
Preskill, J. Quantum Computing in the NISQ era and beyond. Quantum 2, 79 (2018).
Arute, F. et al. Quantum supremacy using a programmable superconducting processor, Nature, Vol. 574, https://doi.org/10.1038/s41586-019-1666-5 (2019).
Harrow, A. W. & Montanaro, A. Quantum Computational Supremacy. Nature 549, 203–209 (2017).
Aaronson, S. & Chen, L. Complexity-theoretic foundations of quantum supremacy experiments. Proceedings of the 32nd Computational Complexity Conference, CCC ’17, pages 22:1-22:67, (2017).
Farhi, E., Goldstone, J., Gutmann, S. & Neven, H. Quantum Algorithms for Fixed Qubit Architectures. arXiv:1703.06199v1 (2017).
Farhi, E. & Neven, H. Classification with Quantum Neural Networks on Near Term Processors, arXiv:1802.06002v1 (2018).
Alexeev, Y. et al. Quantum Computer Systems for Scientific Discovery, arXiv:1912.07577 (2019).
Loncar, M. et al. Development of Quantum InterConnects for Next-Generation Information Technologies, arXiv:1912.06642 (2019).
Shor, P. W. Algorithms for quantum computation: discrete logarithms and factoring. In: Proceedings of the 35th Annual Symposium on Foundations of Computer Science (1994).
IBM. A new way of thinking: The IBM quantum experience. http://www.research.ibm.com/quantum (2017).
Ajagekar, A., Humble, T. and You, F. Quantum Computing based Hybrid Solution Strategies for Large-scale Discrete-Continuous Optimization Problems. Computers and Chemical Engineering 132, 106630 (2020).
Foxen, B. et al. Demonstrating a Continuous Set of Two-qubit Gates for Near-term Quantum Algorithms, arXiv:2001.08343 (2020).
Gyongyosi, L. & Imre, S. Decentralized Base-Graph Routing for the Quantum Internet, Physical Review A, https://doi.org/10.1103/PhysRevA.98.022310, https://link.aps.org/doi/10.1103/PhysRevA.98.022310 (2018).
Gyongyosi, L. & Imre, S. Dynamic topology resilience for quantum networks, Proc. SPIE 10547, Advances in Photonics of Quantum Computing, Memory, and Communication XI, 105470Z; https://doi.org/10.1117/12.2288707 (2018).
Gyongyosi, L. & Imre, S. Topology Adaption for the Quantum Internet, Quantum Information Processing, Springer Nature, https://doi.org/10.1007/s11128-018-2064-x, (2018).
Gyongyosi, L. & Imre, S. Entanglement Access Control for the Quantum Internet, Quantum Information Processing, Springer Nature, https://doi.org/10.1007/s11128-019-2226-5, (2019).
Gyongyosi, L. & Imre, S. Opportunistic Entanglement Distribution for the Quantum Internet, Scientific Reports, Nature, https://doi.org/10.1038/s41598-019-38495-w, (2019).
Gyongyosi, L. & Imre, S. Adaptive Routing for Quantum Memory Failures in the Quantum Internet, Quantum Information Processing, Springer Nature, https://doi.org/10.1007/s11128-018-2153-x, (2018).
Quantum Internet Research Group (QIRG), web: https://datatracker.ietf.org/rg/qirg/about/ (2018).
Gyongyosi, L. & Imre, S. Multilayer Optimization for the Quantum Internet, Scientific Reports, Nature, https://doi.org/10.1038/s41598-018-30957-x, (2018).
Gyongyosi, L. & Imre, S. Entanglement Availability Differentiation Service for the Quantum Internet, Scientific Reports, Nature, (10.1038/s41598-018-28801-3), https://www.nature.com/articles/s41598-018-28801-3 (2018).
Gyongyosi, L. & Imre, S. Entanglement-Gradient Routing for Quantum Networks, Scientific Reports, Nature, (10.1038/s41598-017-14394-w), https://www.nature.com/articles/s41598-017-14394-w (2017).
Gyongyosi, L. & Imre, S. A Survey on Quantum Computing Technology. Computer Science Review, https://doi.org/10.1016/j.cosrev.2018.11.002 (2018).
Rozpedek, F. et al. Optimizing practical entanglement distillation. Phys. Rev. A 97, 062333 (2018).
Humphreys, P. et al. Deterministic delivery of remote entanglement on a quantum network, Nature 558 (2018).
Liao, S.-K. et al. Satellite-to-ground quantum key distribution. Nature 549, 43–47 (2017).
Ren, J.-G. et al. Ground-to-satellite quantum teleportation. Nature 549, 70–73 (2017).
Hensen, B. et al. Loophole-free Bell inequality violation using electron spins separated by 1.3 kilometres, Nature 526 (2015).
Hucul, D. et al. Modular entanglement of atomic qubits using photons and phonons, Nature Physics 11(1) (2015).
Noelleke, C. et al. Efficient Teleportation Between Remote Single-Atom Quantum Memories. Physical Review Letters 110, 140403 (2013).
Sangouard, N. et al. Quantum repeaters based on atomic ensembles and linear optics. Reviews of Modern Physics 83, 33 (2011).
Chakraborty, K., Rozpedeky, F., Dahlbergz, A. & Wehner, S. Distributed Routing in a Quantum Internet, arXiv:1907.11630v1 (2019).
Khatri, S., Matyas, C. T., Siddiqui, A. U. & Dowling, J. P. Practical figures of merit and thresholds for entanglement distribution in quantum networks. Phys. Rev. Research 1, 023032 (2019).
Kozlowski, W. & Wehner, S. Towards Large-Scale Quantum Networks, Proc. of the Sixth Annual ACM International Conference on Nanoscale Computing and Communication, Dublin, Ireland, arXiv:1909.08396 (2019).
Pathumsoot, P. et al. Modeling of Measurement-based Quantum Network Coding on IBMQ Devices, arXiv:1910.00815v1 (2019).
Pal, S., Batra, P., Paterek, T. & Mahesh, T. S. Experimental localisation of quantum entanglement through monitored classical mediator, arXiv:1909.11030v1 (2019).
Caleffi, M. End-to-End Entanglement Rate: Toward a Quantum Route Metric, 2017 IEEE Globecom, https://doi.org/10.1109/GLOCOMW.2017.8269080, (2018).
Caleffi, M. Optimal Routing for Quantum Networks, IEEE Access, Vol. 5, https://doi.org/10.1109/ACCESS.2017.2763325 (2017).
Caleffi, M., Cacciapuoti, A. S. & Bianchi, G. Quantum Internet: from Communication to Distributed Computing, aXiv:1805.04360 (2018).
Castelvecchi, D. The quantum internet has arrived, Nature, News and Comment, https://www.nature.com/articles/d41586-018-01835-3 (2018).
Cacciapuoti, A. S. et al. Quantum Internet: Networking Challenges in Distributed Quantum Computing, arXiv:1810.08421 (2018).
Miguel-Ramiro, J. & Dur, W. Delocalized information in quantum networks. arXiv:1912.12935v1 (2019).
Tanjung, K. et al. Probing quantum features of photosynthetic organisms. npj Quantum Information, 4 2056–6387 (2018).
Tanjung, K. et al. Revealing Nonclassicality of Inaccessible Objects. Phys. Rev. Lett., 119(12), 1079–7114 (2017).
Lloyd, S. & Weedbrook, C. Quantum generative adversarial learning. Phys. Rev. Lett., 121, arXiv:1804.09139 (2018).
Gisin, N. & Thew, R. Quantum Communication. Nature Photon. 1, 165–171 (2007).
Xiao, Y. F. & Gong, Q. Optical microcavity: from fundamental physics to functional photonics devices. Science Bulletin 61, 185–186 (2016).
Zhang, W. et al. Quantum Secure Direct Communication with Quantum Memory. Phys. Rev. Lett. 118, 220501 (2017).
Enk, S. J., Cirac, J. I. & Zoller, P. Photonic channels for quantum communication. Science 279, 205–208 (1998).
Briegel, H. J., Dur, W., Cirac, J. I. & Zoller, P. Quantum repeaters: the role of imperfect local operations in quantum communication. Phys. Rev. Lett. 81, 5932–5935 (1998).
Dur, W., Briegel, H. J., Cirac, J. I. & Zoller, P. Quantum repeaters based on entanglement purification. Phys. Rev. A 59, 169–181 (1999).
Dur, W. & Briegel, H. J. Entanglement purification and quantum error correction. Rep. Prog. Phys. 70, 1381 (2007).
Duan, L. M., Lukin, M. D., Cirac, J. I. & Zoller, P. Long-distance quantum communication with atomic ensembles and linear optics. Nature 414, 413–418 (2001).
Van Loock, P. et al. Hybrid quantum repeater using bright coherent light. Phys. Rev. Lett. 96, 240501 (2006).
Zhao, B., Chen, Z. B., Chen, Y. A., Schmiedmayer, J. & Pan, J. W. Robust creation of entanglement between remote memory qubits. Phys. Rev. Lett. 98, 240502 (2007).
Goebel, A. M. et al. Multistage Entanglement Swapping. Phys. Rev. Lett. 101, 080403 (2008).
Simon, C. et al. Quantum Repeaters with Photon Pair Sources and Multimode Memories. Phys. Rev. Lett. 98, 190503 (2007).
Tittel, W. et al. Photon-echo quantum memory in solid state systems. Laser Photon. Rev. 4, 244–267 (2009).
Sangouard, N., Dubessy, R. & Simon, C. Quantum repeaters based on single trapped ions. Phys. Rev. A 79, 042340 (2009).
Sheng, Y. B. & Zhou, L. Distributed secure quantum machine learning. Science Bulletin 62, 1025–1019 (2017).
Leung, D., Oppenheim, J. & Winter, A. IEEE Trans. Inf. Theory 56, 3478–90 (2010).
Kobayashi, H., Le Gall, F., Nishimura, H. & Rotteler, M. Perfect quantum network communication protocol based on classical network coding, Proceedings of 2010 IEEE International Symposium on Information Theory (ISIT) pp 2686–90. (2010).
Shor, P. W. Scheme for reducing decoherence in quantum computer memory. Phys. Rev. A 52, R2493–R2496 (1995).
Chou, C. et al. Functional quantum nodes for entanglement distribution over scalable quantum networks. Science 316(5829), 1316–1320 (2007).
Muralidharan, S., Kim, J., Lutkenhaus, N., Lukin, M. D. & Jiang, L. Ultrafast and Fault-Tolerant Quantum Communication across Long Distances. Phys. Rev. Lett. 112, 250501 (2014).
Yuan, Z. et al. Nature 454, 1098–1101 (2008).
Kobayashi, H., Le Gall, F., Nishimura, H. & Rotteler, M. General scheme for perfect quantum network coding with free classical communication, Lecture Notes in Computer Science (Automata, Languages and Programming SE-52 vol. 5555), Springer) 622–633 (2009).
Hayashi, M. Prior entanglement between senders enables perfect quantum network coding with modification. Physical Review A 76, 040301(R) (2007).
Hayashi, M., Iwama, K., Nishimura, H., Raymond, R. & Yamashita, S. Quantum network coding, Lecture Notes in Computer Science (STACS 2007 SE52 vol. 4393) ed Thomas, W. and Weil, P. (Berlin Heidelberg: Springer) (2007).
Chen, L. & Hayashi, M. Multicopy and stochastic transformation of multipartite pure states. Physical Review A 83(No. 2), 022331 (2011).
Schoute, E., Mancinska, L., Islam, T., Kerenidis, I. & Wehner, S. Shortcuts to quantum network routing, arXiv:1610.05238 (2016).
Bianco, A., Giaccone, P., Leonardi, E., Mellia, M. & Neri, F. Theoretical Performance of Input-queued Switches using Lyapunov Methodology, In: Elhanany, I., and Hamdi, M. (eds.) High-performance Packet Switching Architectures, Springer (2007).
Leonardi, E., Mellia, M., Neri, F. & Ajmone Marsan, M. On the stability of Input-Queued Switches with Speed-up. IEEE/ACM Transactions on Networking 9, 104–118 (2001).
Leonardi, E., Mellia, M., Neri F., Ajmone Marsan, M. Bounds on delays and queue lengths in input-queued cell switches. Journal of the ACM 50:520-550 23 (2003).
Ajmone Marsan, M., Bianco, A., Giaccone, P., Leonardi, E. & Neri, F. Packet-mode scheduling in input-queued cell-based switches. IEEE/ACM Transactions on Networking 10, 666–678 (2002).
Shah, D. & Kopikare, M. Delay Bounds for Approximate Maximum Weight Matching Algorithms for Input Queued Switches. Proc. of IEEE INFOCOM 2002 (2002).
Mitzenmacher, N. & Upfal, E. Probability and computing: Randomized algorithms and probabilistic analysis. Cambridge University Press (2005).
Imre, S. & Gyongyosi, L. Advanced Quantum Communications - An Engineering Approach. New Jersey, Wiley-IEEE Press (2013).
Kok, P. et al. Linear optical quantum computing with photonic qubits. Rev. Mod. Phys. 79, 135–174 (2007).
Petz, D. Quantum Information Theory and Quantum Statistics, Springer-Verlag, Heidelberg, Hiv: 6 (2008).
Bacsardi, L. On the Way to Quantum-Based Satellite Communication, IEEE Comm. Mag. 51(08), 50–55. (2013).
Biamonte, J. et al. Quantum Machine Learning. Nature 549, 195–202 (2017).
Lloyd, S., Mohseni, M. & Rebentrost, P. Quantum algorithms for supervised and unsupervised machine learning. arXiv:1307.0411 (2013).
Lloyd, S., Mohseni, M. & Rebentrost, P. Quantum principal component analysis. Nature Physics 10, 631 (2014).
Lloyd, S. Capacity of the noisy quantum channel. Physical Rev. A 55, 1613–1622 (1997).
Lloyd, S. The Universe as Quantum Computer, A Computable Universe: Understanding and exploring Nature as computation, Zenil, H. ed., World Scientific, Singapore, arXiv:1312.4455v1 (2013).
Dehaene, J. et al. Local permutations of products of Bell states and entanglement distillation. Phys. Rev. A 67, 022310 (2003).
Munro, W. J. et al. Inside quantum repeaters. IEEE Journal of Selected Topics in Quantum Electronics 78–90 (2015).
Acknowledgements
The research reported in this paper has been supported by the Hungarian Academy of Sciences (MTA Premium Postdoctoral Research Program 2019), by the National Research, Development and Innovation Fund (TUDFO/51757/2019-ITM, Thematic Excellence Program), by the National Research Development and Innovation Office of Hungary (Project No. 2017-1.2.1-NKP-2017-00001), by the Hungarian Scientific Research Fund - OTKA K-112125 and in part by the BME Artificial Intelligence FIKP grant of EMMI (Budapest University of Technology, BME FIKP-MI/SC).
Author information
Authors and Affiliations
Contributions
L.G.Y. designed the protocol and wrote the manuscript. L.G.Y. and S.I. analyzed the results. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gyongyosi, L., Imre, S. Theory of Noise-Scaled Stability Bounds and Entanglement Rate Maximization in the Quantum Internet. Sci Rep 10, 2745 (2020). https://doi.org/10.1038/s41598-020-58200-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-58200-6
This article is cited by
-
Generation of Werner-like states via a two-qubit system plunged in a thermal reservoir and their application in solving binary classification problems
Scientific Reports (2021)
-
Scalable distributed gate-model quantum computers
Scientific Reports (2021)
-
Routing space exploration for scalable routing in the quantum Internet
Scientific Reports (2020)
-
Quantum State Optimization and Computational Pathway Evaluation for Gate-Model Quantum Computers
Scientific Reports (2020)
-
Dynamics of entangled networks of the quantum Internet
Scientific Reports (2020)