Prerpints.org logo
Preprint
Article

A Best-Path Approach to the Design of a Hybrid Space-Ground Quantum Network with Dynamic Constraints

Altmetrics

Downloads

169

Views

167

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

26 February 2024

Posted:

27 February 2024

You are already at the latest version

Alerts
Abstract
The design and operation of quantum networks are both decisive in the current push towards a global quantum internet. Although space-enabled quantum connectivity has already been identified as a beneficial candidate for long-range quantum channels for over two decades, the architecture of a hybrid space-ground network is still a work in progress. Here, we propose an analysis of such a network based on a best-path approach, where either fiber- or satellite-based elementary links can be concatenated to form a repeater chain. The network consisting of quantum information processing nodes, equipped with both ground and space connections, is mapped into a graph structure, where edge weights represent the achievable secret key rates, chosen as the figure of merit for the network analysis. A weight minimization algorithm allows identifying the best path dynamically, i.e., as the weather conditions, stray light radiance, and satellite orbital position changes. From the results we conclude that satellite links will play a significant role on the future large-scale quantum internet, and both a constellation of satellites and significant advances in filtering technology are required to fulfill continuous coverage.
Keywords: 
Subject: Physical Sciences  -   Quantum Science and Technology

1. Introduction

The future quantum internet will enable connectivity between quantum information processing (QIP) nodes through the continuous distribution of entanglement [1,2]. In stricter terms, this translates into high-rate long-distance transmission of multi-partite quantum states exhibiting high fidelity with respect to maximally entangled states of the same dimension – for the simpler case of two QIP nodes, the fidelity with respect to a Bell state [3]. In general, the requirements of high rate and high fidelity play opposing roles as far as entanglement generation processes go, as exemplified by the widely used sources of entangled photon-pairs based on non-linear media [4]: by increasing the pump power, more photon-pairs are generated at the expense of multi-photon-pair states that diminish the overall state fidelity. On top of that, preservation of both rate and fidelity after long-distance transmission represents an even bigger challenge, albeit one that many research groups around the world are willing to undertake [5,6,7]. The core idea is to break down the long distances into smaller stretches, called elementary links, where entanglement can be efficiently distributed in a heralded manner, such that quantum memories preserve the successful results. A chain of elementary links acts as a quantum repeater, where the intrinsic losses of the channel are mitigated by the rounds of entanglement swapping preserved by the quantum memories [8]. Finally, the entanglement is distributed to the end nodes [9], thereby achieving a network of interconnected quantum processors. Although some speculate that the so-called "killer application" of a quantum internet has not yet been identified, the prospects of secure communication relying on quantum key distribution and increased quantum processing power due to distributed quantum computing, both on a global scale, drive the progress and push for innovative solutions.
Satellite-based quantum communications opens the door to long-distance entanglement distribution links with a reduced number of elementary links in the quantum repeater chain. Apart from the tens of kilometers of atmosphere, light propagates through vacuum on its path towards (or from) the satellite, so the dominating loss term is the so-called free-space path loss (FSPL), which scales quadratically with distance instead of exponentially, as the loss in an optical fiber [10]. In ground-based networks, quantum repeater nodes must be located at an accessible location and the fiber can rarely be installed in a straight line between two consecutive nodes, exemplified in the metropolitan experiments of [11,12]. The added benefit of independence concerning the topographical conditions of a certain node could further position satellites as an enabling quantum repeater node even for shorter distances. Furthermore, when a quantum internet that spans the entire globe is considered, the requirement of transoceanic links strongly favors the use of non-ground-based repeater nodes, achievable with space-based links. The expensive and complex endeavor of flying a space-qualified quantum payload on a satellite, coupled to the fact that current quantum memory platforms are not yet capable of supporting the long-distance links that make satellites beneficial, pose significant challenges for space-based quantum communications which, currently, attracts a small number of early adopters: the Micius satellite figures as the pioneer quantum satellite mission [13]; but it is not the only one [14]. Therefore, the main focus remains on ground-based quantum networks [12,15,16,17], and a shift in perspective requires both significant advances in the enabling technologies (sources of entangled photon pairs and quantum memories for light) as well as a clear definition of the benefits of a satellite-based link and the requirements such that these benefits can materialize [18,19,20,21].
Here, a contribution to the latter is put forth: an analysis of entanglement distribution rates between two arbitrary nodes in a quantum network that allows both ground- and satellite-based links. We start by breaking down the building blocks of a quantum repeater and rearranging them to fit a space segment all the while considering its practical constraints. This allows us to write entanglement distribution rate equations for both satellite- and ground-based repeater channels, which are then mapped onto the weights of edges in a graph representing the network. A best-path algorithm is then used to evaluate the beneficial architecture. The analysis is capable of incorporating dynamic conditions of the satellite link, such as the weather, stray light and the satellite orbit, for which realistic data is employed. We show that satellites present significant advantages for longer distance links, as expected. Also note that significant effort towards modeling and analyzing the achievable entanglement distribution rates in a practical quantum network has been made [22,23,24,25]. The present work is complementary, allowing for quick and broad analysis of practical quantum repeater architectures. The paper discusses the thresholds of satellite-based quantum communications, on the identified limits of the analysis as well as of the current technology, with conclusions on the potential avenues towards improvements.

2. Quantum Repeaters: Ground and Space Architectures

