Multilayer Optimization for the Quantum Internet.

We define a multilayer optimization method for the quantum Internet. Multilayer optimization integrates separate procedures for the optimization of the quantum layer and the classical layer of the quantum Internet. The multilayer optimization procedure defines advanced techniques for the optimization of the layers. The optimization of the quantum layer covers the minimization of total usage time of quantum memories in the quantum nodes, the maximization of the entanglement throughput over the entangled links, and the reduction of the number of entangled links between the arbitrary source and target quantum nodes. The objective of the optimization of the classical layer is the cost minimization of any auxiliary classical communications. The multilayer optimization framework provides a practically implementable tool for quantum network communications, or long-distance quantum communications.

Introduction

Quantum Internet is a communication network with quantum nodes and quantum links[1-16] that allows to the parties to perform efficient quantum communications[17-19]. The aim of quantum Internet[5,8] and quantum repeater networks[1,17-23] is to distribute quantum entanglement between distant nodes through a chain of intermediate quantum repeater nodes[24-33]. In the quantum Internet, the quantum nodes share entangled connections that formulate entangled links[1-6]. The quantum nodes store the quantum states in their local quantum memory for path selection and path recovery purposes[1-3,8,24-33]. Since setting several attributes must be optimized in parallel in an arbitrary quantum network, the optimization problem formulates a multi-objective procedure. Formally, multi-objective optimization covers the minimization of quantum memory usage time (storage time), the maximization of entanglement throughput (number of transmitted entangled states per second of a particular fidelity) of entangled links[1-6], and the reduction of the number of entangled links between a source and a target quantum node[8,24-33]. However, the problem does not end here since a quantum repeater network can be approached on the quantum transmission (quantum layer) level and on the auxiliary classical communication (classical layer) level that is required for the dynamic functioning of the quantum layer. Therefore, the problem is not just a multi-objective optimization problem in the quantum layer but also a multilayer optimization issue that covers the development of both the quantum and classical layers of a quantum repeater network.

Here, we define a multilayer optimization method for quantum repeater networks. This covers both the quantum layer and the classical layer of a quantum repeater network.

By utilizing the tools of quantum Shannon theory[1-7], the optimization of the quantum layer includes minimizing the usage of quantum memories in the nodes to reduce the storage time of entangled states, the maximization of entanglement throughput of the entangled links, and also these conditions have to be satisfied for the shortest path between a given source node and target quantum node (i.e., a multi-objective optimization of the quantum layer).

The aim of classical-layer optimization is to curtail the cost of auxiliary classical communications, which is required for such optimization. The cost of the classical communication covers all communication costs required to achieve the optimal quantum network state including the classical communication steps for overall quantum storage time minimization, entanglement throughput maximization, and the selection of a shortest path.

The multilayer optimization employs advanced methods to solve the multi-objective optimization of the quantum layer. We define the structures of the quantum memory utilization graph and the entanglement throughput tree for the multi-objective optimization of the quantum layer of a quantum repeater network. The quantum memory utilization graph models the quantum memory usage for entanglement storage. The entanglement throughput tree shows the entanglement throughput of entangled links with respect to the number of transmittable entangled states at a particular fidelity. Using these advanced constructions, we also define a method for the optimal assignment of entangled states in the repeater nodes. The input of the quantum layer optimization procedure is the quantum memory utilization graph, while the output of the method is a set of entanglement throughput trees. The output identifies the optimal states of the quantum network with respect to the multi-objective optimization function.

Classical-layer optimization focuses on the minimization of the total cost of classical communications by utilizing swarm intelligence[34-38]. This also defines a multi-objective problem since the cost has to be reduced with respect to the classical communication cost required for the minimization of quantum memory usage, the classical cost of entanglement throughput maximization of entangled links, and for the selection of the shortest path in the quantum layer. Classical-layer optimization uses some fundamentals of bacteria foraging models[26,30-35] and probabilistic multi-objective uncertainty characterization[31,39,40,44-46].

The optimization framework requires no changes in the physical layer, so the framework is directly implementable by the current physical devices[1-6,17-19] and quantum networking elements[8,24-33]. The method is useful in quantum networking environments with diverse physical attributes (different quantum memory characteristics, quantum error correction, physical quantum nodes attributes, and transmission capabilities of noisy quantum links).

The novel contributions of our paper are as follows:

We conceive a complex optimization framework for quantum networks. It integrates the development of the quantum and classical layers of quantum repeater networks.

Quantum-layer optimization utilizes the attributes of physical-layer quantum transmissions, quantum memory usage, and entanglement distribution via the framework of quantum Shannon theory.

Classical-layer optimization focuses on minimizing any auxiliary communications related to the quantum layer and optimization.

Multilayer optimization is applicable by current physical devices and quantum networking elements providing a solution for the optimization of arbitrary quantum networking scenarios with diverse physical attributes and environments.

This paper is organized as follows. In Section 2, the system model is proposed. In Section 3, the optimization procedure of the quantum layer of quantum repeater networks is defined. Section 4 studies the optimization of the classical layer. Section 5 provides a performance evaluation. Finally, Section 6 concludes the paper. Supplemental material is included in the Appendix.

System Model

Our model assumes that the quantum repeater network consists of a source and target node with intermediate repeater nodes and a quantum switcher node. A quantum switcher node S operates as follows. Node S is a quantum repeater node capable of switching between the entangled connections stored in its local quantum memory and a permit of applying entanglement swapping on the selected connections. While an i-th quantum repeater node establishes only the entangled connections with the neighbor quantum repeater nodes, a switcher node is equipped with an extended knowledge about the quantum network to select between the entangled links. A general repeater node is not allowed to perform any link selection, since it is assumed in the model that a quantum repeater node has only local knowledge about the network. A switcher node based on its network knowledge can also send entanglement swapping commands to the quantum repeater nodes to define new paths in the network.

Let N be a quantum network, , where V is a set of nodes, is a set of entangled links. Without loss of generality, the level L of an entangled link E(x, y) is defined as follows. For an L-level entangled link, the hop distance between quantum nodes x and y is[1,11]with intermediate nodes between the nodes x and y. The probability that an L-level entangled link E(x, y) exists between nodes x, y is , which depends on the actual network.

The S quantum switcher is modeled as a quantum node with the following attributes and permissions:

knowledge about the physical attributes of distant quantum repeater nodes and the entangled connections of N (e.g., entanglement fidelity, quantum memory status, link noise, etc),

internal quantum memory for the storage of entangled states,

quantum functionality:

permission to set new entangled connections between its local quantum system and a selected quantum node of the quantum network,

permission to switch between the stored entangled states to construct new paths,

classical functionality:

permission to command distant quantum nodes of N via classical links (to construct new entangled connections in the network, to perform entanglement swapping between the selected nodes, other).

The network model is illustrated in Fig. 1. The example network in Fig. 1(a) consists of six quantum repeater nodes R, i = 1, …, 6, and a quantum switcher S that switches between the entangled connections using its local quantum memory. The switcher also can perform entanglement swapping in the network. The switcher node S has knowledge about the physical attributes (e.g., entanglement fidelity, quantum memory status, etc) of quantum repeater nodes to make a decision on a path. The knowledge about the repeater nodes can be transmitted over a classical link to the quantum switcher (classical links are not depicted). As it is depicted in Fig. 1(b), the switcher node has a permission to set new entangled connections via its local quantum state with a selected quantum node. The S switcher node decided to set a new entangled connection between its local quantum system and repeater node R2. A standard quantum repeater is not allowed to perform these operations (except with the direct neighbors in the entanglement distribution phase) without a dedicated command from the switcher.

Figure 1

The network model with a quantum switcher S and quantum repeater nodes R, i = 1, …, 6. The L1, L2- and L3-level entangled links of N are depicted by gray, blue and orange, respectively (additional nodes are not shown). (a) For a current shortest path , the active repeater nodes selected by S are R3 and R4. (b) Node S switches the entangled connections in the local quantum memory from R3 to R2 and from R4 to R6. The switching operation defines a new shortest path .

Note, path in Fig. 1(a) provides a shortest path at a particular network situation at an initial network time T1. Since the quantum network N evolves in time (quality of the entangled links, the status of the nodes, internal quantum memories, etc), at a given time T2, by utilizing the functions of the switcher node, the switcher node determined a new shortest path, , as depicted in Fig. 1(b).

Quantum Memory Scheduling

In this section, we define a structure for scheduling quantum memory usage of the quantum nodes called the quantum memory utilization graph . This is a directed graph mapped from the network model, with several abstracted nodes and links.