A quantum repeater leverages the multimodality (temporal, spectral, etc) enabled by both a quantum memory and the heralded swapping of entanglement to decouple the individual probabilities of success within a so-called elementary link and beats the bounds of direct transmission [26]. Although different choices of protocol, encoding, platform, and wavelengths exist [27], the unavoidable functional blocks of a quantum repeater can be narrowed down to four: sources of entanglement, quantum memories, Bell State Measurement units, and quantum channels. The requirements on the sources are such that entanglement is generated at the best compromise between high rate and high fidelity, in as many multiplexed modes as possible, and with an optical interface, so that optical photons can carry the quantum information with minimum attenuation and decoherence over the quantum channels. The quantum memories come in many flavours, and it is common for a single platform to implement both memory and source in the same physical system [5,7,28]. Due to the latter not being the general case [6], the functional blocks are split. The memory is expected to store a quantum state in a large number of modes with high efficiency. The Bell State Measurement (BSM) unit maps the joint state of two non-interacting photons on the base of maximally entangled bipartite states and is at the core of the entanglement swapping protocol [29]. It is expected to efficiently produce projection measurement results that can be classically conveyed to the quantum memories to swap entanglement in a heralded fashion [30]. Here, efficiencies are understood as the probability that the functional block performs as expected. For the source, this means the probability of successfully emitting a single pair of entangled photons for a given attempt (compared to the probability of emitting zero or two or more pairs). For the quantum memory, it means the probability of successfully storing and retrieving a photon from an entangled pair without disturbing its quantum state. In case of the BSM, it means the probability of successfully mapping and detecting the joint photon state.
Based solely on the natural disposition of the functional blocks within an elementary link [26,27,30], with source and memory on either side with a centralized BSM, it would be natural to assign the latter to the satellite. This way, two ground stations, each equipped with a source and a memory, would beam photons up to the satellite. The successful measurement results would be transmitted down to the ground stations, leaving the states stored in the memory in a joint entangled state. For a ground-based architecture, this configuration is extremely attractive, since the direction of propagation of the photonic qubits is irrelevant. However, the thin layer of atmosphere in the ground-to-satellite path close to the ground breaks the transmission direction symmetry for the space segment, as opposed to the fiber-based ground segment. The configuration in which the satellite receives a beam from a ground station is usually referred to as uplink, as opposed to downlink, where the satellite sends beams to the ground stations. In uplink, the beam experiences the effect of atmosphere at the first few kilometers of propagation, so the slightest angular deviation introduced there is translated into a significant beam spot displacement at the satellite, after hundreds of kilometers of propagation. On the other hand, downlink does not suffer from the same problem simply because the beam only interacts with the atmosphere at the very end of the propagation path. Differences in total loss between the two configurations, depending on the atmospheric turbulence, vary between 10 and 100 times [31]. Therefore, the downlink configuration, i.e., photonic qubits transmitted from the satellite to the ground stations, is chosen as a baseline in this work.
Downlink configuration restricts the design of a quantum repeater’s space segment containing a satellite since the source, rather than the BSM, is at the center-point between the two ground stations. Placing the source on the satellite meets practical constraints, since these have been demonstrated to operate in-orbit [13,14], but relays all the complexity to the ground stations. Due to the probabilistic nature of entanglement distribution over lossy channels, successful rounds must be heralded by a classical signal coming from a BSM-like unit so that the quantum memories can hold on to those specific modes. While proposals such as [32], combine the memories and the BSM heralding protocol into a single system, the functional blocks can be easily kept separated, in which case a source and a BSM must be included in both ground stations. In this final configuration, the satellite beams photon-pairs down to earth, which reach two ground stations. There, photons coming from space meet a photon generated on the ground in a BSM, whose pair is stored in the quantum memory. Based on the simultaneous heralding event of both ground-based BSMs, the states in the quantum memories are left in an entangled state following entanglement swapping. The block diagram representation of both space and ground repeater segments introduced so far is presented in Figure 1. Note that these are not the only possible architectures, but offer a practical short- to mid-term architecture aside from intuitive functional roles of each block.
The functional block description of quantum memory, source of entanglement, BSM, and quantum channel, equipped with efficiencies and fidelities, provides a powerful framework to analyze the overall efficiency of a given concatenation of space and ground segments with rounds of entanglement swapping: a so-called automated repeater chain (ARC) [27,33]. The ARC is assumed to perform subsequent attempts of entanglement distribution over each segment at a fixed rate R, which has the role, in a bipartite entanglement distribution scheme as considered here, of translating efficiencies into probabilities, and those into an average number of Bell pairs available at the end nodes. By enforcing minimum decoherence processes as the photonic qubits propagate through the quantum channels, negligible dark counts on the BSM detectors, and storage times long enough to ensure that the heralding signals reach the quantum memories in a timely manner (before the quantum states are lost), the performance of an ARC can be narrowed down to its achievable rate of entanglement distributed towards its end nodes. One can start by assigning efficiencies η S , η M , η B , and η c h ℓ A , η c h ℓ B to source, memory, BSM, and channels, respectively. Here, ℓ A , B denotes the optical path distance between the center point of an elementary link and either one of its nodes. It is worth noting that the efficiency of the channel is a function of its distance, which will take on a quadratic or exponential scaling in case of space or ground segments.
To arrive at the ground-based elementary link (EL) rate equation [27], one starts from the success probability of a single attempt, which, for a ground segment is simply the product η S 2 η B η c h : g ℓ A η c h : g ℓ B . This value can be subtracted from unity to represent the probability of no success, which, when taken to the power of the available number of multiplexing modes ( Γ ) available, represents the probability that no modes were successful within a given attempt. Subtracting the latter term from unity and multiplying by η M 2 equates to the probability of at least one mode being successfully swapped and stored into the quantum memories, the result of which is summarized in the expressions below:
P g n d E L ℓ A , ℓ B = 1 − 1 − η S 2 η B η c h : g ℓ A η c h : g ℓ B Γ η M 2 ,
η c h : g ℓ ∝ e − α ℓ .
The steps towards the space-based EL rate equation is the same, although the different configuration of the functional blocks must be accounted for as well as the difference in channel losses, as expressed below:
P s p c E L ℓ A , ℓ B = 1 − 1 − η S 3 η B 2 η c h : s ℓ A η c h : s ℓ B Γ η M 2 ,
η c h : s ℓ ∝ ( θ ℓ ) − 2 .
The factors α in Eq. 2 and θ in Eq. 4 denote the fiber attenuation coefficient and half-angle divergence of the free-space optical beam, respectively. It should also be noted that, although the distances ℓ A , B of a ground segment correspond directly to the distances between the nodes and the central point, those for space consider the distances towards the satellite. An ARC is formed through the concatenation of its n segments through n − 1 steps of successful entanglement swapping, leading to a final probability of entanglement distribution of the form:
P A R C = η B n − 1 ∏ i = 1 g P g n d E L ℓ A i , ℓ B i ∏ k = 1 s P s p c E L ℓ A k , ℓ B k ,
where g + s = n .
It is interesting to note that the potential benefit of the space segment relies heavily on whether the quadratic loss scaling can beat the exponential scaling of the fiber and the increased complexity of the ground station hardware. This point is made more stringent due to the fact that the distance covered by a ground segment can be approximated by the distance between the two ground stations, while for a space segment the distance depends on the altitude of the satellite. In other words, two ground stations separated by 100 km can be connected by 100 km of optical fibers, but the fact that the satellite is hovering at ≥500 km altitude causes the optical path for the space segment to be much longer. The expression of Eq. 5, however, tells only part of the story: as previously discussed, entanglement distribution in a quantum network relies on the transmission of entangled particles with high rate and high fidelity. The secret key rate (SKR) combines both requirement aspects of the ARCs into a single meaningful parameter [34] and is widely used as a figure of merit for the quality of quantum networks, although new benchmarking procedures have been proposed, e.g., [35]. The SKR allows for the evaluation of relevant practical issues; the most pressing for satellite quantum channels being the impact of stray light.
Stray light refers to unfiltered photons that are mixed in the path of the photonic qubits transmitted by the satellite. These reduce the quantum information capacity of the link. Although any decoherence process during transmission through the channel is neglected in the present analysis, stray light has a severe impact on the results of the BSM units, especially during the day, when solar radiance is maximum. In order to evaluate the impact of the unwanted photons onto the overall fidelity of the states distributed over the ARCs, a decoherence channel model is assumed. The states produced by the source are considered to be Werner states (noisy EPR-pairs) [36] with an associated density matrix of the form:
ρ S = W S | ϕ + ⟩ ⟨ ϕ + | + 1 − W S 4 I 4 ,
where | ϕ + ⟩ is a maximally-entangled Bell state, W S is the Werner parameter of the source – associated with its fidelity with respect to | ϕ + ⟩ –, and I is a 4 × 4 identity matrix. The final Werner parameter of the states available on the ground stations can be estimated by concatenating the Werner parameters of each functional block, which are assumed to be fixed except for the BSM.
The Werner parameter of the swapped states due to the action of the BSM can be estimated from the Coincidence-to-Accidental Ratio (CAR), which approximates the signal-to-noise ratio (SNR) achievable at the BSM. In the regime of low decoherence introduced by the quantum channel, the CAR provides an upper bound of the fringe visibility extracted from the coincidence measurements performed on a joint state ρ E L delivered by the ARC segment. The visibility, in turn, is directly associated to the Werner parameter of the link, where the term introduced by the BSM
W B ≈ C A R − 1 C A R + 1
is dominant. With this expression, one can connect the stray light radiance impinging on the BSM detectors to the overall entanglement quality generated in an ARC segment. First, the stray photon rate can be derived given the solar radiance (which can be approximated to have a black body radiation profile I s t r a y ); the receiver telescope parameters of aperture area ( A R X ) and solid angle ( Ω ); the photonic qubit bandwidth ( B λ ); and the operational wavelength ( λ ):
R s t r a y = I s t r a y Ω A R X B λ / h c λ ,
where h is Planck’s constant and c the speed of light. The rate R s t r a y is then translated into the term p s t r a y (the probability of measuring a stray photon within a time window of the detectors), with the assumption that the stray detection events are uncorrelated and, thus, homogeneously distributed over time. Therefore, p s t r a y = R s t r a y · δ t , where a time resolution δ t – limited by the detector’s jitter, or response time – is imposed. Finally, the photonic qubits are expected at a well-defined time window w d e t , therefore:
C A R = ∫ w d e t P g n d / s p c E L ℓ A , ℓ B d t ∫ w d e t p s t r a y 2 d t .
In Eq. 9, the term P g n d / s p c E L is not squared since it already presumes pairs of photons arriving at the BSM. From it, the estimated entanglement quality delivered by a given ARC segment (directly proportional to the overall Werner parameter W E L ) is directly calculated from the multiplication of the W B and the remaining individual Werner parameters of the other functional blocks of the elementary links ( W S , W M , for source and memory, respectively). Relevant for the SKR is the quantum bit error rate (QBER), approximated as 1 − W E L 2 [37], which can be used to calculate the entropy per distributed entangled pair over the elementary link:
S q E L = max 0 , 1 − 1 + κ H Q B E R ,
where κ is a constant representing the inefficiency of the classical steps of the quantum key distribution protocol, and H is the Shannon entropy function. Note that the SKR, calculated as the product
S K R E L = R S q E L P g n d / s p c E L ,
is heavily influenced by the entanglement generation probability of the elementary link, since it also impacts W B – as indicated in Eq. 9 – making it the most significant term to optimize. This is relevant for the best-path estimation results of Section 3, since the SKR can only be computed at the end nodes at the edges of an ARC, while the optimization takes into account each individual segment.

2.1. Modeling Free-Space Optical Channels

In the space segment, a satellite is the connecting point between two ground stations (A and B) and the entanglement generation probability will depend on the channel efficiencies η c h : s ( ℓ A ) and η c h : s ( ℓ B ) . Since the orbit of the satellite determines its position at any given point in time and the ground stations are static, these two terms may differ greatly. Evaluating η c h : s ( ℓ A , B ) follows from the identification of its four major contributing factors: the geometrical losses due to the beam divergence and limited receiver telescope dimensions, η f s l ; the pointing losses due to vibrations, atmospheric turbulence and relative movement between transmitter and receiver stations, η p l ; the effects of molecular and aerosol scattering in the atmosphere, η s c a t ; and the coupling efficiency of the incoming free-space optical beam into a medium that allows integration into the BSM unit, in general an optical fiber, η c o u p . Given the individual contributions, η c h : s ( ℓ A , B ) = η f s l η p l η s c a t η c o u p .
A detailed practical analysis of the first term can be found in [38], and takes into account the sizes of both receiving and transmitting telescopes, the size of the beam after propagation with a certain divergence, clipping of the beam at the transmitting telescope, and obscurations in the telescopes. Assuming on-axis propagation and the far-field approximation, an optimal ratio ρ = a / ω between the transmitter telescope’s diameter and beam waist can be derived [38]:
ρ ≈ 1.12 − 1.30 γ 2 + 2.12 γ 4 ,
where γ is the obscuration ratio. The design ratio ρ allows one to write η f s l given a few other design parameters, as follows:
η f s l = A T X A R X ( λ ℓ ) 2 2 ρ 2 e − ρ 2 − e − γ 2 ρ 2 2 ,
where A T X and A R X are the aperture areas of the transmitting and receiving telescopes, respectively, λ is the operational wavelength, and ℓ is the distance between the two telescopes. The beam waist ω – defined by the choice of A T X and Eq. 12 – will, in turn, define the half-angle beam divergence θ , in radians. This is then is used to compute η p l through [39]:
η p l = exp − 2 θ t θ 2 ,
where θ t is the half-angle pointing error, also in radians, which can be mitigated with an active pointing tracking system [13]. Differently from η f s l and η p l , the remaining term η s c a t is only empirically modeled by the following expression [40]:
η s c a t = T c exp − β 0 H 0 cos ( Z ) 1 − exp ℓ cos ( Z ) H 0 ,
where T c is the cloud coverage factor, β 0 the extinction coefficient at sea level, which is a function of the weather visibility V and the operational wavelength λ [41], H 0 = 6600 meters is the scale parameter, and Z is the Zenith angle between the ground station and the satellite. Note that, apart from H 0 and λ , all parameters of Eq. 15 are conditional, since they will change depending on the conditions of the channel. That is also the case for η c o u p , since the atmospheric conditions such as turbulence can impact its value differently depending on their intensity [40]. Assuming the ground stations are equipped not only with pointing tracking hardware to minimize the effects of Eq. 14, but also with an adaptive optics system [42], the optimistic assumption of η c o u p = 0.4 is used for the remainder of the analysis.