Proposition 1

The quantum memory utilization graph is a directed graph with abstract nodes and links to schedule the quantum memory usage mapped from the quantum repeater network.

The graph of quantum memory utilization is constructed as follows. Assuming n quantum nodes (excluding the S switcher node) in the network, the graph contains n abstracted transmitter nodes and n abstracted receiver nodes with directed connections. Note if the quantum switcher S is modeled as the n-th node and n = 1, where n is the number of quantum switchers in N, the graph also can be constructed via (n − 1) transmitter and (n − 1) receiver quantum nodes. A given N with an arbitrary number of quantum switcher nodes defines a particular , therefore the graph is a combination of all other possible switcher modes. The graph contains directed edges between quantum node pairs (V|V) and (V|V) of a particular switcher mode S, i = 1, …, N, where N is the total number of switcher modes, depicted via nodes labeled as (x|y) and (y|z).

Let us assume that N has a single switcher node S, and it has two states, S1 and S2. In both modes S1 and S2, repeater node R1 serves as a transmitter node for node R2.

In switcher mode S1, the shared entangled connection defines the following relations. Repeater node R1 serves as a transmitter node for nodes R3 and R4. Node R3 serves as a transmitter node for nodes R4 and R5. Node R4 serves as a transmitter node for nodes R5 and R6. Node R5 serves as a transmitter node for node R6.

In switcher mode S2, the shared entangled connection defines the following relation. Repeater node R2 serves as a transmitter node for node R6.

The graph of quantum memory utilization derived from the quantum network setting of Fig. 1. is illustrated in Fig. 2.

Figure 2

The graph of quantum memory utilization derived from the network setting in Fig. 1. The graph contains n abstracted transmitter nodes and n abstracted receiver nodes with directed entangled connections. An i-level entangled connect is depicted by L. The S quantum switcher has two states, S1 and S2. The green circles represent quantum nodes operating on switcher mode S1. The yellow circles represent nodes operating on switcher mode S2.

Entanglement Throughput Tree

In this section, we define the structure of the entanglement throughput tree, which aims to extract information from the quantum memory utilization graph. The entanglement throughput tree is also the output format of the multi-objective optimization procedure of the quantum layer.

Lemma 1

The entanglement throughput tree is a structure modeling the multi-objective optimization problem of quantum repeater networks. The tree structure is derived from the quantum memory utilization graph.

Proof.

As a given source node is selected, the next nodes are added to the path with a given probability.

Let us assume that is determined. Let us index the nodes by an ID identifier tag, ID = {A, B, …}, and let S be the set of unvisited neighbor nodes of a node I.

Let refer to the entanglement throughput of a given L-level entangled link between nodes (I, J) measured in the number of d-dimensional entangled states per second at a particular fidelity F[1,3,4].

Further, let J be a neighbor node of I with entangled connection and with entanglement throughput .

A cost function, Ω(I, J), between nodes (I, J) is defined aswhere is the cost of entangled link defined aswhile ζ is the cost of quantum storage in node J.

Let be the entanglement utility coefficient of entangled link between nodes I and J, initialized as . This amount is equivalent to the utility of the entangled link that it has taken to arrive at the current node J from I.

At a given , the initial entanglement utility[4] of link is updated to as

Using these cost functions, the Pr(I, J) probability that from node I a node J is selected is as follows:where ω* and δ are weighting coefficients[4].

Using (2), (3) and (5) for all node pairs, the method to build a random entanglement throughput tree using a graph is as follows[36,37].

Let S′ refer to the set of already reached destination nodes, and let be the set of initial nodes, the set of feasible neighboring nodes to node I. Let be the set of destination nodes.

The method is given in Procedure 1.

Random Entanglement Throughput Tree Construction.

The proof is concluded here.■

The structure of a entanglement throughput tree is illustrated in Fig. 3.

Figure 3

The structure of a entanglement throughput tree. A quantum node has an ID identifier tag, ID = {A, B, …, L}, and all incoming and outcoming entangled links are identified by the s source neighbor node and the t target neighbor node. Node A represents a source, while the destination nodes are J, K, L. The entanglement throughput of all (x, y) node pairs are depicted above the directed lines; the link identifier is depicted under the links.

Entanglement Assignment Cycle

In this section, we propose a solution for an optimal assignment (scheduling) of stored entanglement called the entanglement assignment cycle, . The goal of is to achieve a minimal overall storage time at a given entanglement throughput tree.

Lemma 2

An entanglement assignment cycle can be determined by a weighted graph coloring method.

To determine the minimal overall storage time for a entanglement throughput tree, the conflict graph of that is constructed first. In the graph, each vertex corresponds to a directed link of (an entangled connection). An edge exists between two vertices of , if only the vertices (entangled connections) have a conflict. A conflict occurs if two (stored) entangled connects are associated to the same physical link. The problem is therefore to associate each link of a storage schedule (optimal assignment of stored entanglement), which includes the list of time slots when a given link can transmit the stored entangled states such that total number of time units is minimized. Therefore, our goal is to determine what entangled connects of the given should be scheduled in which time unit, such that the total storage time is minimal in .

Let τ ∈ {0, 1} be an indicator variable, defined as

For a periodic scheduling,where T is a period and i is a constant.

For an entangled connection n, let us define as the set of entangled connects n′ that are scheduled in the same time unit t, but the physical link can transmit only n or n′. As follows, for any ,

As it can be concluded, this problem is analogous to the coloring of the conflict graph .

Then, let us assume that each entangled connection (entangled link) has a weight w(n), which is defined aswhere F is the fidelity of entangled connection i and Fmax is the largest fidelity available. As follows, for a lower-fidelity entangled connection, the transmission of a given amount of information requires more time units.

As the weights are determined, the problem is analogous to a weighted graph coloring of the conflict graph , which is the assignment of at least w(n) distinct colors to each entangled connection n such that no two entangled connections sharing the same color interfere with each other on the physical link. For this purpose, a simple distributed weighted coloring algorithm[46] can be straightforwardly applied.

The method for the weighted coloring of conflict graph is summarized in Procedure 2.

Weighted Coloring of a Conflict Graph.

After an entanglement assignment cycle has been determined by the weighted graph coloring method, in the next step the time intervals between each time unit of a given cycle is computed. For this purpose, a linear programming method[44] can be applied.■

Quantum-Layer Optimization

In this section, we define the multi-objective optimization for the quantum layer. The multi-objective function covers a parallel optimization of quantum memory usage and entanglement throughput for a shortest path.

Theorem 1

At a given , the quantum layer of a quantum repeater network N can be optimized by a procedure that achieves a parallel minimization of the quantum memory usage time t, the maximization of entanglement throughput B, and the minimization of the number of entangled links between two arbitrary nodes.

Proof.

The aim of the procedure is to find an optimal entanglement throughput tree for which the overall storage time is minimal, that is,the entanglement throughput for all links is maximal,and the number of entangled links of a path is minimal, thuswhere is the shortest path.

The procedure is based on the fact that for a given with conflict graph , from the knowledge of weighted coloring of conflict graph , time intervals between each time unit of a given cycle, the required objective values of (10), (11), and (12) can be determined.

Let be the set of optimal entanglement throughput trees. Then the aim of the procedure is to determine .

The problem can be rewritten via a solution set X with decision variables[44]where n is the number of all links in a given quantum memory utilization graph , and x ∈ {0, 1} is defined as

Let be a solution i from solution set . Let refer to solution dominating solution , which isand for these relations, there is at least strict inequality[44]. If dominates all other possible solutions, then is a set of nondominated solutions.

For each set X, let κ refer to the set that contains the best nondominated solutions that have been found at a particular iteration. If κ changes, then the entanglement utilities of the links are updated. If κ is stationary, then the elements of set κ are used to update the entanglement utilities. The former serves as an improvement in the exploration process whereas the latter aims to yield more information via the best solutions[44].

At a given graph of quantum memory utilization, the procedure for the optimization of the quantum layer is defined as follows. The output of is an optimal set of entanglement throughput trees that realize the conditions of the multi-objective optimization function. The steps of the quantum layer optimization are summarized in Algorithm 1.

Note if there are no nondominated solutions, then the values of the weighting coefficients ρ, ς, and in Step 5 of Algorithm 1 can be selected according to the actual trade-off requirements.

The proof is therefore concluded here.■

Quantum Layer Optimization.

Classical-Layer Optimization

In this section we characterize the classical layer optimization procedure.

Theorem 2

The cost function of classical communications can be minimized by a procedure .

Letbe the cost function of classical communication of the multilayer optimization procedure, where Θ ∈ R is a p-dimensional real vector of an i-th system state of the quantum network, , , are the number of nodes that require the determination of optimal , and , , , refer to the number of classical steps required to find , and for a particular node g of network N, while , , and are the costs of classical link L used for the determination of , and for a given g.

Then, let introduce indices for Θ aswhere j is the index of a desired optimal system state, k is the index of an optimal system state reproduction step, l is the index of a non-optimal system state event[39].

Let assume that there is a set of S sub-states {Θ1, …,Θ} in the network, then the T total network state is evaluated as

Let be the cost function of classical communication at a given . Then, ifthen the system state evolves from j to j + 1 aswhere c(i) is the number of random system states, while u(j) quantifies a unit cost of system change[39,40].

Therefore, for a set of S sub-states {Θ1, …,Θ} and current vector , the C total cost of classical communication is yielded aswhich can be rewritten[39] aswhere A is the distribution-entity of a current system state, RA is the information transmission rate of A, ν is the distribution-entity of a system state, R is the information transmission rate of ν, while Θ( is the m-th element of a current network state vector Θ, and is the m-th element of Θ.

Using C(Θ, T(j, k, l)) (see (22)), an environment-dependent cost function C is defined aswhere M is a tuning parameter[39].

From (23), the cost function at a given Θ(j, k, l) is defined as

Using the proposed basic model, we define an optimization procedure of the classical-layer to achieve a C(Θ, T(j, k, l)) minimized cost function. The method is summarized in Algorithm 2.

These results conclude the proof.■

Classical Layer Optimization.

Large-Constrained Optimization

The optimization efficiency can be further improved for a large constrained network scenario. This step can be replayed by a different approach, which allows an optimization of the classical layer for an arbitrary constrained setting. The solution is based on the idea of system state merging.

Lemma 3

The optimization of the classical-layer can be extended to a large-constrained optimization by state merging.

An extended model of the classical-layer optimization is as follows.

Let , , be the number of nodes that require the determination of optimal , and . Then, let J be the objective function[34], aswhereandwhere the coefficients , , refer to the number of classical steps required to find , and at a particular network node g and time t, t = 1, …, T, while , , are the costs of classical link L at a particular g and t. Let assume that Θ(j, k, l), Θ(j, k, l) and Θ(j, k, l) are some system states of the classical layer subject of state merging. From these sub-states, the Θ(j, k, l) merged system is defined aswhere Φ is a merging factor, Φ ∈ [0, 1].

Using (29), an i-th system state Θ(j, k, l) is updated aswhereandwhere u is a uniform random number, x, x and x are random numbers, x, x, x ∈ [0, 1], Θ*(j, k, l) is a system state that minimizes the objective function J. Then for Q ∈ (a, b, c), the objective function value is compared with the objective function J(Θ(j, k, l)) of Θ(j, k, l), and is updated by the following rule:

The proof is concluded here.■

Cost Uncertainty of Large-Scaled Optimization

The determination of the minimal cost function (25) is can be approached by an output variable O = f(ψ), where ψ is the set of input variables, and f(·) is a function that transfers the uncertainty from the independent input random variables ψ to the output variable O[34]. The output set O can be rewritten aswhere q is the set of certain variables, while w is an input variable under certainty with probability function .

Then, for a given variable w, two pc(w) probability concentrations[34], pc(1)(w) and pc(2)(w) are defined asandwhere w, is the poth location of w[34], g = 1, 2, while ζ is a weighting factor.

Therefore, at a given (i, g) parameterization, where i = 1, …, z and g = 1, 2, O is expressed by variable aswhere is the mean of w, while w and w are the poth locations of w. Therefore, the problem is reduced to 2z deterministic equations[34], from which the mean and standard deviation of the output random variable can be computed.

Performance Evaluation

In this section we study the performance of the proposed quantum layer and classical layer optimization methods.

Quantum Layer Optimization

To study the convergence of the quantum layer optimization, we characterize the D(·) closeness function of the elements of the solution set from an optimal Pareto front, and the ζ(·) ratio of optimal solutions in the solution set.

Let be an i-th solution from the solution set found at a N finite number of iterations, and let refer to a solution z from at N → ∞. Then, let be the distance on the Pareto front[44-46] between and , defined aswith , and whereandwhile χ is a control parameter, χ > 0.

Then, let be a ratio of the solution sets and found at finite N and N → ∞ aswhere is the cardinality of set , while is the intersection of solution sets and (solutions included by both solution sets). Note, since the χ control parameter in (39) determines the D distance from the optimal set, χ also affects the ratio of the solution sets (43). If the value of χ is high, the D distance in (39) is low, that results in a high cardinality of the intersection set in (43). If χ is low, the D distance in (39) is high, therefore the cardinality of the intersection set is small.

The quantity of (39) for various N and χ are depicted in Fig. 4. From the results it can be concluded that as N increases, the distance significantly decreases, while the speed of convergence is controllable by χ. It also can be verified that the χ control parameter has a significant impact on , and the distance can be made arbitrarily small via a moderate value for N.

Figure 4

The distance in function of N and χ. (a) The values of for 0 < N ≤ 2 · 103 and 1 ≤ χ ≤ 10. (b) The values of for 0 < N ≤ 2 · 103 and 11 ≤ χ ≤ 20.

In Fig. 5, the quantity of (43) is illustrated for different values of N. As N increases the ratio increases significantly, while the χ control parameter has a moderate effect on the ratio of (43). The ratio exceeds 0.5 at N ≈ 0.5 · 103, and can be increased to arbitrarily high via a moderate increment in N.

Figure 5

The ratio in function of N and χ. (a) The values of for 0 < N ≤ 2 · 103 and 1 ≤ χ ≤ 10. (b) The values of for 0 < N ≤ 2 · 103 and 11 ≤ χ ≤ 20.

The performance of the quantum layer optimization method is therefore approachable via the distance function (39) and ratio (43). The analysis revealed that at moderate N values, the precision of the optimization method can be arbitrary high via the selection of the χ control parameter. It also has been concluded, that a high ratio of the solutions at a finite and moderate N, are identical to the solutions at the limit case of N → ∞.

Classical Layer Optimization

Let ϕ(i) be the step-size function defined for an i-th state aswherewhere and are some lower and upper bounds, , while f(i, j, k, l) is defined aswhere ω is a constant (restriction factor[41-43]), J* is the minimal cost function (for J, see (25)) at a particular parameter setting (i, j, k, l), J is the minimal cost function associated to the population.

Then, at an iteration number j let κ(j) be a ratio as

The step-size (44) in function of and κ(j) is depicted in Fig. 6.

Figure 6

(a) The ϕ(i) step-size function for , and κ(j) ∈ [0, 10]. (b) The ϕ(i) step-size function for , and κ(j) ∈ [0, 40].

From (44) and (46), the optimal cost function J* at a particular setting of (i, j, k, l) is yielded as

Then, letfrom which the iteration number can be rewritten asand also fix f(i, j, k, l) as

Then, from (49) and (51), J* can be evaluated at a given x in function J and ω, as

From (50) and (51), the function in (48) can be rewritten astherefore the cost function J* at a particular (49) and (51) can also be evaluated in function of the step size (44).

The J* cost function values associated to the ϕ(i) step size function values of Fig. 6 at J = 1 and J = 1 and ω = 100 are depicted in Fig. 7.

Figure 7

The J* cost function values associated to the ϕ(i) step size function value, (a) J = 1 and ω = 1, (b) J = 1 and ω = 100.

The analysis is therefore revealed that for any ϕ(i), the cost function J* increases with the κ(j) ratio, and the values of J* are tunable via the control parameter ω for arbitrary J and κ(j) values. The proposed classical layer optimization model is therefore flexible, and allows a dynamics adaption for diverse environmental settings. As future work, our aim is to provide a transmission analysis and comparisons with other schemes.

Conclusions

Here we conceived a multilayer optimization method for the quantum Internet. Multilayer optimization defines separate procedures for the optimization of the quantum layer and the classical layer. Quantum-layer optimization defines a multi-objective function by minimizing the total usage of quantum memories in the quantum nodes, maximizing the entanglement throughput over all entangled links, and reducing the number of entangled links between the arbitrary source and target quantum nodes. We defined the structure of the quantum memory utilization graph and entanglement throughput tree. The classical-layer optimization utilizes the fundaments of swarm intelligence for the minimization of the cost function. Since the proposed multilayer optimization method has no physical layer requirements, it can serve as a useful tool for quantum network communications and future quantum Internet.

Supplemental Information