3. Hybrid Ground/Space ARC Network Simulation

Once equipped with both ground and space segment architectures (Figure 1) and their associated performance metrics (Eqs. 5, 9, and 11), one can attempt to optimize the entanglement generation rate between any two end nodes of an ARC embedded within a network. In order to do so, a framework with two major components is introduced: one is a potential future pan-European quantum network, with over 100 selected nodes representing major European cities distributed throughout the territory, which is presented in Figure 2; the second is a method to optimize the connections between any combination of end nodes in this network connected by an ARC. The latter is achieved by evaluating all the connections between the nodes in the network using the performance metrics introduced so far, associating them with edge weights of a graph, and searching for the best path.
The evaluation process of the hybrid ARC network analyzes each unique combination of end-node pairs after a ground/space graph is built, hereby called a "hybrid" graph. First, a static graph is constructed, where each node is connected to each other via a fiber link. This graph is only constructed once, since it is assumed not to be changing over time. The second graph, based on space segment connections, is constructed for each time window considered in the evaluation (minute accuracy for practical purposes), where the weather conditions (V and T c ) for each node of the network are selected from a historical database [43]. For the simulation results, the weather conditions on the first of January, 2024, are chosen without loss of generality. The time of the day is used to determine both the stray light radiance and the relative position of the satellite (in its polar orbit, at an altitude of 700km) with respect to the ground stations. Since each node is assumed to be able to communicate with any other node directly (either via space or via ground segments), the two graphs are merged into a single hybrid graph. To that end, each node-to-node connection is evaluated individually based on its performance metric. Here, one can use the fact that the rate R of the ARC is assumed to be the same for ground and space segments to introduce the secret key probability, S K P = S q E L P g n d / s p a E L . Since S K P ∈ 0 , 1 , setting the edge weights as
w E = − log 10 S K P ∈ [ 0 , ∞ )
ensures they are always positive and increase as the SKP (SKR) worsens. The segments with better performance are assigned to the edges of the hybrid graph: this allows pruning and converging into a hybridized ground/space graph. Then, Dijkstra’s algorithm [44], which minimizes the path weight in a graph, can be used to estimate the best path between the chosen nodes, as pictorially presented in Figure 3. The choice of ground-only, space-only, or hybrid, as well as the final SKR, are recorded and the process moves on to the next round.
The described procedure is not without its own limitations, starting from the fact that the weight estimated by Dijkstra’s algorithm at the end of each round does not represent the actual achievable SKR. The sum of the weights in a given path translates into the multiplication of the individual SKPs achieved by each segment of the path; although this would represent a meaningful value for the probability of entanglement generation (since probabilities can be multiplied), it does not necessarily equate the achievable SKR across a multi-segment path. In order to avoid edge cases that do not accurately represent the quality of an ARC, a tuple is assigned to each edge containing the values of SKP, S q , and P g n d / s p a E L ℓ A , ℓ B , where the latter term is used to tie-break which segment – ground or space – is chosen for the final graph. At the end of a round, the end-to-end SKR is calculated based on the total entanglement generation probability across the ARC (Eq. 5), and the total Werner parameter across the ARC:
W A R C = ∏ i = 1 n W i E L .
Another constraint is associated with the computational efficiency of the procedure: although Dijkstra’s algorithm can converge quickly to the final solution, the step of building the final hybrid graph for each time window represents a bottleneck. By realizing that the fiber channel does not change significantly over time, as opposed to the free-space channel, the ground connections can be computed once and stored a priori. In doing so, one notices that the number of ground-based elementary links that optimally connects two given nodes changes depending on the set of parameters: i.e., functional block efficiencies, Werner parameters, available modes, and total distance. Although assuming equal functional block parameters for ground and space segments and normalizing the efficiencies would allow more focus on the relative performance of the different architectures, ensuring that the results are obtained with realistic physical parameters is paramount. We set η S = 0.25 , η M = 0.75 , W S = W M = 0.99 , Γ = 300 – in keeping with the long-term (∼10-20 years) expectations for potential functional block physical platforms [27]. As will be seen in the next Section, the algorithm finds solutions in which the total distance in the network shown in Figure 2 is divided into several concatenated elementary links of ∼150 km each – in order to increase the overall rate.

4. Results

We start by evaluating the network of Figure 2 assuming a fixed satellite positioned at 700 km above Zurich, such that the satellite is within a constant line-of-sight for most pairs of ground stations. We further assume night-time operation only. With this approach, the performance impact of the space segment is decoupled from issues associated to the management of the satellite. In Figure 4, the results of all unique combinations of end node pairs, evaluated according to the protocol described in Section 3, are presented in a scatter plot on the left hand side and colour-coded, such that: blue represents ground-only; black represents space-only; and red represents hybrid. There are cases where the key material that can be extracted from the link is null (fail), which do not appear in the scatter plot, but are presented on the bottom right, where the total number for each individual channel choices are stratified. Note that the very choice of node placement introduces a bias in the average distance between QIPs, since metropolitan (shorter than 100km) links are not included in the network of Figure 2. This choice follows from the immediate realization, clear from Eqs. 1 and 3, that satellites are unlikely to benefit short-range quantum networks. Therefore, the comparison at short range distances is fruitless. Finally, Table 1 details the experimental parameters utilized in the simulation. Solar radiation, cloud cover and visibility, and city coordinate data are scraped from an open source weather API [43].
The results depicted in Figure 4 indicate that, for the fixed satellite scenario, a clear distinction between different regimes appears. The first one is associated to the fact that, above 150km distance between end nodes, multiple ground segments are chosen instead of a single segment, as expected from Eq. 1 and the parameters detailed in the previous Section. Second, is the regime of ground preference versus space preference, which happens around 500 km. This result is further highlighted in the histogram on the top right corner. To determine the threshold distance d t h r e s of the regime switch, the histogram distributions of choices involving satellites and no-satellites are fitted to a log-normal distribution and the intersection point between the fit curves is calculated and represented by the segmented red vertical line. The fact that the results indicate a clear preference for space segments after a threshold of ∼ 500 km is noteworthy: although the preference is expected, identifying a practical regime switch can be significant. Also noteworthy is the large spread of black dots, associated to the dynamic nature of the weather conditions. Links with similar distance can have prominent differences in performance due to the weather. Finally, the expressive incidence of hybrid solutions indicates that the possibility to switch between space and ground segments is a valuable feature of the future quantum internet since, even for large distance links, the ARC benefits from having both available.
Moving on from the ideal case where the satellite and weather conditions are static requires taking into account the dynamic variations of the free-space optical channel. To do so, the procedure detailed in Section 3 can be employed to evaluate the performance of the network over time. Unfortunately, combining the dynamics of each different unique combination of end node couples into a single figure is not straightforward. Therefore, we select three combinations of end nodes, expressing different length regimes, and evaluate its performance with a minute accuracy. Here, we also evaluate different satellite constellations of up to 30 satellites in a synchronous polar orbit. The results, presented in Figure 5, uncover the complex interplay between end-node distances and connection choices given the satellite availability (set to 25), the end-to-end fidelity – heavily influenced by the stray light –, and the weather conditions. The end-to-end fidelity is calculated taking into account spectral, temporal, and spatial filters available at the ground stations. The spectral and temporal filters are designed based on the ARC rate R = 10 GHz, so the stray light radiance is spectrally filtered accordingly, and the CAR, calculated following Eq. 9, uses a detection window inversely proportional to the rate. The spatial filter performance is calculated based on a 50 μ r a d field-of-view (FOV), typical for modern satellite and ground station telescopes, and is strongly dependent on the achievable pointing accuracy. As the results of Figure 6 demonstrate, improving the pointing accuracy and further filtering the signal spatially has a two-pronged effect: although it reduces the amount of photonic qubits that are coupled into the ground stations, it also significantly reduces the effect of stray light, and is an important step towards enabling day-time operation. For the results described here, the maximum count rate of the detector is assumed to be in the GHz regime (a parameter achievable in the long-term scenario), and an FOV of 50 μ rad is chosen.
Clear from the figures is the fact that, over a day, the connection choices change significantly and the realistic performance cannot rely solely on the ideal case analysis of Figure 4. In particular, one could question how the value of d t h r e s changes based on the analysis over a longer time-frame, such as a full day, as considered in Figure 5. Fortunately, the chosen figure of merit of key generation provides one with a quantity that can be combined: the key material extracted over a day, for instance. Combining this new performance evaluation parameter with the analysis based on d t h r e s , it is possible to analyze the proposed network for each unique end-node combination including the dynamic constraints. Furthermore, the impact of a constellation of satellites can be extracted, which can then be associated to the total cost of deployment of the network’s space segment. In Figure 7, the values of d t h r e s are depicted, indicating that a similar regime switch as identified in Figure 4 (460km) can be achieved with a number of available satellites in the range of 20.

5. Conclusion

Interconnecting quantum processing nodes via entanglement distribution is the task of the quantum internet, which will require all the efforts available to become a reality. This includes opening the door to complex forms of connectivity, such as a space segment. In this paper, a robust analysis methodology, capable of including static and dynamic constraints of the quantum channels of an Automated Repeater Chain has been proposed and used to evaluate performance estimates. A clear regime switch between ground/space segment choices has been identified, and two major components (narrower spatial filtering and a satellite constellation) have been highlighted as necessary to improve the overall performance. Notwithstanding, the possibility to switch connections dynamically stands out as an important feature of the future European quantum internet; for that, weather forecast and a performance prediction tool, such as the one set forth here, will be crucial for network management. Extending the analysis to include more complex space segment architectures (e.g., [45]), identifying physical platforms to implement different functional blocks (e.g., [32]), and improve the computation speed are essential points in the development of this framework, and are left as future research directions.

Author Contributions

Conceptualization, G.C.A; methodology, G.C.A., Y.J., and D.L.B; software, D.L.B. and Y.J.; validation, D.L.B. and G.C.A; writing—original draft preparation, G.C.A.; writing—review and editing, G.C.A., Y.J., B.P.F.D. and D.L.B.; project administration, B.P.F.D; funding acquisition, G.C.A. and B.P.F.D. All authors have read and agreed to the published version of the manuscript.

Funding

This study was conducted under a Samenwerkings Middelen Onderzoek (SMO) research grant.

Data Availability Statement

Publicly available datasets were analyzed in this study, which can be found here: https://www.open-meteo.com. Simulation results can be made available upon reasonable request from the corresponding author.

Acknowledgments

G.C.A. and D.L.B acknowledge fruitful discussions with F. Silvestri on free-space optical link modeling and R. Saathof on atmospheric turbulence analysis.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
QIP Quantum Information Processing
FSPL Free-Space Path Loss
BSM Bell State Measurement
ARC Automated Repeater Chain
EL Elementary Link
SNR Signal-to-Noise Ratio
SKR Secret Key Rate
EPR Einstein, Podolsky, and Rosen
CAR Coincidence-to-Accidental Ratio
QBER Quantum Bit Error Rate
SKP Secret Key Probability
API Application Programming Interface
FOV Field-of-View

References

  1. Wehner, S.; Elkouss, D.; Hanson, R. Quantum internet: A vision for the road ahead. Science 2018, 362, eaam9288. [Google Scholar] [CrossRef]
  2. Kimble, H.J. The quantum internet. Nature 2008, 453, 1023–1030. [Google Scholar] [CrossRef]
  3. Bell, J.S. On the einstein podolsky rosen paradox. Physics Physique Fizika 1964, 1, 195. [Google Scholar] [CrossRef]
  4. Aspect, A.; Grangier, P.; Roger, G. Experimental tests of realistic local theories via Bell’s theorem. Physical review letters 1981, 47, 460. [Google Scholar] [CrossRef]
  5. Pompili, M.; Hermans, S.L.; Baier, S.; Beukers, H.K.; Humphreys, P.C.; Schouten, R.N.; Vermeulen, R.F.; Tiggelman, M.J.; dos Santos Martins, L.; Dirkse, B.; others. Realization of a multinode quantum network of remote solid-state qubits. Science 2021, 372, 259–264. [Google Scholar] [CrossRef]
  6. Lago-Rivera, D.; Grandi, S.; Rakonjac, J.V.; Seri, A.; de Riedmatten, H. Telecom-heralded entanglement between multimode solid-state quantum memories. Nature 2021, 594, 37–40. [Google Scholar] [CrossRef]
  7. Yu, Y.; Ma, F.; Luo, X.Y.; Jing, B.; Sun, P.F.; Fang, R.Z.; Yang, C.W.; Liu, H.; Zheng, M.Y.; Xie, X.P.; others. Entanglement of two quantum memories via fibres over dozens of kilometres. Nature 2020, 578, 240–245. [Google Scholar] [CrossRef]
  8. Briegel, H.J.; Dür, W.; Cirac, J.I.; Zoller, P. Quantum repeaters: the role of imperfect local operations in quantum communication. Physical Review Letters 1998, 81, 5932. [Google Scholar] [CrossRef]
  9. Duan, L.M.; Lukin, M.D.; Cirac, J.I.; Zoller, P. Long-distance quantum communication with atomic ensembles and linear optics. Nature 2001, 414, 413–418. [Google Scholar] [CrossRef]
  10. Liorni, C.; Kampermann, H.; Bruß, D. Quantum repeaters in space. New Journal of Physics 2021, 23, 053021. [Google Scholar] [CrossRef]
  11. Valivarthi, R.; Puigibert, M.l.G.; Zhou, Q.; Aguilar, G.H.; Verma, V.B.; Marsili, F.; Shaw, M.D.; Nam, S.W.; Oblak, D.; Tittel, W. Quantum teleportation across a metropolitan fibre network. Nature Photonics 2016, 10, 676–680. [Google Scholar] [CrossRef]
  12. Liu, J.L.; Luo, X.Y.; Yu, Y.; Wang, C.Y.; Wang, B.; Hu, Y.; Li, J.; Zheng, M.Y.; Yao, B.; Yan, Z. ; others. A multinode quantum network over a metropolitan area. arXiv 2023, arXiv:2309.00221 2023. [Google Scholar]
  13. Lu, C.Y.; Cao, Y.; Peng, C.Z.; Pan, J.W. Micius quantum experiments in space. Reviews of Modern Physics 2022, 94, 035001. [Google Scholar] [CrossRef]
  14. Villar, A.; Lohrmann, A.; Bai, X.; Vergoossen, T.; Bedington, R.; Perumangatt, C.; Lim, H.Y.; Islam, T.; Reezwana, A.; Tang, Z.; others. Entanglement demonstration on board a nano-satellite. Optica 2020, 7, 734–737. [Google Scholar] [CrossRef]
  15. Zhou, Y.; Malik, P.; Fertig, F.; Bock, M.; Bauer, T.; van Leent, T.; Zhang, W.; Becher, C.; Weinfurter, H. Long-lived quantum memory enabling atom-photon entanglement over 101 km telecom fiber. arXiv 2023, arXiv:2308.08892 2023. [Google Scholar] [CrossRef]
  16. Yang, J.; Jiang, Z.; Benthin, F.; Hanel, J.; Fandrich, T.; Joos, R.; Bauer, S.; Kolatschek, S.; Hreibi, A.; Rugeramigabo, E.P. ; others. High-rate intercity quantum key distribution with a semiconductor single-photon source. arXiv 2023, arXiv:2308.15922 2023. [Google Scholar]
  17. Krutyanskiy, V.; Canteri, M.; Meraner, M.; Krcmarsky, V.; Lanyon, B. Multimode ion-photon entanglement over 101 kilometers of optical fiber. arXiv 2023, arXiv:2308.08891 2023. [Google Scholar]
  18. Miller, A.V.; Pismeniuk, L.V.; Duplinsky, A.V.; Merzlinkin, V.E.; Plukchi, A.A.; Tikhonova, K.A.; Nesterov, I.S.; Sevryukov, D.O.; Levashov, S.D.; Fetisov, V.V.; others. Vector—towards quantum key distribution with small satellites. EPJ Quantum Technology 2023, 10, 52. [Google Scholar] [CrossRef]
  19. Abasifard, M.; Cholsuk, C.; Pousa, R.G.; Kumar, A.; Zand, A.; Riel, T.; Oi, D.K.; Vogl, T. The ideal wavelength for daylight free-space quantum key distribution. arXiv 2023, arXiv:2303.02106 2023. [Google Scholar] [CrossRef]
  20. Wallnöfer, J.; Hahn, F.; Gündoğan, M.; Sidhu, J.S.; Wiesner, F.; Walk, N.; Eisert, J.; Wolters, J. Simulating quantum repeater strategies for multiple satellites. Communications Physics 2022, 5, 169. [Google Scholar] [CrossRef]
  21. Khmelev, A.V.; Duplinsky, A.V.; Bakhshaliev, R.M.; Ivchenko, E.I.; Pismeniuk, L.V.; Mayboroda, V.F.; Nesterov, I.S.; Chernov, A.N.; Trushechkin, A.S.; Kiktenko, E.O. ; others. Eurasian-Scale Experimental Satellite-based Quantum Key Distribution with Detector Efficiency Mismatch Analysis. arXiv 2023, arXiv:2310.17476 2023. [Google Scholar]
  22. Li, Z.; Xue, K.; Li, J.; Chen, L.; Li, R.; Wang, Z.; Yu, N.; Wei, D.S.; Sun, Q.; Lu, J. Entanglement-assisted quantum networks: Mechanics, enabling technologies, challenges, and research directions. IEEE Communications Surveys & Tutorials 2023. [Google Scholar]
  23. Van Milligen, E.A.; Jacobson, E.; Patil, A.; Vardoyan, G.; Towsley, D.; Guha, S. Entanglement Routing over Networks with Time Multiplexed Repeaters. arXiv 2023, arXiv:2308.15028 2023. [Google Scholar]
  24. Labay-Mora, A.; da Silva, F.F.; Wehner, S. Reducing hardware requirements for entanglement distribution via joint hardware-protocol optimization. arXiv 2023, arXiv:2309.11448 2023. [Google Scholar]
  25. Oliveira, R.; Oliveira, R.; Bernardes, N.K.; Chaves, R. Statistical properties and repetition rates for a quantum network with geographical distribution of nodes. arXiv 2023, arXiv:2312.09130 2023. [Google Scholar] [CrossRef]
  26. Pietx-Casas, O.; do Amaral, G.C.; Chakraborty, T.; Berrevoets, R.; Middelburg, T.; Slater, J.A.; Tittel, W. Spectrally multiplexed Hong–Ou–Mandel interference with weak coherent states. Applied Optics 2023, 62, 3284–3288. [Google Scholar] [CrossRef]
  27. Askarani, M.F.; Chakraborty, K.; Do Amaral, G.C. Entanglement distribution in multi-platform buffered-router-assisted frequency-multiplexed automated repeater chains. New Journal of Physics 2021, 23, 063078. [Google Scholar] [CrossRef]
  28. Beukers, H.K.; Pasini, M.; Choi, H.; Englund, D.; Hanson, R.; Borregaard, J. Tutorial: Remote entanglement protocols for stationary qubits with photonic interfaces. arXiv 2023, arXiv:2310.19878 2023. [Google Scholar]
  29. Bennett, C.H.; Brassard, G.; Crépeau, C.; Jozsa, R.; Peres, A.; Wootters, W.K. Teleporting an unknown quantum state via dual classical and Einstein-Podolsky-Rosen channels. Physical review letters 1993, 70, 1895. [Google Scholar] [CrossRef]
  30. Sinclair, N.; Saglamyurek, E.; Mallahzadeh, H.; Slater, J.A.; George, M.; Ricken, R.; Hedges, M.P.; Oblak, D.; Simon, C.; Sohler, W.; others. Spectral multiplexing for scalable quantum photonics using an atomic frequency comb quantum memory and feed-forward control. Physical review letters 2014, 113, 053603. [Google Scholar] [CrossRef]
  31. Maharjan, N.; Devkota, N.; Kim, B.W. Atmospheric Effects on Satellite–Ground Free Space Uplink and Downlink Optical Transmissions. Applied Sciences 2022, 12, 10944. [Google Scholar] [CrossRef]
  32. Chen, K.C.; Dhara, P.; Heuck, M.; Lee, Y.; Dai, W.; Guha, S.; Englund, D. Zero-Added-Loss Entangled-Photon Multiplexing for Ground-and Space-Based Quantum Networks. Physical Review Applied 2023, 19, 054029. [Google Scholar] [CrossRef]
  33. Lütkenhaus, N.; Calsamiglia, J.; Suominen, K.A. Bell measurements for teleportation. Physical Review A 1999, 59, 3295. [Google Scholar] [CrossRef]
  34. Rozpędek, F.; Goodenough, K.; Ribeiro, J.; Kalb, N.; Vivoli, V.C.; Reiserer, A.; Hanson, R.; Wehner, S.; Elkouss, D. Parameter regimes for a single sequential quantum repeater. Quantum Science and Technology 2018, 3, 034002. [Google Scholar] [CrossRef]
  35. Helsen, J.; Wehner, S. A benchmarking procedure for quantum networks. npj Quantum Information 2023, 9, 17. [Google Scholar] [CrossRef]
  36. Werner, R.F. Quantum states with Einstein-Podolsky-Rosen correlations admitting a hidden-variable model. Physical Review A 1989, 40, 4277. [Google Scholar] [CrossRef] [PubMed]
  37. Gisin, N.; Ribordy, G.; Tittel, W.; Zbinden, H. Quantum cryptography. Reviews of modern physics 2002, 74, 145. [Google Scholar] [CrossRef]
  38. Klein, B.J.; Degnan, J.J. Optical antenna gain. 1: Transmitting antennas. Applied optics 1974, 13, 2134–2141. [Google Scholar] [CrossRef] [PubMed]
  39. Degnan, J.J. Millimeter accuracy satellite laser ranging: a review. Contributions of space geodesy to geodynamics: technology 1993, 25, 133–162. [Google Scholar]
  40. Vasylyev, D.; Vogel, W.; Moll, F. Satellite-mediated quantum atmospheric links. Physical Review A 2019, 99, 053830. [Google Scholar] [CrossRef]
  41. Kim, I.I.; McArthur, B.; Korevaar, E.J. Comparison of laser beam propagation at 785 nm and 1550 nm in fog and haze for optical wireless communications. Optical wireless communications III. Spie, 2001, Vol. 4214, pp. 26–37.
  42. Broekens, K.; Klop, W.; Moens, T.; Eschen, M.; do Amaral, G.C.; Silverstri, F.; Oosterwijk, A.; Visser, M.; Doelman, N.; Kaffa, L. ; others. Adaptive optics pre-correction demonstrator for terabit optical communication. International Conference on Space Optics—ICSO 2022. SPIE, 2023, Vol. 12777, pp. 758–766.
  43. Free Weather API. https://open-meteo.com/. Accessed: 2024-01-17.
  44. Dijkstra, E.W. A note on two problems in connexion with graphs. In Edsger Wybe Dijkstra: His Life, Work, and Legacy; 2022; pp. 287–290.
  45. Gündoğan, M.; Sidhu, J.S.; Henderson, V.; Mazzarella, L.; Wolters, J.; Oi, D.K.; Krutzik, M. Proposal for space-borne quantum memories for global quantum networking. npj Quantum Information 2021, 7, 128. [Google Scholar] [CrossRef]
Preprints 99894 g001
Figure 1. Basic space (top) and ground (bottom) segments of a quantum repeater. Dotted black lines represent free-space optics quantum channels, solid grey lines represent fiber optics quantum channels, and solid red arrows denote classical communication rounds. EPPS stands for Entangled Photon-Pair Source.
Figure 1. Basic space (top) and ground (bottom) segments of a quantum repeater. Dotted black lines represent free-space optics quantum channels, solid grey lines represent fiber optics quantum channels, and solid red arrows denote classical communication rounds. EPPS stands for Entangled Photon-Pair Source.
Preprints 99894 g001
Preprints 99894 g002
Figure 2. Conceptual future Pan-European quantum network doted with Quantum Information Processing nodes (yellow triangles) connected via ARCs consisting of either ground or space segments (not depicted).
Figure 2. Conceptual future Pan-European quantum network doted with Quantum Information Processing nodes (yellow triangles) connected via ARCs consisting of either ground or space segments (not depicted).
Preprints 99894 g002
Preprints 99894 g003
Figure 3. Pictorial representation of the best path estimation routine, where solid gray lines represent ground segments, dashed black lines represent space segments, the nodes are represented by yellow circles, and the best path by the green arrow. Left and right images represent situations where the same node selection produces different best path results (space-only or hybrid) due to dynamic constraints such as weather and stray light conditions, or the relative position of the satellite.
Figure 3. Pictorial representation of the best path estimation routine, where solid gray lines represent ground segments, dashed black lines represent space segments, the nodes are represented by yellow circles, and the best path by the green arrow. Left and right images represent situations where the same node selection produces different best path results (space-only or hybrid) due to dynamic constraints such as weather and stray light conditions, or the relative position of the satellite.
Preprints 99894 g003
Preprints 99894 g004
Figure 4. Channel performance for a hybrid ground/space European quantum network, as presented in Figure 2; the ARC rate R is set to 10 GHz. On the left, the colour-coded scatter plot presents expected performance in terms of SKR versus node distance. In the excerpt, cases where the ground-based ARC decides on breaking the channel in more than one segment are highlighted in light blue, which happens mostly above 150km, as expected. The performance results are agglutinated into the histogram presented on the top right, where the threshold that distinguishes the choice between ground and space is highlighted ( d t h r e s ). Finally, the statistical distribution of the channel choices is presented on the bottom right corner, where the cases where key material is null (fail) is also presented.
Figure 4. Channel performance for a hybrid ground/space European quantum network, as presented in Figure 2; the ARC rate R is set to 10 GHz. On the left, the colour-coded scatter plot presents expected performance in terms of SKR versus node distance. In the excerpt, cases where the ground-based ARC decides on breaking the channel in more than one segment are highlighted in light blue, which happens mostly above 150km, as expected. The performance results are agglutinated into the histogram presented on the top right, where the threshold that distinguishes the choice between ground and space is highlighted ( d t h r e s ). Finally, the statistical distribution of the channel choices is presented on the bottom right corner, where the cases where key material is null (fail) is also presented.
Preprints 99894 g004
Preprints 99894 g005
Figure 5. Relevant dynamic parameters with minute-accuracy for 3 different combinations of end nodes: secret key generation probability (SKP), BSM fidelity, and weather-induced loss. The connection choices at each time window inherit the colour-coding scheme. ARC distances are chosen to represent different regimes, according to the results of Figure 4.
Figure 5. Relevant dynamic parameters with minute-accuracy for 3 different combinations of end nodes: secret key generation probability (SKP), BSM fidelity, and weather-induced loss. The connection choices at each time window inherit the colour-coding scheme. ARC distances are chosen to represent different regimes, according to the results of Figure 4.
Preprints 99894 g005
Preprints 99894 g006
Figure 6. Dependency of stray-light-induced dark count rate and BSM fidelity during day-time operation (full solar radiance) for different FOVs of the ground station telescope.
Figure 6. Dependency of stray-light-induced dark count rate and BSM fidelity during day-time operation (full solar radiance) for different FOVs of the ground station telescope.
Preprints 99894 g006
Preprints 99894 g007
Figure 7. Trend of d t h r e s as a function of the number of available satellites, with the ideal case depicted as a black horizontal dotted line. The fitted curve is proportional to 1 x and, when extrapolated, indicates that no more than 30 satellites are necessary to operate the space segment optimally.
Figure 7. Trend of d t h r e s as a function of the number of available satellites, with the ideal case depicted as a black horizontal dotted line. The fitted curve is proportional to 1 x and, when extrapolated, indicates that no more than 30 satellites are necessary to operate the space segment optimally.
Preprints 99894 g007
Table 1. Table of parameters
Table 1. Table of parameters
Parameter Symbol Value
Wavelength λ 1550 nm
Source efficiency η S 0.25
Source fidelity F S 0.99
Source rate R 10 GHz
Memory efficiency η M 0.75
Memory fidelity F M 0.99
Detector efficiency η D E T 0.99
Detector jitter δ t 100 ps
Optical fiber attenuation α 0.146 dB/km
BSM efficiency η B = 0.5 × η D E T 2 0.49005
Coupling efficiency η c o u p 0.4
Number of modes Γ 300
Satellite obscuration γ s a t 0
Satellite aperture d s a t 0.3 m
Satellite pointing error θ t 1.2 μ rad
Field-of-view OGS FOV 50 μ rad
OGS aperture d g r o u n d 1.2 m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated